计算化学公社

标题: 求问 Molpro 势能面交叉相关的问题 [打印本页]

作者
Author: chenguangan 时间: 2021-7-23 13:53
标题: 求问 Molpro 势能面交叉相关的问题
各位大佬，    我目前是想用Molpro计算PH2自由基的基态和第一第二激发态的势能面，采用的Cs对称性。
   基态计算没什么问题，但是第一激发态和第二激发态的势能面存在交叉现象，而且这两个态在Cs对称性下都是A‘
   所以在输出能量点的时候，A‘对称性下第一个态输出的点是能量最低的点，也就是说有一部分点本应该属于第二激发态，
   但是由于计算出来的能量低于第一激发态，所以被Molpro当成第一个态了。
   这就导致输出文件里，A’的前两个态有一部分点的能量给弄颠倒了，
   请问遇到这种情况该怎么设置？
   或者有什么办法能把这些能量颠倒的构型找出来，我再手动换过来？

作者
Author: wzkchem5 时间: 2021-7-23 15:09
我记得你之前说你的分子其实是C2v点群，如果是这样的话，最不用动脑子的方法是用C2v重新算一遍，然后根据C2v的结果来手动换Cs的态的顺序

作者
Author: chenguangan 时间: 2021-7-23 16:38

wzkchem5 发表于 2021-7-23 15:09
我记得你之前说你的分子其实是C2v点群，如果是这样的话，最不用动脑子的方法是用C2v重新算一遍，然后根据C2 ...

用C2v算出来的话倒是能做出来3维的面，能找到交叉的区域，但是我应该怎么确认我把所有的点都找到了呢？万一漏下几个咋办啊

作者
Author: wzkchem5 时间: 2021-7-23 16:53

chenguangan 发表于 2021-7-23 09:38
用C2v算出来的话倒是能做出来3维的面，能找到交叉的区域，但是我应该怎么确认我把所有的点都找到了呢？万 ...

你在C2v对称性下计算只是为了能指认Cs的态，之后的研究都用Cs做

作者
Author: chenguangan 时间: 2021-7-23 18:33

wzkchem5 发表于 2021-7-23 16:53
你在C2v对称性下计算只是为了能指认Cs的态，之后的研究都用Cs做

我明白你的意思，但是我用C2v只能确切的指认出两个PH键长相等的那些点啊，那些不相等的真Cs对称性的点还是没法确切的找出来啊

作者
Author: wzkchem5 时间: 2021-7-23 20:03

chenguangan 发表于 2021-7-23 11:33
我明白你的意思，但是我用C2v只能确切的指认出两个PH键长相等的那些点啊，那些不相等的真Cs对称性的点还 ...

首先你可以指认P-H键近似相等的结构的点，然后再指认P-H键差得稍微多一点的点，以此类推

作者
Author: chenguangan 时间: 2021-7-24 12:36

wzkchem5 发表于 2021-7-23 20:03
首先你可以指认P-H键近似相等的结构的点，然后再指认P-H键差得稍微多一点的点，以此类推

要是用C2v算一下的话，可以少一个坐标变量，势能面能从4维变3维，画个图应该是能找出来两个面相交的那个曲线，我就是担心我选点的时候不彻底，漏下几个点，想问问有没有啥更有效的方法。比如molpro能不能输出每个点的轨道信息，我根据轨道不就可以确定哪个点属于哪个态了吗

作者
Author: wzkchem5 时间: 2021-7-24 15:17

chenguangan 发表于 2021-7-24 05:36
要是用C2v算一下的话，可以少一个坐标变量，势能面能从4维变3维，画个图应该是能找出来两个面相交的那个 ...

我觉得我说的办法应该是更简单的。原则上你也可以从每个点的轨道信息、组态信息来判断，也更严格一些，但是很麻烦，尤其是你不仅要知道轨道的形状，还要知道轨道的相位，可能需要把每个坐标下的轨道都作出来看

作者
Author: chenguangan 时间: 2021-7-24 23:40

wzkchem5 发表于 2021-7-24 15:17
我觉得我说的办法应该是更简单的。原则上你也可以从每个点的轨道信息、组态信息来判断，也更严格一些，但 ...

那我先试试吧，我现在用C2v对称性计算老是报错，应该是算到交叉的那部分点不收敛，好像要用什么态平均算法才行，得再学一学

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