| 请注意，在优化激发态结构时，不合适的 ADDED_MOS 关键词可能会报错，如果报错可以考虑取消该关键词 |
sobereva 发表于 2017-12-17 17:23 偶，明白了，谢谢sob老师指点。 |
小范范1989 发表于 2017-12-17 15:43 基本上是对比时候数据读取问题 根据我的经验，g09和g16，数值不相同的情况，除了默认的ultrafine、acc2e产生的影响外，也就是溶剂模型的一些计算细节发生了变化导致结果会有差异。其它方面基本都没有差异 |
灰灰不啦叽 发表于 2017-12-17 10:07 为什么数据会差这么多呢？按说不应该啊。g09和g16要是差距这么大。那还行吗？ 不知道你测试其他的对比结果如何。谢谢 |
zjxitcc 发表于 2017-12-16 20:17 感谢老师的回复，之前提取激发能的时候两种版本的log我是相同的方法找的，我再确认下自己的数据读取位置，nstates的问题我会改正的 |
|
本帖最后由 zjxitcc 于 2017-12-16 20:19 编辑 激发态优化中间每个结构，都会输出一次oscillator strengths。log文件中有大量的这个词。你是否看错了位置，看的不是最后一个结构的振子强度和激发能？ 还有，td里的子关键词是nstates，不是nstate |
| 参与人数Participants 1 | eV +1 | 收起 理由Reason |
|---|---|---|
|
| + 1 | 我很赞同 |
sobereva 发表于 2017-12-15 14:33 谢谢老师，我统一关键词再测试一下 |
灰灰不啦叽 发表于 2017-12-15 14:21 先把所有关键词统一了。 nstates也会影响结果。而且g16默认用了int=ultrafine，比较的时候应当让格点统一。在完全有可比性的情况下，g09和g16结果基本没有差异。 |
bomsaude 发表于 2017-12-15 01:25 谢谢您的提醒，我又找了一下同一个分子的不同优化结果，激发能相差特别多 |
本帖最后由 灰灰不啦叽 于 2017-12-15 14:24 编辑 sobereva 发表于 2017-12-14 22:58 感谢Sob老师的回复，之前只是怀疑g09和g16对激发态进行优化所使用的方法不一样，可能数据提供的不确切，下面的是我用相同的结构相同结构做的分子的S1态的计算（调用路径、关键词、激发能）： Gaussian16 Initial command: /opt/soft/gauss/g16/g16/l1.exe "/tmp/ram01f/gauss16.507612.mu01/Gau-18729.inp" -scrdir="/tmp/ram01f/gauss16.507612.mu01/" Entering Link 1 = /opt/soft/gauss/g16/g16/l1.exe PID= 18730. %NProcShared=12 Will use up to 12 processors via shared memory. %mem=500MW %chk=bmbpy-nidoB-B3lyp-TD-S-opt-freqNM.chk ------------------------------------------------------------------ #p td=(singlets,root=1,nstate=5) opt Frequency=SaveNM b3lyp/genecp Excitation energies and oscillator strengths: Excited State 1: Singlet-A 1.0877 eV 1139.87 nm f=0.0062 <S**2>=0.000 189 ->190 0.70692 This state for optimization and/or second-order correction. Total Energy, E(TD-HF/TD-DFT) = -2685.44800812 Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Singlet-A 2.0383 eV 608.26 nm f=0.0661 <S**2>=0.000 188 ->190 0.70350 Gaussian09 Initial command: /opt/soft/gauss/g09d01/g09/l1.exe "/tmp/ram01f/gauss09d01.446329.mu01/Gau-3624.inp" -scrdir="/tmp/ram01f/gauss09d01.446329.mu01/" Entering Link 1 = /opt/soft/gauss/g09d01/g09/l1.exe PID= 3637. %NProcShared=12 Will use up to 12 processors via shared memory. %mem=500MW %chk=bmbpy-nidoB-B3lyp-TD-S-opt.chk ------------------------------------------------- #p td(singlets,root=1,nstate=15) opt b3lyp/genecp Excitation energies and oscillator strengths: Excited State 1: Singlet-A 2.2311 eV 555.71 nm f=0.0420 <S**2>=0.000 188 ->190 0.23078 189 ->190 0.66650 This state for optimization and/or second-order correction. Total Energy, E(TD-HF/TD-KS) = -2685.42886721 Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Singlet-A 2.4794 eV 500.07 nm f=0.0220 <S**2>=0.000 188 ->190 0.66537 189 ->190 -0.23358 激发能差的特别多，不知道是不是Gaussian不同版本之间采用的优化方法导致的 |
sobereva 发表于 2017-12-14 22:58 感谢Sob老师的回复，之前只是怀疑g09和g16对激发态进行优化所使用的方法不一样，可能数据提供的不确切，下面的是我用相同的结构相同结构做的分子的S1态的计算（调用路径、关键词、激发能）： Gaussian16 Initial command: /opt/soft/gauss/g16/g16/l1.exe "/tmp/ram01f/gauss16.507612.mu01/Gau-18729.inp" -scrdir="/tmp/ram01f/gauss16.507612.mu01/" Entering Link 1 = /opt/soft/gauss/g16/g16/l1.exe PID= 18730. %NProcShared=12 Will use up to 12 processors via shared memory. %mem=500MW %chk=bmbpy-nidoB-B3lyp-TD-S-opt-freqNM.chk ------------------------------------------------------------------ #p td=(singlets,root=1,nstate=5) opt Frequency=SaveNM b3lyp/genecp Excitation energies and oscillator strengths: Excited State 1: Singlet-A 1.0877 eV 1139.87 nm f=0.0062 <S**2>=0.000 189 ->190 0.70692 This state for optimization and/or second-order correction. Total Energy, E(TD-HF/TD-DFT) = -2685.44800812 Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Singlet-A 2.0383 eV 608.26 nm f=0.0661 <S**2>=0.000 188 ->190 0.70350 Gaussian09 Initial command: /opt/soft/gauss/g09d01/g09/l1.exe "/tmp/ram01f/gauss09d01.446329.mu01/Gau-3624.inp" -scrdir="/tmp/ram01f/gauss09d01.446329.mu01/" Entering Link 1 = /opt/soft/gauss/g09d01/g09/l1.exe PID= 3637. %NProcShared=12 Will use up to 12 processors via shared memory. %mem=500MW %chk=bmbpy-nidoB-B3lyp-TD-S-opt.chk ------------------------------------------------- #p td(singlets,root=1,nstate=15) opt b3lyp/genecp Excitation energies and oscillator strengths: Excited State 1: Singlet-A 2.2311 eV 555.71 nm f=0.0420 <S**2>=0.000 188 ->190 0.23078 189 ->190 0.66650 This state for optimization and/or second-order correction. Total Energy, E(TD-HF/TD-KS) = -2685.42886721 Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Singlet-A 2.4794 eV 500.07 nm f=0.0220 <S**2>=0.000 188 ->190 0.66537 189 ->190 -0.23358 激发能差的特别多，不知道Gaussian不同版本之间采用的优化方法是否一致 |
| 从你的文件名看，应该是gaussian09和16比较，差别是挺大，不仅仅是能量。你确信是同一个分子、关键词都一样？ |
| (1)我看不出你是用谁和谁比。这和Gaussian版本没有任何关系 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-21 09:03 , Processed in 0.252132 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites