24小时热门版块排行榜    

查看: 1379  |  回复: 11

lhjxgy

银虫 (正式写手)

[求助] 计算反铁磁

我是学配位化学的,近来合成了五个锰的配合物,有双核,四核,六核,统一锰链,交替锰链SBU的化合物,做了磁性,但是我想通过量化计算计算它们的理论J值,从而比较一下它们磁性传递情况,我应该从哪入手啊?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingznpmg

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2012-08-28 14:39:42
对于线性的构型了解的不多,但对反铁磁的计算,可以参考G09手册以下部分:
Antiferromagnetic coupling. Here is an example of a fragment guess job. The first step generates guesses for each one of the fragments, and
then combines them together into a guess for the SCF on the full molecule. The second step reads in that guess for a subsequent calculation. The
charge and multiplicity line has the overall molecular charge and multiplicity first, followed by those for each fragment.
This example is Fe2S2 + 4 S-R ligands, where R is a phenyl group. It is a singlet with an overall charge of -2. The guess does each bare sulfur as
S(2-) closed-shell (fragments 2 and 4), the two irons as Fe(3+) sextets antiferromagnetically coupled (fragments 1 and 3, with 1 being alpha-spin
and 3 beta-spin), and each of the S-R fragments (5-8) as closed-shell singlet anions.
需要注意的是,
%chk=FragGuess
%mem=64mw
#P UBP86/6-311G(d) Guess=(Fragment=8,Only) Pop=None

第一步只做初猜,第二步在家优化或其他性质计算的关键词。有时候为了得到你想要的构型,你只能选比较小的基组算,收敛后慢慢往上加基组。

如果用G03,高斯网站上介绍过一种方法,觉得比较复杂。一个简单的方法是,先收敛到多重度为3的构型,然后再把多重度改为1继续算,很快能得到反铁磁的构型。如果仍有困难,计算量不大的话,把高自旋的构型的结果拿出来,再把多重度改为1继续算,应该会有好的结果。

没算过J值,J值得计算可以参考附件的文献!
反铁磁.rar
http://kuai.xunlei.com/d/GSJVPILVWETJ?p=130497
2楼2012-08-28 10:58:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingznpmg

铁虫 (小有名气)

【答案】应助回帖

哦,线性的反铁磁的Fragment上的多重度应该在Mn上间隔设置,如,Mn1:5,M2:-5, Mn3:5,Mn4:-5
3楼2012-08-28 11:02:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨寻

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
gaussian03没用过,用gaussian09可以总过组合片段波函数来设置反铁磁,但是关键是一个分子要分成几块,对应的多重度和电荷设置,非磁性部分是不是都必须是单重态?希望高手多多指点。。。
4楼2012-08-28 13:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨寻

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by lingznpmg at 2012-08-28 10:58:48
对于线性的构型了解的不多,但对反铁磁的计算,可以参考G09手册以下部分:
Antiferromagnetic coupling. Here is an example of a fragment guess job. The first step generates guesses for each one of the fra ...

请问下一个分子到底如何分割几块最合适?奇数和偶数块组合的波函数有差别吗?有时候我把一个分子分割的时候非磁性部分的那块自旋多重度不是单重态 计算就出错。
两步需不需要stable=opt这个命令,直接用结构算不优化得到的J值偏离太离谱,刚刚开始摸索,请多指教。。谢谢。。
5楼2012-08-28 13:30:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingznpmg

铁虫 (小有名气)

【答案】应助回帖

引用回帖:
5楼: Originally posted by 杨寻 at 2012-08-28 13:30:17
请问下一个分子到底如何分割几块最合适?奇数和偶数块组合的波函数有差别吗?有时候我把一个分子分割的时候非磁性部分的那块自旋多重度不是单重态 计算就出错。
两步需不需要stable=opt这个命令,直接用结构算不优 ...

这个不太清楚,我练习的时候用的G09Fe4S4的例子,把它分了8个碎片,一个原子一个碎片,用高基组收敛不了,用小基组可以收敛。后来,发现Fe2S2单重态是闭壳层的抗磁性的(Fe上自旋密度是0,0),stable=opt后就变成反铁磁的。
对于好多都比较糊涂,请有经验的虫虫多帮帮忙!
6楼2012-08-28 16:37:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨寻

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by lingznpmg at 2012-08-28 16:37:54
这个不太清楚,我练习的时候用的G09Fe4S4的例子,把它分了8个碎片,一个原子一个碎片,用高基组收敛不了,用小基组可以收敛。后来,发现Fe2S2单重态是闭壳层的抗磁性的(Fe上自旋密度是0,0),stable=opt后就变成反 ...

Fe4S4分割的很巧妙,除了两个Fe-fragment,其它非金属都是单重态,第一步好像只是组合一个初始波函数留作第二步计算能量用。再问下那个文献里的borken symmetry是第二步计算的模型是在反铁磁的模型基础上加上取消对称行命令scf=nosymm,还是直接就用反铁磁模型?
7楼2012-08-28 18:28:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingznpmg

铁虫 (小有名气)

引用回帖:
7楼: Originally posted by 杨寻 at 2012-08-28 18:28:51
Fe4S4分割的很巧妙,除了两个Fe-fragment,其它非金属都是单重态,第一步好像只是组合一个初始波函数留作第二步计算能量用。再问下那个文献里的borken symmetry是第二步计算的模型是在反铁磁的模型基础上加上取消对 ...

“Fe4S4分割的很巧妙,除了两个Fe-fragment,其它非金属都是单重态,第一步好像只是组合一个初始波函数留作第二步计算能量用。“这个是对的。

“再问下那个文献里的borken symmetry是第二步计算的模型是在反铁磁的模型基础上加上取消对称行命令scf=nosymm,还是直接就用反铁磁模型?”这个问题我不知道,我只算过构型和能量,没算过J值。他们用的是Turbmole,听一个老师说Turbmole处理过渡金属很好!

希望更多的前辈参与讨论,指导!

8楼2012-08-28 21:53:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨寻

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
8楼: Originally posted by lingznpmg at 2012-08-28 21:53:54
“Fe4S4分割的很巧妙,除了两个Fe-fragment,其它非金属都是单重态,第一步好像只是组合一个初始波函数留作第二步计算能量用。“这个是对的。

“再问下那个文献里的borken symmetry是第二步计算的模型是在反铁磁 ...

非常感谢。。只有自己组合尝试了。。。
9楼2012-08-28 22:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
10楼2012-08-28 22:24:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lhjxgy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 帮忙看看fileCode +6 wwncly 2026-08-10 11/550 2026-08-11 09:14 by xiaruohan
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见