在上一篇文章中,我们已经得到了电偶极子的等势面和电场线方程,这应该可以让我们对电偶极子的力场情况有个大致的了解了。当然,我们还是希望能够求出在这样的一个受力情况下,一个带电粒子是如何运动的。简单起见,在下面的探讨中,我们假定带电粒子的质量和电荷量均为1,至于电荷的正负,可以通过改变在$U=-\frac{k \cos\theta}{r^2}$中的k值的正负来控制。我们使用的工具依旧是理论力学中的欧拉-拉格朗日方程。

也许不少读者始终对公式感到头疼,更不用说是博大精深的理论力学了。但是请相信我,如果你花一点点心思去弄懂用变分法研究力学(或其他物理系统,但我目前只会用于力学)的基本思路和步骤,那么对你的物理研究是大有裨益的。因为在我眼中,学习了一丁点的理论力学知识后,我看到的只有物理的简洁与和谐。有兴趣的朋友可以看看我的那几篇《自然极值》等相关文章。

首先写出动能的表达式:$T=\frac{1}{2} (\dot{r}^2+r^2 \dot{\theta}^2)$

还有势能:$U=-\frac{k \cos\theta}{r^2}$

拉格朗日量:$L=T-U=\frac{1}{2}(\dot{r}^2+r^2 \dot{\theta}^2)+\frac{k \cos\theta}{r^2}$

求关于r的微分方程:
$$\begin{aligned}\frac{\partial L}{\partial \dot{r}}=\dot{r} \\ \frac{\partial L}{\partial r}=r\dot{\theta}^2-\frac{2kcos\theta}{r^3}\end{aligned}$$

我们得到:
$$\ddot{r}=r\dot{\theta}^2-\frac{2kcos\theta}{r^3}\tag{10}$$
继续对θ求偏导数我们可以求出关于θ的微分方程,但是我们不这样做,因为保守力场中必然存在着能量守恒:
$$T+U=h=\frac{1}{2}(\dot{r}^2+r^2 \dot{\theta}^2)-\frac{k \cos\theta}{r^2}\tag{11}$$
用(10)和(11)联立就可以得出所有的解。原题本质上是关于两个函数的二阶微分方程,所以应该包含4个积分常数,而(11)式已经包含了一个积分常数,所以还剩下三个。

用心摆弄一会儿,我们会发现一件很巧的事情:
$$h=\frac{1}{2}\dot{r}^2+\frac{1}{2}(r^2 \dot{\theta}^2-\frac{2k \cos\theta}{r^2})=\frac{1}{2}(\dot{r}^2+r\ddot{r})\tag{12}$$
即$$h=\frac{1}{2}(\dot{r}^2+r\ddot{r})\tag{12}$$
太巧了,而且这个结果也太简洁了!乘胜追击,得到
$$\begin{aligned}r\dot{r}=2ht+C_1 \\ \frac{1}{2}r^2=ht^2+C_1 t+C_2\end{aligned}\tag{13}$$
至此,我们已经得到了三个积分常数。让我们用点心思分析一下这个解。从(13)式可以知道,$C_2 >0$(一开始不能就在原点)。如果 h > 0 ,那么带电粒子最终将会飞到无穷远处,同时这必然要求了$\Delta=C_1^2-4hC_2 \leq 0$;如果 h < 0 ,该粒子最终会回到原点,撞到电偶极子上。这也符合我们的知识,具有正能量的系统可以脱离束缚,具有负能量的系统无法脱离束缚。当h为0时,又需要进一步讨论。

剩下最后一个积分常数,我们已经求出了r关于时间t的演化,能量守恒也给出了速度的演化,最后一个是什么呢?答案有很多,但我们有兴趣的是——轨道的形状是什么?以下是简单的过程

$$\begin{aligned}\dot{r}^2+r^2 \dot{\theta}^2=2h+\frac{2kcos\theta}{r^2} \\ 1+r^2(\frac{d\theta}{dr})^2=\frac{2h+\frac{2kcos\theta}{r^2}}{\dot{r}^2}\end{aligned}\tag{14}$$
下面想办法将(14)式中的$\dot{r}^2$消去,根据(13),我们有
$$\begin{aligned}r^2 \dot{r}^2=4h^2 t^2+C_1^2+4hC_1 t \\ 2hr^2=4h^2 t^2 +4hC_1 t+4hC_2 \\ r^2 \dot{r}^2-C_1^2=2hr^2-4hC_2 \\ \dot{r}^2=2h+\frac{C_1^2-4hC_2}{r^2}\end{aligned}$$

代入

$$\begin{aligned}1+r^2(\frac{d\theta}{dr})^2=\frac{2hr^2+2kcos\theta}{2hr^2+C_1^2-4hC_2} \\ r^2(\frac{d\theta}{dr})^2=\frac{2kcos\theta-C_1^2+4hC_2}{2hr^2+C_1^2-4hC_2} \\ (\frac{d\theta}{dr})^2=\frac{2kcos\theta-C_1^2+4hC_2}{(2hr^2+C_1^2-4hC_2)r^2} \\ \frac{dr}{r\sqrt{2hr^2+C_1^2-4hC_2}}=\frac{d\theta}{\sqrt{2kcos\theta-C_1^2+4hC_2}}\end{aligned}\tag{15}$$
最后一个式子已经包含了轨道的形状在内了,两边同时积分即可。很遗憾,这不是一个简单的显式解,它只能够用椭圆积分来表达。一个容易想到的法子是有没有一些什么简单的特例,让我们可以快速看看曲线的形状?好吧,我们令$k=0.5,-C_1^2+4hC_2=1,h=0.5$,那么(15)就变成
$$\frac{dr}{r\sqrt{r^2-1}}=\frac{d\theta}{\sqrt{\cos\theta+1}}$$

得到$arccos(\frac{1}{r})=\frac{1}{\sqrt{2}}ln(tan(\frac{\pi}{4}+\frac{\theta}{4}))+Const$

也就是说$r=\frac{\sqrt{2}}{cos(ln(tan(\frac{\pi}{4}+\frac{\theta}{4})))}$

这些曲线的形状让我也非常头疼...由于对相关知识研究还不够深入,不再多说...权当设置一个悬念,希望有兴趣的读者继续完成.

转载到请包括本文地址:https://spaces.ac.cn/archives/1701

更详细的转载事宜请参考:《科学空间FAQ》

如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。

如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!

如果您需要引用本文,请参考:

苏剑林. (Aug. 18, 2012). 《电偶极子浅探(2) 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/1701

@online{kexuefm-1701,
        title={电偶极子浅探(2)},
        author={苏剑林},
        year={2012},
        month={Aug},
        url={\url{https://spaces.ac.cn/archives/1701}},
}