24小时热门版块排行榜    

查看: 2395  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wydzyx2003

金虫 (正式写手)

[求助] 画能带时遇到奇怪的问题---帮忙看下吧

请问你们有没有遇到过这种情况?

在计算能带的时候,不管怎么改变KPOINTS文件里的K点路径,计算出来的结果跟第一次计算出来的是一模一样的。而且第一次计算的也不是按KPOINTS文件里的路径计算的(并且画出来的能带图很难看,明显不对)。仔细对照以前算的(3*3*3)的发现,现在算的结构优化文件夹和静态计算文件夹中都少了tmp.dat和tmp.dat2文件(现在计算的是2*2*2的超胞)。

KPOINTS文件:
k-points along high symmetry lines
10 ! 10
Line-Mode
rec

0 0 0 !G
0.3333 0.3333 0 !K

0.3333 0.3333 0 !K
0.3333 0.3333 0.5 !H

0.3333 0.3333 0.5 !H
0 0 0.5 !A

0 0 0.5 !A
0 0 0 !G

0 0 0 !G
0 0 0.5 !A

0 0 0.5 !A
0.5 0 0.5 !L

0.5 0 0.5 !L
0.5 0 0 ! M

0.5 0 0 !M
0 0 0 !G

能带图如下

画出来不太对劲的能带图



正确的能带图
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wydzyx2003

金虫 (正式写手)

引用回帖:
12楼: Originally posted by seaso at 2012-01-09 23:49:57:
把incar和outcar的信息贴出来。
[ 发自手机版 http://muchong.com/3g ]

INCAR 文件  (OUTCAR文件太大以附件形式上传了)
YSTEM =4
Startparameter for this run:
   NWRITE =      2    write-flag & timer
   PREC   = high      medium, high low
   ISTART =      0    job   : 0-new  1-cont  2-samecut
   ICHARG =      11    charge: 1-file 2-atom 10-const
#   ISPIN  =      1    spin polarized calculation?
#   LNONCOLLINEAR =      F non collinear calculations
#   LSORBIT =      F    spin-orbit coupling
#   INIWAV =      1    electr: 0-lowe 1-rand  2-diag
#   LASPH  =      F    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.  17.77 17.77 28.48*2*pi/ulx,y,z
   ENINI  =  500.0     initial cutoff
   ENAUG  =  605.4 eV  augmentation charge cutoff
   NELM   =     80;   NELMIN=  2; NELMDL=  0     # of ELM steps
   EDIFF  = 0.1E-03   stopping-criterion for ELM
   LREAL  =      T    real-space projection
#   LCOMPAT=      F    compatible to vasp.4.4
#   LREAL_COMPAT= F    compatible to vasp.4.5.1-3
#   GGA_COMPAT  = T    GGA compatible to vasp.4.4-vasp.4.6
#   LMAXPAW     = -100 max onsite density
#   LMAXMIX     =    2 max onsite mixed and CHGCAR
#   VOSKOWN=      0    Vosko Wilk Nusair interpolation
#   ROPT   =   -0.00040  -0.00040
Ionic relaxation
#   EDIFFG = 0.1E-02   stopping-criterion for IOM
#   NSW    =      100    number of steps for IOM
#   NBLOCK =      1;   KBLOCK =      1    inner block; outer block
#   IBRION =     1    ionic relax: 0-MD 1-quasi-New 2-CG
#   NFREE  =      0    steps in history (QN), initial steepest desc. (CG)
#   ISIF   =      2    stress and relaxation
   IWAVPR =     10    prediction:  0-non 1-charg 2-wave 3-comb
   ISYM   =      2    0-nonsym 1-usesym 2-fastsym
   LCORR  =      T    Harris-Foulkes like correction to forces

   POTIM  =   0.50    time-step for ionic-motion
#   TEIN   =    0.0    initial temperature
#   TEBEG  =    0.0;   TEEND  =   0.0 temperature during run
#   SMASS  =  -3.00    Nose mass-parameter (am)
#   estimated Nose-frequenzy (Omega)   =  0.10E-29 period in steps =****** mass=  -0.217E-26a.u.
#   NPACO  =    256;   APACO  = 16.0  distance and # of slots for P.C.
#   PSTRESS=    0.0 pullay stress

   #RWIGS  =  0.74   1.757   1.402
#   NELECT =     972.0000    total number of electrons
#   NUPDOWN=      -1.0000    fix difference up-down

DOS related values:
#   EMIN   =  10.00;   EMAX   =-10.00  energy-range for DOS
   ISMEAR =    0 ;   SIGMA  =   0.05  broadening in eV -4-tet -1-fermi 0-gaus

Electronic relaxation 2 (details)
   IALGO  =     68    algorithm
   LDIAG  =      T    sub-space diagonalisation
   IMIX   =      4    mixing-type and parameters
     AMIX     =   0.40;   BMIX     =  1.00
     AMIX_MAG =   1.60;   BMIX_MAG =  1.00
     AMIN     =   0.10
     WC   =   100.;   INIMIX=   1;  MIXPRE=   1

Intra band minimization:
#   WEIMIN = 0.0000     energy-eigenvalue tresh-hold
#   EBREAK =  0.43E-07  absolut break condition
#   DEPER  =   0.30     relativ break condition

#   TIME   =   0.10     timestep for ELM

Second variation
#   LSECVAR=     F    do a second variation

Write flags
   LWAVE  =      F    write WAVECAR
   LCHARG =      T    write CHGCAR
   LVTOT  =      F    write LOCPOT, local potential
   LELF   =      F    write electronic localiz. function (ELF)
------------------------------------------------------------------------------------------------

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : OUTCAR
  • 2012-01-10 11:49:16, 534.45 K
14楼2012-01-10 11:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

wydzyx2003

金虫 (正式写手)

注:两个能带图的横坐标也不一样。应该是一样的都到 80 才对吧。
2楼2012-01-05 12:16:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WDD880227

木虫 (正式写手)

LZ用什么方法计算能带图的?
smiling&fighting...
3楼2012-01-05 14:54:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinhua1983

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
WDD880227(金币+1): 感谢交流~~ 2012-01-07 08:57:50
wydzyx2003(金币+2): 谢谢,我先试下,问题解决的话再给您加金币。 2012-01-08 12:14:10
Line-Mode可能不太好用,建议手动定义K网格试一下!
4楼2012-01-07 08:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +8 majunge000 2026-08-11 10/500 2026-08-15 18:55 by Equinoxhua
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见