《vasp在计算磁性的实例和讨论》由会员分享，可在线阅读，更多相关《vasp在计算磁性的实例和讨论（9页珍藏版）》请在金锄头文库上搜索。
1、之答禄夫天创作时间：二O二一年七月二十九日兄弟，问3个问题1,vasp在计算磁性的时候,oszicar中获得的磁 矩和outcar中获得各原子磁矩之和纷歧致,在投稿的是否曾碰到 有审稿人质疑，对这个纷歧致你们一般是怎么解释的了？ 2,另外， 磁性计算应该比力负责.你应该还使用另外法式计算过磁性，与 vasp结果比力是否一致,对磁性计算采纳的法式有什么推荐.ps： 由于曾使用vasp和dmol算过非周期体系磁性，结构对磁性影响非 常年夜，因此使用这两个法式计算的磁性要一致很麻烦.还不敢确 定究竟是哪个法式可能不成靠.3,如果采纳vasp计算磁性,对采纳 的方法和设置有什么推荐.1, OSZICA
2、R中获得的磁矩是OUTCAR中最后一步获得的总磁矩是相 等的.总磁矩和各原子的磁矩（RMT球内的磁矩）之和之差就是间 隙区的磁矩.因为有间隙区存在，纷歧致是正常的.2,如果算磁 性,全电子的结果更精确，我的一些计算结果显示磁性原子对在 最近邻的位置时,PAW与FPLAW给出的能量差纷歧致，在长程时 符合的很好.虽然并没有改变定性结论.感觉PAW似乎不能很好 地描述较强耦合.我试图在找出原因，主要使用exciting和 vasp做比力.计算磁性推荐使用FP-LAPW, FP-LMTO, FPLO很吸 引人（不外是商业的），后者是O（N）算法.3,使用vasp计算磁性, 注意分歧的初始磁矩是否收敛
3、为同一个磁矩.倒没有特别要注 意的处所，个人认为.归根结底，需要一个优秀的交换关联形式 时间：二O二一年七月二十九日呈现VASP计算是否也是像计算DOS和能带一样要进行三步（结构优化, 静态自洽计算，非自洽计算），然后看最后一步的出的磁矩呢？ 一直想计算固体中某个原子的磁矩，根据OUTCAR的结果似乎不能 分析，因为它里面总磁矩跟OSZICAR的值有一定的分歧，据说是 OUTCAR中只考虑WS半径内磁矩造成的.最近看到一个帖子说是可 以用bader电荷分析方法分析原子磁矩.如法炮制之后发现给出的 总磁矩与OSZICAR的结果符合的甚好，可是觉得没有根据，有谁知 道这样做的依据吗，欢迎讨论！设置
4、ISPIN=2计算获得的态密度成为自旋态密度.设置ISPIN=2就可以计算磁性，铁磁和反铁磁在MAG里设置.最后 获得的DOS是分up和down的.磁性计算(2006-12-03 21:02)标签:分类：Vasp磁性计算 顺磁，意味进行non-spin polarized的计算，也就是ISPIN=1.铁磁，意味进行spin-polarized的计算，ISPIN=2,而且每个磁性原 子的初始磁矩设置为一样的值，也就是磁性原子的MAGMOM设置为 一样的值.对非磁性原子也可以设置成一样的非零值（与磁性原子 的一样）或零，最后收敛的结果,非磁性原子的local磁矩很小, 快接近0,很小的情况,很可能
5、意味着真的是非磁性原子也会被极 化而呈现很小的local磁矩.反铁磁，也意味着要进行spin-polarized的计算，ISPIN=2,这是需 采纳反铁磁的磁胞来进行计算,意味着此时计算所采纳的晶胞不再 是铁磁计算时的最小原胞.比如对铁晶体的铁磁状态，你可以采纳 bcc的原胞来计算,可是在进行反铁磁的Fe计算，这是你需要采纳 sc的结构来计算，计算的晶胞中包括两个原子，你要设置一个原子 的MAGMOM为正的，另一个原子的MAGMOM设置为负,可是它们的绝 对值一样.因此在进行反铁磁的计算时，应该确定好反铁磁的磁胞, 以及磁序,要判断哪种磁序和磁胞是最可能的反铁磁状态，那只能 是先做好各种可能的
6、排列组合，然后分别计算这些可能组合的情况, 最后比力它们的总能，总能最低的就是可能的磁序.同样也可以与 它们同铁磁或顺磁的进行比力.了解到该资料究竟是铁磁的、还是 顺磁或反铁磁的.亚铁磁，也意味要进行spin-polarized的计算，ISPIN=2,与反铁磁 的计算类似，分歧的是原子正负磁矩的绝对值不是样年夜.非共线的磁性，那需采纳专门的non-collinear的来进行计算，除 要设置ISPIN,MAGMOM的设置还需要指定每个原子在x,y,z方向上 的年夜小.这种情况会复杂一些.参杂Co原子的CdS稀磁半导体的态密度计算1. 2. 3System = CdSENCUT =500;ISTA
7、RT= 0;ICHARG= 2;GGA = 91;ISPIN=2VOSKOWN=1ISMEAR = -5;SIGMA = 0.1;NSW=165IBRION = 2;ISIF = 3;POTIM = 0.2;EDIFF = 1E-4;EDIFFG = -1E-2;LREAL = AutoPREC = Accurate体系的总磁矩（OZICAR最后一 行MAG=）System = CdS-CoENCUT =500;ISTART= 0;ICHARG= 2;GGA = 91;ISMEAR = -5;SIGMA = 0.1;ISPIN= 2VOSKOWN= 1#NSW=165#IBRION = 1;#
8、ISIF = 3;#POTIM = 0.2;EDIFF = 1E-5;EDIFFG = -1E-3;LREAL=AutoPREC = Accurate3.静态非自洽计算获得自旋态密度以及体系的总磁矩（OZICAR最后一行MAG=）磁矩与2中的相同System = CdS-CoENCUT =500;ISTART= 0;ICHARG= 11;GGA = 91;ISMEAR = -5;SIGMA = 0.1;ISPIN= 2VOSKOWN= 1LORBIT= 10 #NSW=165 #IBRION = 2;#ISIF = 3;#POTIM = 0.2;EDIFF = 1E-5;EDIFFG = -1
9、E-3;LREAL=AutoPREC = AccurateNPAR= 1做铁磁计算或者反铁磁计算需加MAGMOM参数晶胞中只有一种原子时:MAGMOM= 2m*3铁磁计算(2m是总的原子数,3是每个原子初始磁矩的年夜小)MAGMOM二m*-3 m*3反铁磁计算Types of spin ordering in perovskite oxides铁磁FM反铁磁AFMA-type:The intra-plane coupling isferromagnetic while inter-plane coupling is antiferromagnetic.C-type:The intra-plan
10、e coupling isantife磁结构的计算rromagneticwhileinterplanecoupli ng isferromagneti c.G-type:Bothintra-plane and inter-plane coupling are antiferromagnetic.题目：转载如何用VASP计算铁磁、反铁磁和顺磁浏览：95评论：0顺磁，意味进行non-spin polarized的计算，也就是ISPIN=1.铁磁，意味进行spin-polarized的计算,ISPIN=2,而且每个磁性原子的初始磁 矩设置为一样的值，也就是磁性原子的MAGMOM设置为一样的值.对非磁
11、性 原子也可以设置成一样的非零值（与磁性原子的一样）或零,最后收敛的结 果,非磁性原子的local磁矩很小,快接近0,很小的情况,很可能意味着真的是 非磁性原子也会被极化而呈现很小的local磁矩.反铁磁，也意味着要进行spin-polarized的计算,ISPIN=2,这是需采纳反铁磁 的磁胞来进行计算,意味着此时计算所采纳的晶胞不再是铁磁计算时的最小 原胞.比如对铁晶体的铁磁状态，你可以采纳bcc的原胞来计算，可是在进行 反铁磁的Fe计算，这是你需要采纳sc的结构来计算，计算的晶胞中包括两个 原子,你要设置一个原子的MAGMOM为正的，另一个原子的MAGMOM设 置为负，可是它们的绝对值一
12、样.因此在进行反铁磁的计算时，应该确定好反铁 磁的磁胞，以及磁序，要判断哪种磁序和磁胞是最可能的反铁磁状态，那只能是 先做好各种可能的排列组合，然后分别计算这些可能组合的情况，最后比力它 们的总能,总能最低的就是可能的磁序.同样也可以与它们同铁磁或顺磁的进 行比力.了解到该资料究竟是铁磁的、还是顺磁或反铁磁的.NUPDOWN设置为一个为非0而且非-1的数，就是相当是将在自洽迭代循环的计算过程加了一个限制,即限制自旋向 上和自旋向下的电子数之差为一个固定的值，对应的就是给体系设置一个固定 的总磁矩.具体在法式里面实现起来，量子化学的法式（比如Gaussian ）和广泛用 于固体物理的第一原理法式
13、（比如WIEN2K）略有分歧.亚铁，结构，原子，磁矩，整数最近遇到一个问题,如何计算分歧的磁性结构？如铁磁结构,反铁磁 结构，亚铁磁结构.vasp教程上有说明，设置ISPIN=2,MAGMOM=,这样来说，计算铁磁结构只需要把各个原子的磁矩方向设为一致，如,都为正值就可以了；计算反铁磁结构只 需要把分歧子子晶格的原子磁矩方向设为相反，如，一正一负就可以 了；计算亚铁磁结构或许还要用到NUPDOWN.可是有个问题，我做了一个亚铁磁结构的计算，初始MAGMOM设置了 有正有负的原子磁矩，没有设置NUPDOWN值，而优化出来的结果里各 个原子的磁矩却是按同一方向排列的！也就是说没有获得我想要的 亚铁
14、磁结构.如果要设置NUPDOWN进行强制总磁矩的话，又不知道这 个值该设成多年夜？非常希望各位多给意见，只要是有相关内容的回帖都送金币.收藏分享评分回复引用TOP界说：自旋多重度(spin multiplicity)=2S+1, S=n*1/2,n 为 单电子数.所以，关键是单电子的数目是几多.当有偶数个电子时， 例如O2,共有16个电子，那么单电子数目可能是0,即8个alpha 和8个beta电子配对,对应单重态,可是也可能是有9个a电子 和7个8电子，那么能成对的是7对,还剩2个a没有配对，于是 n=2,对应的是多重度3.同理还可以有多重度5,7,9, . 一般而言, 是多重度低的能量低,最稳定，所以，一般来说，偶数电子的体系多 重度就是1.可是也有例外，如果02就是一个年夜家都知道的例子, 它的基态是三重态，其单重态反而是激发态.所以，总结一下，就是 电子数目是偶数，未成对电子数目n=0,2,4,6,.自旋多重度是 1,3,5,7,.电子数目是奇数，未成对电子数目n=1,3,5,7,.自 旋多重度是2,4,6,8,.大都情况是多重度低的能量低，有时(特 别是有“磁”性的时候，例如顺磁的02,以及Fe啊什么的)，可能 会高多重度的能量低.时间：二O二一年七月二十九日
