24小时热门版块排行榜    

查看: 3339  |  回复: 5

huang998

木虫 (小有名气)

[求助] HCP结构的晶格常数怎么取?a和c不相同 已有1人参与

比如镁是HCP结构,它的晶格常数应该如何取?a=b3=3.2094   c=5.2105
回复此楼

» 猜你喜欢

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

huang998

木虫 (小有名气)

2楼2017-06-13 14:34:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cl1994861

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
huang998(月只蓝代发): 金币+2, 感谢帮助! 2017-06-15 10:05:29
这个虽然我没用过hcp,但是lammps官网mannual里lattice command里应该会有解释。。。
3楼2017-06-13 16:10:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huang998

木虫 (小有名气)

引用回帖:
3楼: Originally posted by cl1994861 at 2017-06-13 16:10:03
这个虽然我没用过hcp,但是lammps官网mannual里lattice command里应该会有解释。。。

这个,我早看了,但没我想要的答案

发自小木虫IOS客户端
4楼2017-06-13 16:11:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cl1994861

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
月只蓝: 应助指数+1 2017-06-15 10:04:34
huang998(月只蓝代发): 金币+18, 感谢应助! 2017-06-15 10:05:15
引用回帖:
4楼: Originally posted by huang998 at 2017-06-13 16:11:43
这个,我早看了,但没我想要的答案
...

因为六角密堆结构的截面是六边形的,所以当我们在构建模拟盒子的时候,有可能因为觉得边界周期性不好处理,而觉得不知道如何构建。
而实际上,我们的考虑是多余的,因为LAMMPS为我们提供的HCP的原胞已经将这个问题很好的解决了。因此我们在使用时,只需要按着如下定义:
lattice hcp alength
其中的alength就是晶格常数a的长度。
下面是LAMMPS-lattice命令中对HCP的说明:
Style hcp has a1 = 1 0 0, a2 = 0 sqrt(3) 0, and a3 = 0 0 sqrt(8/3)…A hcp lattice has 4 basis atoms, two in the z = 0 plane and 2 in the z = 0.5 plane.
这里可以看到,LAMMPS在选取的是正交的坐标系,而不是通常用来定义HCP的四轴坐标系中的三轴。下图展示了x-y平面上的坐标轴。
description of HCP
红色的长方形盒子即为原胞在该平面的投影,该盒子在x-y平面包含2个原子,整个盒子包含4个原子。这样的盒子的周期性是很容易理解的。
另外,有些人可能会有疑问:HCP的晶格常数有两个:a 和 c,为什么这里提到了a ?
确实如此,LAMMPS已经将c/a给固定了(不然它没办法确定z轴)。对于c/a不为这个固定比值的情况,只需要对结构进行弛豫即可。
5楼2017-06-13 17:33:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huang998

木虫 (小有名气)

引用回帖:
5楼: Originally posted by cl1994861 at 2017-06-13 17:33:45
因为六角密堆结构的截面是六边形的,所以当我们在构建模拟盒子的时候,有可能因为觉得边界周期性不好处理,而觉得不知道如何构建。
而实际上,我们的考虑是多余的,因为LAMMPS为我们提供的HCP的原胞已经将这个问题 ...

是这样的

发自小木虫IOS客户端
6楼2017-06-13 19:07:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huang998 的主题更新
信息提示
请填处理意见