拉盖尔多项式 matlab,类氢原子的定态波函数

&nbsp

&nbsp

&nbsp

&nbsp

&nbsp

&nbsp

&nbsp

预备知识 球坐标系中的定态薛定谔方程,原子单位制

本文使用原子单位制.类氢原子(hydrogen-like atom)被定义为原子核有 $Z$ 个质子(核电荷为 $+Ze$)有一个核外电子的原子/离子,例如氢原子和失去一个电子的氦原子 $\mathrm{He}^+$,失去两个电子的锂离子 $\mathrm{Li}^{++}$.

类氢原子的定态薛定谔方程为

\begin{equation}

-\frac{1}{2m} \boldsymbol{\nabla}^2 \psi( \boldsymbol{\mathbf{r}} ) - \frac{Z}{r} \psi( \boldsymbol{\mathbf{r}} ) = E \psi( \boldsymbol{\mathbf{r}} )

\end{equation}

类氢原子是唯一存在解析解的原子(离子).

我们这里只讨论束缚态,即 $E < 0$ 的解. 从数学上,$E$ 取小于零的任意值时我们都能找到解,但只有当 $E$ 取特定离散值的时候这些波函数才能归一化(否则没有物理意义).由于类氢原子具有球对称性,球坐标下的波函数具有最简单的形式.波函数的表达式为

\begin{equation}

\psi_{nlm} (r,\theta ,\phi) = R_{nl}(r) Y_{l,m}(\theta, \phi)

\end{equation}

其中 $n$ 是主量子数($n = 1, 2, \dots$),$l$ 是角量子数($l = 0, 1, \dots, n - 1$),$m$ 是磁量子数($m = -l, -l+1, \dots, l$).$R_{nl}(r)$ 是归一化的径向波函数,$Y_{l,m}(\theta, \phi)$ 是归一化的球谐函数(见 “球谐函数列表”).

77f47b6a59d662374621a21a49c3600d.png

图 1:氢原子波函数的概率密度 $ \left\lvert \psi_{nlm}( \boldsymbol{\mathbf{r}} ) \right\rvert ^2$ 的 $x$-$z$ 截面,大小成比例.每个图中的三个数字分别是量子数 $n, l, m$,电子出现在白色圆圈内部的概率为 0.95.色条对应的数值是线性的,每个子图中的色条取值范围不相同.另见 Matlab 画图程序.

径向波函数 $R_{nl}(r)$

如果忽略原子核的运动,以下的 $a$ 是玻尔半径(原子单位中 $a=1$),如果不忽略,$a$ 就是约化玻尔半径.注意 $Z$ 和 $a$ 的作用是把径向波函数关于原点收缩 $Z/a$ 倍(并保持波函数归一化).

\begin{equation}

R_{nl}(r) = \sqrt{ \left(\frac{2 Z}{na} \right) ^3 \frac{(n - l - 1)!}{2n (n + l)!}} \left(\frac{2Zr}{na} \right) ^l L_{n-l-1}^{2l+1} \left(\frac{2Zr}{na} \right) \mathrm{e} ^{-Zr/(na)}

\end{equation}

其中 $L_n^l(x)$ 是连带拉盖尔多项式(associated Laguerre polynomial).

$Z = 1$ 时 $r R_{n,l}(r)$ 的函数图图 2

HWF_1.svg?v=d8b407c574892e

图 2:径向波函数函数图(使用原子单位,$a = 1$)

以下给出前几个径向波函数,注意所有径向波函数的值都是实数.

\begin{equation} n = 1 \qquad

R_{10}(r) = 2 \left(\frac{Z}{a} \right) ^{3/2} \exp\left(-Zr/a\right)

\end{equation}

\begin{equation}

n = 2 \qquad

\left\{\begin{aligned}

R_{20}(r) &= \frac{1}{\sqrt 2} \left(\frac{Z}{a} \right) ^{3/2} \left(1 - \frac12 \frac{Zr}{a} \right) \exp\left(-\frac{Zr}{2a}\right) \\

R_{21}(r) &= \frac{1}{2\sqrt{6}} \left(\frac{Z}{a} \right) ^{3/2} \frac{Zr}{a} \exp\left(-\frac{Zr}{2a}\right)

\end{aligned}\right. \end{equation}

\begin{equation}

n = 3 \qquad

\left\{\begin{aligned}

R_{30}(r) &= \frac{2}{3\sqrt {3}} \left(\frac{Z}{a} \right) ^{3/2} \left(1 - \frac23 \frac{Zr}{a} + \frac{2}{27} \frac{Z^2r^2}{a^2} \right) \exp\left(-\frac{Zr}{3a}\right) \\

R_{31}(r) &= \frac{8}{27\sqrt 6} \left(\frac{Z}{a} \right) ^{3/2} \left(1 - \frac16 \frac {Zr}{a} \right) \frac {Zr}{a} \exp\left(-\frac{Zr}{3a}\right) \\

R_{32}(r) &= \frac{4}{81\sqrt {30}} \left(\frac{Z}{a} \right) ^{3/2} \frac{Z^2r^2}{a^2} \exp\left(-\frac{Zr}{3a}\right) \end{aligned}\right.

\end{equation}

更多 $R_{n,l}$ 可以用 Mathematica 或者 Wolfram Alpha 生成,命令如

性质

我们要求氢原子每个束缚态满足归一化条件

\begin{equation}

\int \left\lvert \Psi_{n,l,m}( \boldsymbol{\mathbf{r}} ) \right\rvert \,\mathrm{d}^{3}{r} = \int \int_0^\infty \left\lvert R_{n,l}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} ) \right\rvert ^2 r^2 \,\mathrm{d}{r} \,\mathrm{d}{\Omega} = 1

\end{equation}

先做角向积分,球谐函数已经满足归一化条件

\begin{equation}

\int \left\lvert Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} ) \right\rvert ^2 \,\mathrm{d}{\Omega} = 1

\end{equation}

得到径向波函数的归一化条件为

\begin{equation}

\int [rR_{n,l}(r)]^2 \,\mathrm{d}{r} = 1

\end{equation}

再来看正交性,我们知道哈密顿算符的本征态之间两两正交(式中 $n',l',m'$ 至少有一个与 $n, l, m$ 不同)

\begin{equation}

\int \Psi_{n,l,m}^*( \boldsymbol{\mathbf{r}} ) \Psi_{n',l',m'}( \boldsymbol{\mathbf{r}} ) \,\mathrm{d}^{3}{r}

= \int \int_0^\infty R_{n,l}(r) R_{n',l'}(r) Y_{l,m}^*( \hat{\boldsymbol{\mathbf{r}}} ) Y_{l',m'}( \hat{\boldsymbol{\mathbf{r}}} ) r^2 \,\mathrm{d}{r} \,\mathrm{d}{\Omega} = 0

\end{equation}

同样可以先对角向做积分,若 $l',m'$ 至少有一个与 $l, m$ 不同,那么积分直接为零,径向波函数不需要任何正交条件(也的确不满足).但若两个球谐函数相同,即 $l' = l$,$m' = m$,$n' \ne n$,那么角向积分等于 1,径向波函数满足

\begin{equation}

\int \int_0^\infty R_{n,l}(r) R_{n',l}(r) r^2 \,\mathrm{d}{r} = 0

\end{equation}

径向概率分布

我们来求径向概率分布 $P(r)$.$P(r)$ 的定义为:发现粒子在 $r \in [a, b]$(厚球壳)内的概率为 $\int_a^b P(r) \,\mathrm{d}{r} $.由于波函数的模长平方就是三维的概率密度,有

\begin{equation}

\int_a^b P(r) \,\mathrm{d}{r} = \int_a^b \int \left\lvert \frac1r \psi_{l,m}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} ) \right\rvert ^2 \,\mathrm{d}{\Omega} r^2 \,\mathrm{d}{r}

= \int_a^b \left\lvert \psi_{l,m}(r) \right\rvert ^2 \,\mathrm{d}{r}

\end{equation}

对任意 $a, b > 0$ 都成立,所以有

\begin{equation}

P(r) = \left\lvert \psi_{l,m}(r) \right\rvert ^2

\end{equation}

任意波函数可以表示为所有本征波函数的叠加

\begin{equation}

\Psi( \boldsymbol{\mathbf{r}} ) = \frac1r \sum_{l,m} \psi_{l,m}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} )

\end{equation}

其径向概率分布为

\begin{equation}

\int_a^b P(r) \,\mathrm{d}{r} = \int_a^b \int \left\lvert \frac1r \sum_{l,m}\psi_{l,m}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} ) \right\rvert ^2 \,\mathrm{d}{\Omega} r^2 \,\mathrm{d}{r}

= \sum_{l,m} \int_a^b \left\lvert \psi_{l,m}(r) \right\rvert ^2 \,\mathrm{d}{r}

\end{equation}

对任意 $a, b > 0$ 都成立,所以有

\begin{equation}

P(r) = \sum_{l,m} \left\lvert \psi_{l,m}(r) \right\rvert ^2

\end{equation}

动量表象 动量分布

要求动量表象下的波函数,我们需要将位置表象的波函数投影到归一化的动量的本征矢上,即三维傅里叶变换

\begin{equation}

\psi_{nlm}( \boldsymbol{\mathbf{p}} ) = \left\langle \boldsymbol{\mathbf{p}} \middle| \psi \right\rangle = \frac{1}{\sqrt{2\pi}} \int \exp\left(- \mathrm{i} \boldsymbol{\mathbf{p}} \boldsymbol\cdot \boldsymbol{\mathbf{r}} \right) \psi( \boldsymbol{\mathbf{r}} ) \,\mathrm{d}^{3}{r}

\end{equation}

这个积分在球坐标中完成才是最方便的,具体方法我们将举例子说明(见例 1).

正如位置表象下位置的分布函数是 $ \left\lvert \psi( \boldsymbol{\mathbf{r}} ) \right\rvert ^2$,动量表象下动量的分布函数是 $ \left\lvert \psi( \boldsymbol{\mathbf{p}} ) \right\rvert ^2$(也符合测量理论).

1. ^ 我们以后会看到 $r R_{l,m}(r)$ 比 $R_{l,m}(r)$ 更常用

致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择会员制,大量广告,内容付费等。 因此,我们请求广大读者热心打赏,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 10 元,我们一个星期内就能脱离亏损, 并保证网站能在接下来的一整年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/356751.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

再论EM算法的收敛性和K-Means的收敛性

标签&#xff08;空格分隔&#xff09;&#xff1a; 机器学习 &#xff08;最近被一波波的笔试面试淹没了&#xff0c;但是在有两次面试时被问到了同一个问题&#xff1a;K-Means算法的收敛性。在网上查阅了很多资料&#xff0c;并没有看到很清晰的解释&#xff0c;所以希望可以…

java xml dom getelementbyid,DOM中常见的元素获取方式

1.getElementById获取元素 返回的是一个元素对象var timer document.getElementById(time);console.dir 打印返回元素对象&#xff0c;更好的查看里面的属性和方法console.dir( timer );2.getElementsByTagName 获取某类标签元素 返回的是 获取过来元素对象的集合 以…

杰尔·地狱

什么是JAR地狱&#xff1f; &#xff08;或者是classpath地狱&#xff1f;还是依赖地狱&#xff1f;&#xff09;在考虑使用Maven或OSGi等现代开发工具时&#xff0c;哪些方面仍然有意义&#xff1f; 有趣的是&#xff0c;似乎没有对这些问题的结构化答案&#xff08;即&#…

matlab radsimp,[转载]MATLAB学习笔记(八)

符号运算符号运算分为以下几类&#xff1a;符号表达式和符号矩阵的操作整体定义为符号微积分符号线性方程符号微分方程A、符号变量、符号表达式和符号方程的分解一、 生成符号变量要使用sym和syms&#xff1a;使用sym函数可以定义符号表达式&#xff0c;此时有两种定义方法&…

windows,python3.x下安装pyspider

由于是初学者&#xff0c;业余学习&#xff0c;习惯使用windows&#xff0c;初次了解到pyspider写代码和调试代码非常简便&#xff1b;作者binux是在Ubuntu下部署测试的。在作者的博客看到windows下安装的讨论。windows直接安装失败主要是lxml、pycurl安装失败&#xff0c;需要…

matlab实验符号计算答案,实验五matlab符号计算

实验五matlab符号计算 实验 5 符号计算 教师评分班级 学号 姓名实验日期 2014 年 6 月 17 日 星期 二 第 1 至 2 节课实验地点实验目的1. 掌握定义符号对象的办法2. 掌握符号表达式的运算法则以及符号矩阵运算3. 掌握求符号函数极限及导数的方法4. 掌握求符号函数定积分和不定积…

junit4 单元测试框架_超越JUnit –测试框架的替代方案

junit4 单元测试框架JUnit是事实上的Java单元测试框架&#xff0c;但是可能有一些新的&#xff08;不是那么新的&#xff09;框架可以用于Web开发。 在采用之前可能要问自己的问题&#xff1a; 它们是否快速&#xff0c;容易开发并因此成本低廉&#xff1f; 他们运行快并因此鼓…

Java学习笔记之:Java String类

一、引言 字符串广泛应用在Java编程中&#xff0c;在Java中字符串属于对象&#xff0c;Java提供了String类来创建和操作字符串。 创建字符串最简单的方式如下: String str "Hello world!"; String类型是特殊的引用类型&#xff0c;我们也可以通过实例化的方式来创建 …

java循环输入直到,使用循环接受其他用户输入,直到用户输入结束输入的值

我是Java的新手 . 我需要一些帮助&#xff0c;使用循环接受其他用户输入&#xff0c;直到用户输入结束输入的值 . 我的问题从语句“System.out.println(”你完成了吗&#xff1f;输入大写的Y / N)开始 . 下面是我的代码 .公共类EmployeeData {//declare variablesprivate Strin…

WildFly 10 CR 2发布– Java EE 7,Java 8,Hibernate 5,JavaScript支持热重载

昨天&#xff0c;WildFly团队发布了最新版本的WildFly 10 。 CR2很可能是预计于十月份发布最终版本之前的最后一个版本。 即使主要支持的Java EE规范是7&#xff0c;WildFly 8和WildFly 9仍具有许多新功能&#xff0c;该版本现在具有三个服务器版本&#xff0c;实现了Java EE 7…

php中什么时候用传值,php中传值与传引用的区别。什么时候传值什么时候传引用?...

java中的this与super的区别java中的this与super的区别 1. 子类的构造函数如果要引用super的话,必须把super放在函数的首位 代码如下: class Base { Base() { System.out.pr ...php传值和传引用的区别php传值和传引用的区别所谓值传递,就是说仅将对象的值传递给目标对象,就相当于…

openshift 部署_OpenShift Express:部署Java EE应用程序(支持AS7)

openshift 部署在过去的几年中&#xff0c;我越来越多地听说过“云”服务。 最初&#xff0c;我并不是很想尝试一下。 但是几个月后&#xff08;一年&#xff1f;&#xff09;&#xff0c;我决定看看这是怎么回事。 我从事Java EE开发已有7年以上&#xff0c;因此&#xff0c;我…

迪克逊准则matlab,浙江科技学院学报

引言传统的农业栽培需要遵循季节性及周期性的栽培规律,受外界环境的影响较大,相应的栽培技术无法得到良好的栽培结果[1]。而温室大棚通过环境控制器调节温室环境,为作物提供了合适的生长环境。因此,它在农业领域的应用越来越广泛,温室环境控制系统的研究成为人们关注的热点[2-5…

Linq表达式和Lambda表达式用法对比

什么是Linq表达式&#xff1f;什么是Lambda表达式&#xff1f;前一段时间用到这个只是&#xff0c;在网上也没找到比较简单明了的方法&#xff0c;今天就整理了一下相关知识&#xff0c;有空了再仔细研究研究 public Program() { List<Student> allStudent new List<…

使用FlexDeploy对融合中间件应用程序进行自动化软件测试

自动化软件测试是任何软件组织都应执行的强制性活动之一&#xff0c;以保证其产品质量。 但是&#xff0c;此过程通常变得非常复杂&#xff0c;尤其是涉及由多个不同部分组成的现代复杂系统的自动化测试时。 所有这些部分都基于不同的技术&#xff0c;显然&#xff0c;应该使用…

uploadify插件的使用

插件&#xff1a; uploadify.css jquery.uploadify.js bootstrap html代码&#xff1a; <input type"file" name"uploadify_coverimg" id"uploadify" /> <div id"the_coverimg"></div> js代码&#xff1a; $("…

wxlogin php,wxlogin.php

if($_GET[echostr]){echo $_GET[echostr];die();}/*** 开发规范&#xff1a;* 1.不同版本不同控制器以及模板* 2.不同版本不同数据库&#xff0c;但是对应数据表表结构必须一致* 3.不同版本共用service层&#xff0c;所以修改表结构必须所有版本统一*/// ---------------------…

关于人生倒计时的一个小玩意,纯属业余

人生倒计时 人生倒计时 出生年份&#xff1a; 出生月份&#xff1a; <!DOCTYPE html> <html id"spLianghui"> <head><meta http-equiv"Content-Type" content"text/html; charsetgb2312" /><title>人生倒计时<…

java 反编译项目_Java 7 –反编译项目硬币

java 反编译项目大家好&#xff0c;该是从2012年开始写作的时候了。正如您在其他博客中可能已经看到的那样&#xff0c;有一些更改可以使您使用Java编程时的开发人员生活变得更加轻松&#xff1a;Diamond运算符&#xff0c;Switchs中的Strings&#xff0c;尝试使用资源&#xf…

日期格式无法识别 oracle,Oracle时间统计 - 不以物喜-不以己悲的个人空间 - OSCHINA - 中文开源技术交流社区...

select count(*) from Member_user where to_char(CREATEDATE,yyyy-mm-dd)to_char(sysdate,yyyy-mm-dd); --今天select count(*) from Member_user where to_char(CREATEDATE,dd)to_char(sysdate,dd); --当天&#xff0c;所有月份的这个号select count(*) from Member_user whe…