欢迎访问喜蛋文章网
你的位置：首页 > 经典文章 > 文章正文

电子云的形状原来是这么来的？!

时间: 2021-12-24 11:30:38 | 作者：zdr0 | 来源: 喜蛋文章网 | 编辑: admin | 阅读: 112次

电子云的形状原来是这么来的？!

电子云结构是这个世界上最精妙绝伦的结构之一。
总注：只想看科普的读者朋友可以忽视以下所有标注”（可忽视）”的引用块中的内容。

与电子云相关的一个重要的数学知识是---球谐函数。

首先，我们必须要知道，球谐函数从何而来？球谐函数是球坐标系下的 Laplace 方程的解。

Laplace 方程是一个非常重要的偏微分方程，其表达式为：

bigtriangleup f=0

其中： bigtriangleup 为 Laplace 算符，是一个标量算符。 f 是一个标量函数。 bigtriangleup 在不同的坐标系下具有不同的形式，这里我们研究的是其在球坐标系下作用于一个标量函数的解。

（可忽视）这里我只给出 bigtriangleup 在直角坐标系和球坐标系下的形式：直角坐标系下的 Laplace 算符： bigtriangleup =frac{partial^2 }{partial x^2}+frac{partial^2 }{partial y^2}+frac{partial^2 }{partial z^2} 直角坐标系下的 Laplace 方程： bigtriangleup f =frac{partial^2f }{partial x^2}+frac{partial^2 f}{partial y^2}+frac{partial^2f}{partial z^2}=0 球坐标系下的 Laplace 算符： bigtriangleup ={frac {1}{r ^{2}}}{frac {partial }{partial r }}left(r ^{2}{frac {partial }{partial r }}right)+{frac {1}{r ^{2}sin theta }}{frac {partial }{partial theta }}left(sin theta {frac {partial }{partial theta }}right)+{frac {1}{r ^{2}sin ^{2}theta }}{frac {partial ^{2}}{partial phi ^{2}}} 球坐标系下的 Laplace 方程： bigtriangleup f ={frac {1}{r ^{2}}}{frac {partial }{partial r }}left(r ^{2}{frac {partial f }{partial r }}right)+{frac {1}{r ^{2}sin theta }}{frac {partial }{partial theta }}left(sin theta {frac {partial f}{partial theta }}right)+{frac {1}{r ^{2}sin ^{2}theta }}{frac {partial ^{2}f}{partial phi ^{2}}}=0left( bigstarright) left( bigstarright) 就是我们要研究的方程。

我们假设标量函数 f 是一个三元函数，即： f=fleft( x,y,z right) 。在坐标系发生转换时 f 中变量的个数并不会改变。这就意味着， f 从直角坐标系下的表达形式转换到球坐标系下的表达形式之后变量的个数还是三个，即 f=fleft( x,y,z right)rightarrow fleft( r,phi,theta right) ，形式变换之后的三个变量分别是： r: 球坐标系径向变量（球的半径）， phi,theta: 球坐标系下的角变量，其中： phi 表示经度， theta 表示维度。

直角坐标系下的某点与对应的球坐标下的点。图片来源：维基百科。
（可忽视）直角坐标系和球坐标系之间的转换关系式为： {displaystyle {r}={sqrt {x^{2}+y^{2}+z^{2}}}} {displaystyle {theta }=arccos left({frac {z}{r}}right)=arcsin left({frac {sqrt {x^{2}+y^{2}}}{r}}right)=arctan left({frac {sqrt {x^{2}+y^{2}}}{z}}right)} {displaystyle {phi }=arccos left({frac {x}{rsin theta }}right)=arcsin left({frac {y}{rsin theta }}right)=arctan left({frac {y}{x}}right)}

现在我们假设 f 在球坐标系下的在三个变量 r,phi,theta 上的分量相互独立，这便可以将 fleft( r,phi,theta right) 写成三个一元函数乘积的形式：

fleft ( r,phi,theta right ):=Rleft ( r right )cdot Phileft ( phi right )cdot Thetaleft ( theta right )

其中： Rleft ( r right ) 是 f 的径向分量， Phileft ( phi right ) 是 f 的经度方向分量， Thetaleft ( theta right ) 是 f 的纬度方向分量。这个过程有点像因式分解对不对？我们也可以将后面的经度和维度分量合写写在一起，称为角分量（即仅与角有关的分量），并定义为：

Yleft( phi,theta right):=Phileft ( phi right )cdot Thetaleft ( theta right )

后面我们就会看到，这个二元函数 Yleft( phi,theta right) 就是所谓的球谐函数。定义了 Yleft( phi,theta right) 之后 f 可以写为：

fleft ( r,phi,theta right ):=Rleft ( r right )cdot Yleft( phi,theta right)

下面我们进入正题：通过球坐标下的 Laplace 方程求解球谐函数 Yleft( phi,theta right) 。（先进行数学推导，最后我们来分析一下得到的解）

（可忽视）我们利用分离变量法求解球坐标下的 Laplace 方程。首先，我们先做一些准备工作： f=Rcdot YRightarrow left{begin{matrix} f_{r}=R_{r}cdot Y f_{phi }=Rcdot Y_{phi } f_{theta }=Rcdot Y_{theta } end{matrix}right.Rightarrow left{begin{matrix} f_{rr}=R_{rr}cdot Y f_{phiphi }=Rcdot Y_{phiphi } f_{thetatheta }=Rcdot Y_{theta theta} end{matrix}right. Rightarrowleft (r^{2}cdot f_{r} right )_{r} =2rcdot f_{r}+r^{2}cdot f_{rr}=2rcdot R_{r}cdot Y+r^2cdot R_{rr}cdot Y left (mathrm{sin}theta cdot f_{theta } right ) _{theta}=mathrm{cos}theta cdot f_{theta } +mathrm{sin}theta cdot f_{thetatheta } =mathrm{cos}theta cdot Rcdot Y_{theta }+mathrm{sin}theta cdot Rcdot Y_{theta theta } 将上面的结果全部代入球坐标下的 Laplace 方程可得： bigtriangleup f ={frac {1}{r ^{2}}}left(2rcdot R_{r}cdot Y+r^2cdot R_{rr}cdot Yright)+{frac {1}{r ^{2}sin theta }}left(mathrm{cos}theta cdot Rcdot Y_{theta }+mathrm{sin}theta cdot Rcdot Y_{theta theta }right)+{frac {1}{r ^{2}sin ^{2}theta }}cdot Rcdot Y_{phiphi }=0 Rightarrowleft(2rcdot R_{r}+r^2cdot R_{rr}right)cdot Y=-left({frac {mathrm{cos}theta}{sin theta }} cdot Y_{theta }+Y_{theta theta }+{frac {1}{sin ^{2}theta }}cdot Y_{phiphi }right)cdot R Rightarrowfrac{1}{R}cdotleft(2rcdot R_{r}+r^2cdot R_{rr}right) =-frac{1}{Y}cdotleft({frac {mathrm{cos}theta}{sin theta }} cdot Y_{theta }+Y_{theta theta }+{frac {1}{sin ^{2}theta }}cdot Y_{phiphi }right):=lambda 从而得到一个二阶常微分方程和一个偏微分方程，分别为： 2rcdot R_{r}+r^2cdot R_{rr}-lambdacdot R=0Leftrightarrow2rcdot frac{mathrm{d} R}{mathrm{d} r}+r^{2}cdot frac{mathrm{d^{2}} R}{mathrm{d} r^{2}}-lambda cdot R=0 {frac {mathrm{cos}theta}{sin theta }} cdot Y_{theta }+Y_{theta theta }+{frac {1}{sin ^{2}theta }}cdot Y_{phiphi }+lambdacdot Y=0left(bigstar right) 我们再将 Yleft( phi,theta right):=Phileft ( phi right )cdot Thetaleft ( theta right ) 代入left(bigstar right) 有： {frac {mathrm{cos}theta}{sin theta }} cdotPhi cdotTheta_{theta }+Phi cdotTheta_{thetatheta }+{frac {1}{sin ^{2}theta }}cdot Phi _{phiphi }cdotTheta+lambdacdot Phi cdotTheta=0 Rightarrowmathrm{cos}thetacdot sin thetacdotPhi cdotTheta_{theta }+Phi cdotTheta_{thetatheta }cdotsin ^{2}theta+Phi _{phiphi }cdotTheta+sin ^{2}thetacdot lambdacdot Phi cdotTheta=0 Rightarrowleft(mathrm{cos}thetacdot sin thetacdotTheta_{theta } +Theta_{thetatheta }cdotsin ^{2}theta +sin ^{2}thetacdot lambdacdot Thetaright)cdotPhi =-Phi _{phiphi }cdotTheta Rightarrowleft(mathrm{cos}thetacdot sin thetacdotTheta_{theta } +Theta_{thetatheta }cdotsin ^{2}theta +sin ^{2}thetacdot lambdacdot Thetaright)cdotfrac{1}{Theta} =-Phi _{phiphi }cdotfrac{1}{Phi}:=m^2 从而再次得到两个二阶常微分方程：mathrm{cos}thetacdot sin thetacdotTheta_{theta } +Theta_{thetatheta }cdotsin ^{2}theta +sin ^{2}thetacdot lambdacdot Theta-m^2cdotTheta=0Leftrightarrowmathrm{cos}thetacdot sin thetacdotfrac{mathrm{d}Theta }{mathrm{d} theta }+frac{mathrm{d^{2}}Theta }{mathrm{d} theta^{2} }cdotsin ^{2}theta +sin ^{2}thetacdot lambdacdot Theta-m^2cdotTheta=0 Phi _{phiphi }+m^2cdotPhi=0Leftrightarrowfrac{mathrm{d^{2}}Phi }{mathrm{d} phi^{2} }+m^2cdotPhi=0 现在，我们通过分离变量法得到了三个常微分方程： left{begin{matrix} 2rcdot frac{mathrm{d} R}{mathrm{d} r}+r^{2}cdot frac{mathrm{d^{2}} R}{mathrm{d} r^{2}}-lambdacdot R=0 left( 1 right) mathrm{cos}thetacdot sin thetacdotfrac{mathrm{d}Theta }{mathrm{d} theta }+frac{mathrm{d^{2}}Theta }{mathrm{d} theta^{2} }cdotsin ^{2}theta +sin ^{2}thetacdot lambdacdot Theta-m^2cdotTheta=0 left( 2 right) frac{mathrm{d^{2}}Phi }{mathrm{d} phi^{2} }+m^2cdotPhi=0 left( 3 right) end{matrix}right. 这三个常微分方程求解的困难程度由低到高分别是： left(3 right),left(1 right),left(2 right) 。其中有关 Theta 的方程进行 t=mathrm{cos}theta 的变量代换之后 mathrm{d}t=-mathrm{sin}thetamathrm{d}theta ，可得到有关 t 的伴随 Legendre 方程，其中 tinleft[ -1,1 right] ,该方程需满足在该区间上取有限值，此时必须有： lambda=lleft( l+1 right) ，其中 linmathbb{N} 。这里就不进行求解了,因为求解过程实在是异常复杂。其实是我不会......（捂脸逃）

最终所求得的球谐函数的表达式为：

Y_{ml}left ( theta ,phi right ):=frac{1}{sqrt{2pi }}N_{lm}P_{lm}left ( cos theta right )expleft ( imphi right )

其中， P_{lm} 称为伴随 Legendre 多项式，其表达式为：

P_{lm}left( x right):=frac{left (-1 right )^{m}}{2^{l}cdot l!}left ( 1-x^{2} right )^{frac{m}{2}}frac{mathrm{d} ^{l+m}}{mathrm{d} x^{l+m}}left ( x^{2}-1 right )^{l}

伴随 Legendre 多项式是伴随 Legendre 方程的解。

N_{lm} 称为归一化常数：

N_{lm}:=sqrt{frac{2l+1}{2}cdot frac{left ( l-m right )!}{left ( l+m right )!}}

如果把上面的式子都写到一起，就是一个看起来非常可（e）爱（xin）的式子，这个式子就是所谓的球谐函数了，我们来一睹它的尊容：

Y_{ml}left ( theta ,phi right ):=frac{1}{sqrt{2pi }}cdotsqrt{frac{2l+1}{2}cdot frac{left ( l-m right )!}{left ( l+m right )!}}cdotfrac{left (-1 right )^{m}}{2^{l}cdot l!}left ( 1-cos^{2} theta right )^{frac{m}{2}}frac{mathrm{d} ^{l+m}}{mathrm{d} x^{l+m}}left (cos^{2} theta -1 right )^{l}cdotexpleft ( imphi right )

=frac{1}{sqrt{2pi }}cdotsqrt{frac{2l+1}{2}cdot frac{left ( l-m right )!}{left ( l+m right )!}}cdotfrac{left (-1 right )^{m+l}}{2^{l}cdot l!}sin^{m}theta frac{mathrm{d} ^{l+m}}{mathrm{d} x^{l+m}}sin^{2l} thetacdotexpleft ( imphi right )

部分球谐函数的表达式。图片来源：维基百科。
球谐函数的3D区域图像（m=n，l=m）。图片来源：维基百科。

得到球谐函数的表达式之后，我们还应该明确一个问题，就是单粒子球对称位势问题下的薛定谔方程的本征解 psileft( r,phi,theta right) 可以描述电子云的形状，该本征解恰好可以被拆成径向分量和角分量的乘积的形式，而其角分量正是球谐函数：

psileft( r,phi,theta right):=R_{nl}left( r right)cdot Y_{ml}left ( theta ,phi right )

下面，我们直接给出单粒子球对称位势问题下的薛定谔方程本征解的径向分量表达式：

R_{nl}=sqrt{left ( frac{2}{na_{0}} right )^{3}cdot frac{left ( n-l-1 right )!}{2n left ( n+l right )! }}cdot left ( frac{2r}{na_{0}} right )^{l}cdot L_{n-l-1}^{2l+1}left ( frac{2r}{na_{0}} right )cdot mathrm{exp}left ( -frac{r}{na_{0}} right )

其中： L_{s}^{k} 称为伴随 Laguerre 多项式，其表达式为：

L_{s}^{k}left ( x right )=frac{mathrm{d}^{s}}{mathrm{d}x^{s}}mathrm{exp}left ( x right )frac{mathrm{d}^{k}}{mathrm{d}x^{k}}left ( x^{k}cdot mathrm{exp}left ( -x right ) right )

a_0 是 Bohr 半径，其值是： 5.29cdot10^{-11}mathrm{m} 。

所以，电子云可以使用以下函数进行描述：

psileft( r,phi,theta right)=R_{nl}left( r right)cdot Y_{ml}left ( theta ,phi right )

其中：

径向本征解（径向波函数）：

R_{nl}=sqrt{left ( frac{2}{na_{0}} right )^{3}cdot frac{left ( n-l-1 right )!}{2nleft ( n+l right )! }}cdot left ( frac{2r}{na_{0}} right )^{l}cdot L_{n-l-1}^{2l+1}left ( frac{2r}{na_{0}} right )cdot mathrm{exp}left ( -frac{r}{na_{0}} right )

L_{s}^{k}left ( x right )=frac{mathrm{d}^{s}}{mathrm{d}x^{s}}mathrm{exp}left ( x right )frac{mathrm{d}^{k}}{mathrm{d}x^{k}}left ( x^{k}cdot mathrm{exp}left ( -x right ) right )

角向本征解（球谐函数）：

Y_{ml}left ( theta ,phi right ):=frac{1}{sqrt{2pi }}N_{lm}P_{lm}left ( cos theta right )expleft ( imphi right )

P_{lm}left( x right):=frac{left (-1 right )^{m}}{2^{l}cdot l!}left ( 1-x^{2} right )^{frac{m}{2}}frac{mathrm{d} ^{l+m}}{mathrm{d} x^{l+m}}left ( x^{2}-1 right )^{l}

N_{lm}:=sqrt{frac{2l+1}{2}cdot frac{left ( l-m right )!}{left ( l+m right )!}}

下标 n,l,m 是指三个量子数，分别是主量子数，角量子数和磁量子数，它们的意义分别是：

  • 主量子数 n ：电子层数，取值为 n=1,2,3,...
  • 角量子数 l ：轨道量子数（亚电子层数），取值为 l=0,1,2,...,n-1

    例如： n=3 ，则 l=0,1,2 ,分别 s,p,d 亚电子层（后面我们会看到相应的分量）

  • 磁量子数 m ：轨道角动量沿某指定轴的投影，取值为 m=0,pm1,pm2,...,pm l

    磁量子数的个数描述了 s,p,d,f,... 亚电子层有多少个分量。

    在上面的“部分球谐函数表示”图片中可以很清楚的看到 m 的取值对球谐函数分量个数的影响。

    下图是径向波函数在 n,l 取不同值时候的曲线：

    图片来源：谷歌图片搜索。

    以及瞩目的电子云图形：

    电子云。图片来源：谷歌图片搜索。

    我们知道，电子云是径向函数和球谐函数共同作用的结果。究竟是如何作用的呢？我们以 psi_{300} 为例来进行解释：

    第一图：psi_300
    第二图：n=3, l=0, m=0

    其中：

    psi_{300}=R_{30}cdot Y_{00}=frac{1}{sqrt{pi}}cdot left ( frac{Z}{3a_{0}} right )^{frac{3}{2}}cdot left ( 1-frac{2Zr}{3a_{0}}+frac{2Z^{2}r^{2}}{27a_{0}^{2}} right )cdot expleft (-frac{Zr}{3a_{0}} right )

    psi_{300}=0Leftrightarrowfrac{1}{sqrt{pi}}cdot left ( frac{Z}{3a_{0}} right )^{frac{3}{2}}cdot left ( 1-frac{2Zr}{3a_{0}}+frac{2Z^{2}r^{2}}{27a_{0}^{2}} right )cdot expleft (-frac{Zr}{3a_{0}} right )=0Leftrightarrow left ( 1-frac{2Zr}{3a_{0}}+frac{2Z^{2}r^{2}}{27a_{0}^{2}} right )=0 Rightarrow r=r_{1,2} 显然， r=0 并不是零点。 psi_{300}left( 0 right)=frac{1}{sqrt{pi}}cdot left ( frac{Z}{3a_{0}} right )^{frac{3}{2}}

    Z 为原子序数。

    从上面的 第二图可以看出， psi_{300} 有两个零点，三个极大值，以及在 rrightarrowinfty 时， psi_{300}rightarrow0 。若我们将极大值处想象成为“最亮的区域”，而零点和趋于零的部分想象成为“最暗的区域”，那么 psi_{300} 的“亮度变化”是：亮 rightarrow 暗 rightarrow 亮 rightarrow 暗 rightarrow 亮 rightarrow 无限趋于暗 。这个变化方式可以从上面的第一图中看出来。

    电子云描述的是在电子在空间某处出现的概率大小，上面所说的亮度刻画的就是这个概率。下图是电子云的空间结构：

    电子云的空间结构。从上至下分别是l=0,1,2,3。图片来源：维基百科。

    还没关注专栏《数学及自然科学》的朋友请赶快关注吧！您的支持是我最大的动力！

    文章标题: 电子云的形状原来是这么来的？!
    文章地址: http://www.xdqxjxc.cn/jingdianwenzhang/130870.html
    • 文章标签：科普  数学  物理科普
    Top