| 查看: 987 | 回复: 10 | |||
| 当前主题已经存档。 | |||
dreampursuer金虫 (正式写手)
|
[交流]
【求助】有没有办法求倒格空间的晶体结构？
|
||
|
有没有办法求倒格空间的晶体结构？ 正空间的原胞用material studio画出来了，想求倒格空间的原胞，也就是WS原胞。想问material studio可以做得到吗？谢谢 [ Last edited by dreampursuer on 2009-12-10 at 15:28 ] |
回复此楼» 猜你喜欢
286求调剂
已经有14人回复
-
一志愿211，0703化学305分求调剂
已经有10人回复
-
362求调剂
已经有8人回复
-
0703调剂
已经有17人回复
-
333求调剂
已经有4人回复
-
085600材料与化工专硕329 求调剂
已经有12人回复
-
一志愿吉大化学327求调剂
已经有8人回复
-
0703化学调剂 348分
已经有8人回复
-
一志愿太原理工大学计算机技术专硕348，求调剂指导
已经有3人回复
-
材料考研求调剂总分280
已经有16人回复
huijiany
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1441.6
- 散金: 2507
- 红花: 5
- 帖子: 695
- 在线: 244.6小时
- 虫号: 426431
- 注册: 2007-07-28
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2009-12-10 10:42:37
dreampursuer
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1212.2
- 散金: 199
- 红花: 3
- 帖子: 693
- 在线: 110.7小时
- 虫号: 517668
- 注册: 2008-03-04
- 专业: 核技术及其应用
3楼2009-12-10 10:59:41
aylayl08
荣誉版主 (文坛精英)
至尊天蝎
- 1ST强帖: 1
- 应助: 3 (幼儿园)
- 贵宾: 4.825
- 金币: 19183.7
- 散金: 2180
- 红花: 34
- 沙发: 48
- 帖子: 16357
- 在线: 824.1小时
- 虫号: 533169
- 注册: 2008-03-26
- 专业: 无机非金属基复合材料
- 管辖: 量子化学
4楼2009-12-10 14:02:36
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★ ★ ★ ★ ★ ★ ★
dreampursuer(金币+1):谢谢参与
dreampursuer(金币+3,VIP+0):你都没有说明你这个程序段属于什么语言的呢。。 12-10 14:27
aylayl08(金币+5,VIP+0):very good.thank you 12-10 15:41
dreampursuer(金币+1):谢谢参与
dreampursuer(金币+3,VIP+0):你都没有说明你这个程序段属于什么语言的呢。。 12-10 14:27
aylayl08(金币+5,VIP+0):very good.thank you 12-10 15:41
|
很简单。 我有写了一个小程序。 你只需要编辑一个输入文件，把正空间的晶格矩阵给它就行。 input.ini 的格式如下： 3.5 0.0 0.0 # Ax , Ay, Az 0.0 3.5 0.0 # Bx, By , Bz 0.0 0.0 3.5 # Cx, Cy , Cz 编译我的程序，执行就可以得到倒空间的晶格矩阵了。 有了倒空间的矩阵，你就可以在ms之类的软件里，看倒空间的样子了。 !--------------------------------------------------------------------------------------------- ! program : from dir to reci ! ! subroutine: main ! ! author: H.Y Liu ! ! contact: ouuing@gmail.com ! !--------------------------------------------------------------------------------------------- program main implicit none real(kind=8) d(3,3) integer i logical alive ! read d_lattice write(*,*)"first, you need write direct lattice in input.ini******important****" write(*,*)"--------------------------------------------------------------------" write(*,*)"this is direct lattice:" write(*,*) inquire(file="input.ini",exist=alive) if(alive)then open(101,file="input.ini"  else write(*,*)"no input.ini" write(*,*)"I can give your example of input.ini" write(*,*)"please edit it!" write(*,*)"format of input.ini :" write(*,*)"3.5 0.0 0.0 # Ax Ay Az" write(*,*)"0.0 3.5 0.0 # Bx By Bz" write(*,*)"0.0 0.0 3.5 # Cx Cy Cz" stop end if do i=1,3 read(101,*)d(1,i),d(2,i),d(3,i) write(*,80)d(1,i),d(2,i),d(3,i) end do 80 format(3x,3f15.9) call d_reci(d) ! write reci lattice stop end !------------------------------------------------------------------------------------- !subroutine: reciprocal lattice !------------------------------------------------------------------------------------- subroutine d_reci(d) implicit none real(kind=8) d(3,3),r(3,3),volume,recv integer i,j ! volume=(a*b.c) volume=(d(2,1)*d(3,2)-d(3,1)*d(2,2))*d(1,3)+(d(3,1)*d(1,3)& &-d(1,1)*d(2,2))*d(2,3)+(d(1,1)*d(2,2)-d(2,1)*d(1,2))*d(3,3) write(*,*)"-------------------------------------------------------------" write(*,*)"volume of the cell is :" write(*,84)volume write(*,*) 84 format(3x,f15.9) ! reci from direct and write reci r(1,1)=(d(2,2)*d(3,3)-d(3,2)*d(2,3))/volume r(2,1)=(d(3,2)*d(1,3)-d(1,2)*d(3,3))/volume r(3,1)=(d(1,2)*d(3,3)-d(3,2)*d(1,3))/volume r(1,2)=(d(2,3)*d(3,1)-d(3,3)*d(2,1))/volume r(2,2)=(d(3,3)*d(1,1)-d(1,3)*d(3,1))/volume r(3,2)=(d(1,3)*d(2,1)-d(2,3)*d(1,1))/volume r(1,3)=(d(2,1)*d(3,2)-d(3,1)*d(2,2))/volume r(2,3)=(d(3,1)*d(1,2)-d(1,1)*d(3,2))/volume r(3,3)=(d(1,1)*d(2,2)-d(2,1)*d(1,2))/volume write(*,*)"-------------------------------------------------------------" write(*,*)"this is reciprocal lattice:" write(*,*) do i=1,3 write(*,81)(r(j,i),j=1,3) end do 81 format(3x,3f15.9) ! volume=(a*b.c) recv=(r(2,1)*r(3,2)-r(3,1)*r(2,2))*r(1,3)+(r(3,1)*r(1,3)-& &r(1,1)*r(2,2))*r(2,3)+(r(1,1)*r(2,2)-r(2,1)*r(1,2))*r(3,3) write(*,*)"-------------------------------------------------------------" write(*,*)"volume of the reciprocal-cell is :" write(*,82)recv write(*,*) 82 format(3x,f15.9) return end subroutine |
5楼2009-12-10 14:23:42
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
6楼2009-12-10 14:36:30
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
7楼2009-12-10 14:37:17
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
8楼2009-12-10 15:07:03
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
9楼2009-12-10 16:52:00
netx_ray
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2358.1
- 红花: 2
- 帖子: 146
- 在线: 140小时
- 虫号: 133380
- 注册: 2005-12-14
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
10楼2009-12-10 16:57:26
