lastzealot 发表于 2021-1-10 23:20 这个是我自己设置的,你也可以算3个态或者两个态,因为之前有老师说尽量多算几个态,否则可能会漏掉第一激发态。 |
| 你好,我看你在每个不可约表示下算了四个态,请问这个数量来自于哪?为啥会计算这么多? |
jiangning198511 发表于 2020-12-18 09:38 好的,谢谢 ![]() |
jiangning198511 发表于 2020-12-18 09:35 您这样一说我感觉您的理解是正确的,感谢 |
原子分子物理 发表于 2020-12-17 20:20 可以这样操作,下面链接是MOLCAS计算的一些例子,你可以参考一些里面的一些技巧和思路 http://bbs.keinsci.com/thread-5632-1-1.html |
原子分子物理 发表于 2020-12-17 18:39 所谓激发能都是相对基态(最低能量)而言的,所以对于N2而言,不管是单重激发态或者是三重激发态的激发能都应该以基态为基准。你所说的1就是这样计算的,而根据你所说的2,我理解为你是通过把最低的三重态当做基准来计算激发能了,这是不对的。(不知道我的理解是否正确) |
jiangning198511 发表于 2020-12-17 17:27 我想再问您一个问题,如果我要算N2+离子的激发能,是不是同N2分子的三重态一样,也用N2分子基态的轨道进行计算? |
jiangning198511 发表于 2020-12-17 17:27 我把文献PDF发出来了您可以直接下载下来 ![]() |
jiangning198511 发表于 2020-12-17 17:27 不好意思,最近查的文献太多,可能我不小心搞混了,文献题目是“Identifying and removing intruder states in multireference Mo ”ller–Plesset perturbation theory” 它的DOI: 10.1063/1.1345510 我用您贴出的计算了一下,有两个疑问1用八个对称性下的三重态的各个态的值减去state1.1(经查看1.1是八个对称性下单重态和八个对称性下三重态各个态中的能量最小的态)从而得出对应的各个态的垂直激发能2在八个对称性下的三重态的各个态中找出能量最小的态(经查看,三重态中能量最小的态是1.5),把数据整理了之后,发现第1种更贴近文献中数据(误差在0.5ev左右),第二种数据相差甚大,所以是不是第一种才是正确的? |
原子分子物理 发表于 2020-12-17 16:08 客气了,共同学习 PS:你给的文献信息页码应该是3958 并且文章里面也没有这个表,是不是你记错了 ![]() |
jiangning198511 发表于 2020-12-17 15:56 首先谢谢老师的指导,我待会尝试一下,对比一下计算结果。帖子里的文献是: Finley, J. P.; Witek, H. A. J Chem Phys 2000, 112, 3968. 再次感谢老师,有问题再来请教您。 |
本帖最后由 jiangning198511 于 2020-12-17 15:57 编辑 原子分子物理 发表于 2020-12-17 15:11 这个地方可以讨论一下,1.CASSCF之前的HF只是提供活性空间,所以原则上只要是基态的计算就可以了,手册里的O2分子例子就是因为O2的基态是三重态 所以HF card 里面是三重态,但CASSCF计算激发态的设置和HF没什么关系 就类似于你计算O2的单重激发态没有必要在HF里设置单重态计算而应该是基态(三重态) 2. 输入文件我估计应该可以这样写 *** , N2 print ,basis , orbitals angstrom geometry={ N,, 0.00000000 0.00000000 0.55920600 N,, 0.00000000 0.00000000 -0.55920600 } basis=avdz {rhf wf,14,1,0 } {multi;occ,3,1,1,0,3,1,1,0;closed,1,0,0,0,1,0,0,0 wf ,14 ,1 ,0 ;state ,4 ... 八个对称性下的单重态 wf,14,1,2;state,4 ...八个对称性下的三重态 } 你可以和你之前的计算比较一下 PS: 你帖子里面参考文献能否给一下 谢谢 |
本帖最后由 原子分子物理 于 2020-12-17 15:56 编辑 jiangning198511 发表于 2020-12-17 11:28 第一步的HF中是用高斯做了一个关于N2的0 3的计算,OCC卡和open卡以及活化空间就是根据算出的轨道得出的,我刚开始也觉得应该写基态的,但是基态下对应的三重态的wf卡和occ卡应该怎么写呢?,这样是不是就没有OPEN卡了?怎么让系统知道我算的是三重态而不是单重态?这是我的困惑所在。把单重态和三重态放在一个文件中,那么对应的{HF}和{multi}应该怎么写呢,针对老师提出的这两个问题,老师能否给我举个例子啊,谢谢老师。 |
| 你的输入文件中第一步HF为什么不是N2的基态,而是激发态,虽然只是给后续计算提供初始轨道,但你这个激发态的轨道是否会导致活性空间有问题?,还有你可以把单重态 和三重态写在一个文件中计算,这样是否更方便处理? |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-21 00:14 , Processed in 0.176636 second(s), 25 queries , Gzip On.