固相下优化出的S1，T1态能量低于S0态的问题

阳光总在风雨后

S0态在这里  attach://11244.gjf

阳光总在风雨后

不好意思，刚刚那个操作失误！  请问：固相下，优化出的S1，T1态能量低于S0态的，这是问题出在哪里？下面贴出S0，S1，T1态输入文件

ggdh

用普通的调控就行了。还非要用pcm调控。。。pcm调控算IP EA还行。算激发态容易低估。

阳光总在风雨后

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

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

ggdh

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

你调控的时候不要加scrf

阳光总在风雨后

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

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

小范范1989

阳光总在风雨后 发表于 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里面有文件例子。

小范范1989

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还是怎么样的方法？
望钟老师指点一二，谢谢钟老师。

ggdh

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

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

阳光总在风雨后

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
给您带来困扰实在不好意思！望赐教

ggdh

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

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

ggdh

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

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

ggdh

小范范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。

小范范1989

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

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

ggdh

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

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