所以我需要在python中定义潜力。我具有无限的方阱势,并且在中间具有势垒。

对于障碍我有这个:

    V_x = zeros([Npts],float)
for i in range(0,Npts,1):
    if x[i] > 0 and x[i]<width:
        V_x[i]=v0


Npts是x的长度,x定义为从xmin到xmax的递增值数组(x = arange(Xmin,Xmax + 0.001,dx))。

如何将无限潜力纳入其中?

它具有一定的长度(例如,从-100到100),并且其行为必须像不可穿透的墙(-100和100处的函数必须为零)。

我可以单独将势垒与该势垒组合在一起(例如V_b势垒,而V_p势无限电势,则V_x可以为V_x = V_b + V_p)。

您认为这可行吗?我有定义无限潜力的问题...

最佳答案

我建议仅对您的无限潜力使用一个比其有限潜力大几个数量级(3-4)的数字。将其增大到更大不会对您的状态产生太大的影响。实际上,求解哈密顿量时,所使用的浮点数的精度很快就会给您带来更多的限制。您实际上只需要将端点设置为“无限”,这样就可以进行以下操作:V_x[0] = V_x[-1] = abs(v0)*10**4。当然,您可以根据需要将其设置为更大,但我认为这不会有太大的不同。如果不将端点设置为无限,那么为什么要在您知道波动函数为零的区域的大部分区域中尝试求解波动函数的值?

另外,请注意,您所使用的单位的无限潜力与h-bar^2/2/m/delta_x^2不可比。所有这些使“无限”势比汉密尔顿方程中的有限势和动能项都大几个数量级。

关于python - 结合python中的两种潜力,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/5084295/

10-10 08:30