请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 950|回复: 17

[Gaussian/gview] 关于反铁磁性耦合计算问题

[复制链接]

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

发表于 2017-9-6 10:05:41 | 显示全部楼层 |阅读模式
我想利用片段组合波函数计算反铁磁性耦合配合物,但是在l122出错,请教各位老师是因为片段设置的不对吗?

040630A-1-s.out

147.11 KB, 下载次数: 20

1万

帖子

25

威望

1万

eV
积分
35736

管理员

公社社长

发表于 2017-9-6 14:06:20 | 显示全部楼层
对于产生片段初猜的目的,没必要让每个片段都做SCF到收敛,写guess(fragment,only)就够了。之后读取这个chk里的片段初猜做反铁磁性耦合体系基三
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

2305

帖子

23

威望

7757

eV
积分
10522

Level 6 (一方通行)

发表于 2017-9-6 15:00:54 | 显示全部楼层
你的配合物片段除Co外都是闭壳层的(水、羧酸根负离子、吡啶),Co 是 Quartet,为什么整体是doublet呢?
整体是doublet不用声明U
另外这种情况叫反铁磁性耦合吗?

snap0001189.jpg

snap0001183.png




51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-6 15:09:18 | 显示全部楼层
liyuanhe211 发表于 2017-9-6 15:00
你的配合物片段除Co外都是闭壳层的(水、羧酸根负离子、吡啶),Co 是 Quartet,为什么整体是doublet呢?
...

老师,那个是吡啶氮氧自由基,有一个单电子,所以整体我就设置成了doublet

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-6 15:13:25 | 显示全部楼层
sobereva 发表于 2017-9-6 14:06
对于产生片段初猜的目的,没必要让每个片段都做SCF到收敛,写guess(fragment,only)就够了。之后读取这个chk ...

老师,我写guess(fragment,only)还是出现同样的错误

1万

帖子

25

威望

1万

eV
积分
35736

管理员

公社社长

发表于 2017-9-6 18:01:45 | 显示全部楼层
wangjianmin 发表于 2017-9-6 15:13
老师,我写guess(fragment,only)还是出现同样的错误


片段的自旋多重度设置和整体根本匹配不上,先改了这个。光是片段1有三个单电子,其它片段都是闭壳层的,而整个体系一个单电子,明显矛盾。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-6 19:07:43 | 显示全部楼层
liyuanhe211 发表于 2017-9-6 15:00
你的配合物片段除Co外都是闭壳层的(水、羧酸根负离子、吡啶),Co 是 Quartet,为什么整体是doublet呢?
...

谢谢老师,我发现我的问题了

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-6 19:08:06 | 显示全部楼层
sobereva 发表于 2017-9-6 18:01
片段的自旋多重度设置和整体根本匹配不上,先改了这个。光是片段1有三个单电子,其它片段都是闭壳层的 ...

谢谢老师,我发现我的问题了

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-9 09:50:00 | 显示全部楼层
wangjianmin 发表于 2017-9-6 19:08
谢谢老师,我发现我的问题了

老师,我还有两个问题想请教您:我设置了好几种片段的组合方式,然后进行了初猜,但是用初猜结构计算后的结果都不收敛,而且最后的.out文件里的每个片段的charge都是0,Unpaired Spin都是Alpha,不是反铁磁性的,而且这个配合物实验测得的就是反铁磁性耦合的。
2.我是不是可以直接设置自旋多重度为低自旋态,进行优化,不去设置片段呢?

1万

帖子

25

威望

1万

eV
积分
35736

管理员

公社社长

发表于 2017-9-9 12:46:34 | 显示全部楼层
wangjianmin 发表于 2017-9-9 09:50
老师,我还有两个问题想请教您:我设置了好几种片段的组合方式,然后进行了初猜,但是用初猜结构计算后的 ...


用低自旋闭壳层态,直接做个stable,立刻就知道当前这个波函数是否是稳定的,若不稳定明显暗示这是开壳层态
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-12 08:42:40 | 显示全部楼层
sobereva 发表于 2017-9-9 12:46
用低自旋闭壳层态,直接做个stable,立刻就知道当前这个波函数是否是稳定的,若不稳定明显暗示这是开壳 ...

老师,我按照您说的算了,得到的结果是不稳定的,我的配合物体系是一个中心离子,和两个吡啶氮氧自由基配位,我想计算它在吡啶氮氧自由基和中心离子是反铁磁性耦合时的单点,是不是直接把整体的自旋多重度改了,可以不设置片段呢?我一直搞不清楚这个问题

1万

帖子

25

威望

1万

eV
积分
35736

管理员

公社社长

发表于 2017-9-12 12:57:20 | 显示全部楼层
wangjianmin 发表于 2017-9-12 08:42
老师,我按照您说的算了,得到的结果是不稳定的,我的配合物体系是一个中心离子,和两个吡啶氮氧自由基配 ...


看不懂你的表述。要算基态,既可以尝试不同初猜然后检查stable情况,或者直接用stable=opt让程序自发优化到稳定波函数
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-12 15:37:22 | 显示全部楼层
sobereva 发表于 2017-9-12 12:57
看不懂你的表述。要算基态,既可以尝试不同初猜然后检查stable情况,或者直接用stable=opt让程序自发优 ...

老师,我已经算过基态了,检查了stable,显示是稳定的,但是我想计算磁耦合常数,看文献说可以再计算它的低自旋态能量,然后利用公式得到磁耦合常数,可是我不知道低自旋态的能量应该怎么得到

1万

帖子

25

威望

1万

eV
积分
35736

管理员

公社社长

发表于 2017-9-12 16:17:43 | 显示全部楼层
wangjianmin 发表于 2017-9-12 15:37
老师,我已经算过基态了,检查了stable,显示是稳定的,但是我想计算磁耦合常数,看文献说可以再计算它的 ...

自旋多重度设1
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

51

帖子

0

威望

526

eV
积分
577

Level 4 (黑子)

 楼主| 发表于 2017-9-13 15:12:45 | 显示全部楼层

特别感谢老师的耐心回答
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-11-17 05:14 , Processed in 0.126061 second(s), 26 queries .

快速回复 返回顶部 返回列表