计算化学公社

标题: 请问gaussian和dalton的计算结果如何增加有效数字（数值精确度） [打印本页]

作者
Author: shaowenhao28 时间: 2019-2-11 23:43
标题: 请问gaussian和dalton的计算结果如何增加有效数字（数值精确度）
如题，例如dalton算Sn-Tn 体系内转换旋轨耦合矩阵元的输出：

@ Transition moment <B | A - <A> | C> in a.u. for
@ A operator label, symmetry, spin:    Z1MNF-SO 1 1
@ B excited state no., symmetry, spin:          2 1 0
@ C excited state no., symmetry, spin:          2 1 1

@ B and C excitation energies, moment: 0.11809436 0.11194147  -0.00000018

以原子单位表示的矩阵元为-0.00000018，如何能显示18的后面几位？或如何用科学计数法表示？

似乎S0-Tn体系内转换的输出格式就是默认以科学计数法表示的：
@ Operator type: Z1MNF-SO
@ Spin-orbit coupling constant (Z1MNF-SO) :  4.46302721E-03 cm-1
                                          2.03350484E-08 au

另外，高斯计算TDDFT的输出文件中，振子强度通常只保留到小数点后第四位
Excited State 1:    Triplet-A    2.7147 eV  456.72 nm  f=0.0000  <S**2>=2.000
如何显示后面的几位？

谢谢各位！

作者
Author: liyuanhe211 时间: 2019-2-12 01:10
似乎只能直接改Fortran代码。

Multiwfn在一些激发态计算时会重新输出振子强度，所以振子强度可以用.fch文件和加了IOp(9/40)的.out文件用Multiwfn计算，可见手册的4.18.10节的例子（高斯输出为0.3935，Multiwfn 0.3935306），不过高斯的跃迁系数只输出到小数点后第五位，我不确定这会不会限制输出的振子强度的精度。

作者
Author: sobereva 时间: 2019-2-12 05:05
-0.00000018已经保留到足够的小数位数了，要求再高也没意义，实际用的计算级别、数值层面的误差因素也不可能把数据算到那么精确

作者
Author: beefly 时间: 2019-2-12 08:06
可以根据gaussian打印出来的原子单位的跃迁偶极矩和激发能，根据振子强度的公式自己算。跃迁偶极矩和激发能在大多数情况下都远小于1，平方再相乘以后就会更小，能显示出来的有效数字就少了

作者
Author: shaowenhao28 时间: 2019-2-12 11:42

sobereva 发表于 2019-2-12 05:05
-0.00000018已经保留到足够的小数位数了，要求再高也没意义，实际用的计算级别、数值层面的误差因素也不可 ...

但是实际上dalton在计算Tn到基态的SOCME时是可以显示成科学计数法的： 2.03350484E-08 au，这也超出了-0.00000018的小数位数。而且原子单位换算成cm-1之后数值上会大很多；因为研究需要有时候需要比较两个都比较小的SOCME，所以有了这样的一个疑惑。除了改程序代码外不能直接在计算时像计算Tn到基态的SOCME那样显示科学技术法吗？谢谢

作者
Author: sobereva 时间: 2019-2-12 12:06
改程序很简单，用grep就很容易定位到对应的源文件，这是最简单的解决当前问题的方法
就算乘以219474.6363转成cm-1，更后面的位数产生的影响依然可以忽略不计，本身计算过程各个层面产生的误差都足以影响到那种程度，如果纠结在这种程度的数值差异上并且做出过度讨论，别人自然会质疑结论的化学意义。我认为默认情况给出的位数已经足够了，如果你在dalton论坛上建议开发者把保留的位数增加，我估计开发者也会持有我这样的态度。
本来dalton算得也慢，倒不如用ORCA干这事

作者
Author: shaowenhao28 时间: 2019-2-12 13:31

sobereva 发表于 2019-2-12 12:06
改程序很简单，用grep就很容易定位到对应的源文件，这是最简单的解决当前问题的方法
就算乘以219474.6363 ...

好的，谢谢老师。我试下ORCA。

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