Pd(111)晶面计算INCAR设置求助

Garyall

请问Pd金属计算的时候需要考虑自旋ISPIN=2

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

jiangning198511

VASP中磁矩是自动优化的

卡开发发

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

Garyall

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

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

Garyall

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

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

卡开发发

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

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

Garyall

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

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

granvia

卡开发发 发表于 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)的

jiangning198511

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

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

卡开发发

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更稳定？这个问题很难讲。

huanglei

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

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

卡开发发

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

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

huanglei

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

好嘞，谢谢老师