加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

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

(2012-08-09 12:47:02)
标签:

杂谈

分类: 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-03 21:02) ? 标签: - 分类: Vasp ·磁性计算 顺磁,意味进行 non-spin polarized 的计算,也就是 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 稀磁半导体的态密度计算 1. 2. 3 1. 结构优化 System = CdS ENCUT =500; ISTART= 0; ICHARG= 2; GGA = 91; ISPIN=2 VOSKOWN=1 ISMEAR = -5; SIGMA = 0.1; NSW=165 IBRION = 2; ISIF = 3; POTIM = 0.2; EDIFF = 1E-4; EDIFFG = -1E-2; LREAL = Auto PREC = Accurate 2.静态自洽计算得到体系的总磁矩(OZICAR最后一行MAG=) System = CdS-Co ENCUT =500; ISTART= 0; ICHARG= 2; GGA = 91; ISMEAR = -5; SIGMA = 0.1; ISPIN= 2 VOSKOWN= 1 #NSW=165 #IBRION = 1; #ISIF = 3; #POTIM = 0.2; EDIFF = 1E-5; EDIFFG = -1E-3; LREAL=Auto PREC = Accurate 3.静态非自洽计算得到自旋态密度以及体系的总磁矩(OZICAR 最后一行 MAG=)磁矩与 2 中的相同 System = CdS-Co ENCUT =500; ISTART= 0; ICHARG= 11; GGA = 91; ISMEAR = -5; SIGMA = 0.1; ISPIN= 2 VOSKOWN= 1 LORBIT= 10 #NSW=165 #IBRION = 2; #ISIF = 3; #POTIM = 0.2; EDIFF = 1E-5; EDIFFG = -1E-3; LREAL=Auto PREC = Accurate NPAR= 1 做铁磁计算或者反铁磁计算需加 MAGMOM 参数 晶胞中只有一种原子时: 晶胞中只有一种原子时: MAGMOM= 2m*3 铁磁计算 (2m 是总的原子数,3 是每个原子初始磁矩的大小) MAGMOM=m*-3 m*3 反铁磁计算 Types of spin ordering in perovskite oxides 铁磁 FM 反铁磁 AFM A-type: The intra-plane coupling is ferromagnetic while inter-plane coupling is antiferromagnetic. C-type: The intra-plane coupling is antiferromagnetic while inter-plane coupling is ferromagnetic. G-type: Both intra-plane and inter-plane coupling are antiferromagnetic. 磁结构的计算标题:[转载]如何用 VASP 计算铁磁、反铁磁和顺磁 浏览:95 评论:0 顺磁,意味进行 non-spin polarized 的计算,也就是 ISPIN=1。 铁磁,意味进行 spin-polarized 的计算,ISPIN=2,而且每个磁性原子的初 始磁矩设置为一样的值,也就是磁性原子的 MAGMOM 设置为一样的值。对 非磁性原子也可以设置成一样的非零值(与磁性原子的一样)或零,最后收 敛的结果,非磁性原子的 local 磁矩很小,快接近 0,很小的情况,很可能 意味着真的是非磁性原子也会被极化而出现很小的 local 磁矩。 反铁磁,也意味着要进行 spin-polarized 的计算,ISPIN=2,这是需采用反 铁磁的磁胞来进行计算,意味着此时计算所采用的晶胞不再是铁磁计算时的 最小原胞。比如对铁晶体的铁磁状态,你可以采用 bcc 的原胞来计算,但是 在进行反铁磁的 Fe 计算,这是你需要采用 sc 的结构来计算,计算的晶胞中 包括两个原子,你要设置一个原子的 MAGMOM 为正的,另一个原子的 MA GMOM 设置为负,但是它们的绝对值一样。因此在进行反铁磁的计算时, 应该确定好反铁磁的磁胞,以及磁序,要判断哪种磁序和磁胞是最可能的反 铁磁状态,那只能是先做好各种可能的排列组合,然后分别计算这些可能组 合的情况,最后比较它们的总能,总能最低的就是可能的磁序。同样也可以 与它们同铁磁或顺磁的进行比较。了解到该材料究竟是铁磁的、还是顺磁或 反铁磁的。 NUPDOWN 设置 设置为一个为非 0 而且非-1 的数, 就是相当是将在自洽迭代循环 的计算过程加了一个限制,即限制自旋向上和自旋向下的电子数之差为一个 固定的值,对应的就是给体系设置 设置一个固定的总磁矩。具体在程序里面实现 设置 起来, 量子化学的程序(比如 Gaussian)和广泛用于固体物理的第一原理程序(比 如 WIEN2K)略有不同。 不 亚铁, 结构, 原子, 磁矩, 整数 最近遇到一个问题,如何计算不同的磁性结构?如铁磁结构,反铁 不 磁结构,亚铁磁结构。vasp 教程上有说明,设置 ISPIN=2, MAGMOM=......,这样来说,计算铁磁结构只需要把各个原子的磁矩 方向设为一致,如,都为正值就可以了;计算反铁磁结构只需要把 不同子子晶格的原子磁矩方向设为相反,如,一正一负就可以了; 计算亚铁磁结构或许还要用到 NUPDOWN。 但是有个问题,我做了一个亚铁磁结构的计算,初始 MAGMOM 设置了 有正有负的原子磁矩,没有设置 NUPDOWN 值,而优化出来的结果里 设置 各个原子的磁矩却是按同一方向排列的!也就是说没有得到我想要 的亚铁磁结构。如果要设置 NUPDOWN 进行强制总磁矩的话,又不知 设置 不 道这个值该设成多大? 非常希望各位多给意见,只要是有相关内容的回帖都送金币。 收藏 分享 评分 回复 引用 TOP 定义: 自旋多重度(spin multiplicity)=2S+1, S=n*1/2,n 为单电子数。所以,关键是单电子 的数目是多少。 当有偶数个电子时,例如 O2,共有 16 个电子,那么单电子数目可能是 0,即 8 个 alpha 和 8 个 beta 电子配对,对应单重态,但是也可能是有 9 个 α 电子和 7 个 β 电子,那么能成对的 是 7 对, 还剩 2 个 α 没有配对, 于是 n=2,对应的是多重度 3。 同理还可以有多重度 5, 9, ... 7, 一般而言,是多重度低的能量低,最稳定,所以,一般来说,偶数电子的体系多重度就是 1。 但是也 有例外, 如果 O2 就是一个大家都知道的例子, 它的基态是三重态, 其单重态反而是激发态。 所以,总结一下,就是 电子数目是偶数,未成对电子数目 n=0,2,4,6,...自旋多重度是 1,3,5,7,... 电子数目是奇数,未成对电子数目 n=1,3,5,7,...自旋多重度是 2,4,6,8,... 多数情况是多重度低的能量低,有时(特别是有“磁”性的时候,例如顺磁的 O2,以及 Fe 啊 什么的) ,可能会高多重度的能量低。

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有