版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
散金: 142
红花: 2
帖子: 360
在线: 150.7小时
虫号: 1415189
注册: 2011-09-24
性别: MM
专业: 凝聚态物性 II ：电子结构

[求助] 如何将VASP中的CONCAR导入到MS中？？

我做的是金属间化合物的表面催化裂解，表面取得是111，构建7层原子和真空层后需要用VASP进行弛豫，以取得稳定结构，但遇到的难题是如何将生成的CONCAR导入到MS中，好像可以用Compaq Visual Fortran这个编译器实现，但由于没学过Fortan语言，不会编程，找到的一些编好的程序也看不懂，下面贴出程序以及CONCAR，还请大师们帮我解释下程序的大致内容，以及那些需要修改的地方，不胜感激！！！

MODULE PREC
   INTEGER, PARAMETER :: q =SELECTED_REAL_KIND(10)
   INTEGER, PARAMETER :: qs=SELECTED_REAL_KIND(5)
   END MODULE

program constucture_of_supercell
   USE PREC
real(q)  deta_d(151)
   real(q) pos(151,3),pos_extend(151,3),pos1(151,3)
real(q) a1,a2,a3,b1,b2,b3,c1,c2,c3,a,x,y,z

integer atom_num,supercellnum_a,supercellnum_b,supercellnum_c
   integer pos_num,flag
integer atom_num1,atom_num2,atom_num3,atom_num4,atom_begin,
   & atom_end
   character*10 temp,temp_Coordination
   supercellnum_a=1
supercellnum_b=1
supercellnum_c=1
   pos_num=1
flag=1

   open(1,file='4.txt',status='old')
   open(2,file='visions-2.xyz',status='unknown')
   open(3,file='deta.out',status='unknown')
write(*,*) "dasafa"

   read(1,*) temp

   read(1,*) a

   read(1,*) a1,a2,a3
read(1,*) b1,b2,b3
read(1,*) c1,c2,c3
! read(1,*) temp
read(1,*) atom_num1,atom_num2,atom_num3 !,atom_num4

   read(1,*) temp_Coordination
   write(*,*) "2dasafa"
atom_num=atom_num1+atom_num2+atom_num3+atom_num4

   write(2,*)  (atom_num*supercellnum_a*supercellnum_b*
   & supercellnum_c)
      write(*,*) "3dasafa"
   write(2,*) "nomame"


   do i=1,atom_num
   read(1,*) x,y,z
   write(*,*) x,y,z,i
   pos(i,1)=x
   pos(i,2)=y
   pos(i,3)=z
   pos1(i,1)=x
   pos1(i,2)=y
   pos1(i,3)=z
enddo


   do i=1,atom_num-1
   deta_d(i)=(pos(atom_num-i+1,3)*8-pos(atom_num-i,3)*8-
   &8.574929257125441E-002)/8.574929257125441E-002


enddo
   if((temp_Coordination(1:1).EQ.'D').OR.
   & (temp_Coordination(1:1).EQ.'d'))then
write(*,*) "This points denote by direct coordinates!"
write(*,*) "convert to cartsian:"

   do i=1,atom_num
   pos1(i,1)=pos(i,1)*a1+pos(i,2)*b1+pos(i,3)*c1
   pos1(i,2)=pos(i,1)*a2+pos(i,2)*b2+pos(i,3)*c2
   pos1(i,3)=pos(i,1)*a3+pos(i,2)*b3+pos(i,3)*c3
   write(*,*) pos1(i,1)/1.5,pos1(i,2)/1.5,pos1(i,3)/1.5,0
   write(3,*) pos1(i,3)/1.5,0
write(3,*) pos1(i,3)/1.5,1
   write(3,*) 0.0,0

enddo

endif

! do i=1,atom_num
! write(2,*) "cu",pos(i,1),pos(i,2),pos(i,3),pos_num,flag
!  pos_num=pos_num+1
! enddo

   do jb=0,supercellnum_b-1
   do ja=0,supercellnum_a-1
         do jc=0,supercellnum_c-1

   atom_begin=1
atom_end=atom_num1

do i=atom_begin,atom_end
   pos_extend(i,1)=pos1(i,1)+a1*ja+b1*jb+c1*jc
   pos_extend(i,2)=pos1(i,2)+a2*ja+b2*jb+c2*jc
   pos_extend(i,3)=pos1(i,3)+a3*ja+b3*jb+c3*jc
   write(8,*) pos_extend(i,1),pos_extend(i,2),pos_extend(i,3),"T T T"
   pos_extend(i,1)=pos_extend(i,1)*a
   pos_extend(i,2)=pos_extend(i,2)*a
   pos_extend(i,3)=pos_extend(i,3)*a

   write(2,*) "Au",pos_extend(i,1),pos_extend(i,2),pos_extend(i,3)

write(31,*) pos_extend(i,1)/a,pos_extend(i,2)/a,pos_extend(i,3)/a

pos_num=pos_num+1




   enddo

   atom_begin=atom_num1+1
atom_end=atom_num2+atom_num1

do i=atom_begin,atom_end
   pos_extend(i,1)=pos1(i,1)+a1*ja+b1*jb+c1*jc
   pos_extend(i,2)=pos1(i,2)+a2*ja+b2*jb+c2*jc
   pos_extend(i,3)=pos1(i,3)+a3*ja+b3*jb+c3*jc
   write(8,*) pos_extend(i,1),pos_extend(i,2),pos_extend(i,3),"T T T"
   pos_extend(i,1)=pos_extend(i,1)*a
   pos_extend(i,2)=pos_extend(i,2)*a
   pos_extend(i,3)=pos_extend(i,3)*a

   write(2,*) "C",pos_extend(i,1),pos_extend(i,2),pos_extend(i,3)
write(31,*) pos_extend(i,1)/a,pos_extend(i,2)/a,pos_extend(i,3)/a
pos_num=pos_num+1
   enddo

   atom_begin=atom_num2+atom_num1+1
atom_end=atom_num3+atom_num2+atom_num1

do i=atom_begin,atom_end
   pos_extend(i,1)=pos1(i,1)+a1*ja+b1*jb+c1*jc
   pos_extend(i,2)=pos1(i,2)+a2*ja+b2*jb+c2*jc
   pos_extend(i,3)=pos1(i,3)+a3*ja+b3*jb+c3*jc
   write(8,*) pos_extend(i,1),pos_extend(i,2),pos_extend(i,3),"T T T"
   pos_extend(i,1)=pos_extend(i,1)*a
   pos_extend(i,2)=pos_extend(i,2)*a
   pos_extend(i,3)=pos_extend(i,3)*a

   write(2,*) "O",pos_extend(i,1),pos_extend(i,2),pos_extend(i,3)
write(31,*) pos_extend(i,1)/a,pos_extend(i,2)/a,pos_extend(i,3)/a
pos_num=pos_num+1

   enddo

atom_begin=atom_num3+atom_num2+atom_num1+1
atom_end=atom_num

     do i=atom_begin,atom_end
   pos_extend(i,1)=pos1(i,1)+a1*ja+b1*jb+c1*jc
   pos_extend(i,2)=pos1(i,2)+a2*ja+b2*jb+c2*jc
   pos_extend(i,3)=pos1(i,3)+a3*ja+b3*jb+c3*jc
   write(8,*) pos_extend(i,1),pos_extend(i,2),pos_extend(i,3),"T T T"
   pos_extend(i,1)=pos_extend(i,1)*a
   pos_extend(i,2)=pos_extend(i,2)*a
   pos_extend(i,3)=pos_extend(i,3)*a

   write(2,*) "Ag",pos_extend(i,1),pos_extend(i,2),pos_extend(i,3)
pos_num=pos_num+1
   enddo

      enddo
      enddo
enddo

   write(*,*)  "the task finished"
   do i=1,atom_num-1
      write(*,*)  "deta_d",i,i+1,"=",deta_d(i)
!    write(3,'(A6,I2.2,I2.2,A6,E)')  "deta_d",i,i+1,"=",deta_d(i)
   !  write(*,'(A6,I2.2,I2.2,A6,E)')  "deta_d",i,i+1,"=",deta_d(i)
   enddo
end


CONCAR

  Ni3Fe-111
1.000000000000000
   4.327023630626480    -2.498208257932070    0.000000000000001
   0.000000000000000    4.996416515864140    0.000000000000000
   0.000000000000000    0.000000000000000    24.238671006281699
  21  7  0
Direct
0.0000000000000000 0.5000000000000000 0.0000000000000001  F F F
0.6666666666666670 0.8333333333333330 0.0841538919571259  F F F
0.3333333333333330 0.1666666666666670 0.1683077839142520  F F F
-0.0000000000000000 0.4999999999999999 0.2524616758713780  T T T
0.6666666666666670 0.8333333333333330 0.3366155678285030  T T T
0.3333333333333330 0.1666666666666670 0.4207694597856290  T T T
-0.0000000000000001 0.5000000000000000 0.5049233517427550  T T T
0.5000000000000000 0.5000000000000000 0.0000000000000001  F F F
0.1666666666666670 0.8333333333333329 0.0841538919571259  F F F
0.8333333333333330 0.1666666666666671 0.1683077839142520  F F F
0.5000000000000000 0.5000000000000000 0.2524616758713780  T T T
0.1666666666666670 0.8333333333333329 0.3366155678285030  T T T
0.8333333333333330 0.1666666666666660 0.4207694597856290  T T T
0.5000000000000000 0.5000000000000000 0.5049233517427550  T T T
0.5000000000000000 0.0000000000000001 0.0000000000000001  F F F
0.1666666666666670 0.3333333333333329 0.0841538919571258  F F F
0.8333333333333330 0.6666666666666670 0.1683077839142520  F F F
0.5000000000000000  -0.0000000000000000 0.2524616758713780  T T T
0.1666666666666670 0.3333333333333329 0.3366155678285030  T T T
0.8333333333333330 0.6666666666666670 0.4207694597856290  T T T
0.5000000000000000  -0.0000000000000001 0.5049233517427550  T T T
0.0000000000000001 0.0000000000000002  -0.0000000000000000  F F F
0.6666666666666661 0.3333333333333330 0.0841538919571259  F F F
0.3333333333333330 0.6666666666666670 0.1683077839142520  F F F
-0.0000000000000000 0.0000000000000000 0.2524616758713780  T T T
0.6666666666666670 0.3333333333333330 0.3366155678285030  T T T
0.3333333333333330 0.6666666666666670 0.4207694597856290  T T T
-0.0000000000000002  -0.0000000000000001 0.5049233517427550  T T T

回复此楼

» 本帖已获得的红花（最新10朵）

lijuanshu

» 猜你喜欢

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

VASP并行计算时，提示“超过CPU时限”而退出计算。请问如何解决？已经有11人回复
Vasp结算结果中如何计算结构因子已经有5人回复
vasp计算中如何比较掺杂前后结合能已经有13人回复
如何在vasp中设定离子的价态已经有18人回复
vasp计算中，如何选择正确的K点已经有7人回复
VASP中如何在体系中引入电子和空穴已经有15人回复
Vasp中删除一个键输入文件如何表现出来？已经有8人回复
当体系中包含两种材料时，如何用VASP计算其中一种材料的功函？已经有11人回复
【求助】vasp中COOP如何理解及输出文件如何处理？已经有12人回复
【其他】关于MS和VASP的比较已经有6人回复
【求助】MS中的fractional坐标和VASP中的direct与cart坐标的区别已经有8人回复
【求助成功】VASP中，如何固定磁矩做结构优化已经有14人回复

第一性原理VASPLinux这些都要学啊

1楼 2012-05-07 16:32:27

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

回帖支持 ( 显示支持度最高的前 50 名 )

mengfc

金虫 (正式写手)

1ST强帖: 1
应助: 20 (小学生)
金币: 1301.2
散金: 1807
红花: 5
帖子: 927
在线: 524.4小时
虫号: 818484
注册: 2009-07-30
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
lijuanshu: 金币+3, ★★★很有帮助 2012-05-07 17:07:27
fzx2008: 金币+1, 谢谢指教 2012-05-07 18:34:18

用VESTA打开CONTCAR，然后保存成cif格式，就可以导入MS了

赞一下(1人)

回复此楼

4楼2012-05-07 17:01:42

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

普通回帖

huangll99

木虫 (职业作家)

应助: 43 (小学生)
金币: 775.9
散金: 2040
红花: 20
沙发: 8
帖子: 4051
在线: 1447.4小时
虫号: 1521828
注册: 2011-12-03
性别: GG
专业: 计算机软件

21 7 0
楼主，我非常不解这个0是干嘛的？

赞一下

回复此楼

2楼2012-05-07 16:47:44

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

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
散金: 142
红花: 2
帖子: 360
在线: 150.7小时
虫号: 1415189
注册: 2011-09-24
性别: MM
专业: 凝聚态物性 II ：电子结构

好像程序里是三种元素，就加了个0

赞一下

回复此楼

第一性原理VASPLinux这些都要学啊

3楼2012-05-07 16:53:23

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

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
散金: 142
红花: 2
帖子: 360
在线: 150.7小时
虫号: 1415189
注册: 2011-09-24
性别: MM
专业: 凝聚态物性 II ：电子结构

VESTA怎么用啊，能不能说的具体点，我对这些东西都不太懂,谢谢啦

赞一下

回复此楼

第一性原理VASPLinux这些都要学啊

5楼2012-05-07 17:05:59

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

mengfc

金虫 (正式写手)

1ST强帖: 1
应助: 20 (小学生)
金币: 1301.2
散金: 1807
红花: 5
帖子: 927
在线: 524.4小时
虫号: 818484
注册: 2009-07-30
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
fzx2008: 金币+1, 谢谢回帖！ 2012-05-07 18:33:44

引用回帖:

5楼: Originally posted by lijuanshu at 2012-05-07 17:05:59:
VESTA怎么用啊，能不能说的具体点，我对这些东西都不太懂,谢谢啦

把CONTCAR拖VESTA的界面里就可以了，需要的话可以问我要

赞一下

回复此楼

6楼2012-05-07 17:10:49

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

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
散金: 142
红花: 2
帖子: 360
在线: 150.7小时
虫号: 1415189
注册: 2011-09-24
性别: MM
专业: 凝聚态物性 II ：电子结构

内容已删除

回复此楼

第一性原理VASPLinux这些都要学啊

7楼2012-05-07 17:29:20

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

youzhizhe

荣誉版主 (职业作家)

有志者

1ST强帖: 8
应助: 17 (小学生)
贵宾: 1.696
金币: 8402
散金: 12847
红花: 49
沙发: 14
帖子: 4491
在线: 763.2小时
虫号: 667979
注册: 2008-12-04
性别: GG
专业: 材料物理化学
管辖: 第一性原理

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教 2012-05-07 18:34:08

引用回帖:

4楼: Originally posted by mengfc at 2012-05-07 17:01:42:
用VESTA打开CONTCAR，然后保存成cif格式，就可以导入MS了

正解。VESTA下载可以在本版资源区找，最新版的请直接登录其官方网站http://jp-minerals.org/vesta/en/。

赞一下

回复此楼

微博http://weibo.com/10986069

8楼2012-05-07 18:09:09

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

jixiao0526

金虫 (小有名气)

应助: 12 (小学生)
金币: 808.3
散金: 158
红花: 2
帖子: 151
在线: 235.6小时
虫号: 1407476
注册: 2011-09-19
性别: GG
专业: 电磁场与波

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2012-05-08 08:40:36

VESTA很好用，直接把CONTCAR拖进去，然后save as 选*.cif就行了

赞一下

回复此楼

9楼2012-05-07 23:30:55

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

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
散金: 142
红花: 2
帖子: 360
在线: 150.7小时
虫号: 1415189
注册: 2011-09-24
性别: MM
专业: 凝聚态物性 II ：电子结构

送鲜花一朵

VESTA我学会用了，先将CONTCAR保存为COMTCAR.vasp，在VESTA中打开，再Expose data，保存成CONTCAR.cif，就可以在MS中打开了，不过有人知道上面的程序是怎么回事吗，继续求救啊！！

赞一下

回复此楼

第一性原理VASPLinux这些都要学啊

10楼2012-05-08 09:19:23

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

相关版块跳转我要订阅楼主 lijuanshu 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55	医学老男孩 2026-08-13	123/6150	2026-08-19 23:44 by stan_et
[基金申请] 时间戳变了，能看出什么问题？ +12	基诺咪客 2026-08-17	16/800	2026-08-19 22:04 by kk 的小木虫
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5	瞬息宇宙 2026-08-15	5/250	2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4	yuleib84 2026-08-19	5/250	2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14	布布和一二 2026-08-19	15/750	2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10	wangze12014 2026-08-14	12/600	2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6	wangzilk 2026-08-18	8/400	2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6	Tide man 2026-08-14	7/350	2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6	CrisMessi 2026-08-18	6/300	2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10	yuleib84 2026-08-18	11/550	2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +5	Shxjjxjkx 2026-08-18	5/250	2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4	wwm09 2026-08-17	6/300	2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14	archvillain 2026-08-15	26/1300	2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护祝所有人高中 +8	gjjjzhong 2026-08-18	9/450	2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13	Ldrop2023 2026-08-13	16/800	2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4	且听虎啸 2026-08-17	5/250	2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6	angus9576 2026-08-17	11/550	2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8	zju2000 2026-08-16	9/450	2026-08-17 12:20 by lmz0216
[精细化工] 招聘金属平磨液，抛光液研发工程师 +3	小天0311 2026-08-14	3/150	2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7	Tide man 2026-08-14	8/400	2026-08-15 01:11 by arzu_hma