计算化学公社

 找回密码 Forget password
 注册 Register

有2个铜和2个铁原子片段之间的相互作用能计算不出来,原因是?

查看数: 2739 | 评论数: 16 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-11-8 17:48

正文摘要:

想计算两个片段之间的相互作用能,有个片段有2个铜原子和2个铁原子,计算的不收敛,不知道是什么原因。两个铁原子距离很远,感觉应该不是反铁磁耦合吧

回复 Reply

zyj19831206 发表于 Post on 2022-11-9 15:14:35
Ferrocene 发表于 2022-11-9 15:00
建立合理的模型是做应用性计算的基本功,跟你的实验合作者提什么要求没有关系。这个体系是Fe3O4和CuBTC的 ...

那样做,体系太大了,感觉跑不动,所以不得已只能切簇算。
Ferrocene 发表于 Post on 2022-11-9 15:00:07
zyj19831206 发表于 2022-11-8 18:58
就是从MOF里面切了个簇,就这么算下去了做实验的人非要加的,没有办法。。

建立合理的模型是做应用性计算的基本功,跟你的实验合作者提什么要求没有关系。这个体系是Fe3O4和CuBTC的复合材料,是Fe3O4填到CuBTC的孔里,还是核壳结构,还是两个组分各成一相?如果是核壳或者各成一相的情况,我觉得最好是建立周期性模型用第一性原理计算,而不是用Gaussian切片段来算。另外既然是考虑材料和有机分子的相互作用,BTC配体的贡献也不能不考虑吧。
wzkchem5 发表于 Post on 2022-11-9 14:38:20
zyj19831206 发表于 2022-11-9 02:16
电荷为3 自旋多重度为1

自旋多重度显然不可能为1,Fe3O4里面Fe2+和Fe3+都是高自旋,这是常识
zyj19831206 发表于 Post on 2022-11-9 09:16:52
wzkchem5 发表于 2022-11-9 05:26
什么叫”3 1“?用规范术语说话,什么等于3,什么等于1

电荷为3 自旋多重度为1
wzkchem5 发表于 Post on 2022-11-9 05:26:40
zyj19831206 发表于 2022-11-8 16:08
如果当作三价,是不是两个铁和两个铜都得分别设置atom group,每个group是 3 1?铜每个都是0 1?这就是4个g ...

什么叫”3 1“?用规范术语说话,什么等于3,什么等于1
zyj19831206 发表于 Post on 2022-11-8 23:08:33
wzkchem5 发表于 2022-11-8 22:51
那就更麻烦一些了,得知道你取的这个簇里的Fe是Fe(II)还是Fe(III)才能算,这个最好能从实验上表征。

如果当作三价,是不是两个铁和两个铜都得分别设置atom group,每个group是 3 1?铜每个都是0 1?这就是4个group,其他原子是group那就是五个group?
wzkchem5 发表于 Post on 2022-11-8 22:51:13
zyj19831206 发表于 2022-11-8 14:55
是四氧化三铁与CuBTC复合材料。

那就更麻烦一些了,得知道你取的这个簇里的Fe是Fe(II)还是Fe(III)才能算,这个最好能从实验上表征。
zyj19831206 发表于 Post on 2022-11-8 21:58:39
wzkchem5 发表于 2022-11-8 20:51
首先,这两个Cu(II)应该是反铁磁耦合,因为醋酸铜是反铁磁耦合的。
其次,两个Fe都是高自旋的,但这两个 ...

您能在这个基础上大概帮我调一下看看嘛?具体怎么设置片段还算有点糊涂
zyj19831206 发表于 Post on 2022-11-8 21:55:54
wzkchem5 发表于 2022-11-8 20:51
首先,这两个Cu(II)应该是反铁磁耦合,因为醋酸铜是反铁磁耦合的。
其次,两个Fe都是高自旋的,但这两个 ...

是四氧化三铁与CuBTC复合材料。
wzkchem5 发表于 Post on 2022-11-8 20:51:23
本帖最后由 wzkchem5 于 2022-11-8 13:52 编辑
zyj19831206 发表于 2022-11-8 13:31
按照您的直觉,这个体系如果按照您的那样修改,换成O配体,一般自旋多重度应该是多少?

首先,这两个Cu(II)应该是反铁磁耦合,因为醋酸铜是反铁磁耦合的。
其次,两个Fe都是高自旋的,但这两个Fe和Cu之间是铁磁耦合还是反铁磁耦合不太好判断,如果没办法从实验上确定的话,建议铁磁耦合和反铁磁耦合都尝试一下。注意要用片段初猜的方法算,光把自旋多重度设对了是不够的。我不知道你这里Fe是几价的,准备片段初猜的时候肯定得知道Fe的价态,不过我以上说的这些不管Fe是几价都成立

zyj19831206 发表于 Post on 2022-11-8 20:31:51
wzkchem5 发表于 2022-11-8 20:09
切个簇没问题,但是第一电荷必须得对,MOF整体是中性的不代表切出来的簇也是中性的;第二要考虑背景电荷

按照您的直觉,这个体系如果按照您的那样修改,换成O配体,一般自旋多重度应该是多少?
wzkchem5 发表于 Post on 2022-11-8 20:09:51
zyj19831206 发表于 2022-11-8 11:58
就是从MOF里面切了个簇,就这么算下去了做实验的人非要加的,没有办法。。

切个簇没问题,但是第一电荷必须得对,MOF整体是中性的不代表切出来的簇也是中性的;第二要考虑背景电荷
zyj19831206 发表于 Post on 2022-11-8 18:58:47
Ferrocene 发表于 2022-11-8 18:19
以你这种拍脑袋建模的做法,计算出问题才是正常的,不出问题你应该去买彩票。你的模型里面Cu是几价,Fe又是 ...

就是从MOF里面切了个簇,就这么算下去了做实验的人非要加的,没有办法。。
zyj19831206 发表于 Post on 2022-11-8 18:56:44
wzkchem5 发表于 2022-11-8 18:26
你应该没见过[Fe(OH)6]3-吧?[Fe(OH)4]-已经需要很强的碱性条件才能存在了。如果Fe是3价,你往Fe上面加5个O ...

好的好的,这一点的确没有考虑到,当时只是为了让其饱和。。

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

GMT+8, 2026-2-24 04:21 , Processed in 0.195762 second(s), 25 queries , Gzip On.

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