24小时热门版块排行榜    

CyRhmU.jpeg
查看: 975  |  回复: 1

hustyh0801

铁杆木虫 (著名写手)

[求助] 结合能

我想算一个分子(实际上是一个盐,但我却以小分子形式来处理)的结合能,本来Gaussian对于盐类这种具有无限重复晶胞结构的物质的结合能计算在理论上是否正确?另外就是阴阳离子的结合方式(配位方式也许这里说配位不准确)不同会影响到结合能,比如说下面的两幅图,这是一个LiPF6的小分子图其中Li可以和两个F配位,也可以和三个F配位,但我现在的问题是因为Li和3个F作用是一种更加稳定的状态(文献上所述),我优化得到的结果也只得到3配位的。当然文献上指出该分子应该有三个势能低点,分别是一配位 二配位 三配位 但我前两个怎么也优化不出来。我也考虑到了初始结构的问题但不管我怎么调整比如我尽量把结构调成接近二齿的但是优化后还是往三齿的方向走,我也降低了优化的步长(当然这样不一定有道理)但依旧没有作用。我想请问大家这样的情况应该怎么设置才能让优化过程往我希望的方向走,我也是初学者可能很多的表述都存在问题希望大家能够不吝指教!

回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

【答案】应助回帖

hustyh0801(金币+10): 首先谢谢你的帮主,低精度的PM6的确可以优化二齿的但是我想采用DFT的方法就不行,不知道该如何作限制性的优化。 2011-05-18 23:13:13
引用回帖:
Originally posted by hustyh0801 at 2011-04-30 20:47:14:
我想算一个分子(实际上是一个盐,但我却以小分子形式来处理)的结合能,本来Gaussian对于盐类这种具有无限重复晶胞结构的物质的结合能计算在理论上是否正确?另外就是阴阳离子的结合方式(配位方式也许这里说配位 ...

高斯计算这个是否准确我不清楚,以前也没有研究过这种问题,在这里就不多说了。
至于2齿的我用高斯09低精度的pm6方法帮楼主优化出来了,1齿的优化过程中变为2齿的了。我猜如果要优化1齿的话可能需要使用限制性优化才有可能得到1齿结构,因为1齿结构应该是非常不稳定的。
2楼2011-05-18 22:00:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hustyh0801 的主题更新
信息提示
请填处理意见