计算化学公社

标题: 关于反铁磁性耦合计算问题 [打印本页]

作者
Author: wangjianmin 时间: 2017-9-6 10:05
标题: 关于反铁磁性耦合计算问题
我想利用片段组合波函数计算反铁磁性耦合配合物，但是在l122出错，请教各位老师是因为片段设置的不对吗？

作者
Author: sobereva 时间: 2017-9-6 14:06
对于产生片段初猜的目的，没必要让每个片段都做SCF到收敛，写guess(fragment,only)就够了。之后读取这个chk里的片段初猜做反铁磁性耦合体系基三

作者
Author: liyuanhe211 时间: 2017-9-6 15:00
你的配合物片段除Co外都是闭壳层的（水、羧酸根负离子、吡啶），Co 是 Quartet，为什么整体是doublet呢？
整体是doublet不用声明U
另外这种情况叫反铁磁性耦合吗？

(, 下载次数 Times of downloads: 42)

(, 下载次数 Times of downloads: 37)

作者
Author: wangjianmin 时间: 2017-9-6 15:09

liyuanhe211 发表于 2017-9-6 15:00
你的配合物片段除Co外都是闭壳层的（水、羧酸根负离子、吡啶），Co 是 Quartet，为什么整体是doublet呢？
...

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

作者
Author: wangjianmin 时间: 2017-9-6 15:13

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

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

作者
Author: sobereva 时间: 2017-9-6 18:01

wangjianmin 发表于 2017-9-6 15:13
老师，我写guess(fragment,only)还是出现同样的错误

片段的自旋多重度设置和整体根本匹配不上，先改了这个。光是片段1有三个单电子，其它片段都是闭壳层的，而整个体系一个单电子，明显矛盾。

作者
Author: wangjianmin 时间: 2017-9-6 19:07

liyuanhe211 发表于 2017-9-6 15:00
你的配合物片段除Co外都是闭壳层的（水、羧酸根负离子、吡啶），Co 是 Quartet，为什么整体是doublet呢？
...

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

作者
Author: wangjianmin 时间: 2017-9-6 19:08

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

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

作者
Author: wangjianmin 时间: 2017-9-9 09:50

wangjianmin 发表于 2017-9-6 19:08
谢谢老师，我发现我的问题了

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

作者
Author: sobereva 时间: 2017-9-9 12:46

wangjianmin 发表于 2017-9-9 09:50
老师，我还有两个问题想请教您：我设置了好几种片段的组合方式，然后进行了初猜，但是用初猜结构计算后的 ...

用低自旋闭壳层态，直接做个stable，立刻就知道当前这个波函数是否是稳定的，若不稳定明显暗示这是开壳层态

作者
Author: wangjianmin 时间: 2017-9-12 08:42

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

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

作者
Author: sobereva 时间: 2017-9-12 12:57

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

看不懂你的表述。要算基态，既可以尝试不同初猜然后检查stable情况，或者直接用stable=opt让程序自发优化到稳定波函数

作者
Author: wangjianmin 时间: 2017-9-12 15:37

sobereva 发表于 2017-9-12 12:57
看不懂你的表述。要算基态，既可以尝试不同初猜然后检查stable情况，或者直接用stable=opt让程序自发优 ...

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

作者
Author: sobereva 时间: 2017-9-12 16:17

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

自旋多重度设1

作者
Author: wangjianmin 时间: 2017-9-13 15:12

sobereva 发表于 2017-9-12 16:17
自旋多重度设1

特别感谢老师的耐心回答

作者
Author: lmy 时间: 2018-11-4 10:38

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

您好，麻烦问下您说的计算磁耦合常数公式的参考文献可以分享下吗，先谢谢啦

作者
Author: wangjianmin 时间: 2018-11-4 17:31
看一下这个

作者
Author: lmy 时间: 2018-11-8 17:00

wangjianmin 发表于 2018-11-4 17:31
看一下这个

不好意思，这几天考试，没来得及查看。谢谢您啦。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)