| 查看: 770 | 回复: 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 ] |
回复此楼» 猜你喜欢
求助火焰封管的时候管子炸了
已经有1人回复
-
PbS量子点如何确定浓度
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有259人回复
-
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）
已经有22人回复
-
PbS量子点紫外吸收
已经有0人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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 (初中生)
- 金币: 12837.4
- 散金: 27307
- 红花: 300
- 帖子: 4667
- 在线: 925.7小时
- 虫号: 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