计算化学公社

标题: 使用Dushin分解重组能和计算Huang-Rhys因子 [打印本页]

作者
Author: sobereva 时间: 2016-5-31 07:35
标题: 使用Dushin分解重组能和计算Huang-Rhys因子
2024-Oct-21补充：Gaussian 16计算振动电子光谱的功能直接就能给出Huang-Rhys因子，非常好用，也没有dushin程序那样有时莫名其妙无法载入文件的问题。在北京科音高级量子化学培训班（http://www.keinsci.com/workshop/KAQC_content.html）里讲授振动电子光谱的部分专门给了具体讲解和例子。

使用Dushin分解重组能和计算Huang-Rhys因子

Using Dushin program to decompose reorganization energy and calculate Huang-Rhys factor

文/Sobereva @北京科音 2016-May-31

经常有做电荷转移的人问怎么用Dushin程序、怎么算重组能、怎么把重组能分解为各个振动模式的贡献、怎么计算Huang-Rhys因子，本文就结合实例专门说说这些问题。

1 重组能(reorganization energy)的计算

重组能是基于Marcus理论计算电子转移速率的关键的量，具体分为内重组能和外重组能，前者衡量电子得/失后（或更广义来说，电子态改变后）因几何结构的弛豫导致的体系能量变化，外重组能则对应将周围环境分子重新极化所花费的能量。外重组能不好算，溶液下往往经验性地取0.2~0.6eV，晶体下则会小一个多数量级，也有些文章专门讨论怎么算，如JPCL,1,941(2010)、JACS,130,12377(2008)，本文我们只说内部重组能。

看这张示意图，有中性和带电状态两个势能面

(, 下载次数 Times of downloads: 733)

图中重组能有两部分，λ(I)和λ(II)。做Marcus理论计算时说的重组能是指二者绝对值之和。注：重组能的计算方式和定义可能令一些人混淆，可参考http://bbs.keinsci.com/thread-35003-1-1.html回帖里的讨论。

具体来说，如果charged态是带+1电荷的情况，算出来的重组能用λ+表示，也称空穴重组能λh。如果charged态是带-1电荷，算出来的重组能用λ-表示，也称电子重组能λe。

根据上面的示意图，很容易知道怎么计算λh和λe。例如计算λh，对应的示意图：

(, 下载次数 Times of downloads: 177)

需要以下步骤
(a) 优化中性状态结构
(b) 在(a)的结构下算中性状态的能量E1
(c) 在(a)的结构下算阳离子状态的能量E2
(d) 优化阳离子状态结构
(e) 在(d)的结构下算中性状态的能量E3
(f) 在(d)的结构下算阳离子状态的能量E4
重组能为
λ(I)=|E3-E1|
λ(II)=|E2-E4|
λh= λ(I)+λ(II)

类似地，我们可以算电子态改变过程的重组能。本文我们实际算一个例子，考察吡啶的基态单重态S0和第一三重态激发态T1之间的两种重组能λ(I)和λ(II)，示意图如下：

(, 下载次数 Times of downloads: 170)

可见，S0和T1态优化的结构有明显不同，前者原子是在同一个平面上的，是C2v对称性，后者就弯折了，是Cs对称性。

通常做单点计算比几何优化用的级别要高，但为了省事，此例优化和单点计算都用B3LYP/def2SVP在G09 D.01下做。实际上，也只有几何优化和单点都在同一级别时算出的重组能和后文提到的用dushin做重组能分解的结果才有可比性。算出来的图上四个点的能量：
E1=-248.1053119 Hartree
E2=-247.9505667 Hartree
E3=-248.0570726 Hartree
E4=-247.9706883 Hartree
重组能为
λ(I)=E3-E1=30.27 kcal/mol
λ(II)=E2-E4=12.63 kcal/mol
可见，从S0平衡结构垂直跃迁到T1后，结构弛豫对应的重组能为12.63 kcal/mol，这个值比从T1垂直返回S0之后结构弛豫对应的重组能30.27 kcal/mol小很多。这也暗示了，在S0-T1极小点结构之间的方向上，S0极小点处的势能面曲率大于T1极小点处的，因为以相同方式位移，S0态能量变化λ(I)比T1态能量变化λ(II)大得多。

2 重组能的分解

对于上面吡啶的例子，在S0极小点结构下做freq任务，可以得到S0下的3N-6个正则坐标。S0和T1两个极小点之间的位移ΔQ可以分解为这些正则坐标的贡献，第i个模式的贡献记为ΔQ_i。谐振势公式是E=(1/2)*k*x^2，基于这个公式，我们也可以将λ(I)分解为各个正则模式的贡献，即λ_i=(1/2)*k_i*(ΔQ_i)^2，这里k_i代表相应振动模式的力常数。当然了，谐振模型终究只是对实际势能面的近似，所以∑λ_i并不精确等于λ(I)，但能定性相符。类似地，在T1极小点结构下做freq得到3N-6个正则坐标后，也可以得到它们对ΔQ和λ(II)的贡献。

S0和T1极小点下分别得到的3N-6个正则坐标是不同的，彼此间是线性变换关系，这个关系也叫Duschinsky旋转或者称Duschinsky混合效应，可表达为Q''=J*Q'+ΔQ，这里Q'和Q''分别代表两个电子态极小点下的正则模式，J称为Duschinsky矩阵，记录了Q'与Q''间线性变换的组合系数。如果两个态极小点下振动模式完全一致，则J是单位矩阵，说明没有混合；J偏离单位矩阵越大说明混合越强烈。

提醒一下，对重组能分解时候都是基于末态极小点结构的正则模式来做分解，所以分解λ(I)要用S0极小点处的正则模式，分解λ(II)要用T1极小点处的正则模式。而不能用比如T1极小点处的正则模式分解λ(I)，因为原理上说不通。

做重组能的分解和计算Duschinsky矩阵并不复杂，也有现成的程序可以做，下面介绍的Dushin就是其一。

3 Dushin程序的使用

Dushin程序由Reimers开发的，关于程序的一些原理细节可以看JCP,115,9103(2001)，引用Dushin程序也应当引这篇文章。Dushin程序源代码包可以在这里下载：http://bbs.keinsci.com/forum.php?mod=viewthread&tid=319。Dushin程序有自带的说明文档README，但写得比较抽象，这里用人话说一下：

3.1 编译方法

Dushin是Fortran写的，机子里得有Fortran编译器才能编译，Dushin默认的是ifort编译器（gfortran应该也行，我没试）。

在Linux下编译的过程是：下载后解压，把compile里的ifort -g plot-modes.for subs.o recalc-freq.o proj0freq.o bmatred.o ddiag.o dmpower.o dmatinv.o -o ~/bin/plot-modes这一行前头加上#给注释掉，因为压缩包里并没有带plot-modes.for文件，这个本身也用不上。然后执行./compile，就会调用ifort进行编译，编译好的可执行文件会产生在用户主目录的bin目录下面，包括dushin、displace和compare-geom三个文件。然后可以在控制台直接输入dushin看是否能启动，如果不能则需要把这个bin目录添加到$PATH环境变量里。

为了便于那些不会编译或不会Linux的人使用，笔者编译了一份Windows版，下载链接： (, 下载次数 Times of downloads: 751) 。ifort编译好的Linux版可以在这里下载： (, 下载次数 Times of downloads: 352) 。

3.2 使用方法

Linux版：进入含有输入文件的目录，直接输入dushin即可运行。
Windows版：把输入文件拷入dushin程序目录下，双击dushin.exe即可运行。运行完毕后会自动关闭窗口，如果不想让窗口自动关闭就自己进DOS然后输入dushin来执行（在dushin.exe所在目录窗口下按住shift点右键选“在此处打开命令行窗口”）。

3.3 输入文件

dushin的运行需要提供原子坐标、量化程序计算出的Hessian等信息，宣称支持许多量化程序，本文只考虑Gaussian的情况。G94、G98、G03、G09的输出文件dushin都支持。

dushin.dat是dushin程序的主输入文件，dushin启动时就会载入当前目录下的dushin.dat。下面是输入文件的例子，详细的参数介绍见README：
1 2
.

2 1 'S0 freq' 'S0_freq.out'
0 1 'T1 freq' 'T1_freq.out'
第一行第一项设定使用什么坐标进行分析，一般设1就行；第二项设定怎么匹配多个输入文件里的原子顺序，一般设2就行。
第二行是指定输入文件所在目录，如果就放在了当前目录下就写.就行。
第三行是空行。
第四行开始定义输入文件。第一项是参考类型，对于第一个文件此项必须设2，之后的文件如果此项是1，代表Duschinsky矩阵是相对于第一个文件计算的，如果是0，代表是对后一个参考类型为1或2的分子计算的。第二项就设1就行了，第三项是这个文件的标签可以随意设，第四项是Gaussian freq任务输出文件名。

一般来说，写dushin.dat的时候就把上面这个例子里面的标签名、文件名改一下就行，其它不用管。

上面dushin.dat中S0_freq.out、T1_freq.out是Gaussian的freq任务的输出文件，都放在当前目录下。计算时用的关键词是#P freq b3lyp/def2svp，注意必须用#P，而且P必须大写，否则dushin无法正确识别其中的信息。算完后把.chk用formchk转换为同名的.fch文件，也放在当前目录下。（对于Linux版，后缀用默认的.log和.fchk也可以被dushin识别）

注意如果用的是Windows版Gaussian，一定要手动在输出文件最开头插入这么一行： Entering Gaussian System, Link 0=g09，否则dushin无法识别文件。

3.4 输出文件

dushin启动后会在屏幕上会输出大量信息，其中大部分都是用于调试的，用户不用管。程序会在当前目录产生一大堆输出文件，主要有用的就这两个：

dushin.out：包含主要输出信息的文件。其中Dushinsky matrix, ncs 2 in terms of ncs 1下面输出的是第二个输入文件的正则坐标(normal coordinates, ncs)是怎么由第一个输入文件的正则坐标组合而成的，包括组合系数和贡献百分比，只有贡献较大的会被输出。下面还会输出Dushinsky matrix, ncs 1 in terms of ncs 2，是描述第一个输入文件的正则坐标怎么由第二个输入文件的正则坐标组合而成的。再往下是Displacement: in terms of nc of 1 THEN of nc of 2，前几列是把位移和重组能按照第一个输入文件中的正则模式分解的结果，后几列是按照第二个输入文件的正则模式分解的结果。Q是指位移在此正则坐标上的分量，lam是每个正则模式对重组能的贡献量λ_i(cm^-1)。末尾total reorg energy (cm**-1, kcal/mol)就是∑λ_i，前两个值单位是cm^-1，分别是第一个和第二个输入文件的正则坐标对重组能贡献的加和，后两个值只不过是把单位换成了kcal/mol。

supplem.dat：包含了被计算的体系的正则坐标、频率，Duschinsky矩阵，各正则模式对位移和重组能的贡献（eV）。其实和dushin.out差不多，只不过输出格式、单位变了变。

4 实例：对吡啶S0-T1的重组能进行分解

我们在B3LYP/def2SVP级别下将电荷和自旋多重度设为0 1和0 3来分别优化吡啶的S0和T1态结构，之后用得到的结构分别做freq计算，得到Gaussian输出文件和fch文件。然后按3.3节的介绍恰当地写个dushin.dat文件。之后把dushin.dat、两个Gaussian输出文件和两个.fch文件都放在当前目录下，启动dushin，由于体系小计算量很低，马上就运行完毕。本例用到的dushin输入输出文件在这里都提供了，请自行查看： (, 下载次数 Times of downloads: 606) 。

输出文件很容易理解。比如我们这里考察对λ(I)的分解，前面提到过λ(I)是T1垂直返回S0态后因结构弛豫造成的能量变化，末态是S0平衡结构，所以要看对S0正则模式的分解结果。从输出文件中看到S0的3N-6个正则模式中对重组能贡献量最大的两个模式是427和721 cm^-1，贡献值分别高达3607.6和2456.7 cm^-1，对位移贡献分别达到4.112和2.610埃，对应的输出信息为：
8 B1 freq= 427. Q= 4.112 lam= 3607.6 A" freq= 342. Q= 0.000 lam= 0.0
11 B1 freq= 721. Q= 2.610 lam= 2456.7 A' freq= 592. Q= 2.058 lam= 1252.4

为什么这两个模式贡献会这么大，用gview看一下振动矢量就知道了，如下所示

(, 下载次数 Times of downloads: 760)

在本文前面的图片我们看到过，吡啶T1极小点结构相对于S0极小点结构，氮原子往上翘了很多，而周围的氢原子则下移了许多，导致分子就像是沿着中间稍微对折了一下，而从上面的振动矢量可见，按照这样方式振动，正是会导致分子以这种方式发生扭曲，所以这两个振动模式对重组能和位移的贡献都很大。此分子的386 cm^-1的振动模式也是偏离平面的扭曲模式，但从振动模式上明显会发现它并不会对S0-T1间结构变化产生任何贡献，所以对重组能的贡献也为0。

从输出文件中的Duschinsky矩阵中我们可以了解S0极小点和T1极小点下正则模式之间的联系。比如有这么一行输出：
8 B1 427. 0.8697 75.6 100.0 7 A' 289.
意思是S0极小点下427 cm^-1这个模式有75.6%都是由T1极小点下289 cm^-1那个模式贡献的，而100.0是代表T1的所有正则模式对S0这个模式的贡献和是100%，程序为了避免信息量太大就没输出那些贡献很小的T1正则模式。我们来看一下T1这个289 cm^-1模式的振动矢量

(, 下载次数 Times of downloads: 166)

可见，这个模式的振动矢量和前面S0的427 cm^-1的振动矢量非常相似，这也解释了为什么贡献值能高达75.6%。

输出文件末尾显示S0极小点的各个正则模式对λ(I)以及T1极小点的各个正则模式对λ(II)的贡献总和分别为32.87和18.23 kcal/mol，这和我们第一节直接按照定义算出来的值30.27和12.63 kcal/mol有一定偏差，但差异还算是可以接受范围。

自行整理一下数据格式，就可以在Origin等程序里绘制不同频率的振动模式对重组能的贡献，使贡献量一目了然，这种图经常出现在文献中。下图是S0振动模式对λ(I)的贡献图

(, 下载次数 Times of downloads: 155)

顺带一提，本例的dushin.dat里的输入文件部分是这么写的，
2 1 'S0 freq' 'S0_freq.out'
0 1 'T1 freq' 'T1_freq.out'
如果调换次序写成
2 1 'T1 freq' 'T1_freq.out'
0 1 'S0 freq' 'S0_freq.out'
实际上结果还是一样，只不过输出的顺序改了一下而已。所以写的次序无所谓。

5 计算Huang-Rhys因子

Dushin程序并不直接输出Huang-Rhys因子，自己简单算一下即可。第i个振动模式对应的Huang-Rhys因子为S_i=λ_i/(h*ν_i)，这个量是无量纲的。计算时应先把振动模式对重组能的贡献λ_i转换为以J为单位，振动频率ν_i转换为以s^-1为单位。比如S0的427 cm^-1模式的λ_i=3607.6 cm^-1，它对应的Huang-Rhys因子即为3607.6/219474.6363*2625500/6.02214179E23 / (427*2.99792458E10*6.6260696E-34)=8.45。这里分子部分先从cm^-1转为Hartree再转为J/mol再转为J。

顺带一提，利用Huang-Rhys因子，可以计算振动态跃迁对应的Franck Condon因子，结合振动能级的改变量，做Lorentzian展宽，就可以获得振动分辨的电子光谱。在JCP,120,7490(2004)当中的式8给出了谐振模型下计算相对强度值的具体公式，是同时考虑p个正则模式的情况：

(, 下载次数 Times of downloads: 150)

式中Franck-Condon积分（FCI）的平方就是FC因子，L是Laguerre多项式，m_i和n_i是分别是第i个振动模式在初态和末态时的振动量子数，S是Huang-Rhys因子，指数项是通过Boltzmann分布计算体系处在各振动态的比例来考虑温度效应。这个式子一个很大的局限性是假定没有Duschinsky混合，即J为单位矩阵，实际中也只能用在Duschinsky混合很轻微的情况，否则初态和末态下也根本没法将振动模式一一对应。上面吡咯的例子S0的721 cm^-1模式对重组能有很大贡献，但是它的Duschinsky混合很强，T1的正则模式对它贡献最大的两个是45.3%和39.0%，此时明显不能用上面的公式，对这种Duschinsky混合不可忽略的情况绘制振动分辨谱可以按此文的做法：《振动分辨的电子光谱的计算》（http://sobereva.com/223），或者自己用FCclasses程序。另外，上面那篇文章中10式是8式在只考虑一个正则模式从振动基态发生跃迁的特例，也就是Dushin程序README里提到的exp(-S) * S^n / n!，但这个用处不大，毕竟只有一个正则模式主导重组能且Duschinsky混合可忽略不计的情况极少。

作者
Author: 初初 时间: 2016-5-31 07:55
本帖最后由初初于 2016-5-31 07:57 编辑

感谢sob老师！

作者
Author: 我本是个娃娃 时间: 2016-5-31 08:19
板凳！虽然现阶段用不上，但是抢占一个好席位强势围观

作者
Author: 暖空 时间: 2016-5-31 09:56
sob老师太赞了

作者
Author: stecue 时间: 2016-5-31 11:03
赞！最近用不到，收藏之先。

作者
Author: yjmaxpayne 时间: 2016-6-2 01:29
Dushin最近我玩了不少，确实是很好用的一个程序。 Sob讲的很仔细到位，赞！

作者
Author: 小苹果 时间: 2016-6-14 06:47
再学习，谢谢群主分享！

作者
Author: aqhuangry 时间: 2016-6-19 17:26
谢谢！

作者
Author: 赵云跳槽 时间: 2016-6-22 17:57
sob老师，帖子中Huang-Rhys因子计算公式为S_i=λ_i/(h*ν_i)，但我看到林圣贤的书籍说的公式，这两个公式怎样转化？
(, 下载次数 Times of downloads: 151)
而且，以前实验室同学说也存在这样的关系：λ_i=(1/2)*v_i*(Q_i)^2；S_i=(1/2)*(Q_i)^2
使用这种简单公式计算的Huang-Rhys因子=0.5*4.112^2=8.45
这样做是否合理？

作者
Author: sobereva 时间: 2016-6-22 18:53

赵云跳槽发表于 2016-6-22 17:57
sob老师，帖子中Huang-Rhys因子计算公式为S_i=λ_i/(h*ν_i)，但我看到林圣贤的书籍说的公式，这两个公式怎 ...

不要看那个书
看这些，都是高引用的：
DOI: 10.1063/1.1687675
Organic Electronics 11 (2010) 1701–1712

作者
Author: 赵云跳槽 时间: 2016-6-22 21:24
谢谢sob老师

作者
Author: shalene 时间: 2016-6-25 18:22
非常感谢社长，momap不算好用的软件

作者
Author: ZHANGZY 时间: 2017-6-9 20:26
sob 老师，我用dushin的输出结果，按照λ=∑ћSiωi计算得到的值和 30.27和12.63 kcal/mol 相差太多，
能请教一下怎么计算出30.27和12.63 kcal/mol的吗？（程序可以直接输出，但是用λ=∑ћSiωi手动算不出这个结果来）

作者
Author: sobereva 时间: 2017-6-9 21:56

ZHANGZY 发表于 2017-6-9 20:26
sob 老师，我用dushin的输出结果，按照λ=∑ћSiωi计算得到的值和 30.27和12.63 kcal/mol 相差太多 ...

30.27和12.63 kcal/mol不是程序输出的，是第一节里直接根据单点能算的

作者
Author: ZHANGZY 时间: 2017-8-4 19:17
sobereva老师，我注意到，您在这个例子中算T1的frequency，不是按照激发态频率分析计算的，而是利用T1态的优化结构，按照T1态的优化结构对应的基态频率计算。这样做是忽略了电子结构的影响，只考虑原子核结构对频率的影响，这样做的依据是什么？这样做的误差是可以接受的吗？

作者
Author: sobereva 时间: 2017-8-4 20:52

ZHANGZY 发表于 2017-8-4 19:17
sobereva老师，我注意到，您在这个例子中算T1的frequency，不是按照激发态频率分析计算的，而是利用T1态的 ...

如果你是指为什么不用0 1结合TD=triplet来优化和算频率而用0 3结合UKS来做，是因为后者便宜得多，而且原理上结果更好，因为轨道是对三重态优化的，而前者只是基于单重态轨道产生的三重态信息

作者
Author: ZHANGZY 时间: 2018-4-18 10:45
sobereva老师，您好。ΔQi 的值能大道10-33吗？

In dushin, the vibrational modes C:\123.png

ΔQi represents the displacement along the i-th normal mode coordinate between the equilibrium positions of the charged state and the neutral state

作者
Author: ZHANGZY 时间: 2018-4-18 20:29
sobereva老师，您好。
我坐tdopt+freq 写入opt=z-matrix。发现计算最后，振动相关内坐标定义发生改变与最初的定义相差很远，导致我利用dushin计算重组能时，没法计算（supplem.geom显示内坐标不匹配）。
请问有什么好的办法固定opt和freq的内坐标形式？

作者
Author: 胖嘟嘟123 时间: 2019-4-30 19:26
老师，您好，我按照您说的方法进行计算，在计算过渡金属的时候都没有问题，但是在计算有机分子的时候，输出数据都是0，这是什么原因导致的呢

作者
Author: sobereva 时间: 2019-5-1 06:23

胖嘟嘟123 发表于 2019-4-30 19:26
老师，您好，我按照您说的方法进行计算，在计算过渡金属的时候都没有问题，但是在计算有机分子的时候，输出 ...

没细节没法说，自行检查输出信息搞清楚，尝试变化不同计算条件
这和是过渡金属体系还是有机体系没有直接关系

作者
Author: 凉白开 时间: 2019-7-23 11:37
卢老师，又没有可能算出来的重组能是0呢，假如算的是S1-T1之间的重组能。

作者
Author: 凉白开 时间: 2019-7-23 11:39

凉白开发表于 2019-7-23 11:37
卢老师，又没有可能算出来的重组能是0呢，假如算的是S1-T1之间的重组能。

补充照片

作者
Author: sobereva 时间: 2019-7-23 18:59

凉白开发表于 2019-7-23 11:37
卢老师，又没有可能算出来的重组能是0呢，假如算的是S1-T1之间的重组能。

在有人回复前，补充图片的时候请直接编辑原先帖子，而非单独发新帖子。

原理上不可能，毕竟S1和T1极小点不同

作者
Author: 凉白开 时间: 2019-7-23 21:24

sobereva 发表于 2019-7-23 18:59
在有人回复前，补充图片的时候请直接编辑原先帖子，而非单独发新帖子。

原理上不可能，毕竟S1和T1极小 ...

谢谢卢老师的回答，我以后注意，这波操作有点菜鸟，哈哈

作者
Author: 凉白开 时间: 2019-7-23 21:26

sobereva 发表于 2019-7-23 18:59
在有人回复前，补充图片的时候请直接编辑原先帖子，而非单独发新帖子。

原理上不可能，毕竟S1和T1极小 ...

所以，这有可能是个计算操作问题，是吗，卢老师

作者
Author: sobereva 时间: 2019-7-24 22:23

凉白开发表于 2019-7-23 21:26
所以，这有可能是个计算操作问题，是吗，卢老师

对

作者
Author: 凉白开 时间: 2019-7-25 10:26

sobereva 发表于 2019-7-24 22:23
对

好的，谢谢卢老师！

作者
Author: ZHANGZY 时间: 2019-11-15 07:08
本帖最后由 ZHANGZY 于 2019-11-15 08:54 编辑

sobereva老师,
QM/mm FREQ算出的fch 和 log ，在应用dushin的时候，说是 “too many atoms, increase M”。
把m调高到10000
重新编译DUSHIN，编译通过
结果计算时，segment 错误，该如何处理？

作者
Author: sobereva 时间: 2019-11-16 13:12

ZHANGZY 发表于 2019-11-15 07:08
sobereva老师,
QM/mm FREQ算出的fch 和 log ，在应用dushin的时候，说是 “too many atoms, increase M” ...

我没试过，可能其它有些地方由于是固定长度的数组，也需要相应地修改
得调试代码

作者
Author: yvette 时间: 2020-4-3 00:13
本帖最后由 yvette 于 2020-4-3 11:43 编辑

sob老师，dushin提取基态S0和激发态S1的重组能使，DOS界面运行dushin.exe报错：int defn changed for coord 50
S1的优化构型是基于S0的优化构型做的TD优化，可还是出现这种错误，该怎么修改呢？谢谢

作者
Author: lianghuashou 时间: 2021-8-12 18:02

凉白开发表于 2019-7-25 10:26
好的，谢谢卢老师！

我用dushin计算重组能的时候也遇到这个问题了，不知道你是怎么解决的

作者
Author: 凉白开 时间: 2022-1-27 20:51

lianghuashou 发表于 2021-8-12 18:02
我用dushin计算重组能的时候也遇到这个问题了，不知道你是怎么解决的

下载最新的版本即可，发邮件询问作者

作者
Author: 天天121 时间: 2022-3-29 20:06
Sob老师你真的是太及时了

作者
Author: Novice 时间: 2022-8-19 18:09
(, 下载次数 Times of downloads: 54)
请问社长，如上图中的“Distribution of bond stretch, bond angle and dihedral angle in total recombination energy”，这个信息从哪里获取呢

作者
Author: sobereva 时间: 2022-8-20 05:34

Novice 发表于 2022-8-19 18:09
请问社长，如上图中的“Distribution of bond stretch, bond angle and dihedral angle in total recombi ...

用下文的功能可以得到各个振动模式中不同类型项的贡献，再结合各个振动模式对重组能的贡献，就能得到不同类型项对重组能的总贡献。需要写点程序算
Gaussian中分析振动模式成份的关键词freq=intmodes
http://sobereva.com/106

作者
Author: 小虎兔 时间: 2022-11-8 21:47

sobereva 发表于 2016-6-22 18:53
不要看那个书
看这些，都是高引用的：
DOI: 10.1063/1.1687675

sob老师，教程中的Huang-Rhys因子计算公式是S_i=λ_i/(h*ν_i)
Organic Electronics 11 (2010) 1701–1712 这篇论文中公式3是S_i=λ_i/(ћ*w_i)
相差一个2π，计算结果是一样的吗？

作者
Author: vverawei 时间: 2022-11-19 22:49

sobereva 发表于 2017-6-9 21:56
30.27和12.63 kcal/mol不是程序输出的，是第一节里直接根据单点能算的

sob老师您好，我看到您上面说优化T1态的结构设置电荷和自旋多重度分别为0和3，但是T1不是带电态吗，假如说带正电的话不应该设置为1，2吗？

作者
Author: sobereva 时间: 2022-11-19 22:55

vverawei 发表于 2022-11-19 22:49
sob老师您好，我看到您上面说优化T1态的结构设置电荷和自旋多重度分别为0和3，但是T1不是带电态吗，假如 ...

T1怎么可能带电，又不是电离态

作者
Author: vverawei 时间: 2022-11-19 23:06

sobereva 发表于 2022-11-19 22:55
T1怎么可能带电，又不是电离态

不好意思老师，我没理解清楚，所以说是对于同一个结构分别设置0，1和0，3就是计算了它的单重态和第三激发态吗？

作者
Author: sobereva 时间: 2022-11-19 23:10

vverawei 发表于 2022-11-19 23:06
不好意思老师，我没理解清楚，所以说是对于同一个结构分别设置0，1和0，3就是计算了它的单重态和第三激发 ...

T1那叫能量最低的三重态
好好把此文看了
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

作者
Author: vverawei 时间: 2022-11-19 23:12

sobereva 发表于 2022-11-19 22:55
T1怎么可能带电，又不是电离态

不好意思老师，我还有一个问题，就是您之前教学帖子说的单重态的E1和E3，第一三重激发态的E2和E4。E1是设置电荷和多重度为0，1进行结构优化得到，E2是在E1结构的基础上改成0，3进行单点能计算，E4是在E2的结构基础上保持电荷和多重度为0，3进行结构优化计算，E3是在E4的结构基础上改成0，1进行单点能计算得到的呀？

作者
Author: vverawei 时间: 2022-11-19 23:32

sobereva 发表于 2022-11-19 22:55
T1怎么可能带电，又不是电离态

老师，请问S0态和T1态该如何确定呢？还麻烦老师多指教几句，拜托了

作者
Author: wzkchem5 时间: 2022-11-20 03:22

vverawei 发表于 2022-11-19 16:32
老师，请问S0态和T1态该如何确定呢？还麻烦老师多指教几句，拜托了

什么叫如何确定？把你要问的问题描述清楚，确保别人可以确切无误地理解，没有歧义

作者
Author: vverawei 时间: 2022-11-20 09:50

wzkchem5 发表于 2022-11-20 03:22
什么叫如何确定？把你要问的问题描述清楚，确保别人可以确切无误地理解，没有歧义

老师好，我想问的就是如何计算一个分子的单重态和第一激发三重态的结构优化，直接对该分子分别改变电荷和多重度为0，1和0，3就可以吗？sob老师之前在dushin 软件的教学帖子说的单重态的E1和E3，第一三重激发态的E2和E4。E1是设置电荷和多重度为0，1进行结构优化得到，E2是在E1结构的基础上改成0，3进行单点能计算，E4是在E2的结构基础上保持电荷和多重度为0，3进行结构优化计算，E3是在E4的结构基础上改成0，1进行单点能计算得到的吗？

作者
Author: I10140317 时间: 2023-11-3 18:29
重组能算法跟势能面相对位置是否相关呢？看论坛里面的帖子，如果是四点法的话好像适用于左图，如果是两点法基本上适用于右图。这两张图的势能面两个态的相对位置会发生反转，在左图中，E4>E3,按四点法算两部分重组能加起来没问题，但是右图中，E4<E3，这时候重组能的两部分应该是怎么标注呢？为什么这个时候就不能用lamda=(E3-E1)+(E2-E4)了

作者
Author: sobereva 时间: 2023-11-4 05:31

I10140317 发表于 2023-11-3 18:29
重组能算法跟势能面相对位置是否相关呢？看论坛里面的帖子，如果是四点法的话好像适用于左图，如果是两点法 ...

右图情况依然是(E2-E4) + (E3-E1)

作者
Author: Abhijit@123 时间: 2024-4-1 23:27
I got the following error after running the dushin.exe from windows- 'Microsoft Windows [Version 10.0.19045.4170]
(c) Microsoft Corporation. All rights reserved.

C:\Users\HP>cd C:\Dushin_win32_bin\Dushin_win

C:\Dushin_win32_bin\Dushin_win>dushin
Enter ifcurve (0=recti, 1=orth curvi, 2=apprx curvi, -1 = skip
And atom ordering flag, 0= never changes, 1= use the first order thenceforth,
  2= deduce it (does read on fail), 3= always read it
Reading file:
./S0_freq.out

  Entering Gaussian System, Link 0=g09

Gaussian format, type: g94, g98, g03= F F F
  Normal termination of Gaussian 09 at Sun Mar 31 02:12:06 2024.
cant determine energy type

C:\Dushin_win32_bin\Dushin_win>'

Can u pls help?

作者
Author: sobereva 时间: 2024-4-3 02:19

Abhijit@123 发表于 2024-4-1 23:27
I got the following error after running the dushin.exe from windows- 'Microsoft Windows [Version 10. ...

The old version of Dushin has poor compatibility with certain output file of Gaussian. There are newer and modified versions of Dushin in this forum, which may work, but I don't have relevant experience.

作者
Author: Abhijit 时间: 2024-4-21 00:35
I am getting the following error when running DUSHIN- "  UNMAPPED #2 atom          1          6
TRYING LOCATION ONLY--- READ IN ORDER
forrtl: severe (24): end-of-file during read, unit 5, file C:\Dushin_win32_bin\Dushin_win\dushin.dat
Image             PC       Routine          Line       Source
dushin.exe       004C3330  Unknown             Unknown  Unknown
dushin.exe       0048A516  Unknown             Unknown  Unknown
dushin.exe       00478722  Unknown             Unknown  Unknown
dushin.exe       00477A64  Unknown             Unknown  Unknown
dushin.exe       00458360  Unknown             Unknown  Unknown
dushin.exe       0045627A  Unknown             Unknown  Unknown
dushin.exe       0043B495  Unknown             Unknown  Unknown
dushin.exe       00401F33  Unknown             Unknown  Unknown
dushin.exe       004C45F3  Unknown             Unknown  Unknown
dushin.exe       004ADF6E  Unknown             Unknown  Unknown
KERNEL32.DLL    7775FCC9  Unknown             Unknown  Unknown
ntdll.dll       77A47C5E  Unknown             Unknown  Unknown
ntdll.dll       77A47C2E  Unknown             Unknown  Unknown"

Can anyone help me?

作者
Author: hu520130 时间: 2024-5-10 01:25
sob老师，您好，我照着您写的吡啶的教程，算了一遍，最后使用Dushin输出重组能这些信息的时候，得到的out文件里，直输出到“The point group of the molecule is C2V
”这一行，后面就没有了，是什么原因呢？我拿您的吡啶文件，能输出全部信息。为了判断是我的fch文件还是out文件有问题，我挨个和您文件里的fch和out文件替换了，发现是我的out有问题，我的fch文件搭配您的out文件，也能输出全部信息。这是什么原因导致的？盼复

作者
Author: heidou 时间: 2024-7-10 18:57
请问一下各位老师，如果window系统没有Fortran编译器，那么就会出现dushin闪退--打不开的现象呢

作者
Author: sobereva 时间: 2024-7-10 23:19

heidou 发表于 2024-7-10 18:57
请问一下各位老师，如果window系统没有Fortran编译器，那么就会出现dushin闪退--打不开的现象呢

命令行下运行

作者
Author: heidou 时间: 2024-7-11 11:07
本帖最后由 heidou 于 2024-7-11 11:09 编辑

sobereva 发表于 2024-7-10 23:19
命令行下运行

sob老师，非常感谢您！
我在网上查询了命令行下运行的方法：以下是我的操作
1.将dushin的文件移入program Files--C 盘（之前是在D盘直接复制到C 盘的）
2.将在命令提示符下输入start c:\program files\dushin\dushin.exe
结果：显示系统找不到文件C:\Users\manli\Desktop\图片1

作者
Author: sobereva 时间: 2024-7-11 15:40

heidou 发表于 2024-7-11 11:07
sob老师，非常感谢您！
我在网上查询了命令行下运行的方法：以下是我的操作
1.将dushin的文件移入progr ...

别放在有空格的路径下，自找麻烦
看https://www.youtube.com/watch?v= ... xT4JbMu&index=1这套视频了解点命令行运行的基本常识

作者
Author: heidou 时间: 2024-7-11 16:03

sobereva 发表于 2024-7-11 15:40
别放在有空格的路径下，自找麻烦
看https://www.youtube.com/watch?v=MBBWVgE0ewk&list=PL6gx4Cwl9DGDV6 ...

非常感谢sob老师

作者
Author: heidou 时间: 2024-7-21 18:25
各位老师，请问我在命令提示符下运行dushin时，输入 start dushin.exe为什么依旧出现闪退的情况呢

作者
Author: sobereva 时间: 2024-7-23 03:39

heidou 发表于 2024-7-21 18:25
各位老师，请问我在命令提示符下运行dushin时，输入 start dushin.exe为什么依旧出现闪退的情况呢

仔细看3.2节

另外，写start完全多余

作者
Author: heidou 时间: 2024-7-24 09:11

sobereva 发表于 2024-7-23 03:39
仔细看3.2节

另外，写start完全多余

收到，谢谢sob老师

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