第十部分_波动方程和薛定谔方程_郑大昉.ppt

第十章　波动方程和薛定谔方程 波动方程 一维波动方程的差分方法 薛定谔方程与氢原子能级 * * (10-1) 以一维弦线横振动为例，波动方程为： 其中: ---为弦线的 线质量密度. (10-2) 若考虑均匀弦线，即：　　　　，则得一维非齐次波动方程： ---为弦线中的张力. ---为外加的强迫力. 其中: （波速）， 其相应的齐次方程为： (10-3) 上式可推广至二维及三维情形，即： (10-4) 与热传导方程比较，波动方程中对时间的偏微商是二阶的． 及： (10-5) 对其波动方程： (10-6) 其解法思路仍然是差分格式． 选： (10-7) 则： (10-8) 差分格式 其二阶中心差分近似为（参看热传导方程做法）： (10-9) 其中： (10-10) 代入 (10-6) 便得波动方程的差分格式： (10-11) (10-12) 上述差分格式的几何意义： a. 初始条件 对 (10-14) 一般采用如下两种差分格式： (10-13) (10-15) 定解条件 (10-14) (1) 向前一阶差商 即： (10-16) 令：　　　　则 (10-11) 变为： 上述出现虚设的　　　，但可联立差分格式(10-11)求出． (10-18) (10-17) (2) 中心一阶差商 即： (10-19) 联立(10-18)及(10-19) 解得： (10-20) b. 边界条件 综合上述定解条件，得以下两种差分格式解： (10-21) 第一种差分格式： (10-22) 应留意以上各 的取值范围. 第二种差分格式： (10-23) 同样应留意以上各 的取值范围. 计算步骤 1. 给定: 2. 计算: 3. 计算: 4. 计算初值(有两种格式供选)和边值. 5. 计算: 例: 用第一种差分格式计算波动方程混合问题: 取: , 计算 的值. 解: 参数值: 初值: 按第一种差分格式有: 初值: 边值: 利用: 可得: 以上已全部算出所需要的值. 以氢原子为例,下面讨论定态薛定谔方程的数值解及能级的确定. 定态薛定谔方程为: (10-24) 对氢原子: (10-25) 则: (10-26) 其中: 其定态波函数: (10-27) ---为径向波函数 ---为球谐函数 对于束缚态,其能量是量子化的,且仅由径向波函数满足的径向方程确定,即: (10-28) 事实上, 能量是 度简并的, 与角动量量子数 无关. 故仅就 对上述方程做数值解, 同时可确定其玻尔能级. 整理得: 其方程进一步简化为: (10-29) 令: 显然, 方程 (10-32) 可用龙格-库塔(R-K)方法求解. (10-30) (10-31) 则: (10-32) 其中: 定态解及能级确定思路: 波函数应满足标准条件,即应是单值、有限且连续的. 例: 选初值: (以上初值的确定可参看量子力学书的有关讨论) 1. 取 (1s态) 则可得其收敛解,因此 该能级即为对应的 n=1 的玻尔能级