24小时热门版块排行榜    

查看: 4567  |  回复: 9

cnmgtb

银虫 (小有名气)

[求助] POSCAR中的晶格矢量的两种方式 已有5人参与

在做体系的优化过程中,同样的原子和坐标,为什么导出来的POSCAR中的晶格矢量是不一样的。
     Ni-mono-fcc                             
     1.00000000000000     
     1.7325000022546899   -1.7325000022546899    0.0000000000000000
     0.0000000000000000    1.7325000022546899   -1.7325000022546899
    17.3205080756887995   17.3205080756887995   17.3205080756887995

      Ni-111-fcc
      1.000
       2.121866162320330      -1.225060000000000       0.000000000000000
       0.000000000000000       2.450120000000000       0.000000000000000
       0.000000000000000       0.000000000000000      30.000000000000000
求大神解释啊·····
回复此楼

» 猜你喜欢

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

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

wmy8802217

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
估计做了一个坐标变换,把坐标轴转了下
凝聚态理论与计算物理
2楼2015-09-23 15:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jiansh

新虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
不明白,一般不会出现这种情况
请确认是否是一个结构

或者你把POSCAR导出。然后用vesta看看结构是否一致。
学术网(www.scholarnet.cn)秉承一切为学术工作者的理念,竭诚为广大科研工作者提供优质的服务。
3楼2015-09-23 16:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cnmgtb

银虫 (小有名气)

引用回帖:
2楼: Originally posted by wmy8802217 at 2015-09-23 15:53:57
估计做了一个坐标变换,把坐标轴转了下

可是坐标轴转了 但是相对的位置和晶格参数是不会变得呀···为什么如今的C方向变化那么大
4楼2015-09-24 11:09:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cnmgtb

银虫 (小有名气)

引用回帖:
3楼: Originally posted by Jiansh at 2015-09-23 16:22:37
不明白,一般不会出现这种情况
请确认是否是一个结构

或者你把POSCAR导出。然后用vesta看看结构是否一致。

确定是一个结构 只是之前的结构是我的导师导出来的POSCAR 然后我把它导入了MS 重新导出的POSCAR 就出现这样的差别····求大神解释
5楼2015-09-24 11:10:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by cnmgtb at 2015-09-24 11:09:42
可是坐标轴转了 但是相对的位置和晶格参数是不会变得呀···为什么如今的C方向变化那么大...

你把C方向的17.3*根号三不还是三十。没问题的,是一个体系
凝聚态理论与计算物理
6楼2015-09-24 11:34:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

get-it

新虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
形状没变,abc的长度没变,两两点乘的值也没变,只是坐标变换成了一个上三角矩阵。
请问你是怎么操作的,我现在想把晶格变成三角矩阵
7楼2015-09-24 16:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
?@??M??????r????????????L???????g??A????]????????M?????????????Q??????M??

?O????M??????????A1=[     1.7325000022546899   -1.7325000022546899    0.0000000000000000
     0.0000000000000000    1.7325000022546899   -1.7325000022546899
    17.3205080756887995   17.3205080756887995   17.3205080756887995]

????M??????????A2=[       2.121866162320330      -1.225060000000000       0.000000000000000
       0.000000000000000       2.450120000000000       0.000000000000000
       0.000000000000000       0.000000000000000      30.000000000000000]

???N????l?FA1?????????A1'*A1??A2?????????A2'*A2???????????????
    6.0031   -3.0015         0
   -3.0015    6.0031         0
         0         0  900.0000

????M????????Q??Q?????????????????????????????????????????????????L????????????????t??????????????????????e??????????@???????????????????c??a,b?????????????????????????a??b?????????????????????????????A???acos( -3.0015/sqrt(6.0031* 6.0031) ) = 120??

???N?????A1??Q??A2??????A2????????????O
A2=[ x y 0
        0 z 0
        0 0 w]
???NA2??????????
[  x^2+y^2  yz      0
   yz             z^2   0
   0              0       w^2
]
????????A1??????????????????????????????????
w=sqrt(900.0000) = 30
z=sqrt(6.0031) = 2.4501
???y=-3.0016/z = -1.2251
?????x=sqrt(6.0031-y^2) = 2.1218
?????
A2 =[ 2.1218  -1.2251 0
        -1.2251 2.4501   0
         0          0           30]
Simpler, stronger, more tolerant
8楼2015-09-24 21:14:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

引用回帖:
8楼: Originally posted by neweroica at 2015-09-24 21:14:44
這兩組矢量是等價的,即它們的長度和彼此之間的夾角都沒變。?牡谝唤M矢量可以如下變換到第二組:

設第一組矢量對應矩陣A1=

第二組矢量對應矩陣A2=

那麼你會發現A1的?埩烤仃嘇1'*A1跟A2的?埩烤仃嘇2'*A2是完 ...

糾正筆誤:A1的?埩烤仃囀茿1*A1',不是A1'*A1。A2同理
Simpler, stronger, more tolerant
9楼2015-09-24 21:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lpf525

新虫 (初入文坛)

【答案】应助回帖

我也遇到这种情况,请问第一种的晶格矢量是怎么导出来的
10楼2017-10-20 09:39:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cnmgtb 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +3 Tide man 2026-08-14 4/200 2026-08-14 16:33 by Yeuchan
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +4 Ldrop2023 2026-08-13 4/200 2026-08-14 14:50 by Tide man
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] filecode +5 cratir 2026-08-14 9/450 2026-08-14 12:25 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青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 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
请填处理意见