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

综合讨论区 Julia入门
1

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行求解释？

2

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

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

定义一个局域势。

 V = potentialoperator(b_position, potential)

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

Hkin = p^2/2m

哈密顿量动能部分。

H = Hkin + dense(V)

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

3

他这个库应该是QuantumOptics.jl…