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

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,4},

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

AxesLabel->{"ξ","|ψₙ(ξ)|²"},PlotRange->{0,0.8}]

(*数值验证归一化积分∫_{-∞}^∞|ψₙ|²dx=1*)

Table[N[Integrate[Rho[n,x],{x,-Infinity,Infinity}]],{n,0,5}]数值计算点符号积分验证归一化，数值浮点输出检验精度；批量生成多能级波函数、概率密度；自动调用高精度特殊函数HermiteH，无需手动编程厄米多项式。三、核心应用2：数值求解薛定谔方程（有限差分/NDSolve）不依赖解析解，直接数值求解二阶常微分本征问题，适合任意势变形谐振子（非标准、含微扰）。1.直接数值求解定态方程mathematica

(*无量纲薛定谔方程，求基态n=0*)

eq=ψ''[x]+(λ-x^2)ψ[x]==0;

bc1=ψ[-5]==0;bc2=ψ[5]==0;

sol=NDSolve[{eq,bc1,bc2},ψ,{x,-5,5},

Eigenvalues->{0,10},Method->{"Shooting","StartingInitialConditions"->{0,1}}]

(*提取最低本征值，理论λ=1,3,5...*)

λval=λ/.sol[[1]]2.含微扰谐振子数值求解（如V(x)=12x标准谐振子无解析解，只能数值计算能级偏移：mathematica

eps=0.1;

V[x_]=1/2x^2+epsx^4;

(*定态方程：-ψ''/2+V[x]ψ=Eψ*)

pertEq=-1/2ψ''[x]+V[x]ψ[x]==Eψ[x];

pertSol=NDSolve[{pertEq,ψ[-6]==0,ψ[6]==0},ψ,{x,-6,6},Eigenvalues->{0,4}]

(*输出微扰后前5个能级，对比无扰E=0.5,1.5,2.5...*)

Elist=E/.pertSol[[1]]优势：Mathematica内置打靶法、有限元本征求解器，自动自适应网格，高精度浮点数（WorkingPrecision可调）。四、核心应用3：含时量子演化（数值传播波包）相干态是谐振子经典对应，数值模拟波包振荡、无扩散；可数值计算任意初态随时间演化。相干态含时演化代码mathematica

(*相干态解析形式，无量纲*)

α=1.5;(*相干态复振幅*)

PsiCoherent[x_,t_]:=Exp[-1/2(x-Sqrt[2]αCos[t])^2-I(α^2Sin[2t]/2-Sqrt[2]αxSin[t])]/Pi^(1/4)

(*概率密度随时间动画*)

Animate[Plot[Abs[PsiCoherent[x,t]]^2,{x,-5,5},

PlotRange->{0,0.6},AxesLabel->{"ξ","|ψ(x,t)|²"}],{t,0,2Pi}]通用数值含时演化（NDSolve求解含时薛定谔PDE）mathematica

(*含时薛定谔i∂ψ/∂t=-1/2∂²ψ/∂x²+x²/2ψ*)

tdse=ID[ψ[x,t],t]==-1/2D[ψ[x,t],{x,2}]+1/2x^2ψ[x,t];

(*初态：基态波包*)

init=ψ[x,0]==Exp[-x^2/2]/Pi^(1/4);

(*边界远场为0*)

bc={ψ[-6,t]==0,ψ[6,t]==0};

(*数值求解PDE*)

tdsol=NDSolve[{tdse,init,bc},ψ,{x,-6,6},{t,0,2Pi}]

(*取解绘图*)

ψnum[x_,t_]=ψ[x,t]/.tdsol[[1]];

Plot[Abs[ψnum[x,Pi/2]]^2,{x,-6,6}]数值特点：自动时空自适应网格，复值波函数高精度传播，可观测量子复兴、波包压缩。五、核心应用4：矩阵数值对角化（离散表象：位置/福克空间）1.福克表象（粒子数表象）矩阵对角化升降算符：哈密顿H^mathematica

dim=10;(*截断希尔伯特空间维度*)

(*产生湮灭算符矩阵*)

a=SparseArray[{{i_,j_}/;j==i+1->Sqrt[i]},{dim,dim}];

adag=Transpose[a];

Hmat=adag.a+1/2IdentityMatrix[dim];

(*数值对角化求能级*)

Eigensystem[N[Hmat]]输出本征值精确等于n+1/2，可扩展到多模谐振子、耦合谐振子系统。2.位置空间离散（有限差分矩阵）将区间离散为格点，二阶差分构造哈密顿矩阵，数值对角化得到波函数与能级，适合教学直观演示离散化数值方法。六、核心应用5：量子期望值、跃迁矩阵元数值积分计算坐标、动量、能量期望值，微扰矩阵元、跃迁概率：mathematica

(*计算<ψ_m|x|ψ_n>矩阵元*)

MatrixEl[m_,n_]:=N[Integrate[Conjugate[Psi[m,x]]*x*Psi[n,x],{x,-Infinity,Infinity}]]

(*生成5×5矩阵元表*)

Table[MatrixEl[m,n],{m,0,4},{n,0,4}]

(*基态坐标期望值<0|x|0>*)

N[Integrate[Conjugate[Psi[0,x]]xPsi[0,x],{x,-∞,∞}]]

(*坐标平方涨落Δx=Sqrt[<x²>-<x>²]*)

ExpX2=N[Integrate[Conjugate[Psi[0,x]]x^2Psi[0,x],{x,-∞,∞}]]

Sqrt[ExpX2]Mathematica混合符号-数值积分，对无穷区间振荡积分自动采用高精度数值算法，避免手动分段近似。七、拓展应用：多维谐振子、耦合谐振子、含时微扰数值模拟二维/三维各向同性谐振子：分离变量+高维数值积分、等高概率密度图；耦合双谐振子：构造高维哈密顿矩阵数值对角化，数值求解简并分裂、正常模；含时微扰（光场驱动谐振子）：数值求解含时薛定谔，计算跃迁概率随时间变化；热谐振子（正则系综）：数值计算配分函数、平均能量、热概率分布。八、Mathematica数值计算相比Python/MATLAB的独特优势符号-数值无缝衔接：先符号推导哈密顿、矩阵元，再一键转高精度数值，无需手动离散；内置特殊量子函数：HermiteH、LaguerreL、相干态、厄米函数原生支持；自适应高精度求解器：NDSolve自动处理ODE/PDE本征问题，可自定义Wor