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

计算化学公社

 找回密码
 现在注册!
查看: 450|回复: 15

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

[复制链接]

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

发表于 2018-6-16 20:12:06 | 显示全部楼层 |阅读模式
我想利用片段组合波函数计算反铁磁性耦合配合物,计算结果两个Fe的Mulliken自旋布居都不高,请教各位老师是因为片段设置的不对吗?麻烦老师帮我看看片段设置对吗。

2Fe_Y122_W_OH_-1_2_r_bs.gjf

6.2 KB, 下载次数: 12

129

帖子

1

威望

1245

eV
积分
1394

Level 4 (黑子)

那个考了三年没考上的

发表于 2018-6-16 20:21:43 | 显示全部楼层
iop(5/13=1,1/8=5)??????????????
http://sobereva.com/331
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

175

帖子

0

威望

1447

eV
积分
1622

Level 5 (御坂)

发表于 2018-6-16 20:38:33 | 显示全部楼层
这关键词过于可怕,我忍不住贴一个sob的另一篇帖子http://sobereva.com/61
最后那张图真实太形象了···

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

 楼主| 发表于 2018-6-16 21:15:23 | 显示全部楼层
不好意思,看了这两篇帖子,iop(5/13=1)不用了,各位老师可以说说怎么设置片段吗?

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-6-17 12:45:28 | 显示全部楼层
根本没必要设那么多片段。
倆金属各设一个片段,其余所有原子都设为另一个片段,总共三个片段足矣。别学Gaussian官方的一些网页和书籍,纯属把简单问题搞复杂。

把垃圾关键词全去掉,把片段按我说的方式设,如果自旋布居还那样,可能当前片段设定还不合理,也可能实际分布就那样(用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!

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

 楼主| 发表于 2018-6-17 12:50:47 | 显示全部楼层
好的,谢谢老师,我马上去试试

93

帖子

0

威望

660

eV
积分
753

Level 4 (黑子)

发表于 2018-6-21 05:31:18 | 显示全部楼层
5/13=1,拖出去打靶,哒哒哒哒哒

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

 楼主| 发表于 2018-6-23 17:15:38 | 显示全部楼层
sobereva 发表于 2018-6-17 12:45
根本没必要设那么多片段。
倆金属各设一个片段,其余所有原子都设为另一个片段,总共三个片段足矣。别学Ga ...

sobereva老师,打扰了,片段算完后的chk文件可以转换成fchk吗?我的转换不了,是我的计算有问题吗?

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-6-24 11:49:58 | 显示全部楼层
wgfy000 发表于 2018-6-23 17:15
sobereva老师,打扰了,片段算完后的chk文件可以转换成fchk吗?我的转换不了,是我的计算有问题吗?

“片段算完后”指代不明,贴关键词

转换时确认此文里提到的formchk转换问题都已经注意了
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379
北京科音自然科学研究中心: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!

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

 楼主| 发表于 2018-7-4 16:40:21 | 显示全部楼层
sobereva 发表于 2018-6-24 11:49
“片段算完后”指代不明,贴关键词

转换时确认此文里提到的formchk转换问题都已经注意了

老师,我的片段组合波函数算完了,可以在保持自旋极化单重态的状态下,算单点能吗?要怎么才能保持自旋极化单重态呢?关键词这样可以吗?   # bp86/6-311+g(d,p) geom=check guess=read (保留电荷、自旋行以上的项)

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-7-4 23:55:42 | 显示全部楼层
wgfy000 发表于 2018-7-4 16:40
老师,我的片段组合波函数算完了,可以在保持自旋极化单重态的状态下,算单点能吗?要怎么才能保持自旋极 ...

guess=read即可
北京科音自然科学研究中心: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!

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

 楼主| 发表于 2018-7-5 14:42:49 | 显示全部楼层

老师,我用 guess=read 算单点能,但是log结果没有保持自旋极化单重态,这是什么原因呢?

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

 楼主| 发表于 2018-7-5 20:47:26 | 显示全部楼层

请问老师,用 guess=read 计算的单点能的两个铁的 spin densities 都只有1.2多 ,没有富集在铁上,但是计算完成的自旋极化单重态两个铁的 spin densities 分别为 3.919132和-3.023736,这是什么原因呢?结果可靠吗?

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-7-6 02:13:29 | 显示全部楼层
wgfy000 发表于 2018-7-5 20:47
请问老师,用 guess=read 计算的单点能的两个铁的 spin densities 都只有1.2多 ,没有富集在铁上,但是计 ...

“计算完成的自旋极化单重态”指什么?产生片段初猜的过程的结果?看这个没用

guess=read完的结果如果你觉得异常,加上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!

9

帖子

0

威望

33

eV
积分
42

Level 2 能力者

 楼主| 发表于 2018-7-6 09:07:57 | 显示全部楼层
sobereva 发表于 2018-7-6 02:13
“计算完成的自旋极化单重态”指什么?产生片段初猜的过程的结果?看这个没用

guess=read完的结果如果 ...

两个铁是高自旋反铁磁耦合,计算完成的自旋极化多重态指用片段组合波函数优化完成后,不同自旋的单电子富集在了两个Fe上,spin densities 分别为 3.919132和-3.023736。但是guess=read算完单点能后,两个铁上的spin densities 不高,这样的单点能计算结果可靠吗?我以为两个Fe会保持单电子富集状态,那我要怎么才能保持高自旋反铁磁的单重态来算单点能或者光谱呢?
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2018-7-18 03:03 , Processed in 0.177834 second(s), 27 queries .

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