gaussian 和gaussview_谈谈Gaussian软件中的guess=mix

本文主要讨论了在Gaussian软件中,guess=mix关键词的使用场景和作用,特别是在涉及自由基的反应计算时。guess=mix主要用于解决自旋极化单重态计算中的自旋污染问题,通过混合HOMO与LUMO轨道来寻找更低的能量解。文章通过实例展示了guess=mix的工作原理,并提醒用户在使用时结合stable=opt以确保波函数的稳定性。对于非自旋极化单重态或二重态、三重态计算,无需使用guess=mix。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

笔者经常碰到小伙伴在用Gaussian软件计算涉及自由基的反应时,不清楚何时该加关键词guess=mix,何时不该加;也可能会有师兄/老师这样告诉新手:碰到自由基一律用guess(mix,always)。前者可能量化基础不扎实,碰到这类问题不懂;后者则可能缺乏实际计算经验。趁假期有空,正好写上一篇,详细解释一下。当然,笔者写的绝对不是标准答案,只能力求合理性和正确性,仅供对这个问题不清楚的小伙伴们参考。

为方便起见,本文仅讨论RHF和UHF,所有内容对DFT同样适用。对UHF计算不熟悉的新手可以阅读本公众号发过的软件教程《用Gaussian做UHF计算》。本文中所涉及计算皆使用G16 A.03。

首先给出简单的结论:(1)guess=mix只在自旋极化单重态(即使用UHF方法在单重态下做计算,发现有严重自旋污染)时需要考虑加,其他情况(如二重态、三重态)等无需考虑这个问题;(2)always表示在结构优化的每一步中都执行guess=mix。顾名思义,这只在结构优化中可能有用,而在单点计算中无需加、加了也没用。

接下来我们一一解释。关于自旋极化单重态和guess=mix的含义,Sob老师的博文《谈谈片段组合波函数与自旋极化单重态》(http://sobereva.com/82)讲得十分详细,强烈推荐经常做此类计算、但又还没看过博文的小伙伴仔细阅读。这里再举一个简单的单点计算实例:在UHF/STO-3G水平下计算键长为2.0 Å的H2分子,对于如下三种关键词写法:

(1) #p RHF/STO-3G nosymm(2) #p UHF/STO-3G nosymm(3) #p UHF/STO-3G nosymm guess=mix
其中,写法(1)和(2)得到的电子能量是一样的,均为−0.783792 a.u., = 0,SCF迭代圈数均为1圈。打开(
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值