计算化学公社

标题: Pd(111)晶面计算INCAR设置求助 [打印本页]

作者
Author: Garyall 时间: 2022-4-1 15:55
标题: Pd(111)晶面计算INCAR设置求助
请问Pd金属计算的时候需要考虑自旋ISPIN=2

MAGMOM=原子个数*磁矩
如果要考虑磁矩是多少

哪里可以查到磁矩信息的相关资料
还有个问题就是ISPIN=2是铁磁材料（是否可以通俗理解为：用吸铁石吸起来的材料就是铁磁材料）
附件是我INCAR设置是不是还需要修改

作者
Author: jiangning198511 时间: 2022-4-1 16:26
VASP中磁矩是自动优化的

作者
Author: 卡开发发 时间: 2022-4-1 16:48
Pd应该没磁性吧，不用考虑ISPIN=2，如果表面上还吸附有自由基等，那就给那些带不成对电子的原子设置磁矩就可以。吸铁石能吸的不一定是铁磁结构，例如磁铁是Fe3O4，本身能被磁铁吸起来，但它是亚铁磁的。

作者
Author: Garyall 时间: 2022-4-1 21:02

卡开发发发表于 2022-4-1 16:48
Pd应该没磁性吧，不用考虑ISPIN=2，如果表面上还吸附有自由基等，那就给那些带不成对电子的原子设置磁矩就 ...

嗯 thanks 我懂了

就是单纯算Pd时不需要考虑，但是后面我算的是Pd表面亚胺物质加氢反应涉及到H原子（H自由基）是要引入磁矩的。这里涉及到一个我的未知领域了怎么样设置（给那些带不成对电子的原子设置磁矩）

希望您有空的时候帮忙回答下~~

作者
Author: Garyall 时间: 2022-4-1 21:13

jiangning198511 发表于 2022-4-1 16:26
VASP中磁矩是自动优化的

就是说当MAGMOM不用设置也可以吗

作者
Author: 卡开发发 时间: 2022-4-1 21:23

Garyall 发表于 2022-4-1 21:02
嗯 thanks 我懂了就是单纯算Pd时不需要考虑，但是后面我算的是Pd表面亚胺物质加氢反应 ...

对吸附的情况而言，那一般加在H上就行，主要看两个要素，一个是要具备不成对电子（这个可能需要结构化学的一些知识），另外就是这个单电子必须局域化。对固体的情况可能比较复杂，这种也只能通过实验和文献以及其他数据综合判断。磁矩必须得猜测，VASP中开ISPIN=2默认每个原子全按1猜，如果是高自旋态会可能会出现难收敛，对于反铁磁或亚铁磁这么猜可能会被算成铁磁结果。

作者
Author: Garyall 时间: 2022-4-2 09:27

卡开发发发表于 2022-4-1 21:23
对吸附的情况而言，那一般加在H上就行，主要看两个要素，一个是要具备不成对电子（这个可能需要结构化学 ...

谢谢了

关于这ISPIN=2的知识我知道的太少了。。我得先再去了解一二

作者
Author: granvia 时间: 2022-4-2 11:00

卡开发发发表于 2022-4-1 16:48
Pd应该没磁性吧，不用考虑ISPIN=2，如果表面上还吸附有自由基等，那就给那些带不成对电子的原子设置磁矩就 ...

Pd(111)不清楚。但只有几层的Pd(100)是有磁性的（10.1063/1.4757409，10.1063/1.4757409），所以通常使用的4-5层slab模型是不能模拟Pd(100)的

作者
Author: jiangning198511 时间: 2022-4-2 11:30

granvia 发表于 2022-4-2 11:00
Pd(111)不清楚。但只有几层的Pd(100)是有磁性的（10.1063/1.4757409，10.1063/1.4757409），所以通常使用 ...

貌似催化计算对表面的磁性不是很关注，不知道这个对催化所关心的吸附能有什么影响

作者
Author: 卡开发发 时间: 2022-4-2 12:25
本帖最后由卡开发发于 2022-4-2 12:26 编辑

granvia 发表于 2022-4-2 11:00
Pd(111)不清楚。但只有几层的Pd(100)是有磁性的（10.1063/1.4757409，10.1063/1.4757409），所以通常使用 ...

建议很好，所以从避雷的角度说最好别用那么少的层数去模拟Pd(100)，可能5~6层也是巧合（10.1103/PhysRevB.75.172402），虽然模拟的目标是NM的（也就是宏观表面）。同样来说，如果slab层高影响如此之大，那么如果对于簇模型而言团簇尺寸的选取影响将会更加巨大，此其一。其二来说，虽然单点能量看薄层FM稳定，但Curie温度多高？是否到了室温其实应该是NM更稳定？这个问题很难讲。

作者
Author: huanglei 时间: 2022-6-29 15:52

卡开发发发表于 2022-4-1 21:23
对吸附的情况而言，那一般加在H上就行，主要看两个要素，一个是要具备不成对电子（这个可能需要结构化学 ...

老师您好，那么对于Fe3O4这种亚磁铁该怎么设置呢

作者
Author: 卡开发发 时间: 2022-6-29 19:24

huanglei 发表于 2022-6-29 15:52
老师您好，那么对于Fe3O4这种亚磁铁该怎么设置呢

你可以找找文献，有文献给出Fe3O4亚铁磁性的设置方式，Materials Project可能也有相应的结构对应的输入。

作者
Author: huanglei 时间: 2022-6-30 19:51

卡开发发发表于 2022-6-29 19:24
你可以找找文献，有文献给出Fe3O4亚铁磁性的设置方式，Materials Project可能也有相应的结构对应的输入。

好嘞，谢谢老师

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