奶牛的高级多极子

# 奶牛的高阶多极矩
来源：https://arxiv.org/html/2504.00506
本杰明·V·莱曼\ affiliation：美国马萨诸塞州剑桥市马萨诸塞理工学院理论物理中心——莱因韦伯研究所，邮编 02139

2026 年 4 月 1 日

###### 摘要

球形奶牛近似法在文献中被广泛使用，但很少有正当的理由。本文提出几种将球形奶牛近似法扩展至完整多极矩展开的方案，其中球形奶牛仅作为首项。这使得我们可以计算超越球对称性的奶牛势场和相互作用，同时提供了一种定义奶牛自身几何形状高阶多极矩的方案。这对于处理被球对称性抑制的物理过程尤为重要，例如旋转奶牛因引力波发射而导致的自转减慢。本文对一个基准奶牛计算了多极矩系数，并阐述了多极矩奶牛在若干重要问题中的应用性。

## I 引言

球形奶牛近似法（SCA）的各种形式是大量子领域关键成果的基础。谷歌学术搜索显示，球形奶牛出现在超过 1000 篇已发表的论文中\[1 (https://arxiv.org/html/2504.00506#bib.bib1)\]。SCA 如此普遍，以至于它已成为我们定期向大一本科生教授的少数近似方法之一。鉴于 SCA 在物理学中的重要性，文献中从未出现对其有效性的定量评估，这令人震惊和担忧。

此外，SCA 作为一种近似技术是贫乏的，因为其通常形式无法系统改进。在大多数近似方案中，虽然人们可能只计算领头阶，但原则上可以考虑次领头阶或更高阶的贡献，以改进结果或评估近似的收敛性。但球形奶牛只是一个球体，没有额外的参数空间来改进其对真实奶牛的匹配。这应该让我们都感到相当惭愧。

然而，这并非 SCA 结构本身固有的。如果存在合适的高阶修正方案，球形奶牛可以被视为一个表述良好的展开式的领头阶项。确实，有一种标准的近似技术非常适合此目的：球谐多极矩展开。多极矩展开的第一项——单极矩——是球对称的，更高阶项则编码了越来越小尺度上对球对称性的偏离。如果球形奶牛可以被重新解释为多极矩展开的单极矩，那么 SCA 就成为一个系统可改进近似的第一项，从而可以在各种情况下定量检验其可靠性。

在本文中，我介绍了几种不同的方案来实现这种扩展，将 SCA 推广到包含奶牛的高阶多极矩。我直接比较了高阶多极矩与单极矩（即球形奶牛）的贡献，首次在各种情况下对 SCA 的有效性进行了真正检验。我利用这一框架确定了 SCA 显然不足的情况，即使在数量级上，也必须包含偶极子或更高阶多极矩的贡献才能获得物理上现实的结果。正如我将解释的，即使是一些最经典的奶牛问题，如奶牛倾倒问题，其主要贡献也来自于高阶多极矩。

由于 SCA 本身有几种不同的形式，扩展到高阶多极矩也会因用例而异。本文重点讨论两种情况。在第 II 节 (https://arxiv.org/html/2504.00506#S2)，我研究了在奶牛源势场中包含高阶多极矩的问题，并讨论了这对被球对称性抑制的过程（包括引力辐射）的一些重要影响。在第 III 节 (https://arxiv.org/html/2504.00506#S3)，我定义了奶牛自身几何形状的多极矩展开，这允许对 SCA 进行更一般的扩展。后一种情况对于奶牛刚体动力学问题（例如奶牛倾倒）尤为重要，这些问题将在第 IV 节 (https://arxiv.org/html/2504.00506#S4) 中处理。我在第 V 节 (https://arxiv.org/html/2504.00506#S5) 中讨论了这些结果的含义并得出结论。

在本研究中，对于数值计算，我使用了来自 libigltutorial 数据的基准奶牛 \[2 (https://arxiv.org/html/2504.00506#bib.bib2),3 (https://arxiv.org/html/2504.00506#bib.bib3)\]。我用 \( \mathcal{C} \) 表示奶牛的内部，用 \( \partial\mathcal{C} \) 表示奶牛的二维表面。我遵循 IUPAC 关于乘数前缀的约定 \[4 (https://arxiv.org/html/2504.00506#bib.bib4)\]，因此，例如，\( \ell=256 \) 的项被称为六百五十四极矩的分量。

## II 奶牛势场的高阶多极矩

让我们首先考虑球形奶牛在引力问题背景下的最简定义。当然，同样的方法也适用于研究带电奶牛的静电势。然而，相对于带电公牛，带电奶牛在自然界中很少被观察到 \[5 (https://arxiv.org/html/2504.00506#bib.bib5),6 (https://arxiv.org/html/2504.00506#bib.bib6)\]。（尽管在实验环境中仍需谨慎，参考文献 \[7 (https://arxiv.org/html/2504.00506#bib.bib7)\]。）

当计算奶牛外部一点 \( \mathbf{x} \) 的引力势 \( \phi(\mathbf{x}) \) 时，拉普拉斯方程的解可以通过不规则固体球谐函数展开为球谐函数 \( Y_{\ell}^{m}(\hat{\mathbf{x}}) \)，得到多极矩展开

\[
\phi(\mathbf{x}) = -G \sum_{\ell=0}^{\infty} \frac{\left[\frac{4\pi}{2\ell+1}\right]^{1/2}}{\|\mathbf{x}\|^{\ell+1}} \sum_{m=-\ell}^{\ell} (-1)^{m} Y_{\ell}^{-m}(\hat{\mathbf{x}}) Q_{\ell}^{m}, \tag{1}
\]

其中球谐多极矩 \( Q_{\ell}^{m} \) 定义为

\[
Q_{\ell}^{m} \equiv \int_{\mathcal{C}} \mathrm{d}^{3}\mathbf{x} \, \rho(\mathbf{x}) \, \|\mathbf{x}\|^{\ell} Y_{\ell}^{m}(\hat{\mathbf{x}}), \tag{2}
\]

给定质量密度 \( \rho(\mathbf{x}) \)。（详细推导见参考文献 \[8 (https://arxiv.org/html/2504.00506#bib.bib8)\]。）我将球谐函数归一化为

\[
Y_{\ell}^{m}(\hat{\mathbf{x}}) = \sqrt{\frac{(2\ell+1)!}{4\pi} \frac{(\ell-m)!}{(\ell+m)!}} \, P_{\ell}^{m}(\cos\theta) \, e^{im\theta}, \tag{3}
\]

其中 \( P_{\ell}^{m} \) 表示关联勒让德多项式，\( \hat{\mathbf{x}} \) 是 \( \mathbf{x} \) 方向的单位向量。我将假设 \( \rho \) 在奶牛上是常数，并在“奶牛坐标系”下工作，其中 \( x \) 轴与奶牛的前进方向对齐，\( y \) 轴垂直于地面，正 \( z \) 轴指向奶牛的右侧。这与基准奶牛的坐标轴一致。我还保留了基准奶牛的尺度，因此长度单位对应于一个尺寸为 \( (1.044, 0.6397, 0.3403) \) 的奶牛包围盒，分别对应 \( x \), \( y \), 和 \( z \) 方向。我称这些单位为“基准单位”，本节后面将提供到国际单位制的转换。基准奶牛的前几个多极矩在这些单位下如表 1 (https://arxiv.org/html/2504.00506#S2.T1) 所示。单极矩对应于球形奶牛，因此更高的多极矩系数表示了对 SCA 的修正大小。

表 1：被视为质量分布的奶牛 \( \mathcal{C} \) 的领先多极矩，以基准单位计算。在质心坐标中，偶极矩为零，因此省略了这些分量。此外，由于 \( Q_{\ell}^{-m} = (-1)^{m} (Q_{\ell}^{m})^{*} \)，仅显示正 \( m \) 的系数。

我们现在可以使用这些系数来研究奶牛势场中的引力现象。有几个潜在应用。例如，四极矩编码了重型奶牛地球赤道隆起的等效效应。此外，具有非零四极矩的奶牛在非均匀引力场中会受到力矩作用，因此从高海拔处落下的奶牛在地球引力场中会倾向于特定的方向。因此，四极矩对于回答下落的奶牛是否倾向于脚着地的问题至关重要。同样，四极矩在轨道上奶牛的潮汐锁定中也起着关键作用。111 对轨道上猪的类似计算也适用，但只有当猪飞起来时才预期能观察到。在这两种情况下，SCA 都显然失效，而且这两种情况都不适合进行实验研究。

然而，目前让我们关注四极矩的另一个应用：计算引力辐射的发射率。真空中自由旋转的奶牛会因向引力波损失能量和角动量而减慢。我们可以使用引力辐射的四极矩公式计算自转减慢的速率。当用球谐四极矩分量表示时，这可读为 \[9 (https://arxiv.org/html/2504.00506#bib.bib9)\]

\[
\dot{E}_{\mathrm{quad}} = \frac{3G}{8\pi c^{5}} \sum_{m=-2}^{2} \left\langle \left\| \dddot{Q}_{2}^{m} \right\|^{2} \right\rangle. \tag{4}
\]

原则上，可以从表 1 (https://arxiv.org/html/2504.00506#S2.T1) 中获取球谐多极矩系数，在旋转下进行变换，并使用这种变换行为来评估绕给定轴旋转时的 \( \dddot{Q}_{2}^{m} \)。然而，为此目的使用笛卡尔多极矩系数更简单，因为笛卡尔四极矩系数的形式是一个对称张量，分量为 \( Q_{ij}^{\mathrm{C}} = \int \mathrm{d}^{3}\mathbf{x} \, \rho(\mathbf{x}) \, (3x_{i}x_{j} - r^{2}\delta_{ij}) \)。基准奶牛在奶牛坐标和基准单位下（其中 \( \rho(\mathbf{x})=1 \)）的笛卡尔四极矩张量为

\[
Q^{\mathrm{C}} = \begin{pmatrix} 4.23004 & \phantom{-}0.83531 & \phantom{-}0.00704 \\ 0.83531 & -1.64753 & \phantom{-}0.00039 \\ 0.00704 & \phantom{-}0.00039 & -2.58252 \end{pmatrix} \times 10^{-3}, \tag{5}
\]

辐射功率由 \( Q_{ij}^{\mathrm{C}} \) 给出

\[
\dot{E}_{\mathrm{quad}} = \frac{G}{45c^{5}} \left\langle \dddot{Q}_{ij}^{\mathrm{C}} \dddot{Q}^{\mathrm{C}\,ij} \right\rangle. \tag{6}
\]

现在，假设奶牛绕轴 \( \hat{\mathbf{a}} \) 以角频率 \( \omega \) 旋转。如果 \( R_{ij}(\theta) \) 表示绕 \( \hat{\mathbf{a}} \) 旋转角度 \( \theta \) 的旋转矩阵，那么在任何时刻 \( t \)，我们可以将四极矩张量写为

\[
Q_{ij}^{\mathrm{C}}(t) = R_{ik}(\omega t) R_{jl}(\omega t) Q_{kl}^{\mathrm{C}}(0), \tag{7}
\]

因此时间的三阶导数为

\[
\dddot{Q}_{ij}^{\mathrm{C}}(t) = Q_{k\ell}^{\mathrm{C}}(0) \, \partial_{t}^{3} \left[ R_{ik}(\omega t) \, R_{jl}(\omega t) \right]. \tag{8}
\]

作为一个简单的例子，考虑绕 \( y \) 轴的恒定旋转，就像一头戴着错位旱冰鞋的奶牛可能经历的那样。这里的旋转矩阵是

\[
R(\theta) = \begin{pmatrix} \cos\theta & 0 & \sin\theta \\ 0 & 1 & 0 \\ -\sin\theta & 0 & \cos\theta \end{pmatrix}. \tag{9}
\]

将其代入公式 8 (https://arxiv.org/html/2504.00506#S2.E8) 并在整个周期上平均结果，得到 \( \left\langle \dddot{Q}_{ij}^{\mathrm{C}} \dddot{Q}^{\mathrm{C}\,ij} \right\rangle \approx 0.00149\omega^{6} \)。

我们现在可以恢复这个结果的物理量纲，它具有 \( M^{2} L^{4} \omega^{6} \) 的单位，其中 \( M \) 和 \( L \) 是基准奶牛的质量和长度单位。一头奶牛的平均长度约为 \( 2.5\,\mathrm{m} \)，基准单位中的长度是 1.044，所以我们可以说 \( L = 2.39\,\mathrm{m} \)。为了确定质量单位，回想一下我们设定了 \( \rho(\mathbf{x})=1 \)，而奶牛的实际密度约为 \( 1\,\mathrm{g} \)。设定 \( M/L^{3} = 1\,\mathrm{g} \) 得到 \( M = 13\,652\,\mathrm{kg} \)。请注意，这不是奶牛本身的质量，而只是我们选择将特定物理值（长度和密度）设为 1 所隐含的质量单位。（奶牛在基准单位中的体积仅为 0.054，对应于约 \( 740\,\mathrm{kg} \) 的物理质量——对于一头奶牛来说相当典型。）代入这些量，在物理单位中，我们得到 \( \left\langle \dddot{Q}_{ij}^{\mathrm{C}} \dddot{Q}^{\mathrm{C}\,ij} \right\rangle \approx 9\times 10^{12}\,\mathrm{kg} \times \omega^{6} \)，或者

\[
\dot{E}_{\mathrm{quad}} \approx 5.5 \times 10^{-41}\,\mathrm{erg} \times \left( \frac{\omega}{1\,\mathrm{Hz}} \right)^{6}. \tag{10}
\]

要理解这个能量损失率的相关性，我们必须首先计算旋转能量。回想一下，在角频率 \( \omega \) 下的旋转能量由 \( E = \frac{1}{2} I \omega^{2} \) 给出，其中 \( I \) 是绕旋转轴的转动惯量。对于绕任意轴的旋转，这推广为一个矢量表达式：转动惯量被惯性张量取代，其分量为 \( I_{ij} = \int \mathrm{d}^{3}\mathbf{x} \, \rho(\mathbf{x}) \, \left( r^{2} \delta_{ij} - x_{i} x_{j} \right) \)。（请注意，这与公式 5 (https://arxiv.org/html/2504.00506#S2.E5) 中的笛卡尔四极矩张量相似但不相同：笛卡尔多极矩张量源于势场的泰勒展开，这在展开的每一阶引入了不同的常数因子乘以笛卡尔坐标。）基准单位下的惯性张量为

\[
I = \begin{pmatrix} \phantom{-}7.95079 & -2.78437 & -0.02348 \\ -2.78437 & \phantom{-}27.5426 & -0.00130 \\ -0.02348 & -0.00130 & \phantom{-}30.6593 \end{pmatrix} \times 10^{-4}, \tag{11}
\]

因此我们可以轻松确定绕 \( y \) 轴旋转相关的动能。（我们还可以看到，在固定 \( \omega \) 的情况下，能量最大的自转将涉及绕 \( z \) 轴旋转，即翻跟头下落。）因此，我们可以获得旋转减慢的时间尺度（忽略角动量辐射的细节）：

\[
\frac{E}{\dot{E}_{\mathrm{quad}}} \approx 1.9 \times 10^{49}\,\mathrm{s} \left( \frac{\omega}{1\,\mathrm{Hz}} \right)^{-4}. \tag{12}
\]