计算化学公社

标题: 开壳层单重态的casscf计算 [打印本页]

作者
Author: xujian 时间: 2017-3-11 11:28
标题: 开壳层单重态的casscf计算
用molpro做开壳层单重态的casscf计算，如何给出电子单占轨道的信息，比如体系有38个电子，闭壳层单重态是19个轨道，开壳层的话是20个轨道，输入文件中的区别体现在哪里，是给出hf轨道的占据数吗？但如果给出20的话，电子数与自旋多重度之和又不是占据数的两倍，运行会报错，想请各位大神指导一下。谢谢！

作者
Author: jiangning198511 时间: 2017-3-11 14:36
报错的是HF部分还是CASSCF部分，CASSCF部分只要指定好对称性和自旋就可以了

作者
Author: xujian 时间: 2017-3-11 15:38

jiangning198511 发表于 2017-3-11 14:36
报错的是HF部分还是CASSCF部分，CASSCF部分只要指定好对称性和自旋就可以了

是HF部分，我应该怎么区分开壳层和闭壳层

作者
Author: jiangning198511 时间: 2017-3-11 20:33
HF部分只是为后面提供初始轨道，闭壳层也可以。只要后面的设置搞对就可以了。

作者
Author: xujian 时间: 2017-3-13 08:43
本帖最后由 xujian 于 2017-3-13 08:55 编辑

jiangning198511 发表于 2017-3-11 20:33
HF部分只是为后面提供初始轨道，闭壳层也可以。只要后面的设置搞对就可以了。

您说的有些抽象

作者
Author: jiangning198511 时间: 2017-3-13 08:50
本帖最后由 jiangning198511 于 2017-3-13 09:12 编辑

xujian 发表于 2017-3-13 08:43
您说的有些抽象，请问加下qq,像您请教吗？

我们以H2为例子，H2的基态 sigma^2 如下图(1)，而它的开壳层单重态为(3),故无法用HF方法描述，需要使用CASSCF方法，对于这个例子，你可以先做一个闭壳层计算，分析一下轨道占据(如下图)，如果你需要计算开壳层单重态(3)，就需要在MCSCF模块中设置这个态，其对称性就是轨道A 乘以轨道B , Spin : 0,最后打印一下CI系数，系数应该是和图中(3)的情况一致。
输入文件可以如下：

作者
Author: xujian 时间: 2017-3-13 08:55
本帖最后由 xujian 于 2017-3-13 09:01 编辑

jiangning198511 发表于 2017-3-13 08:50
我们以H2为例子，H2的基态 sigma^2 如下图(a)

我知道啊，那1和3在molpro输入的时候差别体现在哪里？1和3的hf占据就不一样啊，1的的占据是1,3的占据是2，3写成{hf,
occ,2} ？

作者
Author: jiangning198511 时间: 2017-3-13 09:13

xujian 发表于 2017-3-13 08:55
我知道啊，那1和3在molpro输入的时候差别体现在哪里？1和3的hf占据就不一样啊，1的的占据是1,3的占据是2 ...

请看回复，我好久都没用MOLPRO，只是记得一些流程，请见谅

作者
Author: jiangning198511 时间: 2017-3-13 09:16

xujian 发表于 2017-3-13 08:55
我知道啊，那1和3在molpro输入的时候差别体现在哪里？1和3的hf占据就不一样啊，1的的占据是1,3的占据是2 ...

HF中不用设置开壳层单重态，这个应该是在CASSCF中设置，HF中只是提供猜测轨道，CASSCF中要优化轨道的

作者
Author: xujian 时间: 2017-3-13 09:23

jiangning198511 发表于 2017-3-13 09:16
HF中不用设置开壳层单重态，这个应该是在CASSCF中设置，HF中只是提供猜测轨道，CASSCF中要优化轨道的

奥奥，我的疑惑就在这里，这样的话，hf提供的猜测轨道完全一样，以我的体系为例，一共有38个电子，开壳层单重态的话应该有20个占据，但hf提供的初猜轨道是19个，我现在要六个活性轨道，5个占据，一个非占据，选择的时候还是和闭壳层一样，只能选到4个占据和两个非占据。

作者
Author: jiangning198511 时间: 2017-3-13 09:38

xujian 发表于 2017-3-13 09:23
奥奥，我的疑惑就在这里，这样的话，hf提供的猜测轨道完全一样，以我的体系为例，一共有38个电子，开壳层 ...

把要占据的轨道都选进去就可以了，

，比如说你要单占据第20 和21 ，就将他们都选进来，在CASSCF设置态的对称性就可以了，不过算完后要看看轨道是否跑出去了。CASSCF计算要多多注意，很容易出错。

作者
Author: xujian 时间: 2017-3-13 09:42

jiangning198511 发表于 2017-3-13 09:38
把要占据的轨道都选进去就可以了，，比如说你要单占据第20 和21 ，就将他们都选进来，在CASSCF ...

好的，非常感谢！

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