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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp gk.x文件生成KPOINTS出错
251/3返回列表123下一页
查看: 5607  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
关闭
请填处理意见
关闭