Mathematica 数值计算在量子力学谐振子中的应用.docxVIP

Mathematica数值计算在量子力学谐振子中的应用

一、基础背景：一维量子简谐振子

哈密顿量：

无量纲化（令α=mω/?

本征值λn=2n+1，能级

Mathematica优势：内置HermiteH、高精度数值求解、数值本征值、含时演化、概率密度可视化、微扰数值计算，可同时处理解析推导+数值仿真。

二、核心应用1：解析+数值绘制定态波函数、概率密度

完整代码

mathematica

(*无量纲一维谐振子*)

nMax=5;

(*归一化常数*)

Nn[n_]:=1/Sqrt[2^nn!Sqrt[Pi]]

(*本征波函数*)

Psi[n_,x_]:=Nn[n]Exp[-x^2/2]HermiteH[n,x]

(*概率密度*)

Rho[n_,x_]:=Abs[Psi[n,x]]^2

(*绘制前6个本征波函数*)

Plot[Table[Psi[n,x],{n,0,nMax}],{x,-4,4},

PlotLegends-Table[n=ToString[n],{n,0,nMax}],

AxesLabel-{ξ,ψ?(ξ)},PlotRange-All]

(*概率密度图*)

Plot[Table[Rho[n,x],{n,0,nMax}],{x,-4,