24小时热门版块排行榜    

查看: 879  |  回复: 4

gemucai

木虫 (正式写手)

[交流] 为何我计算出的铁磁性金属在高对称点的能级有跳变?已有2人参与

计算镍的能带,居然发现其L点和X点的能级有突变,比周围点的能量高出2/3,以至于画能带图的时候出现了两个尖峰!

还有我所用的Ni的SCF输入文件是从官网上下载的。

请问这是怎么回事呀?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
sunyang1988(金币+1): 谢谢提示 2011-04-23 19:33:15
贴出来看看是怎么个突变法?我估计是你在计算能带的时候忘记设置某些参数了
2楼2011-04-23 19:09:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phys

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
sunyang1988(金币+1): 谢谢提示 2011-04-23 23:32:44
我估计是对称性的K点取的有问题
JustPhysics
3楼2011-04-23 21:32:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

2楼:
   最后通过处理得出的数据大概就是这样了,图中能量跳变的点是X点,我一开始说错了,能级变化应该只超出1/3。
   我的输入文件:
    &control
    prefix='nickel'
    outdir='temporary directory for large files'
    pseudo_dir = 'directory where pp-files are kept',
/
&system
    ibrav=2, celldm(1) =6.48, nat=1, ntyp=1,
    ecutwfc = 24.0, ecutrho = 288.0,
    occupations='smearing', smearing='marzari-vanderbilt', degauss=0.02
    nspin = 2,  starting_magnetization(1)=0.7,
/
&electrons
/
ATOMIC_SPECIES
Ni 58.69 NiUS.RRKJ3.UPF
ATOMIC_POSITIONS
Ni 0.0 0.0 0.0
……
K点不管是自动布点还是用官网上的布点都会出现跳变。
计算能带时,高对称K点我只选了(0.0 0.0 0.0)、(0.0 0.0 1.0)、(0.5 0.5 0.5)

4楼2011-04-24 08:55:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

对了,我的OUTDIR 和赝势存放地址都已经设置好了的,我贴出来的是下载的原文件
5楼2011-04-24 08:56:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见