科学网

 找回密码
  注册

tag 标签: bands

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

vasp之mBJ计算能带
热度 8 plgongcat 2015-4-14 10:17
1.进行pbe scf calculation,得到IBZKPT #pbe scf calculation #INCAR SYSTEM = BP-00-vdw ENCUT = 500 ISTART = 0 NSW=0 ICHARG = 2 ISMEAR = 0; SIGMA = 0.05 PREC = Accurate #KPOINTS Automatic generation 0 M 20 20 20 0 0 0 2.mBJ band calculation是在PBE基础上的计算,设置ISTART=1,ICHARG=2。(因为hse时哈密顿量变了,也就是CHARG变了,不能从上一步pbe自洽计算读。将pbe得到的IBZKPT和自己需要计算能带Kpoints数据;合二为一,并设置Kpoints的weight是0) ! 注意:mBJ计算需要上一步PBE计算得到的WAVECAR。 #mBJ scf calculation #INCAR SYSTEM = BP-00-vdw ENCUT = 500 ISTART = 1 NSW=0 ICHARG = 2 ISMEAR = 0; SIGMA = 0.05 PREC = Accurate RWIGS = 1.233 LORBIT = 12 METAGGA = MBJ #IBZKPT与能带Kpoints数据合并 Automatically generated mesh 2000 Reciprocal lattice 0.02500000000000 0.02500000000000 0.02500000000000 4 0.07500000000000 0.02500000000000 0.02500000000000 8 0.12500000000000 0.02500000000000 0.02500000000000 8 0.17500000000000 0.02500000000000 0.02500000000000 8 0.22500000000000 0.02500000000000 0.02500000000000 8 0.27500000000000 0.02500000000000 0.02500000000000 8 0.32500000000000 0.02500000000000 0.02500000000000 8 ......... - 3. 上一步完成后得到EIGENVAL,对照KPOINTS删除多余(就是删除那些有weight的数据)并修改 EIGENVAL中K点的数目; 4. 下面就是画能带了,利用小程序 vasp_band.exe 可以得到bands图片(从DOSCAR中看Fermi能级) http://emuch.net/html/201203/4232787.html ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 下面是对单层使用mBJ能带修正,以黑磷为例 1. 需要得到块体黑磷的CMBJ=1.1562(自恰得到,不要在INCAR中人为添加) 2. 将上面的CMBJ数值添加到单层结构的INCAR之中。 如果,不按照上述方法,得到的能带数据发散、古怪!
个人分类: Vasp|28619 次阅读|16 个评论
IDL:删除影像中的波段
lixujeremy 2014-10-19 09:16
问题:删除影像中的若干波段,并生成新影像文件? 方法: IDL 。练习数据中, test.tif 包含 13 个波段,删除其中的 1~N (小于 13 )个波段,得到结果 result.tif 。 附上练习数据及代码( Practice.rar )。
个人分类: ENVI/IDL|5240 次阅读|0 个评论
声子谱的排序(对称性)
plgongcat 2014-9-16 20:36
Ph_by_magnitude.png Ph_by_band.png 对声子谱的重新排序,有两中方式:大小 对称性;(附件图片) 计算方法:pwscf软件 1.进行普通声子谱计算4步骤(pw.x,ph.x,q2r.x,matdyn.x),得到文件matdyn.modes文件(后处理备用) 2.进行投影计算: pw.x(scf),pw.x(bands),projwfc.x,在tmp中得到atomic_proj.xml文件(备用) 3.准备文件input_phon(输入文件) 4.用程序sort_phon.x处理上述备用文件,命令sort_phon.x input_phonoutput_phon 得到phonon_band.dat, phonon_magnitude.dat。 Sort_Phon.tar 5.gnuplot画图:命令 $plot ./phonon_magnitude.dat using 1:2 , u 1:3, u 1:4, u 1:5, u 1:6, u 1:7 6.如果图片中有尖锐峰,可以用插值法平滑处理(fitfunction.f90:命令是./a.outin.datout.dat) fitfunction.f90 in.dat links:http://blog.sciencenet.cn/blog-345795-422819.html
个人分类: pwscf|1 次阅读|0 个评论

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

GMT+8, 2024-6-11 08:53

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部