MATLAB在结构化学中的应用.pdfVIP

2008 NO.36 科技创新导报 Science and Technology Innovation Herald 工 业 技 术 MATLAB在结构化学中的应用 张俊林 张西沙 陈根强 (郑州大学化学系 河南郑州 450001) 摘 要:MATLAB以其精确的数值与符号运算能力,强大的作图与拟合功能,在工程技术领域应用广泛。但其强大计算优势运用于以量 子力学为基础的结构化学中尚很少。本文以休克尔分子轨道的求解,波函数的计算,原子轨道径向和角向分布图的绘制以及节面和最值 的求解等过程为例,将MATLAB应用于结构化学计算。 关键词:MATLAB结构化学 原子轨道 波函数作图计算 中图分类号:Q5 文献标识码:A 文章编号:1674-098X(2008)12(c)-0073-02 MATLAB是矩阵 (matrix)和实验室 (laboratory)两单词前三个字母的组合,最初 由美国New Mexico大学计算机系主任 Cleve Moler设计,经过国外多年的研发,已 发展成为国际市场数值软件的主流。其精 确的数值与符号运算能力,强大的作图与 拟合功能,使之在科研工程技术领域应用 广泛。在欧美,MATLAB已成为自动控制 理论、数理统计、数字信号处理等高级课 [1] 程的基本工具。 以量子力学为主要理论基 础的结构化学,则是通过大量复杂的数学 处理来深入原子和分子内部探究微观世界 的前沿科学。繁杂的数值计算,抽象的图 图1 氢原子3p轨道径向分布图 形绘制成为摆在结构化学科研工作者面前 应用matlab求解编程如下: 图,径向和角向分布图对形象理解波函数 的巨大难题。然而随着计算机技术的飞速 syms x 十分有益。径向分布物理意义在于表示在 发展,随之研发产生的MATLAB能够有效 D=[x 1 0 0;1 x 1 0;0 1 x 1;0 0 两球面间夹层内找到电子的概率。角向分 的解决这些科研难题。但其强大计算优势 1 x]; 布则是重在考虑轨道方向性。 运用于以量子力学为基础的结构化学中尚 d=det(D) 氢原子的 3p径向分布函数为:R=1/ 很少。本文以休克尔分子轨道的计算,波 d =x^4-3*x^2+1 函数的求解,原子轨道径向和角向分布图 solve(x^4-3*x^2+1=0)[4] a0.*exp(-r/a0); 的绘制以及节面和最值的求解等过程为 得四解:1/2*5^(1/2)-1/2 -1/2- 做径向分布函数图:a0=52.93 例,以一斑窥全貌,将以编程和计算机处理 1/2*5^(1/2) 1/2*5^(1/2)+1/2 1/2- r=0:0.001*a0:10*a0; 为