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

导体球或球壳在匀强外电场中的场强分布是电动力学中的经典题目，也是大家所熟悉的题目，但均匀电场中出现两个导体球或球壳(二聚体)时的场强却少有讨论。两导体球电极化所引起的导体球间电磁场的局域增强是表面等离激元增强和针尖等离激元增强最基本的特性之一，也是表面增强光谱的最主要原因。本文通过对单个导体球在均匀外电场中的场强分布和性质的求解分析，考虑在均匀外电场中沿电场线方向引入两个导体球，并对其电势场强进行计算分析，再用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，那么满足