| 查看: 4383 | 回复: 19 | |||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | |||
[交流]
【求助成功】vc-relax优化后的
|
|||
|
昨天试了试vc-relax优化，但是优化后输出的结果只有cell-parameter，和晶胞体积，和原子坐标输出…………请问从中如何看出晶格常数是如何变化的………… [ Last edited by 改王之王 on 2011-4-3 at 10:04 ] |
回复此楼» 猜你喜欢
博士自荐
已经有4人回复
-
317一志愿华南理工电气工程求调剂
已经有9人回复
-
材料化工调剂
已经有3人回复
-
0856化工专硕求调剂
已经有4人回复
-
2025届双非化工硕士毕业，申博
已经有4人回复
-
0856求调剂285
已经有6人回复
-
博士自荐
已经有10人回复
-
高分子化学与物理调剂
已经有8人回复
-
290求调剂
已经有7人回复
-
272求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 晶体模型几何优化后，所得结果反而不如未优化的模型
已经有19人回复
- CASTEP对优化后的晶体进行分态密度分析，应该怎么操作
已经有14人回复
- 优化后分子构型完全转变，求教如何解决
已经有7人回复
- castep优化后得到的xsd文件是倒空间，怎么显示正空间？
已经有7人回复
- Gaussian几何优化后的键长不合理
已经有41人回复
- 【求助】vc-relax得到的晶格常数偏小
已经有7人回复
- 【求助】结构优化后，stress仍然很大，如何解决？
已经有7人回复
- 【求助】只是优化后文件，没有restart文件，能继续跑吗
已经有22人回复
- 【讨论】H-BEA分子筛优化后结果和文献不一样，是怎么回i事？高手指点一下
已经有4人回复
- 【求助】Pwscf计算下的fractional coordinate
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
青岛科技大学JIEQING团队招收材料科学、化学化工相关领域的硕/博研究生数名
+1/90
-
2026年药物化学博士招生
+1/84
-
求助哪里可以做光催化产氢的放大实验
+1/83
-
东北石油大学三亚海洋油气研究院【硕士招生】【地质资源与地质工程】【地质学】
+1/79
-
西安交通大学杜宝吉课题组招收2026统考硕士生
+1/30
-
澳门理工大学人工智能智慧康养2026 年9月入学 博士招生有奖学金
+1/27
-
电子科技大学材料学院SFT创新中心招收博士生及科研助理 长期有效
+1/12
-
加拿大卡尔顿大学【招博后】【心理学人工智能】
+1/12
-
人工智能300分 求调剂，有项目竞赛经历，能写代码
+1/11
-
“教育部长江学者和创新团队发展计划”入选团队招收材料，化学与化工博士研究生
+1/10
-
北京某研究院结构生物学相关专业申请考核制博士招生（春季）
+1/8
-
电子科技大学SFT创新中心（材料、生物医学工程及相关研究方向）诚聘教师和博士后
+1/6
-
蛋白结构质谱实验室202609入学博士生（1名，202603月底报名）
+1/3
-
澳科大诚招2026年秋季药物递送/生物材料方向硕士研究生（倒计时最后一周）
+1/3
-
中国石油大学（华东）生物质绿色化学转化和生物基功能高分子团队诚招26级硕士生
+1/3
-
深圳大学何传新课题组诚聘“百人计划”（副）教授、助理教授
+1/3
-
环境纳米技术@水质风险控制 课题组接收研究生报考/调剂公告
+1/3
-
燕山大学-意大利米兰理工大学国际化科研团队招收硕士研究生
+1/3
-
武汉科技大学黄思煜教授招收博/硕士及师资博士后、特任教授
+1/2
-
PostDoc fellowship: Video-Based 3D Human Kinematics Extraction/Pose Esitimation!
+1/1
4楼2011-04-02 10:07:27
2楼2011-04-02 09:36:39
★ ★ ★ ★
小木虫(金币+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
小木虫(金币+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楼2011-04-02 10:05:40
5楼2011-04-02 16:44:22
6楼2011-04-02 20:56:34
★
zzy870720z(金币+1): 谢谢总结 2011-04-03 10:34:58
zzy870720z(金币+1): 谢谢总结 2011-04-03 10:34:58
|
现在明白了…………大感谢几位的热心解答………… 现在明白了，特别是对于直角坐标更好解释 例如 CELL_PARAMETERS (alat= 6.22300000) 0.997142196 0.000000000 0.000000000 -0.498571098 0.863550473 -0.000000000 0.000000000 -0.000000000 1.611099270 第一排只有第一列有值，就说明晶格常数a只在x轴上有投影，即acosθ=0.997，这时θ=0，所以a=0.977,然后在乘以6.223得6.205bohr就是a的晶格常数了………… 我用的xcrysden是1.4.1好像没这个功能，我下个新的试试…………呵呵 |
7楼2011-04-03 10:18:04
8楼2011-04-03 12:42:38
9楼2011-04-03 21:07:05
10楼2011-04-03 23:38:43
11楼2011-04-03 23:39:44
不知道直接下载原码版，能不能在Windows下编译啊~没用过Windows下的*nix环境，不敢乱说~ 12楼2011-04-03 23:50:45
13楼2011-04-04 08:43:29
14楼2011-04-09 14:59:04
15楼2011-05-10 08:03:26
★
小木虫(金币+0.5):给个红包，谢谢回帖
小木虫(金币+0.5):给个红包，谢谢回帖
|
请问您这个转换方法的来源是哪里？因为我在用另外一个不同的转换方法，即我认为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建立的转化表格。希望论坛里的朋友验证一下。 |
16楼2011-09-14 05:04:58
|
请问您这个转换方法的来源是哪里？因为我在用另外一个不同的转换方法，即我认为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 求夹角alpha要做一个弧度变换，即 alpha = arccos[cos(alpha)*] 这个显然跟你的算法不一样，不知道拿一个是正确的。我还没有检验，应该很简单，在excell里面很容易建立这样的转化。附件是我分别用Microsoft Excel和Linux下的openoffice建立的转化表格。希望论坛里的朋友验证一下。[ Last edited by ChemiAndy on 2011-9-13 at 09:45 ] |
17楼2011-09-14 05:05:30
18楼2011-09-24 03:07:16
19楼2011-09-24 08:13:36
20楼2019-01-02 22:12:12
