请问如何用Gaussian 09用DFT方法进行开壳层的单线态计算(open-shell singlet state),我使用了UB3LYP,但是对于单线态,加U不加U的结果是一样,加了没有作用,请问这种情况该如何设置关键词? 返回小木虫查看更多
不一样吧,至少能量和轨道是不一样的吧
udft,然后交换alpha的homo和lumo
对于开壳和闭壳计算,是有可能结果一样的,不过两种计算的基本原理不一样。 闭壳是假设每个轨道上的两个电子(alpha和beta)除了自旋不一样,其波函数是一样的,计算过程中只需要算一次(alpha),beta的波函数等于alpha(这个可以在计算完结果中看轨道分布得到验证)。 但如果是开壳计算,alpha和beta是分开算的,同一轨道的两个电子的波函数可能不相等。不过开壳计算一般是针对磁性材料 对于你的情况,存在这种可能,如果你计算的分子是没有磁性,有可能几遍alpha和beta分开算,你得到波函数也可能是一样的。 如果你想验证这一点,你可以在计算中前面加p,后面加gfprint pop=full,你可以在out文件中看到分子轨道的参数。如果闭壳应该只有Molecular Orbital Coefficients这个数据块, 如果是开壳则有Alpha Molecular Orbital Coefficients 和Beta Molecular Orbital Coefficients,你的情况应该是这两个数据块是一样的(或者Density Matrix数据块也是如此),
不一样吧,至少能量和轨道是不一样的吧
我比较过了,对于单线态,B3LYP和UB3LYP所有结果都是一样的
udft,然后交换alpha的homo和lumo
对于开壳和闭壳计算,是有可能结果一样的,不过两种计算的基本原理不一样。
闭壳是假设每个轨道上的两个电子(alpha和beta)除了自旋不一样,其波函数是一样的,计算过程中只需要算一次(alpha),beta的波函数等于alpha(这个可以在计算完结果中看轨道分布得到验证)。
但如果是开壳计算,alpha和beta是分开算的,同一轨道的两个电子的波函数可能不相等。不过开壳计算一般是针对磁性材料
对于你的情况,存在这种可能,如果你计算的分子是没有磁性,有可能几遍alpha和beta分开算,你得到波函数也可能是一样的。
如果你想验证这一点,你可以在计算中前面加p,后面加gfprint pop=full,你可以在out文件中看到分子轨道的参数。如果闭壳应该只有Molecular Orbital Coefficients这个数据块, 如果是开壳则有Alpha Molecular Orbital Coefficients 和Beta Molecular Orbital Coefficients,你的情况应该是这两个数据块是一样的(或者Density Matrix数据块也是如此),
请问如何交换alpha的homo和lumo轨道呢?
查guess=alter
你好,我想一下,不论用哪个方法比如M06,我想计算单线态开壳层只要在M06前面加上U就可以了吗,其他关键词都不用改是吗,本人初学高斯还不太懂,感谢!