双势阱、势垒与分界轨道

两个势阱中间隔着势垒，粒子的运动会怎样变？ 看能量能不能越过势垒顶。低于垒顶，粒子在一个阱里来回振荡；高于垒顶，就在两个阱之间来回跑。能量接近垒顶时，粒子在顶上附近停留很久，周期急剧变长。有阻尼的话，耗掉一点能量，就能把本来跨阱的粒子困进一个阱里。

给粒子刚好够越过势垒的能量。比较单阱振动、鞍点附近的缓慢穿越，再看阻尼怎样让它最终留在一个阱里。

质点位置 $0$

速度 $0$

机械能 $0$

高于鞍点的能量 $0$

选中时刻 $0$

保守力 $0$

黏性阻力 $0$

势能 $0$

动能 $0$

耗散能量 $0$

局部小振动角频率 $0$

选定极小值 $0$

初始能量左边界 $0$

初始能量右边界 $0$

保守能量积分周期 $0$

全记录最大能量平衡误差 $0$

半步长位置差 $0$

完整记录的鞍点穿越次数 $0$

周期积分分辨率差 $0$

释放侧实际势垒 $0$

鞍点能量 $0$

鞍点位置 $0$

能量够越过势垒吗？

对称势垒尺度

B

1\,\mathrm{J}

发射动能与势垒尺度之比

K_0/B

1.08

黏性阻尼

c

0\,\mathrm{kg\,s^{-1}}

势阱不对称度

f/f_{\mathrm{crit}}

0

半阱间距

a

1\,\mathrm{m}

发射方向

s_v

1

释放位置

x_0

-1\,\mathrm{m}

质点质量

m

1\,\mathrm{kg}

记录时长

t_{\mathrm{record}}

40\,\mathrm{s}

请求时间步长

h

5\,\mathrm{ms}

选中时刻

t

3.2\,\mathrm{s}

拖动质点设置释放位置，方向键也可操作。能量图和相图共用同一条已求解轨迹。

采用 SI 单位，向右为正。理想质点在给定的四次双势阱中运动，可加入线性黏性阻尼；导轨与球体是示意，不模拟滚动或执行机构。能量曲线是图表，不是真实轨道。灰色相曲线是释放能量下的保守参考；数值轨迹有阻尼时会失能。势垒设置是对称尺度，倾斜后的实际释放侧势垒另算。每步检查机械能加耗散热，半步长位置比较最多覆盖十秒。精确分界轨道需无限时间到达鞍点；能量误差范围内的数值穿越不能当作已分辨结果。

分界轨道分开单阱与跨阱运动

背景四次双势阱有两个稳定的极小值，中间隔着一个鞍点，也就是势垒顶。没有阻尼时，能量低于鞍点的粒子只能待在一侧；能量高于鞍点，它就能到达两个阱。能量恰好等于鞍点时，运动的粒子会无限逼近鞍点，永远走不完一个周期，这条轨道就是分界轨道。

为什么值得理解能量图和相空间图用的是单坐标经典力学的标准方法。势垒高度和允许的最大倾斜，都直接从图中的四次势能推出。加上黏性阻尼后，模型会记下阻尼移走的每一份动能。

核心问题: 只比势垒高一点，为什么穿越会慢那么多？
一句话直觉: 能量越接近势垒顶，粒子在鞍点附近逗留得越久。阻尼移走足够的能量后，原本跨阱的轨道会被困在一侧。

核心数学模型

倾斜双势阱

\begin{aligned}U(x)&=B\left[\left(\frac{x}{a}\right)^2-1\right]^2\\ &\quad+fx\end{aligned}

不倾斜时，势垒尺度就是势垒比任一极小值高出的能量。倾斜后，释放一侧的实际势垒按极小值和鞍点的真实位置重新计算。

两个极小值都保留的条件

\begin{aligned}|f|<f_{\mathrm{crit}}&=\frac{8B}{3\sqrt3\,a},\\ U\prime(x_*)&=0\end{aligned}

倾斜始终小于鞍结分岔的阈值，所以两个极小值和中间的鞍点一直分得开。

有黏性阻尼时的能量平衡

\begin{aligned}m\ddot x&=-U\prime(x)-c\dot x,\\ H(t)+Q(t)&=H(0),\\ Q(t)&=\int_0^t c\dot x^2\,dt\end{aligned}

每个阻尼子步让速度按指数衰减，并把损失的动能精确计入热量。保守力部分的积分误差仍会表现为能量对不上的差额。

超出势垒的能量与无阻尼周期

\begin{aligned}\Delta H&=H-U(x_s),\\ T(H)&=\sqrt{2m}\\ &\quad\times\int_{x_-}^{x_+}\frac{dx}{\sqrt{H-U(x)}}\end{aligned}

临界轨道没有有限的周期，所以图上留了一个缺口，不把势垒上下两支连起来。

鞍点附近的运动

\begin{aligned}U(x)&\simeq U(x_s)\\ &\quad-\frac{m\lambda_s^2}{2}(x-x_s)^2,\\ \lambda_s^2&=-\frac{U\prime\prime(x_s)}{m}\end{aligned}

在鞍点附近，周期按对数增长，微小的能量误差会被放大。粒子如果恰好在鞍点释放，又没有扰动，就会一直静止。

常见难点

临界能量下的穿越来自数值误差

容易误解临界能量下记录到的每次穿越都真实发生了。

正确理解精确的分界轨道在任何有限时间内都不会越过鞍点。有限步长带来的误差可能让它越过去，也可能让它折返，所以要对照能量误差和步长减半后的差别。

倾斜后，势垒尺度不等于实际势垒

容易误解对称势垒参数总等于释放侧实际势垒。

正确理解倾斜会改变两个极小值和鞍点的能量，读数用的是它们之间的实际差值。

阻尼改变轨道

容易误解损失能量以后，测到的周期仍然等于无阻尼参考周期。

正确理解能量一边减少，轨道一边改变。只有记录到的能量比鞍点低出数值不确定度以上，实验才认定粒子已被捕获。

动手实验

01
比较单阱运动与跨阱运动
先预测每次释放能否到达另一侧，再选择“单阱小振动”和“跨过双阱”，比较相图及数值检查中的完整记录穿越次数。
你会看到： 没有阻尼时，只有能量高于势垒的运动能到达两个阱。
02
逼近分界轨道
比较“将到鞍点又折返”和“缓慢越过鞍点”，再在对称、无阻尼且从阱底释放的条件下，把发射动能与势垒尺度之比设为一。解释穿越前先检查能量平衡误差。
你会看到： 两边的周期都变长。到了临界能量，参考周期不再是有限值，结果对数值误差也敏感得多。
03
让阻尼捕获粒子
选择“耗能后留在单阱”，先预测相轨道向哪一侧移动，再展开检查，比较机械能加耗散能量与初始能量。
你会看到： 热量加机械能一直等于初始能量，轨道最后被困在一侧。
04
倾斜但保留两个阱
选择“深浅不同的势阱”，在相图里拖动释放点，展开检查，把实际鞍点能量、释放侧势垒和对称尺度比一比。
你会看到： 鞍点移动了，从两侧看势垒高度不同，但两个极小值都还在。