均匀带电正多边形的静电自能：精确化简、闭式表达式及向圆盘的收敛

《APS Open ScienceNEW》：Electrostatic self-energy of a uniformly charged regular polygon: Exact reduction, closed-form expressions, and convergence to the disk

　　

我们评估了均匀带电正多边形的静电自能，将其视为嵌入三维空间中的平面电荷分布，并通过库仑势进行相互作用。利用边界分段线性的特点，可将四维库仑积分精确简化为有限的初等扇形贡献之和，该结果对任意边数均成立。由此构造得到均匀带电三角形、正方形、五边形和六边形的闭式表达式。五边形的结果虽然涉及四次数域，但形式出人意料地简洁。随着边数增加，能量单调递减，趋近于均匀带电圆盘的对应值。边界微扰分析给出了主导的相对修正项，该修正项随边数的平方反比减小。能量的二阶变分产生了一个非解析修正，其正比于边数的四阶反方与该数对数之积。这一项为多边形角点的直接特征。精确结果的算术结构也解释了为何在均匀带电六边形之后便无法再获得简洁的闭式表达式。所有结果均经高精度数值计算验证，并为二维受限库仑体系的研究提供了有用的基准。

**I. 引言**

静电学在经典场论、势论和边值方法的发展中发挥了核心作用，其中涉及导体上电荷分布的问题为解析技术和数值技术都提供了重要的基准[1–7]。这些问题在当代受限库仑体系的研究中仍然相关，几何形状和边界效应强烈影响平衡构型与能量[8–16]。求导体表面上平衡电荷分布是静电学中的一个基本问题[2]。对于闭合导体，面电荷密度主要由局部几何决定，而对于开放导体，边界和棱边的存在则起主导作用[3]。在平面区域中，圆形圆盘是唯一的例外，其平衡电荷密度可用闭式解析表达[2]。由于已知的精确结果非常有限，因此在有限区域上均匀分布电荷的互补问题常被采用为一个有用的、可处理的模型系统[17]。理解这种均匀分布的静电能如何随约束边界的形状变化具有独立的理论意义，对于解读受限带电粒子的数值模拟也至关重要，其中均匀背景分布经常作为参考态出现[18,19]。

有两点事实使均匀电荷密度能量不仅仅是一个模型量。根据汤姆逊定理，在固定电荷量下，平衡（等势）分布使能量最小，因此均匀电荷密度能量是对应导体能量的严格上界，并校准导体问题的数值解；在蒙特卡罗或分子动力学模拟中，对于被限制在二维（2D）多边形区域（通常为正方形/矩形）内且带有中和均匀背景的电荷，背景自能作为一个必须精确知道的附加常数项，才能可靠地提取相关能。

本文研究一个外接圆半径为 $R$、面积为 $A$ 的正 $n$ 边形，该多边形载有均匀面电荷密度 $\sigma = Q/A$（$Q$ 为总电荷）。该正多边形是嵌入在普通三维空间（3D）中的平面（零厚度）二维电荷分布，因此两个面元之间通过三维库仑势 $k/r$ 相互作用，其中 $r = |\mathbf{r} - \mathbf{r}'|$，$\mathbf{r}$（$\mathbf{r}'$）为位置矢量。这是二维电子体系、量子点、带电胶体单层的设置，必须与沿第三方向具有平移不变性的体系（例如无限长线电荷）的静电学区别开来，后者的相互作用为对数形式。

我们的方法扩展了此前为均匀带电平面区域开发的垂直切片法[20–22]。通过对横向坐标进行积分，库仑相互作用通过一个有效核来表示，该核仅依赖于纵向坐标和局部条带宽度。这一过程将自能的四维原始积分简化为二维形式，大幅简化了计算。对于正多边形，边界分段线性的特点允许将积分区域划分为有限个仿射扇形。在每个扇形内，边界由简单的线性函数描述，使剩余积分可以系统地进行。利用多边形的几何对称性进一步减少了计算量，并得到适合解析研究和数值评估的简洁表达式。

对于正三角形和正方形，我们的分析复现了此前已知的精确结果[17,18]。更重要的是，它还给出了五边形和六边形新的闭式表达式。在这些情形中，五边形尤为引人注目，因为其几何引入了比在简单二次数域内可表达精确能量的多边形更为丰富的代数结构。尽管增加了复杂性，最终结果简化为一个出人意料地简洁的表达式，仅涉及与黄金比例和数字5相关的一小组对数常数。我们还研究了大 $n$ 极限，即正多边形趋近于圆形圆盘的情形。在此区间内，我们恢复了均匀带电圆盘的精确定能，并通过边界微扰分析推导了主导的有限 $n$ 修正。该修正表明，静电自能随边数增加而单调递减，平滑地收敛到圆盘的值。换言之，在外接圆半径固定的正多边形中，边数越多的形状对均匀带电圆盘静电行为的逼近程度越高。

除主导修正外，我们还考察了渐近展开中的高阶项，发现向圆盘的收敛并非纯代数的。相反，四阶贡献中包含对边数的对数依赖，反映了即使边界已近乎圆形，多边形角点仍持续的影响。这一对数项是边界奇点的特征标志，其系数和伴随常数均可由静电能量关于圆盘的二阶变分得到闭式表达（附录 B）。黎曼 zeta 函数的特殊值以透明的原因出现：$n$ 重对称性将边界波纹限制在角调和 $e^{im\theta}$ 上，其中 $m$ 为 $n$ 的倍数，振幅正比于 $n^{-2}$，而库仑能量的二阶变分随调和指标的对数增长；对 $m$ 求和因此产生 $\zeta(2)$ 和 $\zeta(4)$，从而固定了对数振幅及其 $n^{-4}$ 伴随项。

本文所呈现的结果与经典和量子设置下的广泛受限库仑体系相关。例子包括量子点中受限的电子[23]、捕获离子装配体，以及通过光学捕获技术操控的胶体粒子[24]。此外，本文推导的精确表达式和渐近公式为数值方法提供了有价值的基准，包括有限带电体系的蒙特卡罗模拟，其中可靠的解析参考值往往稀缺。

本文的组织如下：我们在第 II 节中阐述问题并推导静电自能的化简核表示。第 III 节建立了精确计算相关积分所依赖的仿射扇形框架。对特定正多边形的应用和推广到任意 $n$ 的讨论分别在第 IV 节和第 V 节中给出。第 VI 节分析了趋向圆形圆盘极限的渐近行为，第 VII 节报告了解析公式的数值检验。第 VIII 节总结了结果，并讨论了该方法向不规则多边形、非均匀密度和等势（导体）问题的扩展。

**II. 形式化与精确核**

我们考虑外接圆半径为 $R$、面积为 $A$ 的正 $n$ 边形，其中包含均匀分布的电荷 $Q$。相应的面电荷密度在整个多边形域上为常数：
$$\sigma = \frac{Q}{A} \tag{1}$$
其中
$$A = \frac{1}{2} n R^2 \sin\left(\frac{2\pi}{n}\right) \tag{2}$$
储存在该物体中的静电自能为
$$E = \frac{\sigma^2}{2} \iint_{\mathcal{D}} \iint_{\mathcal{D}} K(\mathbf{r}, \mathbf{r}') \, d^2\mathbf{r} \, d^2\mathbf{r}' \tag{3}$$
其中 $\mathbf{r}$（$\mathbf{r}'$）为二维矢量，$k$ 为库仑电常数，$\mathcal{D}$ 为积分的多边形域（面积为 $A$）。核 $K$ 是标准库仑势限制在包含电荷的平面 $z=0$ 上的形式。它不是二维拉普拉斯算子的格林函数 $G$，后者描述沿第三 $z$ 方向均匀延伸的电荷。因为 $K$ 在有界平面域上可积，所以公式 (3) 中的自能是有限的。圆盘（$n \to \infty$）是经典的例子：对于均匀面电荷密度，它给出 $E$ [式 (31)]，而对于等势圆盘，面电荷密度 $\sigma(r)$ 和电容 $C$ 可由同一核得到[2]。

我们将多边形定向使一个顶点位于正 $x$ 轴上，图形关于该轴对称，
$$y_{\text{top}}(x) = f(x), \quad x \in [-R, R] \tag{4}$$
其中 $f(x)$ 为分段线性的上边界。

令 $x = s + t$，$x' = s - t$，对 $y$ 和 $y'$ 的内双重积分（其中 $y \in [0, f(x)]$，$y' \in [0, f(x')]$）是初等的。两个区间的卷积具有梯形权重 $w(s) = f(s+t) + f(s-t)$（当 $|t| \leq$ 适当范围时），因此
$$\int_{-R}^{R} dx \int_{-R}^{R} dx' \, w(s) = 2 \int_0^{R} dt \int_{-R}^{R} ds \, [f(s+t) + f(s-t)] \tag{5}$$
执行积分得到精确核
$$K_{\text{eff}}(s) = 2\left[\sqrt{(s+w/2)^2 + (f - f')^2} - \sqrt{(s-w/2)^2 + (f - f')^2}\right] \tag{6}$$
两个平方根项和整体的因子 2 是不可或缺的：它们源于条带的有限宽度，不可省略。利用公式 (6)，自能简化为二重积分
$$E = \frac{\sigma^2}{2} \int \int K_{\text{eff}}(s) \, ds \, dt \tag{7}$$
当 $s \to 0$ 且 $t \to 0$ 时，$K_{\text{eff}} \sim \ln(1/|t|)$，为可积对数奇点，因此公式 (7) 绝对收敛。使用无量纲系数
$$\alpha_n = \frac{E}{\frac{\sigma^2 A}{2k}} \tag{8}$$
该系数仅依赖于 $n$。

**III. 仿射扇形分解**

**A. 几何：必须保留截距**

对于正 $n$ 边形，上边界由 $n/2$ 条直线段组成（图 1）。设顶点为 $V_j$，其中 $j = 0, 1, \ldots, n/2$，按 $x$ 递增排列，$V_0$ 为最左顶点，$V_{n/2}$ 为轴上的顶点。在扇形 $j$ 上，
$$f(x) = a_j + b_j x \tag{9}$$
其中斜率为 $b_j$，截距为 $a_j$。记端点高度为
$$h_j = f(x_{\text{left},j}), \quad h_{j+1} = f(x_{\text{right},j}) \tag{10}$$

图 1. $n=5$ 时的几何结构。上边界 $f$（蓝色）是分段线性的。它仅在顶点 $V_{n/2}$ 处与轴相交；在左端 $V_0$（黑色）处多边形以垂直边闭合；每个内部顶点处于严格正的高度（红色），因此其仿射片段具有非零截距，必须保留。

一个关键的几何事实支配了整个计算：$f(x)$ 在顶点 $V_{n/2}$ 处为零，对于偶数 $n$ 也在 $V_0$ 处为零；对于奇数 $n$，多边形在左侧以垂直边闭合，$f(x_{\text{left}}) > 0$。在每个内部顶点处，边界处于严格正的高度 $h > 0$，因为多边形是凸的且只有位于对称轴上的顶点接触该轴。因此，对于所有内部扇形，截距 $a_j$ 不为零，端点高度 $h_j$ 在扇形积分中与斜率具有同等重要性。

这正是朴素处理失败之处。人们会尝试将每个扇形平移至局部变量，使 $f$ 在扇形的一个端点处为零；这仅对于终止于轴上的扇形是正确的（最右扇形以及对于偶数 $n$ 的最左扇形）。在每个扇形中将 $a_j$ 替换为零会丢弃截距，移除条带宽度的主导贡献，并得到一个不仅不准确还可能为负（因而非物理）的表达式。下面正确的化简在整个过程中保留了 $a_j$（等价地，$h_j$）。

根据积分对 $\mathcal{D}$ 划分为扇形矩形的可加性，公式 (8) 变为有限和
$$\alpha_n = \sum_{j} \alpha_j \tag{11}$$

**B. 自扇形积分**

固定一个扇形并省略下标 $j$。设 $f$ 在宽度为 $w$ 的条带上，换到和差变量 $s = (x+x')/2$ 和 $t = (x-x')/2$，其中 $s \in [0, R]$，$t \in [-R, 0]$。由于 $K$ 关于 $t$ 是偶函数，雅可比行列式和两个半平面组合给出，在 $\eta = |t|$ 之后
$$I_{\text{self}} = 2 \int_0^R ds \int_0^{\eta_{\max}(s)} d\eta \, K_{\text{eff}}(s, \eta) \tag{12}$$
其中 $K$ 为通过 $\eta$ 表达的核 (6)。定义
$$\Phi(s, \eta) = \text{中间函数} \tag{13}$$
$\eta$ 的积分是初等的，得到精确结果
$$I_{\text{self}} = \frac{1}{4}\left[\text{表达式}\right] \tag{14}$$
其中
$$F(\text{端点数据}) = \text{由端点数据定义的函数} \tag{15}$$
公式 (15) 中剩余的 $s$ 积分是初等的：将公式 (13) 代入产生形式为 $\arctanh(\cdot)$ 和 $\ln(\cdot)$ 的项，每一项都积分反双曲和端点数据的代数函数。公式 (14) 已通过与直接二维求积比对验证至 25 位有效数字，包括 $a_j > 0$ 的扇形。在公式 (14) 中令 $a_j = 0$（即虚假的"零截距"选择）可恢复朴素处理的错误表达式，并明确展示了端点高度如何进入。

**C. 交叉扇形积分**

对于不同扇形 $j \neq j'$，两个条带占据不相交的 $x$ 区间，因此 $\eta$ 具有固定符号，令 $u = x - x'$，右侧条带上的内积分给出 $\ln(u)$ 和 $1/u$ 项，在固定 $s$ 下常数为 $c_j$。与公式 (13)–(15) 相同的原函数族适用，因此 $I_{\text{cross}}$ 同样是初等的。相邻扇形共享一个角点，其中 $\eta \to 0$ 而 $u$ 有限；由此产生的 $\ln(\eta)$ 行为是可积的，因此 $I_{\text{cross}}$ 对所有扇形对都是有限的。

## 方法概述

总之，在对横向坐标进行积分消去之后，任意正 $n$ 边形的自能即为扇区对积分的精确有限和[式 (11)–(15)]，每一项都是扇区数据的反双曲函数与平方根的基本组合。该表示对每一个 $n$ 均成立，且不涉及任何近似。该和式能否化为紧凑的闭式表达式是一个独立的算术问题，由多边形几何常数之间的代数关系决定；该问题将在第四、五节及附录 C 中讨论。

---

## IV. $n=3$ 和 $n=4$ 的闭式表达式

### A. 等边三角形（$n=3$）

以一个顶点位于坐标轴上时，上边界为单一线段
$$y = -\sqrt{3}\,x + \sqrt{3} \qquad (16)$$
因此左顶点处 $y=0$，而位于坐标轴上的顶点处 $y=\sqrt{3}$。（注意斜率为 $-\sqrt{3}$，而非 $-1$。）计算式 (14) 并代入 $A$ 的值，可得闭式表达式
$$S_3 = \frac{3\sqrt{3}}{2}\left(\ln 2 + \frac{2\pi}{3}\right) \qquad (17)$$
用边长 $a$ 表示，即得 $S_3 = \frac{3\sqrt{3}}{2}\left(\ln 2 + \frac{2\pi}{3}\right)\frac{R}{a}$，与文献 [18] 一致。

### B. 正方形（$n=4$）

对于以一个顶点位于坐标轴上的正方形，
$$y = \begin{cases} -x + 2, & 0 \le x \le 1, \\ \quad x + 1, & 1 \le x \le 2, \end{cases} \qquad (18)$$
这是两段斜率分别为 $\pm 1$ 的仿射扇区在顶点 $y=2$ 处交汇，其中 $y=2$——同样是一个具有非零高度的内部顶点。计算式 (14) 及面积项 $A$ 的交叉项，可得
$$S_4 = 4\left(\ln 2 + \frac{\pi}{4} + \frac{1}{2}\right) \qquad (19)$$
以边长 $a$ 为单位表示，
$$S_4 = \frac{2\left(\ln 2 + \frac{\pi}{4} + \frac{1}{2}\right)}{1} \qquad (20)$$
与文献 [17,18] 一致。

两种闭式表达式均得到独立验证：式 (17) 的数值积分值与 $S_3$ 在小数点后 18 位吻合，式 (19) 则通过整数关系检测从高精度值中还原得到。

---

## V. 一般 $n$ 的结果：精确约化与基准值

对于一般的 $n$，令 $\omega = 2\pi/n$。有序的上顶点为 $x_k = \cos\left(\frac{2k-1}{n}\pi\right)$（取满足条件的 $k$ 值），由此可直接得到扇区数据 $m_k$、$c_k$，并依据式 (8)、(11) 和 (14) 组装 $S_n$。

由于每一个 $m_k$ 均为初等量，对每一个 $n$ 都存在 $S_n$ 的精确闭式表达式；问题仅在于它能否化为紧凑的表达式。这取决于顶点的算术性质，而非仅由 $n$ 决定。正六边形的常数域为 $\mathbb{Q}(\sqrt{3})$，其自扇区贡献（以 $R^2$ 为单位）简洁明了：
$$\text{自扇区积分} = \frac{3\sqrt{3}}{4}\left(\ln\frac{2+\sqrt{3}}{2} + \frac{2\pi}{3}\right) \qquad (21)$$
$$\text{边界贡献} = \frac{\pi\sqrt{3}}{4} \qquad (22)$$
两个交叉扇区积分同样为初等积分。其内层积分返回第三节 C 中的显式原函数；外层积分经过一次分部积分后，可化为如下类型的积分：
$$\int \frac{dx}{\sqrt{ax^2+bx+c}} \qquad (a \text{ 为 } x \text{ 的二次多项式}),$$
该积分仍为初等积分——不出现双曲对数函数。两个交叉扇区积分均以此方式求解：斜率-斜率积分在计算机代数系统下直接封闭求解；斜率-平段积分则由同一显式原函数得到（经求导验证），在共享顶点处显式取可去极限。对独立计算的高精度值进行整数关系检测，返回相同的闭式表达式（附录 A），构成独立验证。结果如下：
$$\text{斜率-平段积分} = \frac{3\sqrt{3}}{8}\ln\frac{2+\sqrt{3}}{2} \qquad (23)$$
$$\text{斜率-斜率积分} = \frac{\sqrt{3}}{4}\left(\ln\frac{2+\sqrt{3}}{2} - \frac{\pi}{6}\right) \qquad (24)$$
对应斜率-平段（$k=1$）和斜率-斜率（$k=2$）对。所有常数均位于 $\mathbb{Q}(\sqrt{3})$ 中；特别地，$\sqrt{3}$ 出现。组装之后，取 $A = 3\sqrt{3}R^2/2$，
$$S_6 = \frac{3\sqrt{3}}{2}\left(\ln\frac{2+\sqrt{3}}{2} + \frac{2\pi}{3}\right) + \frac{3\sqrt{3}}{8}\ln\frac{2+\sqrt{3}}{2} + \frac{\sqrt{3}}{4}\left(\ln\frac{2+\sqrt{3}}{2} - \frac{\pi}{6}\right) \qquad (25)$$
在精确算术下化简，各项合并，得到紧凑的闭式表达式
$$S_6 = \frac{3\sqrt{3}}{2}\left(\ln 2 + \ln\frac{2+\sqrt{3}}{2} + \frac{\pi}{4}\right) = \frac{3\sqrt{3}}{2}\ln 2 + \frac{3\sqrt{3}}{2}\cdot\frac{\pi}{4} + \frac{3\sqrt{3}}{2}\ln\frac{2+\sqrt{3}}{2} \qquad (26)$$
这是逐精确项求值，而非数值拟合。作为独立验证，对独立计算的 $S_6$ 的 90 位精度值施加整数关系检测，返回相同表达式，残差低于 $10^{-88}$（数值上 $S_6 \approx 0.931002512419454$，与基准值在小数点后全部位数吻合）。据我们所知，此正六边形的值此前未曾出现；结合 $S_3$ 和 $S_4$，表明紧凑闭式表达式在正方形之后依然存在。

正五边形是初等情况中最复杂的。其顶点高度为嵌套根式 $\sqrt{(5\pm\sqrt{5})/8}$，因此常数域为四次域 $\mathbb{Q}(\sqrt{5},\sqrt{2-2/\sqrt{5}})$ 而非 $\mathbb{Q}(\sqrt{3})$。然而其超越内容却出奇地少：出现的每一个 $m_k$——斜率常数 $m_1$ 和 $m_2$（$m_2 = \sqrt{(5-\sqrt{5})/(5+\sqrt{5})}$）以及自扇区和交叉扇区中的端点常数——均可化为 $\ln 2$ 和 $\pi$，而每一个代数前因子（斜率、高度及端点根式 $x_1$、$x_2$）均位于 $\mathbb{Q}(\sqrt{5})$ 中。两个扇区给出（以 $R^2$ 为单位），其中 $x_1 = \sqrt{(5+\sqrt{5})/8}$，$x_2 = \sqrt{(5-\sqrt{5})/8}$：
$$\text{扇区 1} = \frac{1}{2}\left(2x_1\,m_1\,\text{arcsinh}\frac{m_1}{\sqrt{1+m_1^2}} + (1+m_1^2)\ln\frac{m_1\sqrt{1+m_1^2}+1}{m_1\sqrt{1+m_1^2}-1}\right) \qquad (27)$$
$$\text{扇区 2} = \frac{1}{2}\left(2x_2\,m_2\,\text{arcsinh}\frac{m_2}{\sqrt{1+m_2^2}} + (1+m_2^2)\ln\frac{m_2\sqrt{1+m_2^2}+1}{m_2\sqrt{1+m_2^2}-1}\right) \qquad (28)$$
$$\text{交叉项} = x_1 x_2\left(\ln\frac{x_1+x_2}{x_1-x_2} - \frac{2\,\text{arcsinh}(x_1/x_2)}{\sqrt{1+x_1^2/x_2^2}}\right) \qquad (29)$$
组装 $S_5$，取 $A = \frac{5}{2}\sin\frac{2\pi}{5}R^2$，不含对数的部分完全对消，结果为纯对数形式：
$$S_5 = \frac{5\sqrt{5}}{8}\left(\frac{4}{5}\ln 2 + \frac{2}{5}\ln\frac{1+\sqrt{5}}{2} + \frac{2\pi}{5}\right) \qquad (30)$$
数值上 $S_5 \approx 0.970766373674870$，与基准值在小数点后全部位数吻合。因此，正五边形并非第一个非紧凑的情形。紧凑闭式表达式在正五边形和正六边形中均得以延续。正五边形仅将常数域从二次提升到四次，这是其奇数顶点数的代价；先前的工作仅报告了直到正方形的闭式表达式 [17,18]。

一般 $n$ 的 $S_n$ 高精度值收集于表 I；所示各位数均为精确值，可供作为基准使用。

### 表 I

| $n$ | 无量纲自能系数 $S_n$（均匀带电正 $n$ 边形，外接圆半径 $R$） |
|:---:|:---:|
| 3 | 1.268 568 201 195 128 |
| 4 | 1.051 188 334 404 826 |
| 5 | 0.970 766 373 674 870 |
| 6 | 0.931 002 512 419 454 |
| 7 | 0.908 177 021 093 881 |
| 8 | 0.893 786 112 507 559 |
| 9 | 0.884 100 930 176 800 |
| 10 | 0.877 259 474 376 232 |
| 11 | 0.872 242 198 882 440 |
| 12 | 0.868 450 779 575 236 |
| $\infty$ | 0.848826363156775 ($= 8/(3\pi)$) |

$n=3, 4, 5, 6$ 的值为闭式表达式 (17)、(19)、(30) 和 (26)；其余值由式 (8) 的高精度数值积分得到。最后一行为圆盘极限 $n\to\infty$。各值由附录 A 的方案生成，作者可提供更高精度的数据。

该序列关于 $n$ 严格递减：在外接圆半径固定的情况下，多边形趋近于圆盘，系数从上方单调趋向 $8/(3\pi)$。

基准值定量地说明了约束区域的形状对能量的影响有多强。相对于同一外接圆半径的圆盘，三角形储存的能量多 $49.4\%$，正方形多 $23.8\%$，五边形多 $14.2\%$，六边形多 $9.3\%$，七边形多 $6.7\%$，八边形多 $5.1\%$；当 $n>8$ 时，超出部分低于 $5\%$，且由第六节的展开式在优于 $0.1\%$ 的精度上复现。在具有均匀中和背景的多边形胞室中的受限库仑系统模拟中 [18,19]，$S_n$ 即为必须扣除的精确背景-背景项，以分离出关联能；闭式表达式和表 I 提供了任意所需精度的该常数。由于导体的平衡分布使电荷固定下的能量最小化，$S_n$ 同时也是电容为 $C_n$ 的导体多边形板的无量纲能量 $E_n$ 的严格上界（见第八节）。

---

## VI. 圆盘极限与渐近展开

### A. 圆盘极限

当 $n\to\infty$ 且 $R$ 固定时，正 $n$ 边形收敛于半径为 $R$ 的圆盘，无论是逐点收敛还是面积收敛。核函数 (6) 被一个与 $n$ 无关的可积优函数所控制，因此由受控收敛定理，式 (8) 中的积分趋于其对应圆盘的取值。对于均匀带电圆盘，自能为经典结果
$$S_\infty = \frac{2\pi R^2}{3}\left(\frac{3\pi}{4} - \frac{8}{3}\right)^{-1} \cdot R^2 = \frac{8R^2}{3\pi} \qquad (31)$$
这意味着
$$\lim_{n\to\infty} S_n = \frac{8}{3\pi} \approx 0.848826363156775 \qquad (32)$$
表 I 复现了这一结果。该结果确立了圆盘作为正多边形序列的渐近极限。剩余的任务是量化有限的 $1/n^2$ 修正项，并确定多边形边界的几何特征在趋于圆盘极限过程中如何延续。

### B. 首项系数：边界微扰推导

$1/n^2$ 系数可以精确获得。将自能的无量纲系数写为
$$S_n = \frac{8}{3\pi}\left(1 + \frac{c_2}{n^2} + \frac{c_4 \ln n}{n^4} + \frac{c_4'}{n^4} + \cdots\right) \qquad (33)$$
其中纯几何四重积分为
$$\mathcal{I} = \frac{3\pi}{8R^2}\iint_A\iint_A \frac{dA_1\,dA_2}{|\mathbf{r}_1-\mathbf{r}_2|} \qquad (33)$$
圆盘对应 $\mathcal{I} = 1$。

在极坐标下，外接圆半径为 $R$ 的正 $n$ 边形的边界为 $r(\theta) = R\cos(\pi/n)/\cos\left(\theta - k\cdot 2\pi/n\right)$，其中在每条边内 $|k\cdot 2\pi/n + \theta| \le \pi/n$，$\theta$ 从边的中点测量；$r < R$ 且 $r \ge R\cos(\pi/n)$（多边形位于其外接圆内）。

对边界位移作一阶展开，
$$S_n \approx S_\infty + \oint_{\partial D} \phi_D(\mathbf{r})\,\delta r(\theta)\,d\theta \qquad (34)$$
其中 $\phi_D(\mathbf{r})$ 是均匀单位密度圆盘 $\mathcal{D}$ 的势，$\delta r$ 是有符号的边界层。在圆盘边缘上，该势为常数 $\phi_D(R) = 1$，因此仅有平均位移存活：
$$\oint \phi_D\,\delta r\,d\theta = \phi_D(R)\oint \delta r\,d\theta \qquad (35)$$
利用 $\oint \delta r\,d\theta = \Delta A = \frac{A_\infty - A_n}{R}$，
$$\frac{A_\infty - A_n}{R} = \frac{\pi R^2 - \frac{n}{2}R^2\sin(2\pi/n)}{R} = R\left(\pi - \frac{n\sin(2\pi/n)}{2}\right) \qquad (36)$$
因此 $S_n \approx S_\infty + 1\cdot\frac{A_\infty-A_n}{R^2} \cdot\text{(适当归一化)}$。利用 $A_\infty - A_n \approx \frac{\pi R^2}{3n^2}$（当 $n$ 大时），
$$\frac{c_2}{n^2} = \frac{1}{3n^2}\left(\frac{3\pi^2}{8} - 1\right)^{-1} \cdot \frac{\pi}{3} \qquad (37)$$
由此得
$$c_2 = \frac{\pi}{8} \qquad (38)$$
由于边缘势不依赖于 $\theta$，$\delta r$ 的振荡部分（即边缘波纹，具有 $n$ 的倍数的角谐波）在一阶不贡献；它在二阶进入并控制下一项。

### C. 四阶项：角对数

下一个修正项并非 $n$ 的纯幂次。其结构来源于能量关于圆盘的二阶变分，在附录 B 中以闭式形式求解；此处陈述结果并说明其数值验证方式。$n$ 重对称性将边界波纹限制在谐波 $k\cdot n$（$k=1,2,\ldots$）上，且库仑能量的边界海森矩阵随 $n$ 对数增长 [式 (B7)]，因此四阶项获得一个 $\ln n$：
$$S_n = \frac{8}{3\pi}\left(1 + \frac{\pi}{8n^2} + \frac{A\ln n + B}{n^4} + O\left(\frac{\ln n}{n^6}\right)\right) \qquad (39)$$
等价地，$n^4$ 的绝对系数为 $8A/(3\pi)$。数值上，重标度残差 $n^4(S_n - S_\infty)(3\pi/8) - \pi n^2/8$，对 $n=3$ 到 $12$ 以超过 20 位精度计算，并不趋向常数而是随 $n$ 线性漂移；纯幂律拟合病态，而允许 $A\ln n$ 自由参数的回归则稳定，返回的 $A$ 系数等于 $8/(3\pi^2)$，精确到 11 位有效数字。表 I 中 $n=7,8,9,10,11,12$ 的数值中已经显现出相同的漂移趋势，它们趋近于式 (39) 和 (40) 预测的直线，无需任何拟合参数（图 2 插图）。注意 $A = \frac{1}{3\pi^2}$，其中 $\zeta$ 是黎曼 ζ 函数。伴随的常数为 [附录 B，式 (B9)]
$$B = \frac{1}{3\pi^2}\left(\ln 2\pi - 2\gamma\right) \qquad (40)$$
其中 $\gamma$ 是欧拉-马歇罗尼常数（故 $B \approx 0.306757$）。附录 A 5 中的回归（令 $B$ 自由）复现式 (40) 到 13 位有效数字，残差由被忽略的 $n^{-6}$ 塔层而非数据精度决定。

---

**图 2.** 自能系数 $S_n$ 与 $n$ 的关系。实心圆为表 I 中的精确值。蓝色曲线为式 (38) 给出的首项 $1/n^2$ 律；青色虚线在式 (39) 的四阶项基础上叠加了角对数 $\ln n$，在 $n\ge 5$ 时即可明显看出吻合度改善（渐近级数预期无法捕捉 $n=3,4$）。红色虚线为圆盘极限。左上：$n=4$ 正方形，取同样的顶点位于轴上的取向；其上边界（蓝色）为单段帐篷形，具有一个高度为 $y=1$ 的内部顶点，对应的闭式表达式为 $S_4 = 4(\ln 2 + \pi/4 + 1/2)$。右下：重标度四阶残差 $n^4(S_n - S_\infty)(3\pi/8) - \pi n^2/8$，$n=7$ 至 $12$，对 $\ln n$ 作图；直线为式 (39) 及其附录 B 的常数，无任何拟合参数。其斜率 $A$ 即为角对数。

**小 $n$ 下的精度。** 式 (39) 的四项截断对其名义有效范围之外的精确多边形能量表征良好。相对于表 I 中的闭式表达式，该展开式的相对误差为：$n=3$ 时约 $0.3\%$，$n=4$ 时约 $0.1\%$，$n=5$ 时约 $0.03\%$，$n=6$ 时仅约 $0.01\%$——因此即便是正六边形也由在 $n\to\infty$ 极限下推导的展开式在优于二百万分之二的精度上复现。对于此处未以闭式处理的更复杂情形，吻合度迅速提升：利用第七节的高精度积分值，截断的式 (39) 在 $n=7$ 时与真实的 $S_7$ 吻合到 $10^{-5}$ 量级，$n=10$ 时到 $10^{-7}$ 量级，$n=15$ 时为 $10^{-9}$ 量级，$n>20$（$n=20$）时降至 $10^{-11}$ 以下，$n=50$ 时为 $10^{-15}$ 量级，残差随被忽略的 $n^{-6}$ 项同步递减。因此，对每一个 $n$，多边形自能已由圆盘值和三个修正项以优于千万分之一的精度给出，而第四、五节中的精确闭式表达式仅在对于面数最少的多边形 $n\le 6$ 时才真正需要。

**结构性起源。** $A$ 的振幅以及 $B$ 的形式均非巧合；两者均来源于 $S_n$ 关于圆盘的二阶变分。将波纹写为 $\theta$ 的傅里叶级数，$n$ 重周期性将其谱限制在谐波 $k\cdot n$ 上，且抛物线齿形给出当 $k\to\infty$ 时：
$$\int_0^{2\pi}|\delta r(\theta)|^2\frac{d\theta}{2\pi} \sim \frac{C}{k^2 n^2} \qquad (41)$$
在边缘上，库仑核为对数奇异的，因此二阶形状导数（海森矩阵）的本征值以 $n$ 增长。取 $\Delta = R\left(\pi - \frac{n}{2}\sin\frac{2\pi}{n}\right)$，振荡二阶贡献对方括号内为
$$\Delta^2\left[\text{Hessian 本征值}\right] \sim \frac{1}{3\pi^2 n^2}\left(\ln n + \text{const.}\right) \qquad (42)$$
利用式 (41) 和 $\sum_{k=1}^\infty 1/k^2 = \pi^2/6$。因此对数项及其 $n$ 的伴随项共享前因子 $A = 1/(3\pi^2)$，而海森矩阵的正规部分——由奇异边缘核的正则化设定的重整化角能量，也是式 (40) 中 $\ln 2\pi$ 和 $2\gamma$ 的来源——设定剩余的 $B$ 部分。这就是为什么 $B$ 是超越数而非代数数。附录 B 的闭式海森矩阵给出 $A = 1/(3\pi^2)$，将式 (42) 的 $\ln n$ 部分与式 (40) 对比，确定通用海森矩阵系数为简洁的有理数
$$\frac{8A}{3\pi} = \frac{8}{9\pi^3} \qquad (43)$$
附录 B 完整执行了此计算：边界海森矩阵以闭式获得 [式 (B7)]，式 (B2)–(B9) 的组装同时给出振幅 $A$ 和式 (40)，无需任何拟合输入；附录 A 5 的回归构成独立验证。等价地，每个顶点为内角 $\pi(1-2/n)$ 的楔形，且完全此类型的对数角项在具有角点的域上的热核、卡西米尔力以及拉普拉斯行列式展开中早已熟知；它们也解释了为什么对 $S_n$ 在 $n$ 的幂次上做朴素理查逊外推是病态的。

首项修正为正：多边形能量位于圆盘值之上，并随 $n\to\infty$ 从上方弛豫至圆盘值，这是一种曲率控制下的趋近。图 2 将表列值与式 (39) 进行了对比。

---

## VII.

**数值验证**

本文提出的所有解析结果都经过了广泛且独立的数值验证。首先，我们将约化的相互作用核与原库仑积分的直接数值计算进行了对比验证，结果在本质上的机器精度范围内完全吻合。例如，公式(6)中的核函数能够以18位有效数字的精度复现暴力计算的双重积分结果。这证实了分析中采用的降维方法完整保留了精确的静电相互作用。

其次，我们将针对各个扇形区域贡献推导出的表达式与对应二维域上的直接数值积分进行了对比。公式(14)中的自扇形公式与在某一扇形区域内对给定被积函数进行直接二维积分的结果吻合到25位有效数字，适用于端点高度非零的扇形。总体而言，对于具有非零端点高度的典型扇形，吻合精度超过20位小数，这为扇形分解及相关积分约化提供了严格的检验。

第三，文中报道的正多边形基准值是通过将高精度双曲正割-双曲正弦（tanh-sinh）积分法应用于约化积分表示计算得出的。在处理可积对角线奇点方面，我们采取了特别谨慎的态度——在数值计算之前先对其进行解析分离。表I中的系数是通过对方程(8)进行tanh-sinh积分获得的，其中对角线奇点在开放端点处被分离，并与闭合形式进行了交叉验证，随后通过与圆形盘极限的渐近方法对比进一步确认。这些独立计算之间的一致性确立了解析公式和数值实现的正确性。

最后，我们对多边数较多时多边形的渐近展开进行了详细考察。通过边界微扰方法解析推导了对圆盘极限的领头阶修正，随后在多种多边形阶数的数值计算中进行了确认。附录B中以闭合形式导出的四阶系数与精确及高精度数值数据进行了核对：对自由振幅的回归分析恢复了11位有效数字的对数系数和13位有效数字的常数项，而纯幂律模型则在数值上表现为奇异的。这种行为反映了多边形的角点即使在近似圆形的区域中仍持续产生影响，从而对渐近理论提供了额外的、高度非平凡的检验。

综合来看，这些验证确立了整个框架的准确性和内在一致性。因此，本文报道的基准能量构成了蒙特卡洛、分子动力学及相关数值研究中用于研究限制在多边形几何形状中的有限库仑系统的可靠参考值。

**VIII. 结论**

我们将均匀带电的正多边形的静电自能约化为一个具有闭合形式相互作用核的精确二维积分，并进一步约化为有限项仿射扇形贡献之和。该约化利用了边界的分段线性结构，得到了一种表示形式，其中每个扇形贡献一个初等积分。构造中的一个关键点是存在与内部顶点非零高度相关的有限截距项。这些项显式地进入每对相互作用中，对于公式的自洽性是必需的；忽略它们会导致在较低多边数时即出现错误的标度关系。

在这一框架内，我们获得了几种低多边数情形的精确闭合表达式，包括三角形、正方形、五边形和六边形。所得表达式在已知结果可得处复现了这些结果，并在统一的框架下将其推广到新的情形。五边形在代数上最为复杂，涉及的四次数域中的中间量，但最终结果简化为包含对数常数的紧凑表达式。尽管代数复杂度随着对称性降低而增加，所有情形均遵循一致的规律：自能随多边形边数的增加单调递减，并平滑地收敛于圆形盘极限。

大边数下的渐近行为通过围绕圆盘几何的边界微扰展开获得。领头阶修正按特定标度律变化，并具有一个与多边形细节（对称性除外）无关的普适系数。高阶修正揭示了更为微妙的结构：除了解析项外，还出现了一个正比于对数函数的非解析贡献。该项源于静电能量二阶变分中由角点引起的奇点，因此即使在近平面极限中也编码了多边形顶点的持续效应。相关的系数——对数振幅和常数——均从圆盘的边界Hessian矩阵以闭合形式获得，并分别以11位和13位有效数字进行了数值确认。

仿射扇形约化并不局限于正多边形。由于只使用了边界函数的分段线性性质，因此任何关于轴对称的多边形都包含在所述框架内。对于上、下边界并非镜像对称的一般简单多边形，横向积分的结果为：

当边界函数取特定值时，该式退化为公式(6)；此时扇形积分涉及与公式(13)–(15)相同的原函数族。非凸多边形通过在每个顶点处将范围切割来处理，使得每一切片都是梯形条带的并集。以下两个进一步的推广仅需很小的代价即可实现：分段多项式而非均匀密度将积分被乘子乘以坐标的多项式，使得横向积分仍为初等积分；两个不相交多边形的相互作用能则是交叉扇形项的简单求和。对于不规则形状，丧失的不是精确性而是紧凑性：顶点坐标不再位于小的数域中，因此附录C所述的坍缩不会发生。

均匀密度问题是引言中提及的等势（导体）问题的自然对偶问题，两者在定量上相关联。由于平衡分布在固定总电荷下使能量最小化，均匀密度能量是导体能量的严格上界，即对每种形状均成立。对于圆盘，相应关系给出特定比例，即均匀密度值的某一百分比；该差值衡量了电荷向边缘弛豫时所释放的能量。对于多边形板，尚未知平衡密度的闭合表达式：平衡密度随到边缘距离的负二分之一次幂发散，并带有额外的角点奇点，其指数取决于内角[25]，因此该问题目前仅通过数值方法求解的开放型积分方程[26–28]。本文结果从两个方面与之相关。首先，表I为这类数值导体能量提供了上界和校准。其次，由于扇形约化可推广到相同划分上的分段多项式密度（上一段），它为导体问题的里兹处理提供了精确积分的矩阵元，其中密度在扇形上展开为多项式，并在固定电荷下最小化能量。本文中发现的角点机制及其对数特征表明，规则多边形板的电容以带有相同类型非解析修正的方式趋近于圆盘的电容。沿此方向将分析推广至等势多边形是本工作的自然延续。

总体而言，这些结果提供了对均匀带电正多边形静电自能的统一且系统可控的描述。它们阐明了边界几何在确定有限尺寸修正中的作用，并确立了平滑的圆形极限如何从离散的多边形域中涌现出来，同时在高阶项中保留了可辨识的角点奇点特征。