计算化学公社

标题: 关于构建分子势能函数（正氢orthohydrogen） [打印本页]

作者
Author: animao 时间: 2023-4-24 14:00
标题: 关于构建分子势能函数（正氢orthohydrogen）
求助怎么用gaussian或molpro势能面扫描构建正氢的势能函数啊

找了很多文学，是区分正氢和仲氢与别的分子的intermolecular相互作用，但是没有正氢-正氢这种研究啊。也没有看明白具体怎么设置正氢或仲氢的，是用特定的基组吗？还是调量子数之类的？
gaussian只能设置自旋多重度，但这是电子的自旋，不是核自旋。

作者
Author: wzkchem5 时间: 2023-4-25 04:14
正氢和仲氢主要是核自旋和转动量子数的区别，但是后者对分子间相互作用的影响远大于前者，所以可以只考虑转动量子数的区别。对于不太精确的计算可以把氢分子近似成刚性转子，如果刚性转子对你的精度要求来说不够，再讨论怎么处理的问题

作者
Author: animao 时间: 2023-4-25 15:59

wzkchem5 发表于 2023-4-25 04:14
正氢和仲氢主要是核自旋和转动量子数的区别，但是后者对分子间相互作用的影响远大于前者，所以可以只考虑转 ...

多谢您的回答！很专业。请问这个转动量子数在molpro软件可以定义的吗？是“para-H2(j1 = 0) ，ortho-H2(j1 = 1) case, the wave function can be (j1 = 1, k = 0) or (j1 = 1, k = ±1)”这样吗？

作者
Author: animao 时间: 2023-4-26 06:49

wzkchem5 发表于 2023-4-25 04:14
正氢和仲氢主要是核自旋和转动量子数的区别，但是后者对分子间相互作用的影响远大于前者，所以可以只考虑转 ...

转动量子数是磁量子数吗？也就是说正氢角量子数取1时，磁量子数可以取取-1 0 1？
但我理解这些不都是核外电子的状态描述吗，正氢仲氢的区别应该不是核外电子的不同

作者
Author: wzkchem5 时间: 2023-4-27 01:00

animao 发表于 2023-4-25 23:49
转动量子数是磁量子数吗？也就是说正氢角量子数取1时，磁量子数可以取取-1 0 1？
但我理解这些不都是核 ...

我说的是分子的转动能级，不是自旋能级。所以不是磁量子数

作者
Author: animao 时间: 2023-4-30 09:44

wzkchem5 发表于 2023-4-27 01:00
我说的是分子的转动能级，不是自旋能级。所以不是磁量子数

请问这个转动量子数可以在molpro或gaussian中设置吗？

作者
Author: wzkchem5 时间: 2023-4-30 16:09

animao 发表于 2023-4-30 02:44
请问这个转动量子数可以在molpro或gaussian中设置吗？

当使用刚性转子近似的时候，转动量子数一般不是在量化程序里设的，而是你做完计算以后得到转动常数，再用转动量子数套公式自己算出来转动能级。
当不做刚性转子近似的时候，可能需要一开始就把转动量子数输入到程序里，但高斯应该做不了超出刚性转子近似的转动能级计算，molpro我不确定但是估计也不行

作者
Author: animao 时间: 2023-5-1 11:49

wzkchem5 发表于 2023-4-30 16:09
当使用刚性转子近似的时候，转动量子数一般不是在量化程序里设的，而是你做完计算以后得到转动常数，再用 ...

好的我按您说的去查一下。因为我是想区分仲氢和正氢的势能函数，想分别做势能面扫描，这样的话我是不是必须先区分二者的转动量子数设置？

作者
Author: wzkchem5 时间: 2023-5-1 14:24

animao 发表于 2023-5-1 04:49
好的我按您说的去查一下。因为我是想区分仲氢和正氢的势能函数，想分别做势能面扫描，这样的话我是不是必 ...

对，如果程序支持就在程序里设，如果程序不支持就自己按刚性转子的公式，手算不同转动量子数的转动能级

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