MATLAB在推导任意晶系晶面间夹角公式中的应用.docVIP

基于MATLAB的通用晶面间夹角公式之推导与求解 刘 健 （广东科技学院，东莞，523083） 摘 要 充分利用倒易点阵基矢与正点阵基矢互为倒易的特点，推导出通用晶面间夹角求解公式。同时利用MATLAB语言的核心-矩阵与数值计算编制了求解通用晶面间夹角的程序。并介绍了如何利用该程序求解任意类型晶体的各晶面间夹角,再对程序求解的晶体的晶面间夹角进行了验证,证明该程序在晶体学研究中有较好的应用价值。 关键词 晶面夹角，倒易点阵，MATLAB，向量 Derivation and solution of General Formula solving the angle betweentwo crystal planes of arbitrary crystal systems based on MATLAB Liu jian (Guangdong university of sciencetechnology, Dongguan , 52308) Abstract Make full use of the reciprocal lattice basis vectors and punctuality array vector mutually reciprocal characteristics, deduces the general solution formula of crystal plane angle, and, the use of MATLAB language core matrix and numerical calculation for solving general crystal plane angle program, and introduces how to use the program to solve the arbitrary type crystals crystal plane angle, then the procedure of solving crystal plane angle was verified, it is shown that the program is in crystallographic studies have good application value. Kewords crystal plane angle , reciprocal lattice, MATLAB , vector 晶面与晶面在空间的几何关系是材料科学和晶体学等研究领域中的重要参数，其中晶面间夹在许多实际应用，如晶体光学、晶体结构分析以及多晶体择优取向等方面都有广泛应用。对于晶面夹角问题，几乎所有教材均只给出正交、立方、六方三种简单晶系的最后求解公式，没有推导过程，对于其他晶系则几乎从未提及。因各大晶系晶格参数不同，其晶面间夹角公式也不尽相同，如果每遇到一种晶系都去用不同的公式求解的话，则非常麻烦而且难于记忆。如果有一个通用公式适合于任意晶系晶面间夹角的计算，显然会为晶体学研究领域带来前所未有的好处。对于晶体学问题的求解通常局限于晶体投影法和用晶面与晶向指数表示的解析法[1]。前者因涉及到球面投影、极射投影、乌氏网、极式网、标准投影等多方面晶体投影学知识，当遇到高指数晶面时，不但操作过程比较复杂，而且其结果准确性较差。绝大多数教材与资料都采取解析法，并且都引进倒易点阵的概念，但因其涉及到倒易矢量相乘等计算，实际求解很不方便而且过程特别复杂[2]，所以往往只给出最后结果，并无完整的推导过程。对此，文献[3]用C语言编写了任意晶系晶面间夹角的计算程序，但并没有给出具体公式，也没有推导过程，而且该程序比较复杂。有鉴于此，笔者充分利用倒易点阵基矢与正点阵基矢互为倒易的特点推导出该通用公式，并用MATLAB语言编写了相应的计算程序。 1． 任意晶系晶面间夹角通用公式的推导 晶面间夹角指的是两晶面矢间的交角，在晶体学中通常用两晶面法向（一般用晶向代表）间的交角来表示。设任意两晶面（h1 k1 l1） (h2 k2 l2) ，正点阵基矢为 a 、b、 c ，倒易点阵基矢为 a﹡、 b﹡、 c﹡,两晶面对应的倒易点阵为g1﹡ 、g2﹡ 。 [1] 建立晶轴矢量 设有任一晶体坐标系，其晶格参数为a、b、c、α、β、γ，三晶轴与直角坐标系按图1所示安置。令三轴矢量依次为a=(Xa Ya Za) b=(Xb Yb Zb) c=(Xc Yc Zc)　显然有=(a 0 0) b=(bcosγ bsinγ 0) Xc=ccosβ 对于Yc由=bccosα以及b