基态和激发态没有使用同一个方法优化和计算频率,后续的计算结果还可以使用吗

xm.liu · 发表于 Post on 2023-8-16 15:02:58

本帖最后由 xm.liu 于 2023-8-16 15:07 编辑

老师们好，用gaussian16优化了一系列荧光分子的基态和激发态S1,基态优化和频率计算使用的是b3lyp/6-311G(d,p),激发态优化和计算频率使用的是cam-b3lyp/6-311G(d,p),后续又使用这两个freq的计算结果计算了RMSD,重组能和HR因子及非辐射速率Kic,想请教大家这个结果可以使用吗，做这几个分子的横向比较可以吗，如果不行，现在肯定更容易的是把基态也用cam-b3lyp优化，但是这样就会造成现有的S1优化结果的起始点不再是S0,另一条路，如果重新优化激发态时间也要很久（好几个分子，结构也挺大），如果结果不能用想求助老师们有什么方法可以帮我解决这问题，十分感激

wzkchem5 · 发表于 Post on 前天 11:41

njtech_wk 发表于 2026-1-20 17:48
谢谢老师，我再多问一下，下面orca的关键词输入，基于基态优化的结构，来优化S1算发射，并且同时计算5个 ...

一般不要在同一个输入文件里写opt和dosoc true。那样会在带soc校正的情况下优化激发态，而且我不知道会优化第几个态，未必优化的还是S1，也可能程序直接报错。应当先用一个输入文件优化S1态，再另写一个输入文件，在收敛的S1结构下算SOC

njtech_wk · 发表于 Post on 3 day ago

wzkchem5 发表于 2026-1-20 16:04
可以，但建议用ORCA，以高斯的基态结构为初猜，重新做基态结构优化。否则发射虽然可以算，但发射结果和吸 ...

谢谢老师，我再多问一下，下面orca的关键词输入，基于基态优化的结构，来优化S1算发射，并且同时计算5个单重态和5个三重态以及旋轨耦合耦合矩阵元这样输入会默认优化S1的结构吗？
! OPT LibXC(B3LYP) def2-SVP def2/J def2-SVP/C RIJCOSX tightSCF noautostart miniprint
%maxcore 100000
%pal nprocs  36 end

%tddft
nroots    5
DoSOC    true
PrintLevel  3
TPrint    1E-8
TDA       false
triplets true
end

wzkchem5 · 发表于 Post on 3 day ago

njtech_wk 发表于 2026-1-19 15:37
请问一下老师我能用Gaussian优化基态结构然后转去ORCA做激发态的计算吗（做了发射和吸收）用Gaussian做 ...

可以，但建议用ORCA，以高斯的基态结构为初猜，重新做基态结构优化。否则发射虽然可以算，但发射结果和吸收结果不可比。
同一个软件，LibXC的泛函和非LibXC的同一个泛函结果基本可比，问题不大。如果怀疑不可比，自行从你的体系里挑一个分子，用LibXC的泛函和非LibXC的同一个泛函算单点能，看看差多少即可

njtech_wk · 发表于 Post on 3 day ago

本帖最后由 njtech_wk 于 2026-1-20 11:36 编辑

cokie 发表于 2026-1-19 21:06
没看懂你的描述，按说你只用优化好的基态几何结构是不能得到发射的，这需要优化好的激发态几何结构。
如 ...

你好我是用Gaussian只优化了基态结构然后用ORCA优化S1结构（模拟荧光发射，单重态和三重态各算5个输入关键词如下所示），激发的话我就把OPT去掉，可以帮忙看一下这个输入对不对吗

！OPT LibXC(B3LYP) def2-SVP def2/J def2-SVP/C RIJCOSX tightSCF noautostart miniprint
%maxcore 100000
%pal nprocs  36 end

%tddft
nroots    5
DoSOC    true
PrintLevel  3
TPrint    1E-8
TDA       false
triplets true
end

cokie · 发表于 Post on 4 day ago

njtech_wk 发表于 2026-1-19 15:37
请问一下老师我能用Gaussian优化基态结构然后转去ORCA做激发态的计算吗（做了发射和吸收）用Gaussian做 ...

没看懂你的描述，按说你只用优化好的基态几何结构是不能得到发射的，这需要优化好的激发态几何结构。
如果你描述有误，实际指的是用Gaussian做基态、激发态的几何优化，ORCA做单点的话，那没有任何影响。ORCA有很多很好用的高精度的适合算激发态单点的方法，又快又好用。只要你确保横向做对比的几个分子用的方法一致即可。

njtech_wk · 发表于 Post on 4 day ago

wzkchem5 发表于 2023-8-16 15:43
不能，尤其是重组能、HR因子、RMSD的基态、激发态理论级别必须统一，否则就会发现即使激发过程对结构毫无影 ...

请问一下老师我能用Gaussian优化基态结构然后转去ORCA做激发态的计算吗（做了发射和吸收）用Gaussian做发射太慢了基态和激发态用的都是B3LYP/def2-SVP 不过ORCA6.1.1版本里面只能用LibXC版本的B3LYP 我主要研究激发态的能级和振子强度会不会有影响

wzkchem5 · 发表于 Post on 2023-9-25 16:01:25

身在盘丝洞发表于 2023-9-24 16:13
谢谢老师指导，如果考虑这些因素去计算该如何做呢？有没有相关的博文。
另外老师无论什么溶剂我用cam-不 ...

参考https://mattermodeling.stackexch ... e-scf-a/10923#10923

身在盘丝洞 · 发表于 Post on 2023-9-24 23:13:12

wzkchem5 发表于 2023-9-23 19:10
会不会是低极性（或高极性）溶剂里分子的优势构象或聚集行为不同，导致发光的差异？

谢谢老师指导，如果考虑这些因素去计算该如何做呢？有没有相关的博文。
另外老师无论什么溶剂我用cam-不b3lyp都是LE，pbe1pbe都是CT这是怎么回事呢，会不会选一个在这二者之间的泛函能解决？（例如HF成分在二者之间的或者其他因素）

wzkchem5 · 发表于 Post on 2023-9-23 19:10:00

身在盘丝洞发表于 2023-9-23 10:16
谢谢老师指导，我的这个只有弱相互作用，可是溶剂的影响一直体现不出来

会不会是低极性（或高极性）溶剂里分子的优势构象或聚集行为不同，导致发光的差异？

身在盘丝洞 · 发表于 Post on 2023-9-23 17:16:32

wzkchem5 发表于 2023-9-22 23:14
做科研不是这么个做法，如果没有彻底理解原理，导致走了弯路，更加得不偿失。
LM方程假设溶剂可以看作均 ...

谢谢老师指导，我的这个只有弱相互作用，可是溶剂的影响一直体现不出来

wzkchem5 · 发表于 Post on 2023-9-22 23:14:05

身在盘丝洞发表于 2023-9-22 15:42
谢谢老师，我是做TADF的，一般都是CT或者LE，这个分子特殊了点就想着模拟看看，要是跟着推导一遍对我来讲 ...

做科研不是这么个做法，如果没有彻底理解原理，导致走了弯路，更加得不偿失。
LM方程假设溶剂可以看作均匀电介质，如果溶剂和溶质之间有比较强的相互作用（比如比较强的氢键），就会违反LM方程的基本假设，有可能就会导致偏离线性

身在盘丝洞 · 发表于 Post on 2023-9-22 22:42:09

wzkchem5 发表于 2023-9-22 22:21
找LM方程的原始文献，看文献是怎么推导的，逐一分析每步推导用到了什么假设，然后看这些假设你的体系是否 ...

谢谢老师，我是做TADF的，一般都是CT或者LE，这个分子特殊了点就想着模拟看看，要是跟着推导一遍对我来讲有点得不偿失，还是跟着顺一遍并不困难

wzkchem5 · 发表于 Post on 2023-9-22 22:21:54

身在盘丝洞发表于 2023-9-22 15:15
老师我的数据和图片差不多，但是不了解LM方程的使用条件，搜了一下也没有太多介绍的

找LM方程的原始文献，看文献是怎么推导的，逐一分析每步推导用到了什么假设，然后看这些假设你的体系是否满足。
LM方程这种东西，我觉得还是得自己推一遍才能理解

身在盘丝洞 · 发表于 Post on 2023-9-22 22:15:15

老师我的数据和图片差不多，但是不了解LM方程的使用条件，搜了一下也没有太多介绍的

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 基态和激发态没有使用同一个方法优化和计算频率,后续的计算结果还可以使用吗

浏览过的版块