请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 114|回复: 6

[Dalton] 请问gaussian和dalton的计算结果如何增加有效数字(数值精确度)

[复制链接]

8

帖子

0

威望

35

eV
积分
43

Level 2 能力者

发表于 2019-2-11 23:43:05 | 显示全部楼层 |阅读模式
如题,例如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
如何显示后面的几位?

谢谢各位!

2398

帖子

23

威望

8052

eV
积分
10910

Level 6 (一方通行)

发表于 7 天前 | 显示全部楼层
似乎只能直接改Fortran代码。

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

评分

参与人数 1eV +2 收起 理由
sobereva + 2

查看全部评分

1万

帖子

25

威望

2万

eV
积分
38157

管理员

公社社长

发表于 7 天前 | 显示全部楼层
-0.00000018已经保留到足够的小数位数了,要求再高也没意义,实际用的计算级别、数值层面的误差因素也不可能把数据算到那么精确
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

347

帖子

16

威望

2468

eV
积分
3135

Level 5 (御坂)

发表于 7 天前 | 显示全部楼层
可以根据gaussian打印出来的原子单位的跃迁偶极矩和激发能,根据振子强度的公式自己算。跃迁偶极矩和激发能在大多数情况下都远小于1,平方再相乘以后就会更小,能显示出来的有效数字就少了

8

帖子

0

威望

35

eV
积分
43

Level 2 能力者

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

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

1万

帖子

25

威望

2万

eV
积分
38157

管理员

公社社长

发表于 7 天前 | 显示全部楼层
改程序很简单,用grep就很容易定位到对应的源文件,这是最简单的解决当前问题的方法
就算乘以219474.6363转成cm-1,更后面的位数产生的影响依然可以忽略不计,本身计算过程各个层面产生的误差都足以影响到那种程度,如果纠结在这种程度的数值差异上并且做出过度讨论,别人自然会质疑结论的化学意义。我认为默认情况给出的位数已经足够了,如果你在dalton论坛上建议开发者把保留的位数增加,我估计开发者也会持有我这样的态度。
本来dalton算得也慢,倒不如用ORCA干这事
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

8

帖子

0

威望

35

eV
积分
43

Level 2 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
sobereva 发表于 2019-2-12 12:06
改程序很简单,用grep就很容易定位到对应的源文件,这是最简单的解决当前问题的方法
就算乘以219474.6363 ...

好的,谢谢老师。我试下ORCA。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-2-19 03:43 , Processed in 0.148523 second(s), 24 queries .

快速回复 返回顶部 返回列表