计算化学公社

标题: 固相下优化出的S1，T1态能量低于S0态的问题 [打印本页]

作者
Author: 阳光总在风雨后 时间: 2017-11-3 08:52
标题: 固相下优化出的S1，T1态能量低于S0态的问题
S0态在这里 attach://11244.gjf

作者
Author: 阳光总在风雨后 时间: 2017-11-3 08:58
标题: 构型优化
不好意思，刚刚那个操作失误！请问：固相下，优化出的S1，T1态能量低于S0态的，这是问题出在哪里？下面贴出S0，S1，T1态输入文件

作者
Author: ggdh 时间: 2017-11-3 10:41
用普通的调控就行了。还非要用pcm调控。。。pcm调控算IP EA还行。算激发态容易低估。

作者
Author: 阳光总在风雨后 时间: 2017-11-3 12:59

ggdh 发表于 2017-11-3 10:41
用普通的调控就行了。还非要用pcm调控。。。pcm调控算IP EA还行。算激发态容易低估。

普通调控是什么意思？怎么改哪？

作者
Author: ggdh 时间: 2017-11-3 13:06

阳光总在风雨后发表于 2017-11-3 12:59
普通调控是什么意思？怎么改哪？

你调控的时候不要加scrf

作者
Author: 阳光总在风雨后 时间: 2017-11-3 15:15

ggdh 发表于 2017-11-3 13:06
你调控的时候不要加scrf

恕我愚钝不是特别明白，不加scrf的话，那么请问固相下介电常数怎么体现呐？

作者
Author: 小范范1989 时间: 2017-11-3 15:26

阳光总在风雨后发表于 2017-11-3 15:15
恕我愚钝不是特别明白，不加scrf的话，那么请问固相下介电常数怎么体现呐？

看一下钟老师的大作：
Impact of Dielectric Constant on the Singlet-Triplet Gap in Thermally Activated Delayed Fluorescence (TADF) Materials
还有这个Ionization Energies, Electron Affinities, and Polarization Energies of Organic Molecular Crystals: Quantitative Estimations from a Polarizable Continuum Model (PCM)-Tuned Range-Separated Density Functional Approach
里面supporting里面有文件例子。

作者
Author: 小范范1989 时间: 2017-11-3 15:33

ggdh 发表于 2017-11-3 13:06
你调控的时候不要加scrf

钟老师好，正好咨询一下钟老师几个问题，因为最近也在借助钟老师的PCM调控的方法计算solid phase。
1：我看钟老师的方法在计算TADF的能查非常好，光谱如何呢？比如吸收发射，以及激发态等。cam-B3lyp对激发态描述如何呢？我看钟老师好像在计算能差的时候，用LC-wpbe等做了单点矫正。
2：PCM的介电常数，是通过gromacs模拟盒子，得到的介电常数。还是通过计算单分子体积等，然后通过一个公式转换的到的？我看孙老师主页有几个汉语文章，里面是计算单分子方法得到。但是我测了一下我想计算的几个TADF分子，gromacs得到的和单分子得到的这个介电常数还是有不小的差距的。
谢谢钟老师指点。顺带问一下钟老师。我是想最后计算TADF发光效率等。对激发态是采用cam-B3LYP描述，然后在算能查的时候再用LC-wpbe或者wb97xd还是怎么样的方法？
望钟老师指点一二，谢谢钟老师。

作者
Author: ggdh 时间: 2017-11-3 17:54

阳光总在风雨后发表于 2017-11-3 15:15
恕我愚钝不是特别明白，不加scrf的话，那么请问固相下介电常数怎么体现呐？

调控的时候不加，调控完了算td的时候加

作者
Author: 阳光总在风雨后 时间: 2017-11-4 10:13

ggdh 发表于 2017-11-3 17:54
调控的时候不加，调控完了算td的时候加

不好意思，再次打扰您！之前您说的去掉scrf,但是算出的SI态能量依旧比S0的低。这个是调控过w的，然后进行构型优化。在这里优化构型是因为要算S1和T1态的绝热能差。
S0态输入文件：
%nproc=12
%mem=10GB
#p opt  freq=noraman  wB97XD/cc-PVDZ  IOP(3/107=0042000000,3/108=0042000000)

Title Card Required

0 1
C                -7.51448700 1.26844700 -1.32156000
C                -6.64039200 2.23463300 -1.83789300
C                -5.27544500 2.17720500 -1.55246200
C                -4.75339900 1.15652400 -0.73530700
C                -5.66075200 0.22171800 -0.22064100
.....
H                4.96998300 4.02079200 2.71876200
H                3.47045600 5.48113900 1.37203300
H                2.28589900 4.54557100 -0.60709500
H                2.62746900 2.16996300 -1.25374900

eps=3

S1输入文件：
%nproc=12
%mem=10GB
#p opt  wB97XD/cc-PVDZ  td(nstates=3,root=1) IOP(3/107=0042000000,3/108=0042000000)

Title Card Required

0 1
C                -7.53430600 1.04321600 -1.34576400
C                -6.69490700 2.00446500 -1.94316900
C                -5.33001400 2.01403000 -1.68507000
C                -4.72908400 1.06936600 -0.81316300
C                -5.60922300 0.13855100 -0.19853800
。。。
H                2.26249400 4.47242400 -0.67733600
H                2.64166500 2.09013900 -1.28965200

eps=3
给您带来困扰实在不好意思！望赐教

作者
Author: ggdh 时间: 2017-11-5 11:00

阳光总在风雨后发表于 2017-11-4 10:13
不好意思，再次打扰您！之前您说的去掉scrf,但是算出的SI态能量依旧比S0的低。这个是调控过w的，然后进行 ...

我的意思是调控的时候不要加scrf，但是你的w显然是在调控的时候加了scrf后获得的。
因此你要在调控的时候去掉scrf重新调控。
然后你在s1 输入文件中把 scrf加回去。
你好像是弄反了。你的w没有变，反而把s1里面的scrf去掉了。。

作者
Author: ggdh 时间: 2017-11-5 11:00

小范范1989 发表于 2017-11-3 15:26
看一下钟老师的大作：
Impact of Dielectric Constant on the Singlet-Triplet Gap in Thermally Activa ...

这是孙海涛的大作，我只是写了个调控脚本。。。

作者
Author: ggdh 时间: 2017-11-5 11:17

小范范1989 发表于 2017-11-3 15:33
钟老师好，正好咨询一下钟老师几个问题，因为最近也在借助钟老师的PCM调控的方法计算solid phase。
1： ...

1：我看钟老师的方法在计算TADF的能查非常好，光谱如何呢，比如吸收发射，以及激发态等？
答：pcm调控对于CT态容易低估激发态能级，有全局HF的情况稍微好点，真空中调控的还不错。和b3lyp 以及pbe0达到相同水平，但体系适用性要更广。。毕竟是system-depend HF%

2：cam-B3lyp对激发态描述如何呢？我看钟老师好像在计算能差的时候，用LC-wpbe等做了单点矫正。
cam-b3lyp 对于激发态一般高估，但是优化激发态时结果反而比较可靠，尤其是对于CT强的，b3lyp中CT态权重很高，往往为了形成更强的CT（比如TICT)而失去节操，而且用pbe0 或者b3lyp 优化激发态容易被“fashion”的审稿人bb。此外调控的wb97x, lcblyp,lcpbe都可以用。 lc-wpbe不能优化激发态(gaussian09 不支持）

2：PCM的介电常数，是通过gromacs模拟盒子，得到的介电常数。还是通过计算单分子体积等，然后通过一个公式转换的到的？我看孙老师主页有几个汉语文章，里面是计算单分子方法得到。但是我测了一下我想计算的几个TADF分子，gromacs得到的和单分子得到的这个介电常数还是有不小的差距的。
gromas模拟更好，或者castep计算也行，或者是有实验值。不知道你说的不小的差距具体是多少。。。

谢谢钟老师指点。顺带问一下钟老师。我是想最后计算TADF发光效率等。对激发态是采用cam-B3LYP描述，然后在算能查的时候再用LC-wpbe或者wb97xd还是怎么样的方法？
我推荐wb97xd，lc-wpbe相对容易低估ΔST。

作者
Author: 小范范1989 时间: 2017-11-5 15:15

ggdh 发表于 2017-11-5 11:17
1：我看钟老师的方法在计算TADF的能查非常好，光谱如何呢，比如吸收发射，以及激发态等？
答：pcm调控对 ...

谢谢钟老师指点，非常感谢钟老师。
当时我文中说的差距比较大，是说的计算出来的eps。
非常感谢钟老师。

作者
Author: ggdh 时间: 2017-11-5 15:52

小范范1989 发表于 2017-11-5 15:15
谢谢钟老师指点，非常感谢钟老师。
当时我文中说的差距比较大，是说的计算出来的eps。
非常感谢钟老师 ...

我知道，你是说gromacs模拟的和单分子计算出来的eps差很大，不知道具体是差多少？

作者
Author: 小范范1989 时间: 2017-11-5 17:54

ggdh 发表于 2017-11-5 15:52
我知道，你是说gromacs模拟的和单分子计算出来的eps差很大，不知道具体是差多少？

钟老师好，gromacs计算的图一，单分子计算的是图二。
我感觉gromacs计算的还行。不知道对不对。或者是我在计算单分子得到的时候，里面求体积或者polar的时候，有细节没有注意好。
谢谢钟老师指点。

作者
Author: ggdh 时间: 2017-11-5 19:27

小范范1989 发表于 2017-11-5 17:54
钟老师好，gromacs计算的图一，单分子计算的是图二。
我感觉gromacs计算的还行。不知道对不对。或者是我 ...

gmx怎么算的才1点多，要知道真空也才1
是不是算错了。
倒是单分子算的值还比较靠谱

作者
Author: 小范范1989 时间: 2017-11-6 07:26
本帖最后由小范范1989 于 2017-11-6 07:32 编辑

ggdh 发表于 2017-11-5 19:27
gmx怎么算的才1点多，要知道真空也才1
是不是算错了。
倒是单分子算的值还比较靠谱

钟老师好，gmx我就是通过当时sob老师的那个培训班。P117静态介电常数。那个公式求得的。（当时上课的时候，碰到过钟老师和孙老师，两位老师坐在我前面，我就在您后面，打过招呼。哈哈哈）

作者
Author: 阳光总在风雨后 时间: 2017-11-6 10:24

ggdh 发表于 2017-11-5 11:00
我的意思是调控的时候不要加scrf，但是你的w显然是在调控的时候加了scrf后获得的。
因此你要在调控的时 ...

非常感谢老师的指点!

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