计算化学公社

标题: DALTON 老问题求助 [打印本页]

作者
Author:
ZHANGZY    时间: 2017-5-11 16:28
标题: DALTON 老问题求助
我也遇到了“小范范1989”的问题,再发布求助
同时,借助您的帖子,我咨询一个问题,希望高手指点一下。
就是.dal文件中:
.PROPRT
X1MNF-SO
.PROPRT
Y1MNF-SO
.PROPRT
Z1MNF-SO

计算之后,我们得到的有
@ Transition moment <B | A - <A> | C> in a.u. for
@ A operator label,    symmetry, spin:      X1MNF-SO    1    1
@ B excited state no., symmetry, spin:             1    1    0
@ C excited state no., symmetry, spin:             2    1    1

@ B and C excitation energies, moment:   0.09025588   0.08802946   0.00000007


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

@ B and C excitation energies, moment:   0.09025588   0.08802946   0.00000005


@ 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:             1    1    0
@ C excited state no., symmetry, spin:             2    1    1

@ B and C excitation energies, moment:   0.09025588   0.08802946  -0.00000003
上面的数据,最后选取数值的时候:
1:三个方向,这个最后是平方加和,再开方? 把X1MNF-SO,Y1MNF-SO,Z1MNF-SO理解为三个方向?
2:还是说:平方加和,除以3.再开方? 把X1MNF-SO,Y1MNF-SO,Z1MNF-SO理解为是三重态的三个简并态?
以上两种方案,采用哪种呢?
谢谢指点。
作者
Author:
stephenmessi    时间: 2017-11-14 11:44
楼主您好,请问您的问题解决了吗?我也有同样的疑问,谢谢。
作者
Author:
yezhonghua    时间: 2017-11-14 14:13
文献上说的平均值是平方加和,除以3.再开方
作者
Author:
小范范1989    时间: 2017-11-14 15:04
楼主好,咱们做OLED效率计算的,有很多文章都计算了SOC。有的可能用不同的软件计算的,具体的我也不是很清楚计算的原理。只能给几个文章你看看。里面有的是平方和开根号,有的除以3.我记得还有其他的文章,除以根3.你在搜搜查查
Evaluation of spin-orbit couplings with linear-response TDDFT, TDA, and TD-DFTB

Up-Conversion Intersystem Crossing Rates in Organic Emitters for Thermally Activated Delayed Fluorescence: Impact of the Nature of Singlet vs Triplet Excited States


作者
Author:
哇哇吐    时间: 2018-1-11 14:31
楼主您好,各位算发光的老师好,我最近在用BDF (http://182.92.69.169:7226)计算旋轨耦合常数,结果文件中有如下数据不是特别明白,我也看了下小范老师在楼上贴的论文,一知半解。以下是BDF网站上对旋轨耦合结果文件的一个介绍。这是旋轨耦合矩阵元么?请问怎么样得到旋轨耦合常数呢?
“SOC matrix elements
Print selected matrix elements of [Hsoc]
  <  0  0  0 |Hso|  2  1  1 >
  mi/mj        ReHso(au)           cm^-1               ImHso(au)           cm^-1
   1  1        0.0003219734       70.6650036601        0.0009582030      210.3012602778
   1  2        0.0000000000        0.0000000000       -0.0006544171     -143.6279497862
   1  3        0.0003219734       70.6650036601       -0.0009582030     -210.3012602778
这里计算<S0|Hso|T1>分别给出其实部ReHso和虚部ImHso。因为S0只有一个分量,mi为1。T1(spin S=1)有3个分量(Ms=-1,0,1), mj编号这3个分量。

Warning: 在不同程序结果对比时需要注意:这里给出的时所谓spherical tensor,而不是cartesian tensor,即T1是T_{-1},T_{0},T_{1},不是Tx,Ty,Tz,两者之间存在酉变换。”
作者
Author:
柒月小鱼    时间: 2018-1-11 15:32
哇哇吐 发表于 2018-1-11 14:31
楼主您好,各位算发光的老师好,我最近在用BDF (http://182.92.69.169:7226)计算旋轨耦合常数,结果文件 ...

您好 请问用的BDF是自己定义的基组么 我摸了几天的BDF出现了好多问题都没法解决 希望能和您交流交流
作者
Author:
哇哇吐    时间: 2018-1-11 15:42
柒月小鱼 发表于 2018-1-11 15:32
您好 请问用的BDF是自己定义的基组么 我摸了几天的BDF出现了好多问题都没法解决 希望能和您交流交流

我也刚学没几天,我一开始也是因为基组出错,后来根据basis_library文件里的basis.map 选定的基组,然后给老师看了一下说可以,运行也没有出错。
作者
Author:
柒月小鱼    时间: 2018-1-11 15:53
哇哇吐 发表于 2018-1-11 15:42
我也刚学没几天,我一开始也是因为基组出错,后来根据basis_library文件里的basis.map 选定的基组,然后 ...

方便加个qq么 我的451303522 交流交流
作者
Author:
赵云跳槽    时间: 2018-1-13 16:52
哇哇吐 发表于 2018-1-11 14:31
楼主您好,各位算发光的老师好,我最近在用BDF (http://182.92.69.169:7226)计算旋轨耦合常数,结果文件 ...

这里不是都已经写了吗,将实部与虚部平方和再开方,三项加和后平均即可
作者
Author:
赵云跳槽    时间: 2018-1-13 16:57
哇哇吐 发表于 2018-1-11 15:42
我也刚学没几天,我一开始也是因为基组出错,后来根据basis_library文件里的basis.map 选定的基组,然后 ...

关于自定义基组,使用EMSL,导出类型选择Molcas,再略加修改

如果对BDF敢兴趣,可以加群   BDF使用交流 452912262,直接询问索兵兵老师

作者
Author:
哆啦哆啦    时间: 2018-3-3 11:42
楼主,请问一下,B and C excitation energies, moment:   0.09025588   0.08802946   0.00000005 这一行的三个数值分别代表什么
作者
Author:
sobereva    时间: 2018-3-20 05:50
哆啦哆啦 发表于 2018-3-3 11:42
楼主,请问一下,B and C excitation energies, moment:   0.09025588   0.08802946   0.00000005 这一行的 ...

计算的是<B|H_SO|C>,前两个值分别是B和C态的激发能,第三个值就是这个矩阵元数值
作者
Author:
abin    时间: 2018-6-12 15:08
哆啦哆啦 发表于 2018-3-3 11:42
楼主,请问一下,B and C excitation energies, moment:   0.09025588   0.08802946   0.00000005 这一行的 ...

原文写的是"B and C excitation energies, moment:"
所以后面的三个数字分别是, B singlet激发态激发能; C triplet激发态激发能;最后一个是moment, 在当前算例中,指的是SOCME。 Dalton论坛上有讲, 他们的输出文件, 把SOCME叫做某种moment. 当然,如无说明, 默然采用原子单位制。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3