| 查看: 1756 | 回复: 11 | ||||
| 当前主题已经存档。 | ||||
[交流]
【求助】轨道对称性部分出现很多问号！
|
||||
|
请教一个关于对称性的问题： 我用Gaussian优化富勒烯的一个结构（其最高对称性为D3h，但软件刚生成时的初始对称性为C1）。文献报道该结构的五重态能量最低且电子态为5A1'，于是我把自旋多重度设置为5，对其进行几何优化并计算频率以及测试波函数稳定性（不限制其对称性，仍为C1，命令行未加nosymm）。计算完成后，得到其电子态为5A，对称性仍为C1，没有虚频且波函数也是稳定的。但用Gview看其对称性时发现：当前的点群为C1,但更高的点群的框里显示为D3h。于是我用Gview将其对称性提高到D3h，同样进行五重态的优化以及频率和波函数稳定性计算。结果这次算出的能量仅比上次的高0.003kcal/mol；高斯输出中，初猜的电子态为5A1'，但下面的轨道对称性（即Orbital symmetries:下面的轨道对称性输出部分）里却有很多问号（？），接着是：Unable to determine electronic state: an orbital has unidentified symmetry；这次仍没有虚频且波函数也稳定，高斯输出中对称性为D3h。 请问接下来我应该如何处理呢？该结构的对称性到底是什么呢？D3h吗？那它的电子态又是什么呢？恳请大家帮忙！非常感谢！！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 理论化学基本知识 |
» 猜你喜欢
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有295人回复
-
苏州科技大学 物理学专业 接收调剂
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有20人回复
-
河南农业大学，化学学硕，三轮调剂
已经有1人回复
-
招材料相关专业调剂生
已经有10人回复
冬天里的骄阳
铁杆木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 7179.7
- 散金: 10
- 帖子: 387
- 在线: 318小时
- 虫号: 439463
- 注册: 2007-10-13
- 性别: GG
- 专业: 理论和计算化学
2楼2010-01-03 17:31:05
llhxaustin
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4494.5
- 红花: 7
- 帖子: 2749
- 在线: 174小时
- 虫号: 143679
- 注册: 2005-12-22
- 性别: GG
- 专业: 理论和计算化学
3楼2010-01-03 18:31:00
wolfzhong
木虫 (正式写手)
- 应助: 23 (小学生)
- 金币: 3429.5
- 散金: 274
- 红花: 9
- 帖子: 986
- 在线: 751.4小时
- 虫号: 619989
- 注册: 2008-10-08
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★ ★
fly9810(金币+5,VIP+0):多谢！ 1-3 20:13
yjcmwgk(金币+2,VIP+0): 1-8 13:51
fly9810(金币+5,VIP+0):多谢！ 1-3 20:13
yjcmwgk(金币+2,VIP+0): 1-8 13:51
|
电子态出现问号的情况非常多尤其是对含有3重轴的体系。 考虑这个问题基本出发点有以下2种： 1.你的体系是否真的具有这个对称性？ 例如：Oh对称的某个分子因为John-Tayler畸变得到最终的D4h对称性。 2.体系确实有这个构型，而计算过程中出现“对称性破却”。 有以下解决方式： 1） 使用小basis sets，看是否破却，如果没有，读小基组的轨道；原因：基组越大，轨道越多，轨道分部越紧密，导致破却的可能性越大。 2） 强制保持轨道及密度矩阵的对称性，dsymm； 3） 改变初始结构，键长越长，energy gaps越小，所以保持对称性的基础上缩短键长。 4）其他？ 对于你的体系，可以用1）和2）。 |
4楼2010-01-03 18:49:01
5楼2010-01-03 20:12:10
6楼2010-01-03 20:20:28
544577578
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3593.7
- 散金: 22
- 红花: 1
- 帖子: 136
- 在线: 52.7小时
- 虫号: 819733
- 注册: 2009-07-31
- 性别: MM
- 专业: 理论和计算化学
7楼2010-01-03 20:33:36
wolfzhong
木虫 (正式写手)
- 应助: 23 (小学生)
- 金币: 3429.5
- 散金: 274
- 红花: 9
- 帖子: 986
- 在线: 751.4小时
- 虫号: 619989
- 注册: 2008-10-08
- 专业: 理论和计算化学
8楼2010-01-03 22:25:03
wolfzhong
木虫 (正式写手)
- 应助: 23 (小学生)
- 金币: 3429.5
- 散金: 274
- 红花: 9
- 帖子: 986
- 在线: 751.4小时
- 虫号: 619989
- 注册: 2008-10-08
- 专业: 理论和计算化学
9楼2010-01-03 22:30:01
测试结果！
|
按照您所说的方法，我重新进行了一些测试计算，结果如下： 1）正如您所说，“对称性破却”确实是由所用基组增大引起，我用6-31G*和6-311G* 优化该结构时，都能得到电子态为5A1'，且轨道对称性部分不出现问号，但当提高到6-311+G*时，又会出现问号，且Unable to determine electronic state: an orbital has unidentified symmetry。 2）接下来我按您说的，我先用6-31G* 计算该结构的单点，然后用6-311+G* 读取以上的chk，继续优化（命令：geom＝check guess=read opt），命令中不论加不加“scf=dsymm” ，最后结果均为：电子态：5A1；对称性（高斯输出中的）为C2v；轨道对称性部分不出现问号。 3）然后，我直接用6-31G* 优化该结构，接着用6-311+G* 读取其chk，继续优化（命令：geom＝check guess=read opt scf=dsymm），或者直接用6-311+G* 优化该结构（opt scf=dsymm ）。最后结果均为：电子态：5A1'；对称性（高斯输出中的）为D3h；但轨道对称性部分都出现了问号。 4）所有上述用6-311+G* 基组优化得到的能量值以及上面发的帖子里的能量值都很接近，相差不到0.038kcal/mol，频率计算正在进行中，估计不会有虚频。 现在的问题是：我能说该结构的对称性就是D3h吗？尽管轨道对称性部分的问号仍然没有消除。 |
10楼2010-01-05 11:24:42
5