计算化学公社

标题: VDEs计算问题求助 [打印本页]

作者
Author: VASP-DFT 时间: 2019-6-11 12:09
标题: VDEs计算问题求助
大家好，想请教大家一个问题，我想算Al7Na2-分子的VDEs，VDE我用的是阴离子构型得到的中性分子，并用中性分子能量减去阴离子能量，而计算其他的VDE，也就是阴离子减去中性分子激发态能量，我用的guess=alter来得到，这是我的输入文件


算阴离子与中性分子的第一个激发态的VDE，我用的是

%chk=8141.chk
%mem=6GB
%nprocshared=24
# b3lyp/6-311++g(3df,3pd) guess=alter

opt

0 2
Al                0.16131300 -1.76741600 0.00434600
Al                1.66206500 -0.14153300 -1.52205500
Al             -2.66449300 -1.26393000 0.00293500
Al             -0.96251800 -0.08488800 -1.64739500
Al             -0.96279600 -0.07692300 1.64777300
Al                0.11393200 1.68430300 -0.00405600
Al                1.66178500 -0.13389600 1.52269500
Na                4.13179200 0.38500400 -0.00079900
Na             -2.96095100 1.72369200 -0.00421500

56 57

算阴离子与中性分子的第二个激发态的VDE，我用的输入文件是

%chk=Al7Na2-.chk
%mem=6GB
%nprocshared=24
# b3lyp/6-311++g(3df,3pd) guess=alter SCF=Symm

opt

0 2
Al                0.16131300 -1.76741600 0.00434600
Al                1.66206500 -0.14153300 -1.52205500
Al             -2.66449300 -1.26393000 0.00293500
Al             -0.96251800 -0.08488800 -1.64739500
Al             -0.96279600 -0.07692300 1.64777300
Al                0.11393200 1.68430300 -0.00405600
Al                1.66178500 -0.13389600 1.52269500
Na                4.13179200 0.38500400 -0.00079900
Na             -2.96095100 1.72369200 -0.00421500

55 57

文献图片所示：HOMO-1与下面的轨道的VDE可以用guess=alter得到，但是MOMO与HOMO-1的VDE是怎么得到呢，他们之间不是通过交换轨道就可以得到，而是从两个电子占据一个轨道变成分占两个轨道，感觉不是通过guess=alter就可以得到。

中性分子的轨道图以及文献参考图已上传在附件。

希望大家可以帮忙解答一下，非常感谢！

作者
Author: sobereva 时间: 2019-6-11 12:46
MOMO指代不明

至于“从两个电子占据一个轨道变成分占两个轨道”，显然限制性闭壳层占据轨道都是双占据的，非限制性开壳层都是单占据的

作者
Author: VASP-DFT 时间: 2019-6-11 13:46

sobereva 发表于 2019-6-11 12:46
MOMO指代不明

至于“从两个电子占据一个轨道变成分占两个轨道”，显然限制性闭壳层占据轨道都是双占据的 ...

老师，您好，那个图片我又重新弄上去了。我在算

阴离子与中性分子的第一个激发态的VDE，我用的是 56还有57换了轨道；算阴离子与中性分子的第二个激发态的VDE，我首先是用的中性分子基态算的55和57交换了轨道，然后我又用第一个激发态的构型做了55和56换了轨道，这两种算法都又回到了中性分子状态，能量也是中性分子基态能量。那请问怎么计算阴离子与中性分子的第二个激发态的VDE呢？希望老师可以帮忙解答一下，谢谢老师！

作者
Author: sobereva 时间: 2019-6-12 22:23

VASP-DFT 发表于 2019-6-11 13:46
老师，您好，那个图片我又重新弄上去了。我在算阴离子与中性分子的第一个激发态的VDE，我用的是 56还有57 ...

仔细看此文弄懂ΔSCF的原理
乱谈激发态的计算方法
http://sobereva.com/265（http://bbs.keinsci.com/thread-415-1-1.html）

作者
Author: VASP-DFT 时间: 2019-6-14 10:25
本帖最后由 VASP-DFT 于 2019-6-14 10:29 编辑

sobereva 发表于 2019-6-12 22:23
仔细看此文弄懂ΔSCF的原理
乱谈激发态的计算方法
http://sobereva.com/265（http://bbs.keinsci.com/t ...

老师，您好。

我用guess=alter交换轨道，加了SCF=qc/dm，Log写了贝塔电荷交换，

No Alpha orbitals switched.

Pairs of Beta  orbitals switched:
19  20

但是能量还是没变，貌似轨道并没有换过去，请问是怎么回事呢？望老师帮忙解答。

输入文件如下
%chk=8141.chk
%mem=6GB
%nprocshared=24
# b3lyp/6-311++g(3df,3pd)
# guess=(alter,read)  geom=allcheck SCF=qc


55 57

作者
Author: sobereva 时间: 2019-6-14 23:59

VASP-DFT 发表于 2019-6-14 10:25
老师，您好。

两个轨道的不可约表示必须不同，而且一个是占据的一个是非占据的，这样的交换做SCF才有意义，否则结果和不交换必然相同

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