24小时热门版块排行榜    

查看: 4581  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +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能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +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
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
信息提示
请填处理意见