| 查看: 1226 | 回复: 7 | |||
chembetsey木虫 (小有名气)
|
[求助]
QCISD(T) 频率计算：稳定点有多个虚频问题
|
|
在计算CH3O这样一个简单的分子的时候，QCISD（T)的优化成功了，但是频率计算却得出数个很大的虚频。有没有人也碰到同样的事情过？ 一个可能性是，优化虽然成功，但是得到的并不是一个稳定的结构。我用QCISD算基于QCISD（T)优化结构的频率却没有虚频。（当然这样做出来的结果是meanless的）这是否说明QCISD（T)下优化的结构并没有错？ 另，QCISD(T)的结构优化和频率计算都只能是numerical 的，初始结构输入似乎只能是用z-matrix。根据输入的z-matrix的不同，似乎优化出来的结果也略有差异。这是否正常？ |
回复此楼» 猜你喜欢
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有191人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
26博士申请
已经有0人回复
-
调剂招生
已经有7人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】Gaussian优化和频率计算结果出现虚频怎么办呢？
已经有12人回复
- 【求助】castep的IR计算输出频率的前三个虚频是什么？
已经有4人回复
- 【原创】GAUSSION计算常见错误及解决方案
已经有32人回复
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
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
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
5楼2012-11-05 00:51:21
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
6楼2012-11-05 01:09:16
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13046.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
7楼2012-11-05 02:53:43
virtualzx
木虫 (著名写手)
- QC强帖: 2
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
8楼2012-11-05 23:48:26
