径向分布函数..docVIP

实验一 径向分布函数、角度分布函数电子云图形的绘制 一、实验目的 1．绘制波函数及其各种分布以及电子云的图像，观察各种函数的分 2．了解计算机绘图方法。 二、实验 1．程序原理：本程序可绘制类氢原子的径向分布函数，角度分布函 数及原子轨道、 杂化轨道和分子轨道等电子几率密度图， 绘制过程中 的各函数形式列于下列各表中。式中 2 Zr ，n a为0主量子数， na 0 =0.0529nm，为波尔半径， Z 是有效核电荷，由 Slater规则计算得到 的周期表中前四个周期元素的有效核电荷列于表 1.1 中，下面简要叙 述对各类图形的处理方案。 ①径向分布函数图： 径向分布函数 D(r)=r 2R2(r) 反映了电子的几率随半径 r 的分布情况， D(r)dr 代表半径 r 到 r+dr 两个球壳夹层内找到电子的几率。其中 R(r)为类氢原子的径向函数， 本程序所采用的径向函数 R(r)分别列于表 2－2 中。 ②角度分布函数图： nlm (r , , )的角度部分 lm ( , ) 以及角度分布函数 2 lm ( , ) 表示同 一球面不同方向上 所采用的角度函数  nlm (r , , ) 或 nlm ( r , , ) 的相对大小，本程序 2 lm ( , ) 分别列于表 3－3 pz , pz2 , f z3 , f xz2 , ( f yz2 ),Ysp ,Yd 2sp3 角度分布图是画的 X-Z 平面的截面图， 其余角度分布图都是画的 X-Y 平面的截面图。角度分布函数图中， ③等电子几率密度图： (r , , ) 2 称为电子几率密度函数， 它描述在该 轨道中的电子在三维空间的分布情况， 为了在平面上表示出这种分布 往往采用某一切面上的等值面图， 程序按指定的轨道在该切面上逐点 计 算 2 的值 ， 及 找 出 max 2的最大值，求出相对几率密度 P 2 / max 2 ，该值在 X-Y 平面上是位置坐标 (x,y)的函数 (对于 3d z2 轨 道是在 X-Z 平面 )，绘图时不是将取值相同的点连成曲线，而是打印 一系列符号表示相对几率密度的分布区域。当 P＜0.01时为空白， 0.01≤P＜0.02时用“：”，0.02≤P＜0.1 时用“／”，0.1≤P＜0.25 时用 “O”，0.25≤P＜0.5 时用“ ”和 P＞0.5 时用“ # ”符号表示。根据这些符号可以粗略看出几 在 X-Y 平面内，坐标变化范围为 -2.4≤x≤2.4(步长 =0.08) -1.42≤y≤1.42(步长 =0.133) 所有距离的长度单位都是 10-10m 原子轨道使用的波函数如表 1－4 所示。对 3d z2 ,4d z2 ,4 f z3 , 和 4 f xz2 轨道采 用 X-Z 平面做截面，所有其它原子轨道都画在 X-Y 平面上，程序使用原子轨道的四重轴对称性， 首先计算第三象限内， 即-2.4≤x≤ 0,-1.42 ≤y≤0 的Ψ值，随后被 P 2 / max 2 代替，在其它三个象限内的相应 值由对称性得到，用 P(x,y)代表电子在坐标 (x，y)点的几率密度，则： P(-x,-y)=P(-x,y)=P(x,-y)=P(x,y) 表 1－1 Slater轨道中的 Z * 参量值 Z 原 有效核电荷 Z * Z 原 有效核电荷 Z * ls 2s,2p 3s,3p 1s 2s,2p 3s,3p 3d 4s,4p 子 子 1 H 1.00 19 K 18.70 14.85 7.75 2.20 2 He 1.70 20 Ca 19.70 15.85 8.75 2.85 3 Li 2.70 1.30 21 Sc 20.70 16.85 9.75 3.00 3.00 4 Be 3.70 1.95 22 Ti 21.70 17.85 10.75 3.65 3.15 5 B 4.70 2.60 23 V 22.70 18.85 11.75 4.30 3.30 6 C 5.70 3.25 24 Cr 23.70 19.85 12.75 5.60 2.95 7 N 6.70 3.90 25 Mn 24.70 20.85 13.75 5.60 3.60 8 O 7.70 4.55 26 Fe 25.70 21.85 14.75 6.25 3.75 9 F 8.70 5.20 27 Co 26.70 22.85 15.75 6.90 3.90 10 Ne 9.70 5.85 28 Ni 27.70 23.85 16.75 7.55 4.05 11 Na 10.70 6.85 2.20 29 Cu 28.70 24.85 17.75 8.85 3.70 12 Mg 11.70 7.85