24小时热门版块排行榜

返回列表

本帖产生 1 个 1ST强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
帖子: 817
在线: 113.5小时
虫号: 968169

[交流] 【求助成功】vc-relax优化后的

昨天试了试vc-relax优化，但是优化后输出的结果只有cell-parameter，和晶胞体积，和原子坐标输出…………请问从中如何看出晶格常数是如何变化的…………

[ Last edited by 改王之王 on 2011-4-3 at 10:04 ]

回复此楼

» 猜你喜欢

表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有12人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有14人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有11人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复

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

晶体模型几何优化后，所得结果反而不如未优化的模型已经有19人回复
CASTEP对优化后的晶体进行分态密度分析，应该怎么操作已经有14人回复
优化后分子构型完全转变，求教如何解决已经有7人回复
castep优化后得到的xsd文件是倒空间，怎么显示正空间？已经有7人回复
Gaussian几何优化后的键长不合理已经有41人回复
【求助】vc-relax得到的晶格常数偏小已经有7人回复
【求助】结构优化后，stress仍然很大，如何解决？已经有7人回复
【求助】只是优化后文件，没有restart文件，能继续跑吗已经有22人回复
【讨论】H-BEA分子筛优化后结果和文献不一样，是怎么回i事？高手指点一下已经有4人回复
【求助】Pwscf计算下的fractional coordinate 已经有5人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-04-02 08:58:03

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

ChemiAndy

木虫 (正式写手)

1ST强帖: 4
应助: 52 (初中生)
金币: 2573.3
帖子: 590
在线: 553.8小时
虫号: 836907

★
小木虫(金币+0.5):给个红包，谢谢回帖

引用回帖:

3楼: Originally posted by souledge at 2011-04-01 14:05:40:
cell-parameter不就是晶格常数么？比如cell-parameter的单位为alat=x，则x就是晶格常数的基数，其单位为Bohr。
然后，比如cell-parameter的第一行，意义为a方向的矢量，从左到右三个数分别代表：
1->在迪卡尔 ...

请问您这个转换方法的来源是哪里？因为我在用另外一个不同的转换方法，即我认为CELL PARAMETERS下面的三行并不是cell parameters, 而是cell vectors, 比如：
CELL_PARAMETERS (alat= 19.80400009)
0.997610862 0.000000000  -0.003022554
0.000000000 0.596367719 0.000000000
  -0.257837471 0.000000000 1.127252205
假如上面的矩阵元素依次为
x1, y1, z1
x2, y2, z2
x3, y3, z3

那么，先做一个转化，即上面矩阵的所有元素都乘以alat，得到绝对值，然后，

A = x1^2 + y1^2 + z1^2
B = x2^2 + y2^2 + z2^2
C = x3^2 + y3^2 + z3^2

cos(alpha) = (B^2 +C^2 - |BC|^2) / 2*B*C
cos(beta)  = (A^2 +C^2 - |AC|^2) / 2*A*C
cos(gama)  = (A^2 +B^2 - |AB|^2) / 2*A*B

其中，
|AB|^2 = (x2-x1)^2 + (y2-y1)^2 + (z2-z1)^2
|BC|^2 = (x3-x2)^2 + (y3-y2)^2 + (z3-z2)^2
|AC|^2 = (x3-x1)^2 + (y3-y1)^2 + (z3-z1)^2

这个显然跟你的算法不一样，不知道拿一个是正确的。我还没有检验，应该很简单，在excell里面很容易建立这样的转化。附件是我用Linux下的openoffice建立的转化表格。希望论坛里的朋友验证一下。

赞一下(1人)

回复此楼

16楼2011-09-14 05:04:58

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

查看全部 20 个回答

souledge

专家顾问 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+3): 谢谢指点 2011-04-02 10:07:45
zzy870720z(1ST强帖+1): 很详细，呵呵 2011-04-02 14:15:51
改王之王(金币+3): 感谢…………这下理解了……呵呵 2011-04-03 10:03:35

cell-parameter不就是晶格常数么？比如cell-parameter的单位为alat=x，则x就是晶格常数的基数，其单位为Bohr。
然后，比如cell-parameter的第一行，意义为a方向的矢量，从左到右三个数分别代表：
1->在迪卡尔坐标系中x轴上投影长度，单位是alat。
2->与迪卡尔坐标中y轴夹角的余弦。
3->与迪卡尔坐标中z轴夹角的余弦。
这里一定要分清楚作为参考系的迪卡尔坐标和晶向a,b,c，是不一定重合的，所以才有了夹角和投影长度。所以，通过该cell-parameter可以回推晶格的变化。
其实CASTEP里也用了同样的方法表示晶格，只不过CASTEP最后还把这个表示方法做成了xsd模型图，看起来一目了然了而已~

赞一下(3人)

回复此楼

3楼2011-04-02 10:05:40

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