变分法数值求解资料.ppt

理论探讨 双参数变分模型 四参数变分模型 实验值 -2.84766 -2.87566 -2.9033 ①后者试探波函数更符合实际情况，因此精度高 ②若用同样方法考虑类氦原子，如 等，精度将更高！ ③目前，有人利用39个变分参数，计算出-2.9037， 与实验值-2.9033更接近 Q A 我们可以直接由薛定谔方程求解势场中存在束缚态的精确条件，但是当定态方程不能求解析解的时候，用变分法相对简单些。 下面我们比较变分法求解势场中存在束缚态问题与精确解法的区别。 用变分法研究势场 中存在束缚态的条件 研究对象 1. 2. 3. 4. 5. 6. 首先变分法登场 1.用变分法求中心力场存在束缚态的条件 设中心力场V(r)满足条件V(r)0,及当r→∞时,V(r)→0。如果在此势场中只存在一个束缚态,则必定是l=0的球对称基态Ψ(r)。取归一化基态试探波函数 （1）计算基态近似能量： （2） 由E(α)取极小值条件 ，求出α，代 入式（1）与式（2），得基态近似波函数 与能量 ，令 ，式（2）变 为 (3) 由于在具有负值的势场 中束缚态能量 ，便有 (4) 或 (5) 由式（5）可得势场V(r)存在束缚态的条件，不难看出由极值条件 （6） 可以求得，因此，将式（6）求出的不为0的α代入式（5），可得势场存在束缚态的条件。当α=0时，E=0（见式3），束缚态不能形成，如果不存在 的极大值，则α要由 （7） 决定。 为验证上述结论，将 代入式（3）， 求出 （8） 将式（8）代入极值条件（7），给不出使E(α)取极小值的α的解析式，只能通过数值计算求Emin，为此令势场参数 （9） f是表示势场 强度的量纲参数。 再令 ，式（8）变为 （10） 其中 ，选定f，即确定势场强度， 由式（10）计算不同x值的E,找出E的最小值 Emin,由表可以看出， 是势场 存在束缚态的条件，这个条件正是将 取极大值 得到的，这就表明，由方程（6）求出α代 入（5）得到的束缚态条件是正确的。 1.对