BOMD能算出绝热和非绝热两条解离途径吗？

算一算 · 发表于 Post on 2021-3-12 19:38:53

本帖最后由算一算于 2021-3-12 21:37 编辑

正打算研究一下高斯的BOMD的输入文件编写方式，但是不清楚这个方法到底能不能解决我现在这个问题：

就是碳碘在一起 (CI)3+这个离子，算它解离方式，其中一个是绝热的，是解离成为 C+ 和 I2+，另一个是非绝热的，解离成不带电的 C 和 I3+。

我刚试了一下高斯的BOMD，发现它只给我绝热的那个，我是不是要在输入文件里写些特殊的东西才能让它把非绝热的那个也算出来？还是说高斯的BOMD压根就不能算这个？
如果不能，那用什么软件的什么方法能把这两个解离方式都模拟出来呢？想比较一下这两条解离路线的动能势能变化
=======
刚才下面有人有疑问，然后我发现是没上角标，我这些数字其实全是电荷数，应该都是上角标，但是论坛里没有上下角标的格式……是碳碘一共带三个正电荷，然后两条路线通向不同的电荷分配方式

biogon · 发表于 Post on 2021-4-27 15:09:57

wzkchem5 发表于 2021-4-27 14:31
我就是说不+Q的结果应该不能发表吧

发表应该是没问题的吧，毕竟合理性是没问题的，看了看这年头还是有用MRCISD发表的，发好点的刊物那就另当别论

wzkchem5 · 发表于 Post on 2021-4-27 14:31:14

biogon 发表于 2021-4-27 10:01
应该是不能实现的，你可以问问@zjxitcc
另外不+Q只是没有大小一致性吧，可信度还是有的，能和其他程序的 ...

我就是说不+Q的结果应该不能发表吧

zjxitcc · 发表于 Post on 2021-4-27 12:10:42

wzkchem5 发表于 2021-4-27 00:14
这样做MRCISD支持+Q吗？不+Q没有可信度吧

不支持+Q。+Q还是得用更专业的软件（如OpenMolcas, Molpro, ORCA）

biogon · 发表于 Post on 2021-4-27 10:01:40

wzkchem5 发表于 2021-4-27 00:14
这样做MRCISD支持+Q吗？不+Q没有可信度吧

应该是不能实现的，你可以问问@zjxitcc
另外不+Q只是没有大小一致性吧，可信度还是有的，能和其他程序的MRCISD结果对的上

wzkchem5 · 发表于 Post on 2021-4-27 00:14:16

biogon 发表于 2021-4-21 20:36
jxzou-MOKIT作者 2021/2/23 18:07:30
RASCI某种特殊情况=MRCISD
jxzou-MOKIT作者 2021/2/23 18:07:53

这样做MRCISD支持+Q吗？不+Q没有可信度吧

biogon · 发表于 Post on 2021-4-26 21:21:38

喵星大佬发表于 2021-4-26 19:59
话说我一直没试过高斯里cas(x,x) mp2，这玩意和caspt2有木有区别。。。。

这是非主流的MRPT2:OVB-MP2，没啥用，大部分基本都是用的CASPT2和NEVPT2

zjxitcc · 发表于 Post on 2021-4-26 20:04:30

喵星大佬发表于 2021-4-26 19:59
话说我一直没试过高斯里cas(x,x) mp2，这玩意和caspt2有木有区别。。。。

有区别。OVB-MP2几乎没人用，CASPT2仍是主流多参考二阶微扰

喵星大佬 · 发表于 Post on 2021-4-26 19:59:10

biogon 发表于 2021-4-21 20:36
jxzou-MOKIT作者 2021/2/23 18:07:30
RASCI某种特殊情况=MRCISD
jxzou-MOKIT作者 2021/2/23 18:07:53

话说我一直没试过高斯里cas(x,x) mp2，这玩意和caspt2有木有区别。。。。

biogon · 发表于 Post on 2021-4-21 20:36:43

喵星大佬发表于 2021-4-21 12:25
高斯怎么实现的

jxzou-MOKIT作者 2021/2/23 18:07:30
RASCI某种特殊情况=MRCISD
jxzou-MOKIT作者 2021/2/23 18:07:53
RASCI更为广义，通过恰当定义活性空间，可以退化为MRCISD
下面是算某体系用mokit自动生成的睾斯算MRCISD的关键词

#p CAS(14,28,ras(2,5,2,19))/chkbasis nosymm guess=read geom=allcheck scf(maxcycle=-1)

复制代码

喵星大佬 · 发表于 Post on 2021-4-21 12:25:41

biogon 发表于 2021-4-21 09:35
那有多半是TDDFT算的能量不准，不过用了MRCI也不能说一定有希望
ps：高斯是可以实现MRCISD的

高斯怎么实现的

算一算 · 发表于 Post on 2021-4-21 10:41:51

biogon 发表于 2021-4-21 09:35
那有多半是TDDFT算的能量不准，不过用了MRCI也不能说一定有希望
ps：高斯是可以实现MRCISD的

除了TDDFT以外高斯里的那些复杂的方法我也都试了，都算不出来，已经放弃了

biogon · 发表于 Post on 2021-4-21 09:35:54

算一算发表于 2021-4-2 20:28
我的两条势能曲线都是单调下降到解离态的，都没有极小点。我在研究的是能不能从下面那条曲线（C带1+，I带 ...

那有多半是TDDFT算的能量不准，不过用了MRCI也不能说一定有希望
ps：高斯是可以实现MRCISD的

算一算 · 发表于 Post on 2021-4-3 18:32:56

wzkchem5 发表于 2021-4-2 20:48
假如TDDFT计算结果准确的话，那基本没戏了。但是保险起见还是要用比如说MRCISD+Q重新算一下

你好，我还想继续问一个问题：我刚做完的这个是双原子体系，所以扫势能曲线是rigid scan，那如果是三原子或更多原子的体系呢？是要rigid scan 还是 relaxed scan？我猜还是要rigid scan？
还是这个实验，我今天跟同事讨论了一下，觉得有可能一开始是从 H-C-I 离子解离，并不一定是从 C-I 离子开始的。所以想再扫一次这两个势能曲线 (绝热和非绝热两条曲线），扫 H-C-I 离子。还是能量 vs C-I键长。我在想是不是只能做 rigid scan？如果是 relaxed scan，随着 C-I 的解离，H 和 C-I 的相对位置也会逐渐变化，那么就算是这两条势能曲线相交了，可是在交点，虽然C-I键长和总能量是一样的，但是这个H-C-I的几何结构不同，所以这种charge transfer也是无法实现？所以必须 rigid scan，我这个思路对吗？

算一算 · 发表于 Post on 2021-4-2 21:00:21

wzkchem5 发表于 2021-4-2 20:48
假如TDDFT计算结果准确的话，那基本没戏了。但是保险起见还是要用比如说MRCISD+Q重新算一下

好的，多谢了

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

[综合交流] BOMD能算出绝热和非绝热两条解离途径吗？

浏览过的版块