vasp在计算磁性的实例和讨论

创作时间：二零二一年六月三十天之邯郸勺丸创作创作时间：二零二一年六月三十天兄弟,问3个问题1,vasp在计算磁性的时候,oszicar中获取的磁矩和outcar中获取各原子磁矩之和纷歧致,在投稿的能否曾碰到有审稿人怀疑,对这个纷歧致你们一般是怎么解说的了？2,此外,磁性计算应当比力负责.你应当还使用此外法式计算过磁性,与vasp结果比力能否一致,对磁性计算采用的法式有什么推荐.ps：因为曾使用vasp和dmol算过非周期系统磁性,构造对磁性影响特别大年夜,所以使用这两个法式计算的磁性要一致很麻烦.还不敢确立终究是哪个法式可能不可靠.3,假如采用vasp计算磁性,对采用的方法和设置有什么介绍.1，OSZICAR中获取的磁矩是OUTCAR中最后一步获取的总磁矩是相等的.总磁矩和各原子的磁矩(RMT球内的磁矩)之和之差就是空隙区的磁矩.因为有空隙区存在,纷歧致是正常的.2,假如算磁性,全电子的结果更精准,我的一些计算结果显示磁性原子对在近来邻的地点时,PAW与FPLAW给出的能量差纷歧致,在长程时切合的很好.固然并无改变定性结论.感觉PAW仿佛不可以很好地描绘较强耦合.我试图在找出原由,主要使用exciting和vasp做比力.计算磁性介绍使用FP-LAPW,FP-LMTO,FPLO很吸引人(不外是商业的),后者是O(N)算法.3,使用vasp计算磁性,注意分歧的初始磁矩能否收敛为同一个磁矩.倒没有特别要注意的地方,个人以为.归根结底,需要一个优异的交创作时间：二零二一年六月三十天创作时间：二零二一年六月三十天换关系形式体现VASP计算能否也是像计算DOS和能带相同要进行三步（构造优化,静态自洽计算,非自洽计算）,而后看最后一步的出的磁矩呢？向来想计算固体中某个原子的磁矩,依据OUTCAR的结果仿佛不可以剖析,因为它里面总磁矩跟OSZICAR的值有必定的分歧,听说是OUTCAR中只考虑WS半径内磁矩造成的.近来看到一个帖子说是可以用bader电荷剖析方法剖析原子磁矩.依样画葫芦以后发现给出的总磁矩与OSZICAR的结果切合的甚好,但是感觉没有依据,有谁知道这样做的依照吗,欢迎议论！设置ISPIN=2计算获取的态密度成为自旋态密度.设置ISPIN=2就能够计算磁性,铁磁和反铁磁在MAG里设置.最后获取的DOS是分up和down的.磁性计算(2006-12-0321:02)?标签：-分类：Vasp·磁性计算顺磁,意味进行non-spinpolarized的计算,也就是ISPIN=1.铁磁,意味进行spin-polarized的计算,ISPIN=2,并且每个磁性原子的初始磁矩设置为相同的值,也就是磁性原子的MAGMOM设置为相同的值.对非磁性原子也能够设置成相同的非零值（与磁性原子的相同）或零,最后收敛的结果,非磁性原子的local磁矩创作时间：二零二一年六月三十天创作时间：二零二一年六月三十天很小,快靠近0,很小的状况,很可能意味着真的是非磁性原子也会被极化而体现很小的local磁矩.反铁磁,也意味着要进行spin-polarized的计算,ISPIN=2,这是需采用反铁磁的磁胞来进行计算,意味着此时计算所采用的晶胞不再是铁磁计算时的最小原胞.比方对铁晶体的铁磁状态,你可以采用bcc的原胞来计算,但是在进行反铁磁的Fe计算,这是你需要采用sc的构造来计算,计算的晶胞中包含两个原子,你要设置一个原子的MAGMOM为正的,另一个原子的MAGMOM设置为负,但是它们的绝对值相同.所以在进行反铁磁的计算时,应当确立好反铁磁的磁胞,以及磁序,要判断哪一种磁序和磁胞是最可能的反铁磁状态,那只好是先做好各样可能的摆列组合,而后分别计算这些可能组合的状况,最后比力它们的总能,总能最低的就是可能的磁序.相同也能够与它们同铁磁或顺磁的进行比力.认识到该资料终究是铁磁的、仍是顺磁或反铁磁的.亚铁磁,也意味要进行spin-polarized的计算,ISPIN=2,与反铁磁的计算近似,分歧的是原子正负磁矩的绝对值不是样大年夜.非共线的磁性,那需采用特意的non-collinear的来进行计算,除要设置ISPIN,MAGMOM的设置还需要指定每个原子在x,y,z方向上的大年夜小.这类状况会复杂一些.参杂Co原子的CdS稀磁半导体的态密度计算System=CdSENCUT=500;创作时间：二零二一年六月三十天创作时间：二零二一年六月三十天ISTART=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;#ISIF=3;#POTIM=0.2;EDIFF=1E-5;EDIFFG=-1E-3;LREAL=AutoPREC=Accurate静态非自洽计算获取自旋态密度以及系统的总磁矩（OZICAR最后一行MAG=）磁矩与2中的相同System=CdS-CoENCUT=500;ISTART=0;ICHARG=11;GGA=91;ISMEAR=-5;SIGMA=0.1;ISPIN=2VOSKOWN=1创作时间：二零二一年六月三十天创作时间：二零二一年六月三十天LORBIT=10#NSW=165#IBRION=2;#ISIF=3;#POTIM=0.2;EDIFF=1E-5;EDIFFG=-1E-3;LREAL=AutoPREC=AccurateNPAR=1做铁磁计算或许反铁磁计算需加MAGMOM参数晶胞中只有一种原子时：MAGMOM=2m*3铁磁计算（2m是总的原子数,3是每个原子初始磁矩的大年夜小）MAGMOM=m*-3m*3反铁磁计算Typesofspinorderinginperovskiteoxides铁磁FM反铁磁AFMA-type:Theintra-planecouplingisferromagneticwhileinter-planecouplingisantiferromagnetic.创作时间：二零二一年六月三十天创作时间：二零二一年六月三十天C-type:Theintra-planecouplingis磁构造的计算intra-planeandinter-planecouplingareantiferromagnetic.

antiferromagneticwhileinter-planecouplingisferromagnetic.G-type:Both创作时间：二零二一年六月三十天创作时间：二零二一年六月三十天题目：[转载]怎样用VASP计算铁磁、反铁磁温顺磁阅读：95议论：0顺磁,意味进行non-spinpolarized的计算,也就是ISPIN=1.铁磁,意味进行spin-polarized的计算,ISPIN=2,并且每个磁性原子的初始磁矩设置为相同的值,也就是磁性原子的MAGMOM设置为相同的值.对非磁性原子也能够设置成相同的非零值（与磁性原子的相同）或零,最后收敛的结果,非磁性原子的local磁矩很小,快靠近0,很小的状况,很可能意味着真的是非磁性原子也会被极化而体现很小的local磁矩.反铁磁,也意味着要进行spin-polarized的计算,ISPIN=2,这是需采用反铁磁的磁胞来进行计算,意味着此时计算所采用的晶胞不再是铁磁计算时的最小原胞.比方对铁晶体的铁磁状态,你能够采用bcc的原胞来计算,但是在进行反铁磁的Fe计算,这是你需要采用sc的构造来计算,计算的晶胞中包含两个原子,你要设置一个原子的MAGMOM为正的,另一个原子的MAGMOM设置为负,但是它们的绝对值相同.所以在进行反铁磁的计算时,应当确立好反铁磁的磁胞,以及磁序,要判断哪一种磁序和磁胞是最可能的反铁磁状态,那只好是先做好各样可能的摆列组合,而后分别计算这些可能组合的状况,最后比力它们的总能,总能最低的就是可能的磁序.相同也能够与它们同铁磁或顺磁的进行比力.认识到该资料终究是铁磁的、仍是顺磁或反铁磁的.NUPDOWN设置为一个为非0并且非-1的数,,就是相当是将在自洽迭代循环的计算过程加了一个限制即限制自旋向上和自旋向下的电子数之差为一个固定的值,对应的就是给系统设置一个固定的总磁矩.详细在法式里面实现起来,量子化学的法式(比方Gaussian)和宽泛用于固体物理的第一原理法式(比方WIEN2K)略有分歧.亚铁,构造,原子,磁矩,整数近来碰到一个问题,怎样计算分歧的磁性构造？如铁磁构造,反铁磁构造,亚铁磁构造.vasp教程上有说明,设置ISPIN=2,MAGMOM=......,这样来说,计算铁磁构造只需要把各个原子的磁矩方向设为一致,如,都为正当就能够了；计算反铁磁构造只需要把分歧子子晶格的原子磁矩方向设为相反,如,一正一负就能够了；计算亚铁磁构造也许还要用到NUPDOWN.但是有个问题,我做了一个亚铁磁构造的计算,初始MAGMOM设置了有正有负的原子磁矩,没有设置NUPDOWN值,而优化出来的结果里各个原子的磁矩倒是按同一方向摆列的！也就是说没有获取我想要的亚铁磁构造.假如要设置NUPDOWN进行强迫总磁矩的话,又不知道这个值该设成多大年夜？特别希望各位多给建议,只假如有有关内容的回帖都送金币.珍藏分享评分创作时间：二零二一年六月三十天创作时间：二零二一年六月三十天答复引用TOP界说:自旋多重度（spinmultiplicity）＝2S+1,S=n*1/2,n为单电子数.所以,重点是单电子的数量是若干.当有偶数个电子时,比如O2,共有16个电子,那么单电子数量可能是0,即8个alpha和8个beta电子配对,对应单重态,但是也可能是有9个α电子和7个β电子,那么能成对的是7对,还剩2个α没有配对,于是n=2,对应的是多重度3.同理还能够有多重度5,7,9,...一般而言,是多重度低的能量低,最稳固,所以,一般来说,偶数电子的系统多重度就是1.但是也有例外,假如