| 查看: 1563 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
chembetsey木虫 (小有名气)
|
[求助]
QCISD(T) 频率计算：稳定点有多个虚频问题
|
||
|
在计算CH3O这样一个简单的分子的时候，QCISD（T)的优化成功了，但是频率计算却得出数个很大的虚频。有没有人也碰到同样的事情过？ 一个可能性是，优化虽然成功，但是得到的并不是一个稳定的结构。我用QCISD算基于QCISD（T)优化结构的频率却没有虚频。（当然这样做出来的结果是meanless的）这是否说明QCISD（T)下优化的结构并没有错？ 另，QCISD(T)的结构优化和频率计算都只能是numerical 的，初始结构输入似乎只能是用z-matrix。根据输入的z-matrix的不同，似乎优化出来的结果也略有差异。这是否正常？ |
回复此楼» 猜你喜欢
元素硫的化合价和氧化态
已经有1人回复
-
江苏科技大学2026年招聘启事
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有78人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）
已经有78人回复
-
北京工业大学2026年公开招聘公告
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】Gaussian优化和频率计算结果出现虚频怎么办呢？
已经有12人回复
- 【求助】castep的IR计算输出频率的前三个虚频是什么？
已经有4人回复
- 【原创】GAUSSION计算常见错误及解决方案
已经有32人回复
virtualzx
木虫 (著名写手)
- QC强帖: 2
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★
小红豆: 金币+5, 3ks 2012-11-01 18:10:29
小红豆: 金币+5, 3ks 2012-11-01 18:10:29
|
波函数在交叉点附近不稳定。建议限制波函数对称性，分别指定电子态的对称性为A'和A''进行两次优化。限制对称性后的波函数是稳定的。 最好用有cs可约表示的内坐标，然后冻结所有A''内坐标(计算数值导数时也不考虑这些内坐标)，这样应该可约比较顺畅地得到这两个限制优化的极小点。a'和a''中一个能量更高的应该是鞍点，低的那个是极小点。 在姜泰勒系统使用qcisd和ccsd需要额外小心。ccsd和qcisd的哈密顿矩阵都不是赫米矩阵，所以ccsd和qcisd势能面间的交叉点的拓扑性质都不正确。 |
4楼2012-11-01 11:10:05
virtualzx
木虫 (著名写手)
- QC强帖: 2
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
chembetsey: 金币+5, ★★★很有帮助 2012-11-01 09:05:22
小红豆: 金币+3, 3ks 2012-11-01 18:10:17
感谢参与，应助指数 +1
chembetsey: 金币+5, ★★★很有帮助 2012-11-01 09:05:22
小红豆: 金币+3, 3ks 2012-11-01 18:10:17
|
这个分子的计算并不是很容易!它的基态在c3v下是简并的，因为jahn teller效应而分裂，产生三个对称性破缺的极小点和三个鞍点。 有两种可能性，一个是这里所有c3v结构的点都是圆锥交叉点，导数不连续，而且会有两个负无穷的虚频，如果的jahn teller distortion比较小，计算数值导数时可能会在圆锥交叉附近取了点，那么得到的结果就会不准确，可以观察一下目前的结构，考虑降低数值导数的步长，步长一定要远小于到c3v结构的距离 还有可能就是你可能取到了那三个鞍点中的一个 另外对于这个问题，由于接近简并，最好还是用多组态方法。在有，为什么选择用qcisd而不是ccsd呢？现在很少见到qcisd计算了 |
2楼2012-11-01 07:36:11
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
3楼2012-11-01 09:03:43
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
5楼2012-11-05 00:51:21
