导体球或球壳在匀强外电场中的场强分布是电动力学中的经典题目,也是大家所熟悉的题目,但均匀电场中出现两个导体球或球壳(二聚体)时的场强却少有讨论。两导体球电极化所引起的导体球间电磁场的局域增强是表面等离激元增强和针尖等离激元增强最基本的特性之一,也是表面增强光谱的最主要原因。本文通过对单个导体球在均匀外电场中的场强分布和性质的求解分析,考虑在均匀外电场中沿电场线方向引入两个导体球,并对其电势场强进行计算分析,再用Matlab软件画出直观的电势场强图像,进一步讨论了在匀强电场中两导体球之间间距和电场局域增强的关系。通过对这一电动力学经典习题的分析和扩展,增进对表面增强光谱的初步认识和了解,同时锻炼对计算物理和Matlab计算软件的运用能力。

1 问题引出

首先,简要讨论均匀电场中导体球的外电场分布。

设有一半径为R0的接地导体球置于均匀外电场E0中。导体球使空间分为两均匀区域——球外区域和球内区域。球外区域没有自由电荷,因此电势φ满足拉普拉斯方程。由球坐标系中拉普拉斯方程的通解,代入后可得球外电势为

(1)

其中an,bn为待定常数,由边界条件确定。最后可得到球外电势为

(2)

再对电势求负梯度,得到球外的电场强度为

(3)

这是电动力学教材中的典型例题[1]。由式(3)可以看出,导体球外靠近导体球表面的电场由于导体球在外电场的作用而重新分布。球外电场线在球表面处比原有均匀场密集,即电场线向导体端面会聚。分析发现沿外场方向的导体球端点处的电场强度是均匀电场强度的3倍,而垂直外场方向的导体球端点处的电场强度接近零。

考虑到导体球端点处的电场强度是外加均匀电场强度的3倍,电场强度极大地增强了。可以引入一个有趣的问题:如果在均匀外电场中,沿电场线方向放置两个相近的导体球,那么球外电势和电场强度的分布会如何,电场强度能增大到多少倍?这一问题可以看作表面等离激元增强和针尖等离激元增强等科研前沿的简化模型,即电子显微镜针尖和纳米颗粒之间的相互作用。课堂上通过从科研前沿和实验技术引入这一问题,可以增加学生学习电动力学的兴趣和锻炼学生对计算物理和Matlab计算软件的运用能力[2]。

2 问题求解

如果采用求解拉普拉斯方程的办法求解双球问题,无论在双球坐标系还是切球坐标系中进行,通常会因繁琐的分离变量过程而束手无策。而镜像法是电动力学中一种常用的数学方法,可以避免对拉普拉斯方程复杂繁琐的求解。下面使用镜像法来求解均匀电场中双导体球的空间电势。

图1 (a) 导体球和一个点电荷的镜像法示意图;(b) 匀强电场中单个导体球的镜像法示意图

首先讨论一个导体球和一个点电荷的静电感应现象,如图1(a)所示,处于真空的导体球半径为R0,右端点电荷Q放在距离球心a处。由于静电感应现象,球体表面将会产生感应电荷。如果导体球接地,由镜像法和对称性可看作在连线上有一点电荷Q1,并且在边界P点满足

(4)

由几何关系可以得到,只要选取Q1的位置使△OQ1P∽△OPQ,那么满足

(5)

由三角形相似关系式(4)和式(5)可以得出,

式(6)、式(7)分别确定了由Q1像电荷在导体球内部的位置和电荷大小。如果导体球不接地,要同时满足上述边界条件和导体球总电荷为零条件,那么在球心O添加一电荷量为-Q1的点电荷即可。

外加匀强电场可以看成是两个相距无限远的等量异号点电荷在其连线中点处产生的电场[3],单个导体球的静电感应现象,如图1(b)所示,由对称性条件与式(6)和式(7)可以得出,由±Q产生的两组像电荷在球心处的叠加为零,球内像电荷大小和距球心距离为

(8)

a趋向于无限大,则b趋向于零,因此在匀强电场中的像电荷可以看作一对电偶极子,其偶极矩p的大小为

(9)

由±Q在导体球表面产生的电场大小为

(10)

式(10)和式(9)联立消去Q/a2得,

(11)

其中p的方向和E0方向相同。这与教材和其他文献中单介质球的结果一致,其中的导体球可以看作介电常数为无限大的介质球[1-4]。

下面看电偶极子在导体球中产生的像电荷,其模型如图2所示。

图2 (a) 电偶极子在空间一个导体球时产生的像电荷;(b) 两个相同导体球时产生的像电荷

取+q和-q的间距为Δl,由式(3)、式 (4)可以得出电荷-q的像电荷q1和-q1,电荷+q的像电荷q2和-q2的大小和位置(其中-q1和q2在球心处):

(12)

对于电偶极子p0,Δl应为一无穷小量,满足:

因此电偶极子在球心处的像电荷为零,另外两个像电荷组成一对电偶极子,其偶极矩为

(15)

如图2所示,当空间电场方向有两个相同导体球时,由于空间中电场作用,会在两球心位置各产生一对像电偶极子p1,其偶极矩大小可由式(11)给出。像电偶极子p1会分别在另一个导体球中产生次级的电偶极子p2,接着产生p3…pn就像两面镜子一样无穷反射下去[5]。

由式(11)、式(14)、式(15)得出下表:

由上可以得到两相同导体球在均匀电场中的镜像电荷分布情况,空间电势情况只需要将所有像偶极子的电势和匀强电场的电势叠加即可。

由于电场关于球心连线旋转对称,所以以任意一个经过球心连线的切面建立平面直角坐标系,其中球心连线的中点为原点,以电场方向为x轴,其垂直方向为y轴。可由电偶极子空间电势公式得出在空间中导体球外任一点(x,y)处的电势为

(16)

导体球表面上电势处处相等,其电势可由边界电势得出。由静电场唯一性定理及此处级数表达式可得出此电势是收敛的。

当y=0时,得到x轴上的电势分布为

(17)

对电势求负梯度,便得到x轴上的场强分布:

(18)

3 结果讨论

3.1 空间中电势和电场强度分布

上面通过计算得到了均匀电场中沿电场线两导体球外的电势和电场强度的分布。下面通过Matlab作图来直观地分析空间中球外的电势和电场。

令两球心距离r=10m,E0=5N/C,金属球半径R0=2.5m,并取多项式(16)的前100项作空间电势方程的近似解。得到xOy平面空间等势线图像,如图3所示。图中A点、B点为两金属球沿x轴线方向在表面上的点,C点为沿x轴线两球心的中点。

图3 均匀外电场中两导体球球外空间的电势示意图

图3中两个圆为金属球边界,可以看出边界处于一条等势线上,因此满足上述边界条件。由静电场的唯一性定理可知,上述函数φx,y满足金属外部空间电势分布,是正确的解。

使用Matlab对函数φx,y进行梯度运算,可以得到电场的分布及其大小。将球内电场强度置零,可以得到空间中各位置电场强度的大小,如图4所示。图4中可以看出,在无穷远处电场强度恒定为E0,在沿电场线方向的金属球边缘处的电场强度会明显增大。这正是由于空间电场的电极化作用和两球之间的库仑力相互作用使球体产生感应电荷,其作用可用无限对电偶极子代替。

图4 均匀外电场中不同距离两导体球外的电场强度大小分布(图4(a)球心距离为10m;图4(b)球心距离为6m)。色带从0.0到1.0表示强度从零到最大

在图4中,两球心连线上电场强度的分布关于中心点C对称(其中A、B、C点位置同图3)。在图4(a)中,趋于两导体球外端表面A点处的电场强度明显增强。靠近导体球表面的B点处的电场强度强同样明显增强,而在两导体球中心处C点的电场强度稍弱。电场在导体球靠近中心处B点的场强稍稍大于该导体球另一端的A点处的场强(如图5)。这说明在均匀电场中导体球的旁边再放一个导体球时,会使原来那个导体球端点处的电场强度进一步增强。并且随着两球心距离的缩小,B点和C点的电场强度将会进一步增强。

图5 均匀电场中不同距离两导体球x轴上的场强分布(图5(a)球心距离为10m;图5(b)球心距离为6m)

3.2 x轴上ABC点处电场强度

当移动两个导体球,改变两导体球的相对位置时,会发现x轴上导体球附近的场强发生明显变化(如图4,图5)。固定A、B和C 3点为球心连线上金属球表面上两点和球心连线中点,分别求出这3个点的电场强度随两导体球的相对位置R0/r的值,如图6,得到如下结果:

图6 (a) x轴上A、B和C 3点处电场强度随R/r的变化;(b) 电场强度以10为底的对数

(1) x轴上第一个导体球的左端A点处的电场强度:当两导体球相距无穷远时,A处的电场强度为均匀场强的3倍,这与单个导体球在均匀电场处的情形相同;而当两个导体球相互靠近时,A点处的电场强度有所增加,但变化不明显。

(2) x轴上第一个导体球的右端B点处的电场强度:当两导体球相距无穷远时,B处的电场强度为均匀场强的3倍,这与单个导体球在均匀电场处的情形相同;当两个导体球相互接近时,一开始B点处场强的增强不明显(如图4(a));但当两导体球心间距离小于2.22倍球半径后,B点处的电场强度随相对位置R0/r的增加发生指数形式的增长。特别是当两导体球无限接近时,B点处的电场强度达到极值。其中,当两导体球球心间距离为2.0002倍导体球半径时,B点处的电场强度为均匀电场强度的100倍。

(3) x轴上两导体球中心C点处的电场强度:当两导体球相距无穷远时,C点处的电场强度为均匀场强,这说明此时该点处的电场强度基本不受两个导体球的影响。当两个导体球相互接近时,一开始C点处场强的增强不明显,但当两导体球心间距离约为2.6倍导体球球半径时,中心C点处的电场强度开始大于x轴上第一个导体球的左端A点处的电场强度。在这以后,C点处的电场强度增加越来越快,最后也呈现指数增长。

(4) 用Matlab作出x轴上A、B和C 3点处电场强度以10为底随R0/r的指数变化(如图6(b)所示)。可以看出:A点处的场强变化不大,B点场强在当两导体球心间距离小于2.22个球半径后的增长速率大于指数增长,C点一开始增长速率与指数增长相似但后来的增长速率也超过指数增长。

4 结论

本文从电动力学的经典题目——均匀电场中一个导体球在端点处的电场强度是外加均匀电场强度的3倍出发,提出均匀电场中存在两个导体球时的问题,通过模型求解和数值计算对问题进行详细的分析讨论。发现当两导体球球心间距离为2.0002倍导体球半径时,两导体球间端点处的电场强度将增加为均匀电场强度的100倍,并用数值计算给出了两球心连线上不同高对称点的场强随不同球心距离时的变化。通过对这一电动力学的经典题目的延伸拓展,使同学增加学习兴趣,并锻炼对计算物理和Matlab计算软件的运用能力。

参考文献

[1] 郭硕鸿. 电动力学[M]. 3版. 北京: 高等教育出版社, 2008:5.

[2] 华玲铃,杨阳. 大学物理教学中培养科学计算能力的研究[J]. 物理与工程, 2013, 23(3):37-41.

HUA L L,YANG Y. Research on the cultivation of scientific computing ability in college physics teaching[J]. Physics and Engineering, 2013, 23(3): 37-41. (in Chinese)

[3] 黄爽. 外电场中导体球体系的电势及荷电相关问题研究[D]. 贵阳: 贵州大学, 2017.6.

[4] 桑芝芳,李振亚.介质球在基质中的等效电偶极矩[J].大学物理, 2014, 33(6):4-6.

SANG Z F, LI Z Y. Equivalent electric dipole moment of dielectric sphere embedded in a matrix[J]. College Physics, 2014, 33(6): 4-6. (in Chinese)

[5] 吕晶. 点偶极子场中两导体球的电极化与相互作用研究[D]. 黑龙江: 哈尔滨工业大学, 2014.6.

基金项目: 教育部高等学校教学研究项目(DWJZW201701zn);河南省教育科学规划课题([2011]-JKGHAD-0241);中国高等教育学会大学素质教育研究课题(CALE201662);郑州大学教育教学改革研究与实践项目(2019ZZUJGLX308)。

作者简介: 王飞,郑州大学副教授,主要从事凝聚态物理和计算材料学研究工作, wfei@zzu.edu.cn。

引文格式: 董延更,赵佳羿,王飞. 双导体球在匀强外电场中的空间电场分布[J]. 物理与工程,2019,29(4):29-34.

END

更多精彩文章请点击下面“蓝字”标题查看:

  • 《物理与工程》2019年第3期目录

  • 顾牡:对于重新制定的《非物理类理工学科大学物理课程教学基本要求》的认识和体会

  • 黄昆:永远怀念吴大猷老师

  • 黄昆:我的研究生涯

  • 葛惟昆:香港科技大学授予荣誉理学博士的颂词

  • 朱邦芬院士:一代宗师黄昆

  • 朱邦芬院士:纪念中国半导体物理及固体物理奠基人黄昆先生

  • 朱邦芬院士:“减负”误区及我国科学教育面临的挑战

  • 朱邦芬院士:从基础科学班到清华学堂物理班

  • 朱邦芬院士:对培养一流拔尖创新人才的思考

  • 王炳燊、葛惟昆:玻恩-黄昆方程的历史回顾——纪念黄昆诞辰100周年

  • 贾惟义、葛惟昆:在黄昆先生指导下研究固体中的晶格振动、光色散及电声子耦合——纪念黄昆先生诞辰100周年

  • 李学潜教授:物理是一种文化

  • 李学潜教授:如何帮助物理系学生迈过从高三到大一这个坎

  • 李焱、杨宏:从超短光到超短超强光的突破——解读获2018年诺贝尔物理学奖的啁啾脉冲放大技术

  • 王雯宇:光和引力波专题 Ⅰ ——广义相对性原理、光速不变原理及引力论

  • 王雯宇:光和引力波专题Ⅱ——电磁及引力介质理论

  • 穆良柱:什么是物理及物理文化?

  • 穆良柱:什么是ETA物理认知模型

  • 穆良柱:什么是ETA物理教学法

  • 罗洪刚等:“课程群”组织体系: 对物理学基础课程教学体系的思考

  • 邓崇林:创新量纲分析重导玻尔模型能级公式

  • 徐劳立等:航班飞行中的物理问题

  • 物理教育的春天来了!

  • 吴国祯教授:我的国外研究生经历印象——应清华大学物理系“基科班20年·学堂班10年纪念活动”而写

  • 《物理与工程》2019年第1期目录

  • 陈佳洱,赵凯华,王殖东:面向21世纪,急待重建我国的工科物理教育

  • 王亚愚教授:清华物理系本科人才培养理念与实践

  • 葛惟昆教授:关于中外人才培养的几点思考

  • 王青教授:小班教学与翻转课堂:《费曼物理学Ⅱ》的10年教学实践——纪念费曼先生百年诞辰

  • 安宇教授:为什么传统的课堂讲授模式需要改变

  • 刘玉鑫教授:关于本科生物理基础课程教学和教材编著的一些思考

  • 沈乾若:重创理科教育的美加课程改革

  • Henderson C:美国研究基金支持下的物理教育研究及其对高等物理教育的影响

《物理与工程》期刊是专注于物理教育教学研究的学术期刊,是中国科技核心期刊,1981年创刊,欢迎踊跃投稿,期刊投审稿采编平台:

http://gkwl.cbpt.cnki.net

欢迎关注

《物理与工程》微信公众号

matlab中表示拉普拉斯分布_双导体球在匀强外电场中的空间电场分布相关推荐

  1. matlab中表示拉普拉斯分布_拉普拉斯分布的随机数

    一.功能 产生拉普拉斯分布的随机数. 二.方法简介 1.产生随机变量的组合法 将分布函数\(F(x)\)分解为若干个较为简单的子分布函数的线性组合 \[F(x)=\sum_{i=1}^{K}p_{i} ...

  2. matlab中表示拉普拉斯分布_神奇的正态分布

    在统计学中有各种各样的分布,称为统计分布,例如有离散型的伯努利分布.二项分布.超几何分布.几何分布.负二项分布.泊松分布,有连续型的均匀分布.指数分布.t分布.卡方分布.F分布.正态分布等等,其中正态 ...

  3. matlab中表示拉普拉斯分布_分布拟合——正态/拉普拉斯/对数高斯/瑞利 分布

    作者:桂. 时间:2017-03-16  20:30:20 声明:欢迎被转载,记得注明出处~ 前言 本文为曲线与分布拟合的一部分,主要介绍正态分布.拉普拉斯分布等常用分布拟合的理论推导以及代码实现. ...

  4. matlab中表示拉普拉斯分布_深度优化局部拉普拉斯金字塔滤波器。

    微信公众号:OpenCV学堂关注获取更多计算机视觉与深度学习知识 觉得文章有用,请戳底部[好看]支持 算法概述 基于局部拉普拉斯金字塔的Edge-aware滤波器是在2011年由Adobe 公司的研究 ...

  5. matlab中表示拉普拉斯分布_CHAPT1:场论;电磁学和微波学的基本的数学手段和表示...

    物理学中把某个物理量在空间一个区域内的分布称为场.从各种场的取值性质来看可以分成两大类,一类是每个点对应一个数值,这种场统称为标量场,如温度场.密度场等;另一类是每 个点对应一个向量,这种场称为向量场 ...

  6. 拉普拉斯分布_中心极限定理:从高尔顿板到麦克斯韦分布

    ​神奇的正态分布源于"加". 撰文 | 张和持 时隔多年,或许你早就记不得16岁那年夏天高中闷热的教室,但可能会记得有一天数学老师说着要给大伙看个稀奇--一块祖传的高尔顿板.尽管班 ...

  7. mysql查询两个表中的不同数据_如何实现大数据在多工作表中精准查询

    大家好,今日继续讲解VBA数据库解决方案,今日讲解第40讲:利用ADO实现同一文件的多个工作表精确查询.其实这讲的内容利用我之前讲解的知识点也可以自己完成.这里只是拿出来让大家再熟悉一下EXCEL的A ...

  8. sc在java中是什么意思_基础入门:如何在你的HTML中使用Javasc_html

    JavaScript实现html的动态效果,通常我们称为网页特效,在本站网页特效栏目代码比较多!那么如何把Javsscript加到HTML中呢? javascript是一种与HTML一起工作的客户端脚 ...

  9. 拉普拉斯分布_理解拉普拉斯特征映射中的优化问题的约束条件

    引言:在学习拉普拉斯特征映射(Laplacian Eigenmaps, LE)的过程中,发现大多数参考资料仅列出了其中的最优化问题,然后直接过渡到特征值问题,对于该优化问题,特别是其中的约束条件解释的 ...

最新文章

  1. 「AI初识境」什么是深度学习成功的开始?参数初始化
  2. 智能安防赛道上 国产AI芯片力挫群雄
  3. mxnet基础到提高(24)-C++-符号,变量读写
  4. DBus glib 各数据类型接收与发送详解—C语言(3)
  5. 数据库备份策略 分布式_管理优秀的分布式数据团队的4种基本策略
  6. qtreewidget 获取根节点_详解去中心化信任根dRoT技术
  7. NopCommerce 2.5的部署
  8. c语言一年日历程序,c语言日历程序
  9. linux批量重命名脚本,Linux批量命名文件SH脚本整理
  10. 电子围栏判断_脉冲电子围栏和张力围栏之间的区别
  11. SQL常见错误及处理方法
  12. 几款ReactJS最优秀的UI框架
  13. E-Prime教程 | 打mark标记
  14. opc服务器协议gsd,S7-300如何通过GSD文件实现PROFIBUS DP主从通讯
  15. 科普 | 天才在左,疯子在右
  16. html5摄像头手势识别,基于OpenCV的摄像头动态手势轨迹识别及其应用
  17. 【LeetCode】Algorithms 题集(二)
  18. 微信公众号Java开发记录(二)接收、回复消息
  19. 纯java编写的联机五子棋项目(附带开源链接)
  20. idea 报错Output directory is not specified错误

热门文章

  1. 取出字符串strName前4个汉字
  2. java枚举的线程安全及序列化
  3. 设计模式-策略模式(Strategy)
  4. Leetcode475.Heaters供暖器
  5. 给RabbitMQ发送消息时,设置请求头Header。
  6. 手动安装vue-devtools
  7. 反射(3)—动态、静态代理
  8. AndroidStudio_开发工具的设置_快捷键设置_编辑器设置---Android原生开发工作笔记72
  9. k8s集群资源监控-监控指标和方案---K8S_Google工作笔记0052
  10. SpringCloud工作笔记068---登录堡垒机_下载文件报错_Xshell向Linux上传下载文件_lrzsz_ZModem