sobMECP能计算S0和S1交叉点吗？

ZHANGZY · 发表于 Post on 2016-4-30 13:02:00

我知道求s1/s0 intersection points 可以用 casscf
S1/T1 ISC point 用 MECP，
那么想问一下 s1/s0 IC point 可不可以用 MECP 程序（sobMECP程序）？

sobereva · 发表于 Post on 2016-4-30 17:56:53

原理上可以，只要修改一下代码中收敛标准的判断即可。原先用于系间窜越MECP搜索，收敛条件之一是必须能量差小于一定值，但S0-S1交叉一般是避免交叉会有一定能量裂分，所以能量差必须足够接近0这个条件应去掉。

小范范1989 · 发表于 Post on 2016-6-21 10:30:09

sobereva 发表于 2016-4-30 17:56
原理上可以，只要修改一下代码中收敛标准的判断即可。原先用于系间窜越MECP搜索，收敛条件之一是必须能量差 ...

sob老师好，想在想通过这个sob-mecp来搜索S1-S0的交叉点。我看您说需要修改：
如果要调整收敛判断阈值，应修改temp/MECP.f的第300行。总共有5个标准。TDE是能量变化，TDXMax和TDXRMS是位移最大值和方均根，TGMax和TGRMS是受力最大值和方均根。默认情况下，受力和位移判断标准比Gaussian09几何优化默认判断标准松将近一倍。
我找到程序300行
PARAMETER (TDE=5.d-5,TDXMax=4.d-3,TDXRMS=2.5d-3,TGMax=7.d-4,TGRMS=5.d-4)
修改TDE数值吗？多少？谢谢sob老师。

sobereva · 发表于 Post on 2016-6-21 15:02:27

小范范1989 发表于 2016-6-21 10:30
sob老师好，想在想通过这个sob-mecp来搜索S1-S0的交叉点。我看您说需要修改：
如果要调整收敛判断阈值， ...

不是改这个。要搜索S1-S0交叉点得把“两个态之间能量差必须很小”这个收敛要求去掉或者设得很松，这是在其它地方改的。代码不长，从头到尾看看就能找到。

小范范1989 · 发表于 Post on 2016-6-23 09:32:07

sobereva 发表于 2016-6-21 15:02
不是改这个。要搜索S1-S0交叉点得把“两个态之间能量差必须很小”这个收敛要求去掉或者设得很松，这是在 ...

sob老师好，我比较笨，没有从.f程序中找到老师说的这个地方，不知道改哪个地方，希望sob老师能给指明一下，同时数值改为多少？
顺带问一下sob老师，这个MECP也能搜索T2-S1的交叉点是吧？谢谢sob老师。

sobereva · 发表于 Post on 2016-6-24 03:54:59

小范范1989 发表于 2016-6-23 09:32
sob老师好，我比较笨，没有从.f程序中找到老师说的这个地方，不知道改哪个地方，希望sob老师能给指明一下 ...

MECP.f里
   IF (GMax .lt. TGMax) THEN
      PConv(1) = .true.
      flags(1) = "YES"
   END IF
   IF (GRMS .lt. TGRMS) THEN
      PConv(2) = .true.
      flags(2) = "YES"
   END IF
   IF (DXMax .lt. TDXMax) THEN
      PConv(3) = .true.
      flags(3) = "YES"
   END IF
   IF (DXRMS .lt. TDXRMS) THEN
      PConv(4) = .true.
      flags(4) = "YES"
   END IF
   IF (DE .lt. TDE) THEN
      PConv(5) = .true.
      flags(5) = "YES"
   END IF

就这5个判断收敛的条件，把IF (DE .lt. TDE) THEN改成IF (.true.) THEN就行了，使这个条件判断总是为真。

原理上可以搜索T2-S1，只要让程序正确读取激发态计算给出的激发态能量信息就行了。

小范范1989 · 发表于 Post on 2016-6-24 07:33:03

sobereva 发表于 2016-6-24 03:54
MECP.f里
IF (GMax .lt. TGMax) THEN
PConv(1) = .true.

谢谢sob老师的指点，谢谢老师。
在搜索的时候S1这样写
%mem=2GB
%nproc=12
%chk=S1-gas.chk
#n pbe1pbe/6-311G* force guess(read) td

First State

0 1

T2这样写
%mem=2GB
%nproc=12
%chk=T2-gas.chk
#n pbe1pbe/6-311G* force guess(read) td(triplet,root=2)

Second State

0 1
对吧？谢谢老师

sobereva · 发表于 Post on 2016-6-24 16:26:42

小范范1989 发表于 2016-6-24 07:33
谢谢sob老师的指点，谢谢老师。
在搜索的时候S1这样写
%mem=2GB

对

ZHANGZY · 发表于 Post on 2016-7-10 13:07:40

有没有已发表的文献用这个方法？？

ZHANGZY · 发表于 Post on 2016-7-12 22:21:38

首先要得到激发态S1的PES 用tddft optimization 不合适，因为在conical intersection附近TDDFT失效
A serious drawback for photochemical applications, however, is that TDDFT predicts the wrong topology for any CI that involves the reference state
在此基础上上应用MECP程序的方法求S0/S1交叉点，似乎就不行了
当然求S0/T1就不存在这个问题。

sobereva · 发表于 Post on 2016-7-13 07:25:53

ZHANGZY 发表于 2016-7-12 22:21
首先要得到激发态S1的PES 用tddft optimization 不合适，因为在conical intersection附近TDDFT失效
A seri ...

对于避免交叉，只要能量分裂够大显然是可以用的
当然交叉上的情况TDDFT没法用，只能MCSCF或多参考方法

blueyangliu · 发表于 Post on 2016-12-4 23:40:44

小范范1989 发表于 2016-6-24 07:33
谢谢sob老师的指点，谢谢老师。
在搜索的时候S1这样写
%mem=2GB

你用的什么服务器算这个的？我看你设置的%nproc=12，是单机还是集群

小范范1989 · 发表于 Post on 2016-12-5 08:03:36

blueyangliu 发表于 2016-12-4 23:40
你用的什么服务器算这个的？我看你设置的%nproc=12，是单机还是集群

集群，我们点买了好几百个核。

blueyangliu · 发表于 Post on 2016-12-5 14:51:47

小范范1989 发表于 2016-12-5 08:03
集群，我们点买了好几百个核。

SobMECP程序能放到集群上计算吗？我没成功呢，怎么做到的

小范范1989 · 发表于 Post on 2016-12-5 17:52:24

blueyangliu 发表于 2016-12-5 14:51
SobMECP程序能放到集群上计算吗？我没成功呢，怎么做到的

能，就是在集群上。我之前用过，可以的。教程我的都有啊。你不能用，是因为你没把文件全部放进去，少文件吧

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

[综合交流] sobMECP能计算S0和S1交叉点吗？