简单量子化学问题答疑专帖（旧版，已关帖）

柒月小鱼

请问老师，对于计算m->k跃迁偶极矩时，不是应该先右乘初态波函数么，为什么在您的ppt中是左乘
是否是因为跃迁偶极矩算符的特殊性质哦
麻烦老师给小白解惑了

Scienceme

卢老师 请问一下：
1. 荧光和磷光都有明确的定义，分别对应S1到S0和T1到S0的直接跃迁，但是如果基态是三重态的话（就是基态是T0），各激发态返回到T0态是什么过程？
2. 就刚才的问题，考虑旋轨耦合，S1到T1也应该有跃迁，这个绝热和垂直跃迁能量可以通过S1到S0和T1到S0的差值得到，但是这个跃迁的振子强度怎么获得呢？
3. 不同的激发态上面的激子束缚能不同，这个out文件有输出么？ 怎么可以得到呢？

sobereva

Scienceme 发表于 2019-6-11 11:21
卢老师 请问一下：
1. 荧光和磷光都有明确的定义，分别对应S1到S0和T1到S0的直接跃迁，但是如果基态是三重 ...

1 就说成“发射”就完了
2 利用旋轨耦合矩阵元，可以自己写程序通过微扰方式得到，旋轨耦合矩阵元的计算看
使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱
http://sobereva.com/462（http://bbs.keinsci.com/thread-12114-1-1.html）
3 如果是指exciton binding，这里说了怎么算
使用Multiwfn做空穴-电子分析全面考察电子激发特征
http://sobereva.com/434（http://bbs.keinsci.com/thread-10775-1-1.html）

Scienceme

多谢卢老师指点，看上去比较复杂，估计要琢么琢么~

sobereva

柒月小鱼 发表于 2019-6-11 10:43
请问老师，对于计算m->k跃迁偶极矩时，不是应该先右乘初态波函数么，为什么在您的ppt中是左乘
是否是因为 ...

ij和ji调换结果是一样的

柒月小鱼

sobereva 发表于 2019-6-11 13:07
ij和ji调换结果是一样的

谢谢老师回复

happyknighthawk

老师，您好。手边有一个天然产物二聚体，是两个monomer通过一个C-C单键连接形成的。现在想通过扫描的方法来判断该化合物是否因为该C-C单键的旋转受阻而存在阻旋异构的现象，针对这个问题，有以下几个问题想跟您请教：
（1）对这个C-C单键对应的二面角做扫描的时候应该做柔性扫描，而不是刚性扫描是吧？
（2）这个分子比较大（C43H47N4O8），是不是可以选用比较便宜的方法（比如B3LYP/6-31G(d)）来做扫描，然后对每一帧用更好的方法来做单点能计算？还是后面的单点能计算没有太大必要？
（3）一般判断常温下是否存在阻旋异构的现象的旋转能垒是多大啊？
非常感谢！

CayJ

老师，您好！请问一下：一个带有+1价的无金属元素的有机阳离子与二氧化碳的体系中，结构优化后，输入文件是这样的：

%chk=C:\Users\Administrator\Desktop\Single opt\cation-co2\energy\energy-D3-c309-co2.chk
# b3lyp/6-31+g(d,p) geom=connectivity em=gd3

Title Card Required

1 1
C                 -1.80175238   -1.72358561   -1.02832222
C                 -3.09453656   -2.12811811   -0.33388681
C                 -2.81496550   -1.72690615    1.12143602
C                 -2.08238876   -0.38466707    1.00191535
N                 -1.42203384   -0.37533245   -0.39932871
H                 -1.85841143   -1.59921468   -2.11054066
......

我想计算它们之间的弱相互作用能，于是输入文件改成：

%chk=C:\Users\Administrator\Desktop\Single opt\cation-co2\energy\energy-D3-c309-co2.chk
%nprocshared=2
# b3lyp/6-31+g(d,p)  em=gd3 counterpoise=2

Title Card Required

1 1 1 1 0 1
C(fragment=1)                 -1.80175238   -1.72358561   -1.02832222
C(fragment=1)                 -3.09453656   -2.12811811   -0.33388681
C(fragment=1)                 -2.81496550   -1.72690615    1.12143602
C(fragment=1)                 -2.08238876   -0.38466707    1.00191535
N(fragment=1)                 -1.42203384   -0.37533245   -0.39932871
H(fragment=1)                 -1.85841143   -1.59921468   -2.11054066
.......
C(fragment=2)                 -6.23285332    1.00647298    0.23449179
O(fragment=2)                 -7.31653895    1.40860177    0.35777947
O(fragment=2)                 -5.13634186    0.59901825    0.11025009

提交作业后总是报错，这是为什么？该怎么改呢？

Scienceme

CayJ 发表于 2019-6-11 22:14
老师，您好！请问一下：一个带有+1价的无金属元素的有机阳离子与二氧化碳的体系中，结构优化后，输入文件是 ...

也许你应该把 geom=connectivity 去掉~；对应的连接关系也都去掉

Scienceme

卢老师 请问一下：
1. 对弱相互作用的研究，好像分子动力学和高斯量化程序都可以，一般来说这两种方法哪个更好一些？比如研究纵横交错的两根碳纳米管之间的相互吸附力大小
2. 如果考虑的弱相互作用问题发生在溶液中，那高斯的SMD关键词（隐式溶剂模型）可以与考虑弱相互作用方法（D3校正）联合在一起使用么？ 进一步在溶液中考虑BSSE修正，是不是也可以同时使用 关键词 countpoise ？

CayJ

Scienceme 发表于 2019-6-11 23:19
也许你应该把 geom=connectivity 去掉~；对应的连接关系也都去掉

谢谢你的回复，我的想法是：geom=connectivity  应该没有什么影响，关键是Title下面的0 1 0 1 0 1出错了，一开始结构优化后此处出现的是1 1 ，第一个1应该是表明这个体系所带的电荷，第二个1是自旋多重度（看别人讨论的，具体自旋多重度是什么，怎么计算不晓得），然后我把他们改过来，改成1 1 1 1 1 1 ，再进行作用能计算就没有报错了。但为什么写1 1 1 1 1 1 ，我也是照猫画虎。

CayJ

Scienceme 发表于 2019-6-12 09:41
那就是双自由基情形

一个是阳离子，一个是CO2，为什么是双自由基呢？

CayJ

CayJ 发表于 2019-6-11 22:14
老师，您好！请问一下：一个带有+1价的无金属元素的有机阳离子与二氧化碳的体系中，结构优化后，输入文件是 ...

哪位老师能帮我解答一下疑惑，感激

sobereva

happyknighthawk 发表于 2019-6-11 17:29
老师，您好。手边有一个天然产物二聚体，是两个monomer通过一个C-C单键连接形成的。现在想通过扫描的方法来 ...

1 对
2 可以。扫描时候最好加D3
3 如果按照常温下是否能发生单分子反应的规则来判断，那就是21kcal/mol

sobereva

CayJ 发表于 2019-6-12 08:41
谢谢你的回复，我的想法是：geom=connectivity  应该没有什么影响，关键是Title下面的0 1 0 1 0 1出错了 ...

写成1 1 1 1 1 1明显是错的，哪可能每个片段净电荷加和都不等于整体电荷。不要以为不报错结果就有意义
仔细看此文
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

还搞不明白就把完整输入文件上传