请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 238|回复: 5

[新手求助] VDEs计算问题求助

[复制链接]

6

帖子

0

威望

39

eV
积分
45

Level 2 能力者

发表于 2019-6-11 12:09:25 | 显示全部楼层 |阅读模式
大家好,想请教大家一个问题,我想算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就可以得到。


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


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

图片_20190611.png
图片_20190611-1.png

2万

帖子

25

威望

2万

eV
积分
44986

管理员

公社社长

发表于 2019-6-11 12:46:29 | 显示全部楼层
MOMO指代不明

至于“从两个电子占据一个轨道变成分占两个轨道”,显然限制性闭壳层占据轨道都是双占据的,非限制性开壳层都是单占据的
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

6

帖子

0

威望

39

eV
积分
45

Level 2 能力者

 楼主| 发表于 2019-6-11 13:46:42 | 显示全部楼层
sobereva 发表于 2019-6-11 12:46
MOMO指代不明

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

老师,您好,那个图片我又重新弄上去了。我在算
阴离子与中性分子的第一个激发态的VDE,我用的是 56还有57换了轨道;算阴离子与中性分子的第二个激发态的VDE,我首先是用的中性分子基态算的55和57交换了轨道,然后我又用第一个激发态的构型做了55和56换了轨道,这两种算法都又回到了中性分子状态,能量也是中性分子基态能量。那请问怎么计算阴离子与中性分子的第二个激发态的VDE呢?希望老师可以帮忙解答一下,谢谢老师!
图20190611-2.png

2万

帖子

25

威望

2万

eV
积分
44986

管理员

公社社长

发表于 2019-6-12 22:23:41 | 显示全部楼层
VASP-DFT 发表于 2019-6-11 13:46
老师,您好,那个图片我又重新弄上去了。我在算阴离子与中性分子的第一个激发态的VDE,我用的是 56还有57 ...

仔细看此文弄懂ΔSCF的原理
乱谈激发态的计算方法
http://sobereva.com/265http://bbs.keinsci.com/thread-415-1-1.html
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

6

帖子

0

威望

39

eV
积分
45

Level 2 能力者

 楼主| 发表于 2019-6-14 10:25:22 | 显示全部楼层
本帖最后由 VASP-DFT 于 2019-6-14 10:29 编辑
sobereva 发表于 2019-6-12 22:23
仔细看此文弄懂ΔSCF的原理
乱谈激发态的计算方法
http://sobereva.com/265http://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                                                                  
                                                                        


2万

帖子

25

威望

2万

eV
积分
44986

管理员

公社社长

发表于 2019-6-14 23:59:15 | 显示全部楼层

两个轨道的不可约表示必须不同,而且一个是占据的一个是非占据的,这样的交换做SCF才有意义,否则结果和不交换必然相同
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-8-26 16:39 , Processed in 0.185207 second(s), 26 queries .

快速回复 返回顶部 返回列表