DALTON 老问题求助

ZHANGZY

我也遇到了“小范范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理解为是三重态的三个简并态？
以上两种方案，采用哪种呢？
谢谢指点。

stephenmessi

楼主您好，请问您的问题解决了吗？我也有同样的疑问，谢谢。

yezhonghua

文献上说的平均值是平方加和，除以3.再开方

小范范1989

楼主好，咱们做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

哇哇吐

楼主您好，各位算发光的老师好，我最近在用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，两者之间存在酉变换。”

柒月小鱼

哇哇吐 发表于 2018-1-11 14:31
楼主您好，各位算发光的老师好，我最近在用BDF （http://182.92.69.169:7226）计算旋轨耦合常数，结果文件 ...

您好 请问用的BDF是自己定义的基组么 我摸了几天的BDF出现了好多问题都没法解决 希望能和您交流交流

哇哇吐

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

我也刚学没几天，我一开始也是因为基组出错，后来根据basis_library文件里的basis.map 选定的基组，然后给老师看了一下说可以，运行也没有出错。

柒月小鱼

哇哇吐 发表于 2018-1-11 15:42
我也刚学没几天，我一开始也是因为基组出错，后来根据basis_library文件里的basis.map 选定的基组，然后 ...

方便加个qq么 我的451303522 交流交流

赵云跳槽

哇哇吐 发表于 2018-1-11 14:31
楼主您好，各位算发光的老师好，我最近在用BDF （http://182.92.69.169:7226）计算旋轨耦合常数，结果文件 ...

这里不是都已经写了吗，将实部与虚部平方和再开方，三项加和后平均即可

赵云跳槽

哇哇吐 发表于 2018-1-11 15:42
我也刚学没几天，我一开始也是因为基组出错，后来根据basis_library文件里的basis.map 选定的基组，然后 ...

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

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

哆啦哆啦

楼主，请问一下，B and C excitation energies, moment:   0.09025588   0.08802946   0.00000005 这一行的三个数值分别代表什么

sobereva

哆啦哆啦 发表于 2018-3-3 11:42
楼主，请问一下，B and C excitation energies, moment:   0.09025588   0.08802946   0.00000005 这一行的 ...

计算的是<B|H_SO|C>，前两个值分别是B和C态的激发能，第三个值就是这个矩阵元数值

abin

哆啦哆啦 发表于 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. 当然，如无说明， 默然采用原子单位制。