| 查看: 827 | 回复: 5 | ||||
| 当前主题已经存档。 | ||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 fly9810 的 15 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会 | ||||
[交流]
【求助】电子态问题
|
||||
我用DFT方法优化富勒烯C28的经典异构体,电荷和自旋多重度均设为: 0 1,发现能量最低结构的初始对称性为Td(Gview下得到的结果),不施加对称性限制,只在输入文件中 加入 "symmetry=loose ", 但DFT优化后降为C1(Gview下的结果,近似更高的对称性为S4),但文献报道的最优结构是Td对称的(跟我得到的最优结构相同),而且还说该结构含有4个未成对电子,电子态为5A2.于是我首先将该异构体的初始对称性在Gview中(point group 项)改为Td, 将自旋多重度改为"5",重新优化之,结果发现此次得到的能量值确实比原来的能量低(约十几kcal/mol),但电子态却是5A1,而非文献报道的5A2,请问我该如何处理呢?非常感谢!!! [ Last edited by zeoliters on 2009-11-13 at 18:36 ] |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有2人回复
-
085404（计算机技术）293分 希望调剂到材料相关专业
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有123人回复
-
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gwdavid(金币+5,VIP+0):谢谢 9-29 15:57
fly9810(金币+5,VIP+0): 10-6 19:08
gwdavid(金币+5,VIP+0):谢谢 9-29 15:57
fly9810(金币+5,VIP+0): 10-6 19:08
|
你这个问题没有说清楚 能量高低的判断是根据gaussian计算的结果来判断了，怎么是GaussView下得道的结果呢？ 一般情况下，Gaussian计算是不会改变分子的对称性的，除非你在关键词中设置了允许改变的关键词。所以你的分子不是优化后降为C1，而是一开始就是C1。 根据你的描述，我怀疑你这个体系可能是一个双自由基或者四自由基的开壳层单重态体系，不一定就是传说中的五重态，你试试计算中将电荷与自旋多重度设为0 1，但使用开壳层来计算，也就是UDFT（比如ub3lyp），此时需要加关键词nosymm。有时候可能能量更低。 Gaussian好像没法指定电子态，只能靠软件自己去判断。 想要精确判断体系的电子态需要molpro等更专业的软件用多参考组态或者活性空间方法来处理。 |
2楼2009-09-29 12:35:59
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13087.2
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★ ★ ★ ★ ★ ★
yjcmwgk(金币+2,VIP+0):感谢您耐心的指导 9-30 08:53
yjcmwgk(金币+1,VIP+0):追加1分 9-30 08:54
fly9810(金币+5,VIP+0):感谢！！！ 10-7 15:31
yjcmwgk(金币+2,VIP+0):感谢您耐心的指导 9-30 08:53
yjcmwgk(金币+1,VIP+0):追加1分 9-30 08:54
fly9810(金币+5,VIP+0):感谢！！！ 10-7 15:31
|
关于gaussian指定电子态： 如果组态能产生对称性唯一的电子态，比如(t)6(a1)1(a2)1组态，产生的一定是单重或三重的a2态。如果gaussian的初始猜测不符合要求，可以交换轨道。 如果存在简并的开壳层轨道，并且导致组态能产生多个不同对称性的电子态，除了极个别情况，一般没办法。例如，(t)1(a1)1会产生单重或三重的t态，但像(t)1(t)1就不好判断了，只能用其他方法才能指定对称性，如casscf，eom-ip/ea-cc，多重态dft等。特殊的例子如氧原子2p4，这个组态产生3个态1S，1D，3P。如果指定为三重态，则为3P。若指定为单态，通常是1S和1D中能量较低的态，但也不一定，特别是对复杂体系。 |
3楼2009-09-30 01:04:45
★ ★
fly9810(金币+2,VIP+0): 10-6 10:02
fly9810(金币+2,VIP+0): 10-6 10:02
|
Gaussian里现在有eom-ip吗? http://gaussian.com/g_tech/g_ur/k_eom.htm 没找到 |
4楼2009-09-30 03:04:06
|
这位朋友说的没错，该结构的初始对称性的确是C1，但Approximate higher-order point group 中显示的对称性为Td(不好意思，刚开始开始没说清楚）。按照文献，第二次计算中在用半经验方法优化前,先在Gview中(point group 项)将其对称性改为Td，将自旋多重度改为"5",重新优化之,计算结果如上文所述。 接下来我按你说的重新算了一下： %rwf=1,2000mb,2,2000mb,3,2000mb,4,2000mb %nproc=2 %mem=1000mb %nosave %chk=28-00-01 # ub3lyp/6-31g* opt nosymm freq 28-00-01 0 1 坐标…… 发现，这次得到的能量值较第一次计算结果还略高一点。对称性由开始的Td降为C1，没有输出电子态的信息。频率计算得到2个虚频。所以，可能不是阁下说的这种情况吧。不过，Thank you all the same!!! |
5楼2009-10-06 19:07:29
6楼2009-10-06 19:10:59
5