量子力学中势阱的哈密顿量的算法

xmin = -3
xmax = 3
Nsteps = 200
L = xmax - xmin
m = 0.5
E0 = 20

b_position = PositionBasis(xmin, xmax, Nsteps)
xpoints = samplepoints(b_position)

x = position(b_position)
p = momentum(b_position)

potential = x -> E0 + x^4 - 8*x^2
V = potentialoperator(b_position, potential)
Hkin = p^2/2m
H = Hkin + dense(V)
Julia小白
程序后4行求解释?

这里应该是调用了某个库,建议附上库的链接。以下是我根据函数命名的猜测:


potential = x -> E0 + x^4 - 8*x^2

定义一个局域势。

 V = potentialoperator(b_position, potential)

这里看函数名称是指生成势能算符。

Hkin = p^2/2m

哈密顿量动能部分。

H = Hkin + dense(V)

动能和势能相加,求得全部的哈密顿量。

他这个库应该是QuantumOptics.jl

备案号:京ICP备17009874号-2