24小时热门版块排行榜    

查看: 4195  |  回复: 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 的主题更新
信息提示
请填处理意见