重离子束探针(HIBP)· Proca–Green 平行板能量分析器 · 45分 · 每次求导、每步几何全部写全
几何部分请随时对照原卷图 2a、图 2b
磁约束聚变装置(托卡马克)芯部的等离子体温度高达几千万度,任何实体温度计伸进去都会瞬间气化,还会污染等离子体。但物理学家偏偏想知道芯部每一点的电势 \(\varphi\)——它直接关系到等离子体的约束品质、涨落和输运,是聚变研究的核心诊断量。
怎么办?思路是打一发"子弹"进去,让它在里面"变身",再飞出来接受检查。这就是重离子束探针(Heavy Ion Beam Probe, HIBP),1970–80 年代由苏联学者在 T-10 托卡马克上发明,此后装在美国的 TEXT、西班牙的 TJ-II、中国的 HL-2A 等装置上,至今是测芯部电势的主力手段。
① 加速器把重离子(题中取电荷 \(q_1=e\),动能 \(K_1\),常用 Tl⁺、Cs⁺ 等)射入等离子体;
② 在"测量位置"(那里电势为 \(\varphi\)),离子与电子碰撞再被电离一次,电荷从 \(+e\) 变成 \(+2e\)——这发"子弹"变身了;
③ 磁场对运动电荷不做功,但会使其偏转。二次离子电荷不同,偏转半径与一次离子不同,于是沿着另一条轨道自己飞出等离子体;
④ 在装置外接住它,精确测出它的动能 \(K_2\)。
电荷从 \(e\) 跳到 \(2e\) 是全部机关所在:穿过同一段电势差,二次离子"领到的"电势能是一次离子的两倍,一里一外,\(\varphi\) 就被换算成了动能差 \(K_2-K_1\)(第(2.3)问)。
最后一步"测 \(K_2\)"用的是 1970 年 Proca 和 Green 设计的平行板静电能量分析器(图 2b):两块平行板加电压 \(V_A\),板间匀强电场,离子从入射缝斜射进去走一段抛物线再出来,打到接收板上。不同能量的离子落点不同——本质上是台"弯出去的斜抛测能计"。
它的招牌设计就是本题要证的东西:入射角取 30° 时实现"二级聚焦"——入射角有小误差时,落点误差是三阶小量。这个 30° 在文献里就叫 Proca–Green 角。
全题 45 分,只用三样东西:高中斜抛运动、静电场能量守恒、泰勒展开。区分度全在"聚焦阶数"的概念和"成像"的思想上。
水平方向匀速、竖直方向匀加速(加速度 \(a\) 向下)。初速率 \(v\),仰角 \(\theta\)。取出发点为原点:
$$x=v\cos\theta\cdot t,\qquad y=v\sin\theta\cdot t-\frac12at^2$$"回到同一高度"即 \(y=0\):\(t\left(v\sin\theta-\frac12at\right)=0\),舍去 \(t=0\)(出发时刻),得飞行时间
$$T=\frac{2v\sin\theta}{a}$$代回 \(x\):射程 \(=v\cos\theta\cdot\dfrac{2v\sin\theta}{a}=\dfrac{v^2\cdot 2\sin\theta\cos\theta}{a}\)。用二倍角公式 \(2\sin\theta\cos\theta=\sin2\theta\):
$$X_{\text{射程}}=\frac{v^2\sin2\theta}{a}\quad\blacksquare$$对称性(本题反复用):\(y\) 是开口向下的抛物线,\(y=0\) 的两根 \(0\) 与 \(T\) 对称,故落地速度与起飞速度关于水平面对称——出射角仍是 \(\theta\),只是转为俯角,速率不变。所以离子穿过电场区后,方向"复原",只是整体平移了一段 \(X_{\text{射程}}\)。
① \((\tan x)'=\left(\dfrac{\sin x}{\cos x}\right)'=\dfrac{\cos x\cos x-\sin x(-\sin x)}{\cos^2 x}=\dfrac{1}{\cos^2 x}=\sec^2 x\)
② \((\cot x)'=\left(\dfrac{\cos x}{\sin x}\right)'=\dfrac{-\sin^2 x-\cos^2 x}{\sin^2 x}=-\csc^2 x\)
③ \((\sec^2 x)'\):链式法则,\(\sec^2 x=(\cos x)^{-2}\),导数 \(=-2(\cos x)^{-3}(-\sin x)=2\dfrac{\sin x}{\cos^3 x}=2\sec^2 x\tan x\)
④ \((\csc^2 x)'\):同理,\((\sin x)^{-2}\) 的导数 \(=-2(\sin x)^{-3}\cos x=-2\csc^2 x\cot x\)
再备两条恒等式:和角公式 \(\sin(A+B)=\sin A\cos B+\cos A\sin B\);二倍角 \(\sin 2\theta=2\sin\theta\cos\theta\)。全题三角运算就这些。
光滑函数在某点 \(x_0\) 附近总可以写成
$$f(x_0+\Delta x)=f(x_0)+f'(x_0)\Delta x+\frac{f''(x_0)}{2}(\Delta x)^2+\frac{f'''(x_0)}{6}(\Delta x)^3+\cdots$$现在 \(f\) 是"落点位置",\(\Delta x\) 是某个误差(入射角误差、能量误差)。误差对落点的影响被哪一阶主导,就叫几级聚焦:
举例感受三阶的威力:角度误差 \(\Delta\theta=0.01\) rad(约 0.6°),一阶影响是 \(0.01\),三阶只有 \(10^{-6}\)——差了四个数量级。这就是为什么仪器要追高级聚焦。
判断"扔哪项"的标准还是看量级,本文每个展开都会写明代数依据。
第(2.1)问里落点同时依赖角度误差 \(\delta\theta\) 和能量误差 \(\delta G\),两个都是小量。展开式是一维的直接推广——每个自变量各阶都列出来,按"总阶数"从低到高排:
$$N(\theta_0+\delta\theta,\,G_0+\delta G)=N_0+\underbrace{N_\theta\,\delta\theta+N_G\,\delta G}_{\text{一阶}} +\underbrace{\tfrac12N_{\theta\theta}(\delta\theta)^2+N_{\theta G}\,\delta\theta\,\delta G+\tfrac12N_{GG}(\delta G)^2}_{\text{二阶}} +\underbrace{\tfrac16N_{\theta\theta\theta}(\delta\theta)^3+\cdots}_{\text{三阶}}$$其中 \(N_\theta\) 表示对 \(\theta\) 的偏导数(求导时把 \(G\) 当常数),其余同理。别怕符号,其实就是"一次只对一个字母用工具3"。
一个光学(或粒子光学)系统把"物"上的点映到"像"上的点,若关系是
$$s=s_0+m\,l$$(\(l\):物上坐标,\(s\):像上坐标,\(m\) 常数),就说它成像:\(m\) 叫放大率,\(m<0\) 叫倒立(物上移、像下移),\(|m|=1\) 叫等大。关键判据是:\(m\) 必须不随"从哪个角度进"而变,否则同一点发出的不同角度的粒子会落到不同地方,像就糊了。二级聚焦恰好保证了"角度差异的影响是三阶小量",所以第(1.3)问才能谈成像——这两问是因果关系。
电荷 \(q\) 从电势 \(\varphi_A\) 处移到 \(\varphi_B\) 处,静电力做功 \(W=q(\varphi_A-\varphi_B)\),转化为动能增量:
$$K_B-K_A=q(\varphi_A-\varphi_B)$$取等离子体外 \(\varphi=0\) 为参考。正电荷往高电势处走("爬坡")动能减少 \(q\varphi\),往低电势处走动能增加。就这一句,第(2.3)问全用完。磁场始终不做功(洛伦兹力垂直速度),所以经过磁场不改变动能数值。
取下极板为 \(y=0\),上极板 \(y=d\),板间匀强电场 \(E=V_A/d\) 竖直向下,二次离子(\(q_2=2e\),质量 \(m\))在板间受恒定向下加速度
$$a=\frac{q_2E}{m}=\frac{q_2V_A}{md}=\frac{2eV_A}{md}$$对照图 2b,离子的一生分三段:
式子里 \(\dfrac{v^2}{a}\) 会反复出现,给它起个名字——特征长度 \(G\)(量纲:\(\frac{(\text{m/s})^2}{\text{m/s}^2}=\text{m}\) ✓):
$$G\equiv\frac{v^2}{a}=\frac{\frac{2K_2}{m}}{\frac{q_2V_A}{md}}=\frac{2K_2d}{q_2V_A}=\frac{K_2d}{eV_A}$$最后一步用了 \(q_2=2e\)。\(G\) 的物理含义:速度全转成竖直向上时恰好能爬的高度(的两倍),它是"能量有多难被电场掰弯"的度量。第(2)部分会看到,测 \(K_2\) 等价于测 \(G\)。
此后全程用缩写 \(H\equiv Y_{D1}+Y_{D2}\)(两段漂移的总竖直落差)。
对 \(X_D=H\cot\theta+G\sin2\theta\) 求导(工具2②,外加 \((\sin2\theta)'=2\cos2\theta\)):
$$\frac{\partial X_D}{\partial\theta}=-H\csc^2\theta+2G\cos2\theta=0 \quad\Longrightarrow\quad G=\frac{H}{2\sin^2\theta\cos2\theta}\tag{1}$$对一阶导数再导一次(工具2④用于 \((\csc^2\theta)'\),\((\cos2\theta)'=-2\sin2\theta\)):
$$\frac{\partial^2 X_D}{\partial\theta^2}=2H\csc^2\theta\cot\theta-4G\sin2\theta=0\tag{2}$$把(1)代入(2):
$$2H\csc^2\theta\cot\theta=4\cdot\frac{H}{2\sin^2\theta\cos2\theta}\cdot\sin2\theta$$两边约去 \(2H\),注意 \(\csc^2\theta=\dfrac{1}{\sin^2\theta}\):
$$\frac{\cot\theta}{\sin^2\theta}=\frac{1}{\sin^2\theta}\cdot\frac{\sin2\theta}{\cos2\theta} \quad\Longrightarrow\quad \cot\theta=\tan2\theta$$用 \(\cot\theta=\tan(90^\circ-\theta)\):
$$\tan(90^\circ-\theta)=\tan2\theta\quad\Longrightarrow\quad 90^\circ-\theta=2\theta\quad\Longrightarrow\quad \theta=30^\circ$$物理解围:\(\theta_I\) 实际只取 \(0\sim45^\circ\)(再大抛物线太平、仪器过长),\(\tan\) 在此区间单调,解唯一。这个 30° 就是文献里的 Proca–Green 角。
\(\theta=30^\circ\) 时 \(\sin^2 30^\circ=\frac14\),\(\cos60^\circ=\frac12\),代入(1):
$$G=\frac{H}{2\cdot\frac14\cdot\frac12}=4H$$由 \(G=\dfrac{2K_2d}{q_2V_A}\) 反解电压(\(q_2=2e\)):
光令前两阶导数为零还不够,得看三阶项系数。由(2)左边继续求导(乘积法则拆开,工具2③④):
$$\frac{\partial^3 X_D}{\partial\theta^3} =2H\big[(\csc^2\theta)'\cot\theta+\csc^2\theta(\cot\theta)'\big]-8G\cos2\theta =2H\big[-2\csc^2\theta\cot^2\theta-\csc^4\theta\big]-8G\cos2\theta$$代入 \(\theta=30^\circ,\ G=4H\)(\(\csc^2 30^\circ=4,\ \cot^2 30^\circ=3,\ \csc^4 30^\circ=16,\ \cos60^\circ=\frac12\)):
$$=2H\big[-2\cdot4\cdot3-16\big]-8\cdot4H\cdot\frac12=2H\cdot(-40)-16H=-96H$$所以
$$\Delta X_D=\frac{1}{6}\cdot(-96H)(\Delta\theta)^3+O(\Delta\theta^4)=-16H(\Delta\theta)^3+O(\Delta\theta^4)$$对照图 2b,按官方解答的读法:入射板为"\"形(上端在左、下端在右),\(+l\) 沿板指向右下;接收板为"/"形(上端在右、下端在左),\(+s\) 沿板指向左下。板与水平夹角分别为 \(\theta_{b1},\theta_{b2}\)。于是当入射点从缝中心沿板移动 \(l\):
离子的速度方向不变(仍为 \(\theta_I=30^\circ\))——移动的是出发点,不是角度。接收端同理:接收点偏离缝中心 \(s\) 时坐标为 \(\big(X_D-s\cos\theta_{b2},\;-Y_{D2}-s\sin\theta_{b2}\big)\)。
段一:入射点比原来低了 \(l\sin\theta_{b1}\),要到 \(y=0\) 需上升 \(Y_{D1}+l\sin\theta_{b1}\),故到达下极板时的水平坐标
$$x_1=l\cos\theta_{b1}+\big(Y_{D1}+l\sin\theta_{b1}\big)\cot\theta_I$$段二:进入电场的角度、速率都没变,斜抛射程不变(\(G\) 只依赖速率和加速度):
$$x_2=x_1+G\sin2\theta_I$$段三:从 \(y=0\) 以俯角 \(\theta_I\) 下降 \(Y_{D2}+s\sin\theta_{b2}\) 到接收点,水平走 \(\big(Y_{D2}+s\sin\theta_{b2}\big)\cot\theta_I\)。到达即:
$$x_2+\big(Y_{D2}+s\sin\theta_{b2}\big)\cot\theta_I=X_D-s\cos\theta_{b2}$$把 \(x_1,x_2\) 代入上式左边:
$$l\cos\theta_{b1}+\big(Y_{D1}+l\sin\theta_{b1}\big)\cot\theta_I+G\sin2\theta_I+\big(Y_{D2}+s\sin\theta_{b2}\big)\cot\theta_I=X_D-s\cos\theta_{b2}$$用 (1.1) 的结果 \(X_D=H\cot\theta_I+G\sin2\theta_I\) 消去两边的 \(Y_{D1}\cot\theta_I+Y_{D2}\cot\theta_I+G\sin2\theta_I\),剩下:
$$l\cos\theta_{b1}+l\sin\theta_{b1}\cot\theta_I+s\sin\theta_{b2}\cot\theta_I+s\cos\theta_{b2}=0$$含 \(l\) 的归一边、含 \(s\) 的归一边:
$$l\big(\cos\theta_{b1}+\sin\theta_{b1}\cot\theta_I\big)=-s\big(\cos\theta_{b2}+\sin\theta_{b2}\cot\theta_I\big)$$括号里通分(乘 \(\sin\theta_I\) 后用和角公式,工具2末):\(\cos\theta_b\sin\theta_I+\sin\theta_b\cos\theta_I=\sin(\theta_I+\theta_b)\),即
$$\cos\theta_b+\sin\theta_b\cot\theta_I=\frac{\sin(\theta_I+\theta_b)}{\sin\theta_I}$$代入(两边 \(\sin\theta_I\) 正好约掉):
检验①:\(\theta_{b1}=\theta_{b2}=0\)(两板水平),则 \(s=-l\)——板水平时"沿板移动"就是纯水平平移,整条轨迹平移 \(l\),落点当然平移 \(-l\)(方向相反是坐标定义使然)。✓
检验②:\(\theta_{b1}=\theta_{b2}\)(两板平行放置),放大率 \(=1\)——几何上平行的两板对同一条直线束"截距"变化相同。✓
检验③(物理意义):负号 = 倒立像。物上移、像下移,正是聚焦成像系统的标志(对照工具5)。而且 \(m\) 只含角度常数、与离子从哪个方向来无关——这只有在二级聚焦压住角度误差后才成立。(1.2) 是 (1.3) 的前提。
段一:\(x_1=Y_{D1}\cot\theta_I\);段二:\(x_2=x_1+G\sin2\theta_I\);段三到达接收板(坐标 \(X_D-s_0\cos\theta_{b2}\),深度 \(Y_{D2}+s_0\sin\theta_{b2}\)):
$$Y_{D1}\cot\theta_I+G\sin2\theta_I+\big(Y_{D2}+s_0\sin\theta_{b2}\big)\cot\theta_I=X_D-s_0\cos\theta_{b2}$$含 \(s_0\) 的项归左边,其余归右边,并代入 \(X_D=3\sqrt3H\) 与 \(Y_{D1}+Y_{D2}=H\):
$$s_0\big(\cos\theta_{b2}+\sin\theta_{b2}\cot\theta_I\big)=3\sqrt3H-H\cot\theta_I-G\sin2\theta_I$$右边要化成"以 \(\tan\theta_I\) 为单位"的形式,两边同乘 \(\tan\theta_I\)(\(\cot\theta_I\tan\theta_I=1\)):
其中用了 \(\sin2\theta\tan\theta=2\sin\theta\cos\theta\cdot\dfrac{\sin\theta}{\cos\theta}=2\sin^2\theta\)。另外分母可写成 \(\dfrac{\sin(\theta_I+\theta_{b2})}{\cos\theta_I}\)(通分+和角公式),后文用哪个形式方便就用哪个。
记分子 \(N(\theta,G)=3\sqrt3H\tan\theta-H-2G\sin^2\theta\),在设计点 \((30^\circ,4H)\) 展开(工具4)。逐个算偏导:
继续算需要的非零偏导:
于是分子的领头项(一阶的 \(\delta\theta\) 项、二阶的 \((\delta\theta)^2\) 项都为零):
$$N\approx-\frac12\,\delta G-\sqrt3\,\delta\theta\,\delta G+\frac{1}{6}\cdot32\sqrt3H\,(\delta\theta)^3 =-\frac{\delta G}{2}-\sqrt3\,\delta\theta\,\delta G+\frac{16\sqrt3H}{3}(\delta\theta)^3$$分母取设计点的值即可(它的修正只会乘在已是小量的分子上,升阶,扔):
$$\cos\theta_{b2}\tan30^\circ+\sin\theta_{b2}=\frac{1}{\sqrt3}\cos\theta_{b2}+\sin\theta_{b2} =\frac{2}{\sqrt3}\Big(\frac12\cos\theta_{b2}+\frac{\sqrt3}{2}\sin\theta_{b2}\Big)=\frac{2}{\sqrt3}\sin(30^\circ+\theta_{b2})$$入射缝宽 \(w\),束流均匀地照亮整条缝(\(l\in[-\frac w2,\frac w2]\) 上粒子数均匀)。接收板沿 \(s=0\) 分成上半、下半,分别收集电流 \(i_U,i_L\)。要由这两个读数反推 \(K_2\)。
由 (1.3) 的线性成像加上 (2.1) 的中心偏移,任意 \(l\) 处的落点为
$$s=s_0-\frac{\sin(30^\circ+\theta_{b1})}{\sin(30^\circ+\theta_{b2})}\,l$$斜率为负(倒像)。于是存在一个临界物点 \(l_0\):从 \(l_0\) 进的粒子恰好落在分界 \(s=0\) 上。\(l>l_0\) 的粒子成倒像落到 \(s<0\)(上半侧)被上半板收集;\(l
两式相减、相加再取比(\(I_{\text{总}}\) 消掉——用比值测,免去绝对流强标定,这是实验上极关键的优点):
$$\frac{i_L-i_U}{i_L+i_U}=\frac{2l_0}{w}\quad\Longrightarrow\quad l_0=\frac w2\cdot\frac{i_L-i_U}{i_L+i_U}\tag{3}$$在 5.1 的追踪里把 \(l\) 加回来(段一变成 \(x_1=l\cos\theta_{b1}+(Y_{D1}+l\sin\theta_{b1})\cot\theta_I\)),令落点 \(s=0\):
$$l_0\big(\cos\theta_{b1}+\sin\theta_{b1}\cot\theta_I\big)+H\cot\theta_I+G\sin2\theta_I=3\sqrt3H$$同乘 \(\tan\theta_I\) 整理(与 5.1 完全同构):
$$l_0=\frac{3\sqrt3H\tan\theta_I-H-2G\sin^2\theta_I}{\cos\theta_{b1}\tan\theta_I+\sin\theta_{b1}}$$安装误差的影响是三阶小量(5.2 已证),取 \(\theta_I=30^\circ\):分子 \(=3\sqrt3H\cdot\frac{1}{\sqrt3}-H-2G\cdot\frac14=2H-\dfrac G2\);分母同 5.2 的算法 \(=\dfrac{2}{\sqrt3}\sin(30^\circ+\theta_{b1})\)。所以
$$l_0=\frac{\sqrt3\left(2H-\frac G2\right)}{2\sin(30^\circ+\theta_{b1})}\tag{4}$$两边乘 \(2\sin(30^\circ+\theta_{b1})/\sqrt3\):
$$2H-\frac G2=\frac{w\sin(30^\circ+\theta_{b1})}{\sqrt3}\cdot\frac{i_L-i_U}{i_L+i_U}$$解 \(G\)(两边乘 \(-2\) 再移项,注意 \(\frac{i_L-i_U}{i_L+i_U}=-\frac{i_U-i_L}{i_U+i_L}\)):
$$G=4H+\frac{2w}{\sqrt3}\sin(30^\circ+\theta_{b1})\,\frac{i_U-i_L}{i_U+i_L}$$最后由 \(G=\dfrac{K_2d}{eV_A}\) 即 \(K_2=\dfrac{eV_A}{d}G\):
检验①:\(i_U=i_L\) 时 \(K_2=\dfrac{4eV_AH}{d}\),正是 (1.2) 的设计关系 \(G=4H\)。✓ 检验②:\(K_2\) 偏大时由 5.2 末的符号分析,中心束偏向上半侧,\(i_U>i_L\),公式给出的修正项为正——\(K_2\) 大于设计值,自洽。✓
能量链条三步走(工具6;按题意忽略碰撞衰减、电离能、电子带走的小量):
"二级聚焦条件" = 落点对误差的泰勒展开中一阶、二阶项系数都为零,领头误差是三阶。
"特征长度 \(G\)" = 把 \(v^2/a\) 打包起绰号;本题里 \(G\) 就是能量的化身,测位置 = 测 \(G\) = 测 \(K_2\)。
"上下半平面收集的电流 \(i_U,i_L\)" = 把位置测量换成电流比值测量,消掉绝对流强,实验上才好做。
三板斧复盘:斜抛射程 + 无场漂移分段相加、聚焦条件 = 导数方程组、成像 = 线性映射 + 电流比反演。
① 三段运动相加得 \(X_D=H\cot\theta_I+G\sin2\theta_I\);② 令 \(X_D'(\theta)=X_D''(\theta)=0\) 解出 \(30^\circ\) 与 \(G=4H\),并验证三阶项 \(-16H(\Delta\theta)^3\);③ 偏移 \(l\) 的轨迹追踪得线性成像 \(s=-m l\),\(m=\sin(30^\circ+\theta_{b1})/\sin(30^\circ+\theta_{b2})\);④ 非设计态中心束落点 \(s_0\) 的严格式 + 展开(分子前两阶为零是聚焦条件的化身);⑤ 电流比给出临界物点 \(l_0=\frac w2\frac{i_L-i_U}{i_L+i_U}\),与几何式联立解出 \(K_2\);⑥ 能量链条 \(K_2=K_1+e\varphi\)。
能讲清"为什么二级聚焦是成像的前提""为什么用电流比而不是电流绝对值"这两个为什么,这题你就吃透了——45 分的骨架全在这两个观念上。