计算化学公社

标题: 结构优化完怎么查看点群 [打印本页]

作者
Author: zgdx 时间: 2022-8-26 09:48
标题: 结构优化完怎么查看点群
各位老师，我用b3lyp泛函和lanl2dz基组，对7个原子的过渡金属团簇WGe6进行结构优化后，单重态、三重态、五重态优化得到的结构是不一样的，但是查看点群的时候，当前点群全部显示是C1点群，请问我优化的结果是错的，还是我查看点群的位置不对（在产生的fchk文件中的Edit-Point Group中查看的右上方的当前点群）？应该怎么得到正确的点群呢？

作者
Author: 乐平 时间: 2022-8-26 10:02
首先，请在帖子里说明你用的软件。当然，读完你的帖子之后能看出你用的是 Gaussian……

其次，计算的初始结构是否已经设置好了点群？

再次，计算是否要保持初始结构的点群？是否存在 Jahn-Teller 效应而破坏对称性？

作者
Author: zgdx 时间: 2022-8-26 12:29
是的，所用软件是gaussian，初始结构点群是Cs没有特意设置,不管是单重态、三重态还是五重态的初始结构的笛卡尔坐标都一样，所以点群都是Cs，优化后的结构不同了，但是点群都是C1。是不是初始结构的点群没有改变才导致这种情况出现的？

作者
Author: zgdx 时间: 2022-8-26 12:32

乐平发表于 2022-8-26 10:02
首先，请在帖子里说明你用的软件。当然，读完你的帖子之后能看出你用的是 Gaussian……

其次，计算的初 ...

我只想这到该初始结构的最稳定时的结构的点群，应该不用特意保持与初始结构一致吧

目前所查阅的文献中没有提到存在 Jahn-Teller 效应

作者
Author: zjxitcc 时间: 2022-8-26 12:52
本帖最后由 zjxitcc 于 2022-8-26 12:53 编辑

zgdx 发表于 2022-8-26 12:32
我只想这到该初始结构的最稳定时的结构的点群，应该不用特意保持与初始结构一致吧

目前所查阅的文献中 ...

（1）你认为它应该是xx点群，这点有实验光谱数据或更高级理论计算结果的支撑否？如果都没有，那没法说你认为的点群是对的。

（2）如果初始结构没有对称性，或计算过程中使用了非限制性方法，或写了nosymm，那结果不符合任何点群对称是正常的。但要注意具体偏离多少，例如优化甲烷分子肉眼看上去是Td点群，但程序没识别出，你量了一下键长，有两根键长相差0.001，那这种情况下应该认为符合Td点群。GaussView点群工具里有个Tolerance可以调节，略微放宽，如果符合目标点群就没问题。

（3）UHF/UDFT计算一般不能保证点群对称性。尤其是要检验波函数稳定性的话，通常带着nosymm，就更没法保证了。说到这一点，还建议你检验一下波函数稳定性，如果不稳定，频率分析就没有意义了。

（4）你要想什么点群，最好提供的初始结构应严格满足该点群。

（5）要严格限制结构符合某点群，需要用内坐标书写输入文件，将等价的坐标用同一个变量表示，模板文件见http://qchem.pw/data/coord/coord.html
但要注意，点群限制是一种可能使电子能量偏高的限制。如果当前计算级别在nosymm下算得的结果根本不可能符合某个点群，而你使用内坐标强行限制了那个点群，会导致你算出来的电子能量偏高，且可能有虚频。

（6）若你坚信它符合某点群，那么开壳层过渡金属最好的办法使用CASSCF方法开点群对称性做计算。若做这种计算建议使用Molpro/OpenMolcas。

作者
Author: chands 时间: 2022-8-26 12:52

zgdx 发表于 2022-8-26 12:32
我只想这到该初始结构的最稳定时的结构的点群，应该不用特意保持与初始结构一致吧

目前所查阅的文献中 ...

Gaussian默认的点群判断标准比较严，经常把点群判低了，我试过H2O判成Cs的。Shermo判得要准一些，有些高对称性的只能靠眼睛看。

作者
Author: sobereva 时间: 2022-8-26 14:15
对于优化收敛限足够严，肉眼看上去明显是高对称的，而且gview里把默认阈值稍微放宽就能判断出更高阶的点群，就认为是相应的点群
若肉眼看上去都没有更高阶的点群特征，要么本来就是低对称性的，要么当前计算级别存在问题。
始终记得检验虚频情况

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