xiaoqiugood的个人博客分享 http://blog.sciencenet.cn/u/xiaoqiugood

博文

加SOC对晶格参数的影响

已有 7124 次阅读 2016-12-14 15:54 |个人分类:结构优化|系统分类:科研笔记

关注:

1) 磁性优化对晶格参数的影响

2) 考虑SOC后优化,对晶格参数的影响



SOC优化耗时太长:

1) 或可在正常优化得到结构的基础上,加入SOC做一次自洽,以此次自洽得到的能量为加SOC优化后的能量值;【误差可在一百至几百meV左右】

2) 或可在考虑铁磁/反铁磁优化得到结构的基础上,加入SOC做一次自洽,以此次自洽得到的能量为加SOC优化后的能量值;

3) 考虑SOC时,自旋默认为1。-- ISPIN  =      1




问题:

不同参数下优化,晶格常数变化较大
(1)加自旋轨道耦合(SOC)-(铁磁)优化
Pnma结构:a=7.5965,b=4.9924,c=3.8962

(2)加-(铁磁)优化
Pnma结构:a=6.9756,b=4.5444,c=3.5738

(3)正常优化
Pnma结构:a=6.9723,b=4.5421,c=3.5716



从OUTCAR中看初始参数设置: grep 'Startparameter'


Startparameter for this run   :



OUTCAR摘录:SOC-AFM

Startparameter for this run:
  NWRITE =      2    write-flag & timer
  PREC   = accura    normal or accurate (medium, high low for compatibility)
  ISTART =      0    job   : 0-new  1-cont  2-samecut
  ICHARG =      2    charge: 1-file 2-atom 10-const
  ISPIN  =      1    spin polarized calculation?
  LNONCOLLINEAR =      T non collinear calculations
  LSORBIT =      T    spin-orbit coupling
  INIWAV =      1    electr: 0-lowe 1-rand  2-diag
  LASPH  =      T    aspherical Exc in radial PAW
  METAGGA=      F    non-selfconsistent MetaGGA calc.

Electronic Relaxation 1
  ENCUT  =  500.0 eV  36.75 Ry    6.06 a.u.  14.09  9.21  7.29*2*pi/ulx,y,z
  ENINI  =  500.0     initial cutoff
  ENAUG  =  455.5 eV  augmentation charge cutoff
  NELM   =     60;   NELMIN=  2; NELMDL= -5     # of ELM steps
  EDIFF  = 0.1E-03   stopping-criterion for ELM
  LREAL  =      F    real-space projection
  NLSPLINE    = F    spline interpolate recip. space projectors
  LCOMPAT=      F    compatible to vasp.4.4
  GGA_COMPAT  = T    GGA compatible to vasp.4.4-vasp.4.6
  LMAXPAW     = -100 max onsite density
  LMAXMIX     =    6 max onsite mixed and CHGCAR
  VOSKOWN=      1    Vosko Wilk Nusair interpolation
  ROPT   =    0.00000   0.00000




对XH3的YF3结构check时发现,FM与SOC-FM结构晶格常数不大


FM

XH3                                  
  1.00000000000000    
    6.6126159350751834    0.0000000000000000    0.0000000000000000
    0.0000000000000000    6.8490055081721524    0.0000000000000000
    0.0000000000000000    0.0000000000000000    3.8177718559422864
  H    X
   12     4
Direct
 0.5000174978855227  0.2500000000000000  0.4999984464310501

.....


SOC-FM

XH3                                  
  1.00000000000000    
    6.6181079580488928    0.0000000000000000    0.0000000000000000
    0.0000000000000000    6.8537718526809712    0.0000000000000000
    0.0000000000000000    0.0000000000000000    3.8232081204009312
  H   X
   12     4
Direct
 0.5002032561368708  0.2499999431915754  0.4999995155845411
 0.4997967438631434  0.7499999431915754  0.5000004844154589

摘录:

http://wenku.baidu.com/link?url=JJ4iSFbGgHFdcyPMqVM3fcGnpss5KZKlnt83EaGhtAIUyJ1NecDxoPPtTy0lPJClr0qkyn2YE1HlOzn0fZriON3JL7OZPw6ciSF7L4z2jtO





https://blog.sciencenet.cn/blog-567091-1020684.html

上一篇:绝缘体:Mott-Hubbard insulator and charge-transfer insulator
下一篇:再谈Phonopy的安装:学而时习之 python3.6 setup.py install
收藏 IP: 125.67.185.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-26 22:37

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部