计算化学公社

标题: xtb单点计算体系是否不一定真实存在? [打印本页]

作者
Author:
WaterMirror    时间: 2022-2-28 23:51
标题: xtb单点计算体系是否不一定真实存在?
尊敬的各位老师:
  我用GV6的模板建立一个环己烷,将其转换成.xyz文件,用xtb cyclohexane.xyz --chrg +1,计算环己烷阳离子在环己烷几何下的单点能,xTB正常结束,输出charge 1; spin 0.5。也即在计算阳离子时xTB自动判断了此时体系是二重态,S=0.5。接下来用xtb cyclohexane.xyz --uhf 1计算,xTB也正常结束,输出charge 0; spin 0.5这样一个不可能的状态。同时给出单点能,HOMO-LUMO GAP; 接下来用 xtb cyclohexane.xyz --uhf 1 --chrg 0,计算依然正常输出charge 0; spin 0.5这样一个不可能的状态。同时给出单点能,HOMO-LUMO GAP。这是为什么呀?这些状态应该不可能存在的,这是怎么计算出来的呢?




作者
Author:
exity    时间: 2022-3-1 09:37
没opt直接上的吗?
作者
Author:
WaterMirror    时间: 2022-3-1 09:42
exity 发表于 2022-3-1 09:37
没opt直接上的吗?

没有--opt 直接单点
作者
Author:
RAL    时间: 2022-3-1 09:46
xtb判断电子数的时候,chrg的优先级似乎比uhf高,可以看看计算结果中的轨道占据数,--uhf 1没起作用
作者
Author:
WaterMirror    时间: 2022-3-1 10:36
RAL 发表于 2022-3-1 09:46
xtb判断电子数的时候,chrg的优先级似乎比uhf高,可以看看计算结果中的轨道占据数,--uhf 1没起作用

可是输出结果显示spin=0.5呀
作者
Author:
RAL    时间: 2022-3-1 10:38
WaterMirror 发表于 2022-3-1 10:36
可是输出结果显示spin=0.5呀

那只是它把你输入的参数打印出来了而已
作者
Author:
sobereva    时间: 2022-3-2 03:38
实际该是什么状态就设什么,不报错不代表结果有意义




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3