计算化学公社

 找回密码 Forget password
 注册 Register

如何计算基态是双自由基单重态还是双自由基三重态

查看数: 312 | 评论数: 10 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-11-6 14:08

正文摘要:

本帖最后由 Aocan 于 2025-11-10 15:33 编辑 如题,目标分子不带电荷,是稠环芳烃,但是EPR有信号,并且根据UNO算了双自由基特征y0为0.335。所以我认为这个分子在基态下是开壳状态,但是我想进一步确认此分子在基 ...

回复 Reply

Aocan 发表于 Post on 2025-11-11 11:18:38
wjc404 发表于 2025-11-10 18:05
方便以附件的方式分享分子坐标吗?这样能获得更有针对性的帮助

谢谢回复,我大部分问题已经解决了
wjc404 发表于 Post on 2025-11-10 18:05:46
本帖最后由 wjc404 于 2025-11-10 18:07 编辑

方便以附件的方式分享分子坐标吗?这样能获得更有针对性的帮助
Aocan 发表于 Post on 2025-11-10 16:26:30
zjxitcc 发表于 2025-11-10 14:31
没有gauss-summary这种东西,可能你想说的是GaussView面板上的Results -> Summary,这是给新手做简单计算后 ...

Hartree已更正,刚从log文件读取了SCF DONE后面的值,和gaussview上的值是一样的,感谢指正提问细节
zjxitcc 发表于 Post on 2025-11-10 14:31:46
本帖最后由 zjxitcc 于 2025-11-10 14:35 编辑

没有gauss-summary这种东西,可能你想说的是GaussView面板上的Results -> Summary,这是给新手做简单计算后看的,你做的这种计算已经不是简单计算了,要自己根据Gaussian输出文件内容判断,不能盲目相信GaussView读取log文件的总结,否则GaussView万一读错了,你就受到了误导。

没有Harterr这种东西,可能你想说的是Hartree。

提问时最好先挑两个你认为最重要的核心问题入手,进行准确、详细的描述,如果一上来就提很多问题,有些是次要、细枝末节的,而且描述不严谨,容易干扰阅读和读者回复。待核心问题收到回复、解决问题后,可以再跟帖提问其他小问题。
Aocan 发表于 Post on 2025-11-10 10:57:56
sai77 发表于 2025-11-10 10:35
一般结构优化完都要做一个频率计算,以验证收敛到的是极小点还是过渡态。不是自由基必须,而是基本所有的 ...

OK, 我去试试,感谢
sai77 发表于 Post on 2025-11-10 10:35:51
本帖最后由 sai77 于 2025-11-10 10:36 编辑
Aocan 发表于 2025-11-10 10:22
谢谢,请问为什么要加freq关键词呢,这是开壳单线态构型优化必须的吗

一般结构优化完都要做一个频率计算,以验证收敛到的是极小点还是过渡态。不是自由基必须,而是基本所有的优化都需要做这一步
Aocan 发表于 Post on 2025-11-10 10:22:12
sai77 发表于 2025-11-10 08:45
稠环芳烃体系建议B3LYP优化一定要加色散校正。
第一步关键词是做了单点能计算,然后测试了波函数稳定性 ...

谢谢,请问为什么要加freq关键词呢,这是开壳单线态构型优化必须的吗
sai77 发表于 Post on 2025-11-10 08:45:53
Aocan 发表于 2025-11-7 09:52
好的,谢谢。我依照博文的第三个方法输入了如下关键词进行计算
%nprocshared=16
%mem=10GB

稠环芳烃体系建议B3LYP优化一定要加色散校正。
第一步关键词是做了单点能计算,然后测试了波函数稳定性,结果说明收敛到了稳定的波函数。
然后需要读取该波函数进行结构优化,关键词
  1. %nprocshared=16
  2. %mem=10GB
  3. %oldchk=name.chk
  4. %chk=name_2.chk
  5. #p opt freq UB3LYP/6-31G(d,p) em=gd3bj guess=read geom=allcheck


复制代码
有geom=allcheck下面不需要坐标,注意空行。
这里加了色散校正,跟你第一步的不一样,建议第一步加上色散校正重跑一遍。

优化出的结构再测试一遍波函数稳定性
%nprocshared=16
%mem=10GB
%oldchk=name_2.chk
%chk=name_3.chk
#p UB3LYP/6-31G(d,p) em=gd3bj guess=read geom=allcheck stable







sai77 发表于 Post on 2025-11-6 17:06:46
Singlet里是RB3LYP,那说明是优化的闭壳层而不是双自由基
要优化出自旋极化单重态,参考社长的这篇博文
《谈谈片段组合波函数与自旋极化单重态》
http://sobereva.com/82

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 12:29 , Processed in 0.174828 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list