应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp gk.x文件生成KPOINTS出错
251/3返回列表123下一页
查看: 5977  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mengfc

金虫 (正式写手)

[求助] vasp gk.x文件生成KPOINTS出错 已有1人参与

请教各位，我自己写了一个syml文件。但用编译的gk.x执行时只出现如下内容：
“
k-points along high symmetry lines
-1225672974
Reciprocal
”
请问这是什么问题?
我的syml文件如下：
2
25   
G 0.0    0.0    0.0
M 0.0    0.5    0.0
50.041423917  0.005031325  0.000000000     0.019983444 -0.000002011  0.000000000     
0.000248027  2.465241129  0.000000000    -0.000040784  0.405639837  0.000000000     
0.000000000  0.000000000 29.937787689     0.000000000  0.000000000  0.033402602     
-8.0  8.0
-3.7306
gk.f文件如下，我是使用 ifort -O gk.f命令进行编译的。
c    +---------------------------------------------------------
c    |  For generating k-points along the high-symmetry lines in
c    |  Brillouin zone and for calculate band-structures !
c    |
c    |  OK for "VASP" package
c    |                                Modified by   Hou Zhufeng
c    +----------------------------------------------------------
C    ---------'syml'---------
c    6                     : nhighk
c    20  20  20  10  20    : ndiv(i)
c    X 0.5    0.0    0.5   : labhk(1),phighk(1,1),........
c    G 0.0    0.0    0.0
c    L 0.5    0.5    0.5
c    W 0.5    0.25   0.75
c    K 0.375  0.375  0.75
c    G 0.0    0.0    0.0
cc   0.000  2.449  2.449  -0.204158033  0.204158033  0.204158033
cc   2.449  0.000  2.449   0.204158033 -0.204158033  0.204158033
cc   2.449  2.449  0.000   0.204158033  0.204158033 -0.204158033   
c    -20.0  20.0           :emin, emax
c    0.0                   :efermi
c    direct & reciprocal lattice vectors over 'emin, emax' line
C    -----------------------
c     max k-points = 200
      program  gk
      implicit real*8 (a-h,o-z)
      character*2 labhk
      dimension tkpt(200,3),pk(200,3),phighk(10,3)
      dimension disk(200),dish(10),labhk(10)
      dimension ndiv(10)
c
      open(5,file='syml',status='old')
      open(7,file='inp.kpt')
      open(8,file='KPOINTS')
c
      read(5,*) nhighk
      read(5,*) (ndiv(i),i=1,nhighk-1)

      do i=1,nhighk-1
            ntkp=ntkp+ndiv(i)
      enddo
      ntotkpt=ntkp+1  

      if(nhighk>10)then
        write(*,*)'Number of high-symmetry k points must < 10!'
        STOP
      endif
      if(ntotkpt>200)then
        write(*,*)'Total number of k points must <= 200!'
        STOP
      endif

      do i=1, nhighk
         read(5,*)  labhk(i),(phighk(i,j),j=1,3)
      enddo
      write(*,*) (labhk(i),i=1,nhighk)
c
c----- generating k-points along high symmetric lines --------
c
c
      pk(1,1)=phighk(1,1)
      pk(1,2)=phighk(1,2)
      pk(1,3)=phighk(1,3)
      ii=1
       do i = 2, nhighk
         delx = (phighk(i,1) - phighk(i-1,1))/float(ndiv(i-1))
         dely = (phighk(i,2) - phighk(i-1,2))/float(ndiv(i-1))
         delz = (phighk(i,3) - phighk(i-1,3))/float(ndiv(i-1))
           do j=1, ndiv(i-1)
              ii = ii + 1
              pk(ii,1) = pk(ii-1,1) + delx
              pk(ii,2) = pk(ii-1,2) + dely
              pk(ii,3) = pk(ii-1,3) + delz
           enddo
       enddo
c
      write(8,10)'k-points along high symmetry lines'
  10  format(A34)
      write(8,*)ntotkpt
      write(8,'(A10)')'Reciprocal'
         weight=1.d0
         do i=1,ntotkpt
         write(7,200) pk(i,1),pk(i,2),pk(i,3),weight
         write(8,200) pk(i,1),pk(i,2),pk(i,3),weight
         enddo
200  format(3F10.6,F6.2)
      stop
      end
c----------------------- end ---------------------------
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-04 23:12:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Mr_WuSH

金虫 (小有名气)
不受意思··记错了··使用这个编译
g77 -o gk.x gk.f
赞一下(2人) 回复此楼
15楼2012-06-11 16:23:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr_WuSH

金虫 (小有名气)
★
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-10 10:58:22
重新用gfortran编译就好使了····
赞一下(1人) 回复此楼
7楼2012-06-09 20:41:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
mengfc: 金币+15, ★★★很有帮助, 谢谢您 2012-05-05 00:03:32
liliangfang: 金币+2, 谢谢指教 2012-05-05 09:54:57

两个高对称k点，G 0.0    0.0    0.0    M 0.0    0.5    0.0
中间插入25个点，excel表中一下就能生成~

更简单的，vasp的KPOINTS都已经支持LINE模式了~
赞一下 回复此楼
2楼2012-05-04 23:45:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gfei2008

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
mengfc: 金币+5, ★有帮助, 谢谢您 2012-05-05 00:04:31
liliangfang: 金币+2, 谢谢指教 2012-05-05 09:55:16
你这样编译出之后是不是有个a.out的可执行文件，那个就是用来生产KPOINTS的程序，使用方法是：./a.out
要么你就是ifort -o gk gk.f这样编译。然后执行./gk就可以生成你要的KPOINTS
赞一下 回复此楼
3楼2012-05-04 23:52:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfc

金虫 (正式写手)
引用回帖:
2楼: Originally posted by fzx2008 at 2012-05-04 23:45:55:

两个高对称k点，G 0.0    0.0    0.0    M 0.0    0.5    0.0
中间插入25个点，excel表中一下就能生成~

更简单的，vasp的KPOINTS都已经支持LINE模式了~

唉。。。我也决定用LINE模式了。。。就是用pbnd.x程序时还是要用到syml文件。不知道怎么弄了。还请问您一下，计算能带是非自洽计算，是不是因为读取电荷密度，而不是再重新计算？我在计算能带时的INCAR中ediff和ediffg两个参数跟几何优化时的一样？
赞一下 回复此楼
4楼2012-05-05 00:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfc

金虫 (正式写手)
引用回帖:
3楼: Originally posted by gfei2008 at 2012-05-04 23:52:32:
你这样编译出之后是不是有个a.out的可执行文件，那个就是用来生产KPOINTS的程序，使用方法是：./a.out
要么你就是ifort -o gk gk.f这样编译。然后执行./gk就可以生成你要的KPOINTS

谢谢您，但执行之后KPOINTS里是
“
k-points along high symmetry lines
-1225672974
Reciprocal
”
赞一下 回复此楼
5楼2012-05-05 00:04:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yxz476440843

铁虫 (小有名气)
引用回帖:
5楼: Originally posted by mengfc at 2012-05-05 00:04:22
谢谢您，但执行之后KPOINTS里是
“
k-points along high symmetry lines
-1225672974
Reciprocal
”...

你好，我也遇到了同样的问题，请问该怎么解决啊，谢谢!
赞一下 回复此楼
6楼2012-06-09 19:32:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr_WuSH

金虫 (小有名气)
重新用gfortran编译不好使·？
赞一下 回复此楼
8楼2012-06-11 09:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfc

金虫 (正式写手)
引用回帖:
8楼: Originally posted by Mr_WuSH at 2012-06-11 09:10:07
重新用gfortran编译不好使·？

谢谢您，这个用gfortran重新编译之后也不太行，这是什么原因呢？
赞一下 回复此楼
9楼2012-06-11 10:32:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr_WuSH

金虫 (小有名气)
我之前也遇到一样的问题··
gfortran编译后就好使了···
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : gk.f
    • 2012-06-11 11:34:59, 2.75 K
10楼2012-06-11 11:36:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mengfc 的主题更新
251/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 00:25 by xxniao123
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 国自然面上复盘~欢迎讨论 +4 晴天加油 2026-08-26 5/250 2026-08-26 21:27 by 木风9012
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭