导学：动边界势阱中的绝热演化
量子绝热定理 & 绝热不变量 · 模拟器的设计思路与图表解读
1. 这个模拟器是做什么的？
一个一维无限深势阱，阱壁按指定规律运动（收缩、膨胀、周期振荡等）。 阱中的粒子（一个量子态）随着阱壁运动而演化。 我们想考察的问题是：
- 阱壁运动得多"慢"，粒子才能始终留在基态？
- 阱壁运动得快了会怎样？粒子跃迁到高能级，绝热被破坏。
- 能量 ⟨E⟩ 与阱宽 l 是什么关系？绝热不变量 ⟨E⟩·l² 是否守恒？
2. 物理系统
一维无限深势阱的瞬时本征态能量为：
En(t) = π²ħ² / 2m·l(t)² · n²
系统的完整量子态可展开为瞬时本征态的叠加。系数 cn(t) 满足含时薛定谔方程，阱壁运动引起的能级间耦合为：
⟨Em|∂t|En⟩ = 2mn / m²−n² · ł / l · (−1)m−n
这个耦合是非绝热跃迁的根源：阱壁运动越快（ł 越大）、阱宽越小（l 越小），耦合越强。
3. 三张图怎么看
图1 · 基态 P₁
概率 |c₁|² 随时间变化。P₁ 接近 1 → 留基态；大幅下降 → 非绝热。图2 · 激发态
高能级 |cₙ|² 的演化。观察跃迁到哪些能级，反映耦合的分布。图3 · ⟨E⟩ vs 1/l²
平均能量随阱宽的变化。绝热时呈直线或窄回线，非绝热时偏离。4. 量子绝热定理
核心思想：如果哈密顿量 H(t) 变化得足够慢，初始处于第 n 个瞬时本征态的粒子将始终保持在同一个态：
|ψ(t)⟩ = cn(t) · |ψn(t)⟩ （绝热近似下）
能级间没有净跃迁：Pn = |cn|² ≈ const。
绝热不变量：对无限深势阱，若系统始终处于某一瞬时本征态，则：
I = En · l² = n²π²ħ²/(2m) = const
即 ⟨E⟩·l² 在绝热条件下守恒。图3正是对此的检验——绝热时 ⟨E⟩ vs 1/l² 呈直线，非绝热时偏离。
5. 绝热判据（矩阵元条件）
绝热定理成立的数学条件是能级间耦合远小于能级间距：
|⟨Em|∂t|En⟩| ≪ |Em − En| / ħ (m ≠ n)
代入矩阵元可得各模式的实用绝热带宽：
| 运动模式 | 主要参数 | 绝热条件（近似） | 非绝热效应 |
|---|---|---|---|
| 匀速收缩/膨胀 | 壁速 v | v ≲ 0.05 | v ≥ 0.2 时跃迁显著 |
| 周期振荡 | 频率 ω | ω ≲ 0.1 | ω ≥ 0.5 时大幅激发 |
| 先缩后胀 | 壁速 v | v ≲ 0.05（可逆） | v 大时不可逆 |
以上为经验参考，建议在模拟器中调节参数观察变化。
6. ⟨E⟩ vs 1/l² 图的解读
横轴 1/l²，纵轴 ⟨E⟩ = Σ Pn·n²π²/(2·l²)。不同参数下呈现不同特征：
- 绝热：⟨E⟩ ∝ 1/l²，呈直线（匀速）或窄回线（周期振荡），绝热不变量守恒。
- 轻度非绝热：直线微弯或回线增宽，不变量缓慢漂移。
- 强非绝热：曲线严重偏离线性。
- 数值崩溃（匀速收缩且 r < 0.3 时）：能量暴增到荒谬值，这是 RK4 步长不足的数值产物，非真实物理。
技巧：先跑一组绝热参数作为参照（ω=0.05, b=1 周期振荡，或 v=0.05 匀速收缩），再增大参数观察偏离。
7. 数值方法的说明
使用 RK4 积分 cn(t) 的薛定谔方程，步长自动取全程约 6000 步。
- 精度正常：阱壁速度和阱宽适中时，概率偏差 < 10⁻¹²
- 精度下降：匀速收缩 r < 0.3 时，特征时间尺度 ~ l² 急剧缩短，固定步长跟不上
- 应对：保持 r ≥ 0.3；增大 N 有时改善；右下角"概率和偏差"实时监控
8. 论文对照
本模拟器的物理基础基于顾卓成等（2023）关于动边界势阱绝热不变量的研究。
- 论文采用周期性边界驱动 l(t) = a − b·cos(ωt) 为主要模型
- 本模拟器保留论文核心公式，扩展了匀速收缩/膨胀、先缩后胀等模式
- 所有模式共享同一耦合矩阵元公式，仅阱壁运动学 l(t) 不同
- 矩阵元公式与论文公式(14)一致，完整耦合方程(5)是 RK4 积分的核心
- 论文 Fig.2 的 P₁ 演化对应模拟器的基态图