уравнение Шредингера; проблема в расходящемся решении - PullRequest
Новая
       8

уравнение Шредингера; проблема в расходящемся решении

0 голосов
/ 19 ноября 2018

Я пытаюсь построить волновую функцию до и после ступенчатой ​​функции, например, потенциала.

Я обнаружил, что если начальная энергия больше, чем потенциал, она выглядит хорошо, но если потенциал больше, тоРешение будет расходиться отрицательно.я думаю, что решение в таком случае: 

enter code here

psi(x)=A*exp(x)+B*exp(-x)

В квантовой механике мы обычно отказываемся от первого слагаемого, устанавливая А в ноль.Тем не менее, я также хочу удалить такое нефизическое решение, но я не знаю, как это сделать в python.

, поэтому мой вопрос заключается в том, есть ли способ устранить расходящиеся решения, когда я решил ODE с помощью odeint??Решение после функции шага должно быть exp (-x) как функция.

Большое спасибо 

enter code here

import matplotlib.pyplot as plt
%matplotlib inline
from scipy.integrate import odeint
import numpy as np

def ps(y,t):
    if t<10:
        u=0
    else:
        u=2
    psi=y[0]
    dpsi=y[1]
    dp=dpsi
    dp2=(-1+u)*psi
    return [dp,dp2]
yo=[1,0]
t=np.linspace(0,100,400)
psis=odeint(ps,yo,t)    
plt.ylim(-10,10)
plt.xlim(0,20)
plt.plot(t,psis[:,1],'r')

plt.plot([0,10,10,40],[0,0,2,2],'g')
...