应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】各位大侠！小女子求助钙钛矿结构计算磁性的问题？
371/1返回列表
查看: 5554  |  回复: 36
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hyff0809

银虫 (小有名气)

[交流] 【求助】各位大侠！小女子求助钙钛矿结构计算磁性的问题？

我最近试图用vasp算关于钙钛矿结构物质的磁性，关于铁磁和反铁磁的磁胞不知道到底该怎么选取。  举个简单的例子说吧，如果要算SrMnO3,它应该是G型反铁磁的，算磁性的话磁胞应该怎么取？是原来原胞的8倍吗（2x2x2）?还是什么？


请各位大侠踊跃发言，因为实验室没师兄师姐带，这个问题已经纠结了很久了，才来求助万能的小木虫！
回复此楼

» 收录本帖的淘帖专辑推荐

第一原理资料汇编 第一性原理杂

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-04-11 20:53:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

songdewang

金虫 (小有名气)
★ ★ ★
hyff0809(金币+1):谢谢参与
hyff0809(金币+5): 2011-04-11 22:01:57
zzy870720z(金币+2): 谢谢分享 2011-04-12 08:41:54
我也是做磁性计算的，可以探讨探讨。是这样的，你先取好一个supercell，比如说2X2X2的，然后在INCAR里面设置参数，下面是我的INCAR：
System = **************
    LPLANE = .TRUE.
    NPAR = 10
    LSCALU = .FALSE.
    NSIM = 4

    ISTART = 0   
    ICHARG=2   
    ISPIN = 2
    VOSKOWN = 1
    MAGMOM=22*0  24*0  5  5
    ENCUT  = 500
    ISMEAR = 0  
    SIGMA = 0.2
    EDIFF = 1E-5
    EDIFFG = -0.02
    NSW = 100
    IBRION = 1
    POTIM=0.6
    PREC=Accurate

    ALGO = Fast
    NELMIN = 6
    WEIMIN=0
    MAXMIX = 80

    LREAL = Auto
    LCHARG=F
    LWAVE=F
注意MAGMOM，如果是非磁性元素的话就写原子数*0，如果是磁性原子，比方两个Mn原子，如果计算FM，就5  5；如果是AFM计算，则设置5  -5。计算结果比较FM和AFM的总能量高低，得出哪种态稳定。
赞一下(4人) 回复此楼
4楼2011-04-11 21:33:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
★ ★
hyff0809(金币+1):谢谢参与
hyff0809(金币+3): 2011-04-13 18:57:16
zzy870720z(金币+1): 谢谢分享 2011-04-13 22:39:57
example
G型反铁磁 CaMnO3  
system=CaMnO3
1.00000
5.2790999413         0.0000000000         0.0000000000
0.0000000000         7.4566001892         0.0000000000
0.0000000000         0.0000000000         5.2705998421
Ca Mn O
4  4  12
Direct
0.032299999         0.250000000         0.990700006
0.967700005         0.750000000         0.009300000  
0.467700005         0.750000000         0.490700006  
0.532299995         0.250000000         0.509299994  
0.000000000         0.000000000         0.500000000  
0.500000000         0.000000000         0.000000000  
0.000000000         0.500000000         0.500000000  
0.500000000         0.500000000         0.000000000  
0.493200004         0.250000000         0.068300001  
0.506799996         0.750000000         0.931699991  
0.006799996         0.750000000         0.568300009  
0.993200004         0.250000000         0.431699991  
0.290399998         0.030200001         0.710999966  
0.709599972         0.969799995         0.289000005  
0.209600002         0.969799995         0.210999995  
0.790400028         0.030200001         0.789000034  
0.709599972         0.530200005         0.289000005  
0.290399998         0.469799995         0.710999966  
0.790400028         0.469799995         0.789000034  
0.209600002         0.530200005         0.210999995  

MAGMOM      =    4*0 3 -3  -3  3  12*0
赞一下(3人) 回复此楼
28楼2011-04-13 17:34:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hyff0809

银虫 (小有名气)
自己顶一下，没人理我。坐等大侠指导啊！！
赞一下 回复此楼
2楼2011-04-11 21:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmhyingyu

至尊木虫 (文坛精英)
★
hyff0809(金币+1):谢谢参与
hyff0809(金币+1): 2011-04-11 22:02:10
bless！
回复此楼
5楼2011-04-11 21:45:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by zyxme at 2011-04-11 21:09:58:

非常感谢您的交流!!!!

我有两个问题：
1. Mn 的初始磁矩为什么取5呢？对MAGMOM的取值一直不知道怎么取
2.我看到有的文献上取得磁胞在xy平面内是原胞的2倍，而依照您所说的应该是4倍吧？而且貌似A-AFM和G-AFM取得还不一样吧？
赞一下 回复此楼
6楼2011-04-11 22:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)

回复songdewang

不好意思，刚才上面引错了，是对你的回复
赞一下 回复此楼
7楼2011-04-11 22:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
  没人了吗？
回复此楼
8楼2011-04-11 22:29:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)
★ ★
hyff0809(金币+1):谢谢参与
zzy870720z(金币+1): 谢谢指点 2011-04-12 08:42:12
MAGMOM的取值都是参考实验值或者是文献中的取值。取的太大或者太小会影响计算结果和计算速度。vasp 的handson上面有一节专门说这个的
赞一下(2人) 回复此楼
9楼2011-04-11 22:57:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by hn19870519 at 2011-04-11 22:57:21:
MAGMOM的取值都是参考实验值或者是文献中的取值。取的太大或者太小会影响计算结果和计算速度。vasp 的handson上面有一节专门说这个的

如果没有实验值或查不到相关文献的话，取1或-1能给出较好的结果吗？
赞一下 回复此楼
10楼2011-04-11 23:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

songdewang

金虫 (小有名气)
★
ellsaking(金币+1): 谢谢指点 2011-04-12 09:24:52
引用回帖:
Originally posted by hyff0809 at 2011-04-11 22:09:40:
非常感谢您的交流!!!!

我有两个问题：
1. Mn 的初始磁矩为什么取5呢？对MAGMOM的取值一直不知道怎么取
2.我看到有的文献上取得磁胞在xy平面内是原胞的2倍，而依照您所说的应该是4倍吧？而且貌似A-AFM和G-A ...

Mn的初始磁矩就是5啊，Mn的d电子有5个
赞一下(1人) 回复此楼
11楼2011-04-12 09:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)
★
bingmou(金币+1): 感谢交流！ 2011-04-12 15:01:34
引用回帖:
Originally posted by hyff0809 at 2011-04-11 23:02:45:
如果没有实验值或查不到相关文献的话，取1或-1能给出较好的结果吗？

一般都能查到的吧。。。handson上面讲的取太小可能会得到错误的结果。具体你要自己测试下。
赞一下(1人) 回复此楼
12楼2011-04-12 13:18:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by songdewang at 2011-04-12 09:24:04:
Mn的初始磁矩就是5啊，Mn的d电子有5个

这个我完全不懂，d电子有几个就设几吗？


欢迎大家多多交流哈！
赞一下 回复此楼
13楼2011-04-12 18:37:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxiyue

铁杆木虫 (著名写手)
★
hyff0809(金币+1):谢谢参与
祝福一下，早日解决
回复此楼
14楼2011-04-12 19:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

franch

版主 (著名写手)
★
hyff0809(金币+1):谢谢参与
引用回帖:
Originally posted by hyff0809 at 2011-04-12 18:37:16:
这个我完全不懂，d电子有几个就设几吗？


欢迎大家多多交流哈！

不确定，，， 建议楼主自己测试下，，
赞一下(1人) 回复此楼
15楼2011-04-12 19:28:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by franch at 2011-04-12 19:28:49:
不确定，，， 建议楼主自己测试下，，

请教您，测试的标准是什么呢？是指设定的值和计算出来的值差不多吗？
赞一下 回复此楼
16楼2011-04-12 21:45:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangla

新虫 (正式写手)
★
hyff0809(金币+1):谢谢参与
顺便问下，磁性原子在不同的磁性体系下的磁距设置一样吗？
赞一下(1人) 回复此楼
17楼2011-04-12 21:53:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangla

新虫 (正式写手)
也就是在d电子附近取值测试吗？我看到有说测试标准是看能带的带隙，是这样吗？
赞一下 回复此楼
18楼2011-04-12 21:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangla

新虫 (正式写手)
追问下，同一磁性原子在不同的磁性体系下的初始磁距设置一样吗？设置都是参考文献得来的还是计算得到的？本人也是初学着，就弱弱的问下，谢谢！
赞一下 回复此楼
19楼2011-04-12 22:03:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmhyingyu

至尊木虫 (文坛精英)
bless！
回复此楼
20楼2011-04-12 22:11:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmhyingyu

至尊木虫 (文坛精英)
bless！
回复此楼
21楼2011-04-12 22:11:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)
★ ★ ★
hyff0809(金币+1):谢谢参与
hyff0809(金币+3): 2011-04-13 12:54:48
zzy870720z(金币+2): 谢谢指教 2011-04-13 22:38:30
怎么变成了祝福贴？
这个初试值真的只是“初试值”而已，对于一般简单的磁体系只要符号设置对了就可以找到正确的结果。
只有体系比较复杂，某一磁态附近有多个基态或极小值点时，初始值不对会导致得不到想要的结果，这个时候才需要仔细设置。
至于超胞的选择，你可以这么看：在这种磁状态下，把同一元素不同自旋的原子看成不同元素，那么这个新物质的最小结构是多大的，你需要选择的超胞就是多大的。
赞一下(2人) 回复此楼
23楼2011-04-12 23:22:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by acridine at 2011-04-12 23:22:28:
怎么变成了祝福贴？
这个初试值真的只是“初试值”而已，对于一般简单的磁体系只要符号设置对了就可以找到正确的结果。
只有体系比较复杂，某一磁态附近有多个基态或极小值点时，初始值不对会导致得不到想要的结 ...

谢谢您的热心回复。我也是像您所说那样想的，看成不同的原子选取最小的单胞，只是如果要比较不同磁序的能量的话，是不是应该在相同的体积下进行呢？不同磁序下所选的单胞应该不一样的吧？

盼望您的回复！
赞一下 回复此楼
25楼2011-04-13 13:01:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)
★ ★ ★
hyff0809(金币+1):谢谢参与
zzy870720z(金币+2): 谢谢指点 2011-04-13 22:39:32
G-type的ABO3的AFM可以参考如下文献，是BiFeO3的。你的材料应该也有实验报道的。

http://scripts.iucr.org/cgi-bin/paper?S0108768190006887
赞一下(2人) 回复此楼
26楼2011-04-13 14:58:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)
★ ★
zzy870720z(金币+2): 谢谢指教 2011-04-13 22:39:44
引用回帖:
Originally posted by hyff0809 at 2011-04-13 13:01:23:
谢谢您的热心回复。我也是像您所说那样想的，看成不同的原子选取最小的单胞，只是如果要比较不同磁序的能量的话，是不是应该在相同的体积下进行呢？不同磁序下所选的单胞应该不一样的吧？

盼望您的回复！

这个应该不是问题，因为你的能量可以平均成 eV/Mn,就是把能量平均到每个磁性原子，这样一来就可以比较不同体积了。
体积不同确实是容易造成能量的差别，因为体积的改变可能会导致计算参数，如k点密度、平面波数目的不一致。但如果不同体积的计算精度都足够高(比如，所有计算参数都满足收敛测试)，那么这种不一致是不会影响到计算结果的。即得到的平均每个晶胞的能量应该是可比的。
否则的话，如果增加几个重复单元，得到平均每个晶胞的能量不一致，那么计算岂不是变的没有意义了吗?我们也就没办法讨论结合能、形成能这些量了
赞一下(2人) 回复此楼
27楼2011-04-13 15:54:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by jghe at 2011-04-13 17:34:12:
example
G型反铁磁 CaMnO3  
system=CaMnO3
1.00000
5.2790999413         0.0000000000         0.0000000000
0.0000000000         7.4566001892         0.0000000000
0.0000000000         0.000000 ...

非常感谢您的热心交流。
但我想知道您上面提到的G-AFM CaMnO3 ,是CaMnO3的什么相？是不是不同的相算磁性时取得胞也不一样呢
赞一下 回复此楼
29楼2011-04-13 19:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by acridine at 2011-04-13 15:54:28:
这个应该不是问题，因为你的能量可以平均成 eV/Mn,就是把能量平均到每个磁性原子，这样一来就可以比较不同体积了。
体积不同确实是容易造成能量的差别，因为体积的改变可能会导致计算参数，如k点密度、平面波数 ...

谢谢您的回复！感觉您是这方面的高手，再请教您一个跟这个无关的问题，就是您在做PDOS的时候有没有出现过负的值，但都比较小？
赞一下 回复此楼
30楼2011-04-13 19:11:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
★ ★
zzy870720z(金币+2): 谢谢指教 2011-04-14 08:55:27
hyff0809(金币+2): 2011-04-14 09:28:47
我也是略知一二，只是恰好，我算过类似的体系。
我也碰到过DOS中有比较小的负值，大概在-2 左右，可以不管。如果比较大的负值就是计算的问题了，看有没有完全收敛。算DOS时候用 ISMEAR   =      -5，K点可以取密一些
赞一下(1人) 回复此楼
31楼2011-04-14 04:33:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyff0809

银虫 (小有名气)
引用回帖:
Originally posted by jghe at 2011-04-14 04:33:11:
我也是略知一二，只是恰好，我算过类似的体系。
我也碰到过DOS中有比较小的负值，大概在-2 左右，可以不管。如果比较大的负值就是计算的问题了，看有没有完全收敛。算DOS时候用 ISMEAR   =      -5，K点可以取密 ...

谢谢您！那我就释怀了，我的也不大，最大超不过0.5
赞一下 回复此楼
32楼2011-04-14 09:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

k_now

木虫 (正式写手)
★
hyff0809(金币+1): 谢谢参与
钙钛矿还可以算磁性？
赞一下 回复此楼
33楼2012-07-10 21:37:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

umklapp

银虫 (小有名气)
★
hyff0809(金币+1): 谢谢参与
引用回帖:
599743楼: Originally posted by hyff0809 at 2011-04-11 22:09:40
非常感谢您的交流!!!!

我有两个问题：
1. Mn 的初始磁矩为什么取5呢？对MAGMOM的取值一直不知道怎么取
2.我看到有的文献上取得磁胞在xy平面内是原胞的2倍，而依照您所说的应该是4倍吧？而且貌似A-AFM和G-AFM取 ...

话说这“小女子”的帖子 就这么多人回复，看看我发的帖子真是应者寥寥；
赞一下 回复此楼
34楼2012-07-11 00:27:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amberkaka

木虫 (初入文坛)
★ ★
hyff0809(金币+1): 谢谢参与
uuv2010: 金币+1, 多谢指导 2012-07-12 09:54:05
引用回帖:
602536楼: Originally posted by hyff0809 at 2011-04-11 23:02:45
如果没有实验值或查不到相关文献的话，取1或-1能给出较好的结果吗？...

MAGMOM一般取实验值的1.3~1.5倍，最好不要太小。
这个手册里可以查到的。
赞一下 回复此楼
35楼2012-07-11 17:20:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amberkaka

木虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
1270201楼: Originally posted by umklapp at 2012-07-11 00:27:08
话说这“小女子”的帖子 就这么多人回复，看看我发的帖子真是应者寥寥；...

难道我下次发帖也要加个“小女子”？（摸下巴思索状）
赞一下 回复此楼
37楼2012-07-13 17:25:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
zyxme3楼
2011-04-11 21:09   回复  
hyff0809(金币+1):谢谢参与
zzy870720z(金币-2): 你好，已经警告过你三次了，以后都双倍扣除所得，谢谢 2011-04-11 21:11:58
bingmou(金币-1): 鞭尸处理 2011-04-12 15:00:50
spc0822楼
2011-04-12 22:28   回复  
hyff0809(金币+1):谢谢参与
祝福祝福
wittiger24楼
2011-04-12 23:41   回复  
hyff0809(金币+1):谢谢参与
zzy870720z: 请不要在悬赏贴中纯表回复，警告一次，相信你也不希望在悬赏求助是得到的全是表情吧，谢谢合作，下次没收所得 2011-04-13 22:38:48
uuv201036楼
2012-07-12 09:53   回复  
hyff0809(金币+1): 谢谢参与
引用回帖:
1270201楼: Originally posted by umklapp at 2012-07-11 00:27:08 话说这“小女子”的帖子 就这么多人回复，看看我发的帖子真是应者寥寥；...

相关版块跳转 我要订阅楼主 hyff0809 的主题更新
371/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 5/250 2026-08-24 21:51 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭