第五章-薛定谔方程数值解法-1.pptVIP

第五章 薛定谔方程数值解法;H可表示为;§5.1 一维方势阱的计算机求解;图5.1.1 方势阱 ;由于势能和时间无关，属于定态问题。考虑一维问题，故薛定谔方程(5.0.1)式可简化。为此设;由(5.1.3)式，直接可得;几率密度为;为简单起见，我们令 ，于是(5.1.6)式变为;在方势阱内：;在方势阱外：;其中j = 0，表示势阱左边的波函数，j = 1，表示势阱右边的波函数。 ;因为在 时，波函数要有物理意义必须有 。所以，可以得到： 和 。 ;下面求势阱内的波函数。 在 处（即左边界），波函数是连续的，根据(5.1.9)式和(5.1.11)式，有;可解得;令(5.1.14-1)=(5.1.14-2)，建立 和 满足的一个方程式。;由此解得 和 为;由上面过程可见，C 0和D 0的值完全确定了A1，B1，C1和D1的值，结果完全确定了波函数。一般情况下C1的值不一定为零。为了满足 时， 为零，必须要求C1=0。这就对E的取值进行了限制，不能取任意值，只能取某些确定值，才能保证C1=0的要求。;我们用图解来说明上面的情况。 设E1是能量本征值。对能量E可能有三种情况：E E1，E = E1，E E1。我们画出波函数，如图5.1.2所表示。由图可见，这3个波函数都满足当 时， 。但是，当 时，只有对于 的波函 数， ，即C1 = 0。这个波函数称为方程(5.1.7)式的本征波函数。;图5.2;对E E1的情况。由C1的表达式(5.1.14)和(5.1.15)式可知，C10，因此在 时，波函数向下发散。对E E1的情况，C1 0，当 时，波函数向上发散。;图5.1.3 本征值为 E2时不同能量 E对应的波函数;由上述讨论可知，粒子的能量只能取得某些 分裂的值，如图5.1.4所示。E的值与势阱的参数 V0和W有关，我们的中心问题是求解薛定谔方程 (5.1.7)式的本征波函数和能量本征值。 ;图5-4 一维方势阱内粒子的能级;我们进一步分析图5.1.2和图5.1.3中的波函数。波函数通过 x 轴的交点称为结点。在图5.1.2中能量E小于E1的的波函数没有结点，能量E大于E1的波函数有一个结点。在图5.1.3中能量E小于E2的波函数有一个结点，能量E大于E2的波函数有两个结点。;由此推广到一般规律：若 和 是薛定谔方程的两个解（不一定是本征波函数），而相应的能量是 和 ， 对应着n-1个结点， 对应着n个结点，则第n个能量本征值必然处在 ?? 之间，即 。更一般说，若 和 分别有n-1个和n+k个结点，则能量本征值 必然处在 和 之间。 ;根据上面分析，我们得到结论：E1的本征波函数没有结点E2的本征波函数有一个结点，E3的本征波函数有二个结点等等。对一维方势阱情况，我们用计算机通过确定结点来求解全部本征值和本征波函数。 ;第一步，给定参数。设 V=势阱的深度 W=势阱的宽度 Emax为猜测的能量本征值的上限。在本问题中，其值为零。由图5.1.4可见，若Emax大于零，势阱不起作用，变为自由粒子的运动。 Emin为猜测的能量本征值的下限。在现在的问题中，其值为 。因为由图5.1.4可见，粒子在势阱内运动，最小的可能能量为 。;第二步，将能量下限和能量上限之间分成M个能量。即从能量Emin开始，其增量为;因为 是半波长的数目。如果它小于1，取值为零，没有结点。若大于1而小于2，取值为1，有一个结点。若大于2而小于3，取值为2，有两个结点。如此下去，即可求得全部结点情况。在FORTRAN语言中，实数可以自动转化为整数。在方势阱右边，波函数是否与x轴相交，引起对结点数贡献，这需要计算波函数，由波函数的系数C1和D1决定。若 ，则结点数增加1，若 ，则对结点数没有影响。;第三步，由结点数计算结果，定出能级E1，E2，…。 根据前面分析