UA PHYS515A 电磁理论V 电磁波与辐射7 运动点电荷的辐射

实际问题中辐射的source都是比较复杂的charge density与currency density,但作为比较简单直观易于理解的模型,我们可以先学习point charge与point currency。用ro(t~)\textbf r_o(\tilde t)ro​(t~)表示point charge的位移,这里的ttt加了上标是为了强调这是retarded time,也就是source在过去某个时刻产生电磁波然后在未来某个时刻被我们观察到。用x\textbf xx表示观察者的位移,当观察者观察到电磁波时,point charge的位移为ro(t)\textbf r_o(t)ro​(t),用这个位移作为参考位移。

point charge的density为
ρ(r′,t′)=qδ(r′−ro(t′))J(r′,t′)=qvδ(r′−ro(t′))\rho(\textbf r',t')=q\delta(\textbf r'-\textbf r_o(t')) \\ \textbf J(\textbf r',t')=q\textbf v \delta(\textbf r'-\textbf r_o(t'))ρ(r′,t′)=qδ(r′−ro​(t′))J(r′,t′)=qvδ(r′−ro​(t′))

假设∣v∣<<c|\textbf v|<<c∣v∣<<c,否则需要引入狭义相对论来完成推导,标量势与向量势满足
(∇2−1c2∂2∂t2)A=−4πcJ(∇2−1c2∂2∂t2)Φ=−4πρ(\nabla^2-\frac{1}{c^2}\frac{\partial^2}{\partial t^2})\textbf A = -\frac{4 \pi}{c} \textbf J \\ (\nabla^2-\frac{1}{c^2}\frac{\partial^2}{\partial t^2}) \Phi = -4 \pi \rho(∇2−c21​∂t2∂2​)A=−c4π​J(∇2−c21​∂t2∂2​)Φ=−4πρ

用含时Green函数法写出potential的积分解
Φ(x,t)=∬ρ(r′,t′)δ(t′−t+∣x−r′∣c)∣x−r′∣dr′dt′=∫ρ(r′,t−∣x−r′∣c)∣x−r′∣dr′=q∫δ(r′−ro(t−∣x−r′∣c))∣x−r′∣dr′A(x,t)=1c∫J(r′,t−∣x−r′∣c)∣x−r′∣dr′=q∫β⃗(t−∣x−r′∣c)δ(r′−ro(t−∣x−r′∣c))∣x−r′∣dr′\Phi(\textbf x,t)= \iint \rho(\textbf r',t') \frac{\delta(t'-t+\frac{|\textbf x-\textbf r'|}{c})}{|\textbf x - \textbf r'|}d \textbf r' dt' \\= \int \frac{\rho(\textbf r',t-\frac{|\textbf x-\textbf r'|}{c})}{|\textbf x - \textbf r'|}d \textbf r' = q \int \frac{\delta(\textbf r'-\textbf r_o(t-\frac{|\textbf x - \textbf r'|}{c}))}{|\textbf x - \textbf r'|}d\textbf r' \\ \textbf A(\textbf x ,t)=\frac{1}{c}\int \frac{\textbf J (\textbf r',t-\frac{|\textbf x-\textbf r'|}{c})}{|\textbf x - \textbf r'|}d \textbf r' \\ =q \int \frac{\vec \beta(t-\frac{|\textbf x-\textbf r'|}{c})\delta(\textbf r'-\textbf r_o(t-\frac{|\textbf x - \textbf r'|}{c}))}{|\textbf x - \textbf r'|}d\textbf r' Φ(x,t)=∬ρ(r′,t′)∣x−r′∣δ(t′−t+c∣x−r′∣​)​dr′dt′=∫∣x−r′∣ρ(r′,t−c∣x−r′∣​)​dr′=q∫∣x−r′∣δ(r′−ro​(t−c∣x−r′∣​))​dr′A(x,t)=c1​∫∣x−r′∣J(r′,t−c∣x−r′∣​)​dr′=q∫∣x−r′∣β​(t−c∣x−r′∣​)δ(r′−ro​(t−c∣x−r′∣​))​dr′

接下来就比较难搞了,因为Dirac函数的自变量是函数,所以需要引入一些关于Dirac函数的新的技巧来完成后续积分。

记R=x−r′\textbf R=\textbf x - \textbf r'R=x−r′,n^\hat nn^表示R\textbf RR的方向,考虑Φ\PhiΦ的表达式:
q∫δ(r′−ro(t−Rc))Rdr′q \int \frac{\delta(\textbf r'-\textbf r_o(t-\frac{R}{c}))}{R}d\textbf r' q∫Rδ(r′−ro​(t−cR​))​dr′

定义
r∗=r′−ro(t−Rc)dr∗=[1−n^(t~)⋅β⃗(t~)]d3r′\textbf r^* = \textbf r'-\textbf r_o(t-\frac{R}{c}) \\ d \textbf r^* = [1-\hat n(\tilde t) \cdot \vec \beta(\tilde t)] d^3 \textbf r'r∗=r′−ro​(t−cR​)dr∗=[1−n^(t~)⋅β​(t~)]d3r′

变换积分变量,
Φ=q∫δ(r∗)d3r∗∣x−r∗−r0(t~)∣(1−n^⋅β⃗)=q(1−n^⋅β⃗)R∣t~\Phi = q \int \frac{\delta(\textbf r^*)d^3 \textbf r^*}{|\textbf x-\textbf r^*-\textbf r_0(\tilde t)|(1-\hat n \cdot \vec \beta)} = \frac{q}{(1-\hat n \cdot \vec \beta)R}|_{\tilde t}Φ=q∫∣x−r∗−r0​(t~)∣(1−n^⋅β​)δ(r∗)d3r∗​=(1−n^⋅β​)Rq​∣t~​

同样的方法可以得到
A=qβ⃗(1−n^⋅β⃗)R∣t~\textbf A = \frac{q \vec \beta}{(1-\hat n \cdot \vec \beta)R }|_{\tilde t}A=(1−n^⋅β​)Rqβ​​∣t~​

称这两个解为Lienard-Wiechert potentials。知道了势之后可以写出电磁场
E=−∇Φ−1c∂∂tAB=∇×A\textbf E = -\nabla \Phi-\frac{1}{c}\frac{\partial}{\partial t} \textbf A \\ \textbf B = \nabla \times \textbf AE=−∇Φ−c1​∂t∂​AB=∇×A

一种比较常用的计算方法是将势重新写为关于时间的积分形式
Φ(x,t)=q∫δ(t′−t+R(t′)c)R(t′)dt′A(x,t)=q∫β⃗(t′)δ(t′−t+R(t′)c)R(t′)dt′\Phi(\textbf x,t) = q\int \frac{\delta(t'-t+\frac{R(t')}{c})}{R(t')}dt' \\ \textbf A(\textbf x,t) = q\int \frac{\vec \beta(t')\delta(t'-t+\frac{R(t')}{c})}{R(t')}dt'Φ(x,t)=q∫R(t′)δ(t′−t+cR(t′)​)​dt′A(x,t)=q∫R(t′)β​(t′)δ(t′−t+cR(t′)​)​dt′

先引入一个关于Dirac函数的性质,其实就是对前面换元法的抽象:
∫δ(f(x))dx=∫δ(f(x))dfdf/dx=1dfdx∣x=0\int \delta(f(x))dx = \int \delta(f(x)) \frac{df}{df/dx} = \frac{1}{\frac{df}{dx}|_{x=0}}∫δ(f(x))dx=∫δ(f(x))df/dxdf​=dxdf​∣x=0​1​

如果
f(t′)=t′−t+R(t′)cf′(t′)≈1−2(x−r′)⋅r′˙2Rc=1−n^⋅β⃗f(t')=t'-t+\frac{R(t')}{c} \\ f'(t')\approx 1-\frac{2(\textbf x- \textbf r')\cdot \dot{\textbf r'}}{2Rc} = 1-\hat n \cdot \vec \betaf(t′)=t′−t+cR(t′)​f′(t′)≈1−2Rc2(x−r′)⋅r′˙​=1−n^⋅β​

这是对前文变量替换时多出来的contraction effect的解释。在这个积分形式下计算导数并带入电磁场的表达式(计算过程有点长,就省略了,有需要的话可以翻一下Jackson的书)
E=−q∫∇(δ(t′−t+Rc)R)dt′−qcddt∫β⃗δ(t′−t+Rc)Rdt′=qn^(1−n^⋅β⃗)R2∣t~+qcddtn^−β⃗(1−n^⋅β⃗)R∣t~=q((n^−β⃗)(1−β⃗2)(1−n^⋅β⃗)R2+n^×[n^−β⃗]×β⃗˙c(1−n^⋅β⃗)3R)t~\textbf E = -q \int \nabla \left( \frac{\delta(t'-t+\frac{R}{c})}{R} \right)dt' - \frac{q}{c} \frac{d}{dt} \int \frac{\vec \beta \delta(t'-t+\frac{R}{c})}{R}dt' \\ = q \frac{\hat n}{(1-\hat n \cdot \vec \beta)R^2}|_{\tilde t}+\frac{q}{c} \frac{d}{dt} \frac{\hat n - \vec \beta}{(1-\hat n \cdot \vec \beta)R}|_{\tilde t} \\ = q \left( \frac{(\hat n-\vec \beta)(1-\vec \beta^2)}{(1-\hat n \cdot \vec \beta)R^2}+\frac{\hat n \times [\hat n - \vec \beta] \times \dot{\vec \beta}}{c(1-\hat n \cdot \vec \beta)^3R}\right)_{\tilde t}E=−q∫∇(Rδ(t′−t+cR​)​)dt′−cq​dtd​∫Rβ​δ(t′−t+cR​)​dt′=q(1−n^⋅β​)R2n^​∣t~​+cq​dtd​(1−n^⋅β​)Rn^−β​​∣t~​=q((1−n^⋅β​)R2(n^−β​)(1−β​2)​+c(1−n^⋅β​)3Rn^×[n^−β​]×β​˙​​)t~​

第一项被称为velocity term或者Coulomb term,它只与速度有关且形式与库仑field的形式一致,并且当β⃗→0\vec \beta \to 0β​→0,也就是点电荷运动速度相比光速可以忽略时,这一项退化为库仑field;如果点电荷的运动没有加速度,那么第二项为0;如果点电荷的运动存在加速度,第二项不为0,称其为radiation field

UA PHYS515A 电磁理论V 电磁波与辐射7 运动点电荷的辐射相关推荐

  1. UA PHYS515A 电磁理论V 电磁波与辐射11 简单辐射问题 电偶极子的辐射

    UA PHYS515A 电磁理论V 电磁波与辐射11 简单辐射问题 电偶极子的辐射 一对带等量相反电量的点电荷构成一对电偶极子,假设电量为qqq,两个点电荷的距离为aaa,dipole moment为 ...

  2. UA PHYS515A 电磁理论V 电磁波与辐射10 简单辐射问题 一根通电电线的辐射

    UA PHYS515A 电磁理论V 电磁波与辐射10 简单辐射问题 一根通电电线的辐射 假设zzz轴上放了一根电线,观察者位移的单位向量为r^\hat rr^,用I\textbf II表示电线中的电流 ...

  3. UA PHYS515A 电磁理论V 电磁波与辐射9 简单辐射系统

    UA PHYS515A 电磁理论V 电磁波与辐射9 简单辐射系统 前文讨论了单个直线运动的带电粒子的辐射,但在实践中只有在实验室才能观察到这种现象,应用中遇到的情况会更加复杂,比如多个粒子做直线运动, ...

  4. UA PHYS515A 电磁理论V 电磁波与辐射8 单个粒子的辐射 匀速运动与匀加速运动的情况

    UA PHYS515A 电磁理论V 电磁波与辐射8 单个粒子的辐射 匀速运动与匀加速运动的情况 单个粒子的辐射场满足: E=q((n^−β⃗)(1−β⃗2)(1−n^⋅β⃗)3R2+n^×[n^−β⃗ ...

  5. UA PHYS515A 电磁理论V 电磁波与辐射6 波导

    UA PHYS515A 电磁理论V 电磁波与辐射6 波导 波导(wave guide)是用来定向引导电磁波的结构,它为我们研究边界条件对电磁波传播的影响提供了一个简单的模型.假设电磁波的形式为 E=E ...

  6. UA PHYS515A 电磁理论V 电磁波与辐射4 反射与折射

    UA PHYS515A 电磁理论V 电磁波与辐射4 反射与折射 假设平面z=0z=0z=0是两种不导电介质的交界面,z<0z<0z<0的空间中介质的介电常数与磁导率为ϵ1,μ1\ep ...

  7. UA PHYS515A 电磁理论V 电磁波与辐射3 偏振

    UA PHYS515A 电磁理论V 电磁波与辐射3 偏振 之前谈到过电动力学讨论的就是作为源头的电荷与电流如何产生电磁场,以及由源头激发出的电磁场又如何反作用于这些电荷与电流以及存在于介质中的电荷.因 ...

  8. UA PHYS515A 电磁理论V 电磁波与辐射5 电磁波在介质中的传播

    UA PHYS515A 电磁理论V 电磁波与辐射5 电磁波在介质中的传播 在介绍麦克斯韦方程的时候,我们提到过 D⃗=E⃗+4πP⃗\vec D = \vec E + 4 \pi \vec PD=E+ ...

  9. UA PHYS515A 电磁理论V 电磁波与辐射2 电磁波的能量

    UA PHYS515A 电磁理论V 电磁波与辐射2 电磁波的能量 在讨论电磁场的能量时,我们引入了Poynting矢量,为了描述波动,我们把电场与磁场描述为时空的复变函数,因此Poynting矢量也需 ...

最新文章

  1. AVERAGE函数参数说明
  2. ubuntu java反编译,[原创]Ubuntu下简单配置反编译工具以及其开发环境
  3. 算法与数据结构 -- 排序和查找(五)
  4. 科研实习 | 中国科学院深圳先进技术研究院徐坤老师课题组招聘客座实习生
  5. redis 源码 dict.c 实现
  6. Android多媒体学习五:调用Android自带的播放器播放Audio
  7. python中的所有功能_python – 是否可以列出模块中的所有功能?
  8. Oracle错误:ORA-27121: unable to determine size of shared memory segment
  9. wp8.1 Study11:APP里文件读写和使用XML和Json序列化
  10. 绝地求生7月5日服务器维护,绝地求生7月5日更新到几点 吃鸡更新维护公告
  11. 【概率论】对弈输光模型,ruin model
  12. linux+patches工具_[经验分享] 最简安装: ubuntu+gnome3+常用工具
  13. three.js 使用DragControls.js 拖动元素
  14. C++之explicit构造函数
  15. 学计算机专业选择设计为类需要艺考,想学设计必须参加艺考吗?设计学类专业,新高考选科怎么选?...
  16. 【自动控制原理】渐进稳定与临界稳定的区别与联系
  17. local variable ‘count‘ referenced before assignment
  18. Tensorflow2.* 加载和预处理数据之用 tf.data 加载 Numpy数据(2)
  19. 二维光子晶体禁带的遗传优化算法实现
  20. 短视频风口持续 今日头条再投10亿补贴火山小视频

热门文章

  1. 【IM】关于在线学习(被动攻击学习和适应正则化学习)的理解
  2. 【正一专栏】曼彻斯特德比落幕——曼城独孤求败继续超神
  3. swiper左右切换按钮无效_BMD为ATEM切换台发布新的2 M/E和4M/E控制面板
  4. Python 技术篇-sha256()加密的使用方法,sha1、md5加密方法
  5. PyQt5 图形界面-基础准备:安装Qt Creater、Qt Designer工具,使用Qt Designer设计界面,安装PyQt5工具库
  6. CodeBlacks16\17\20版本下载
  7. 电源系统中纹波和噪声的处理
  8. 三角剖分(delaunay)拓扑结构 高维近邻
  9. This function or variable may be unsafe
  10. 锐捷交换机清除mac-address-table(mac地址表)