求问 Molpro 势能面交叉相关的问题

查看数: 6761 | 评论数: 8 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2021-7-23 13:53

各位大佬，我目前是想用Molpro计算PH2自由基的基态和第一第二激发态的势能面，采用的Cs对称性。基态计算没什么问题，但是第一激发态和第二激发态的势能面 ...

回复 Reply

chenguangan 发表于 Post on 2021-7-24 23:40:53

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

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

wzkchem5 发表于 Post on 2021-7-24 15:17:52

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

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

chenguangan 发表于 Post on 2021-7-24 12:36:59

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

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

wzkchem5 发表于 Post on 2021-7-23 20:03:00

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

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

chenguangan 发表于 Post on 2021-7-23 18:33:09

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

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

wzkchem5 发表于 Post on 2021-7-23 16:53:44

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

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

chenguangan 发表于 Post on 2021-7-23 16:38:30

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

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

wzkchem5 发表于 Post on 2021-7-23 15:09:12

我记得你之前说你的分子其实是C2v点群，如果是这样的话，最不用动脑子的方法是用C2v重新算一遍，然后根据C2v的结果来手动换Cs的态的顺序

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register