文章提出了一种计算非均匀、非球形行星引力势的数值方法。该方法通过将物理行星映射到一个几何球形的参考行星上，将复杂的几何问题转化为参考域中的泊松方程求解。文章详细介绍了前向计算和伴随理论，用于计算引力势及其灵敏度核。以火卫一为例，展示了该方法在均匀和非均匀密度模型下的应用。结果表明，该方法能够高效准确地计算复杂形状天体的引力势，为行星内部结构研究提供了有力工具。
精确计算牛顿引力势在行星科学、地球物理学和天体物理学中具有重要意义。对于几何球形的天体，引力势可以通过球谐展开来表示。然而，大多数天体并非完美的球形，这使得传统的球谐展开方法难以直接应用。
文章采用参考域方法，将物理行星映射到一个几何球形的参考行星上。这种方法的关键在于将几何复杂性转化为泊松方程的复杂性。通过定义一个微分同胚映射，可以将物理行星的内部和外部边界转换为参考行星中的同心球壳。
在参考域中，泊松方程被重新表述为一个包含变形张量的弱形式。这种方法允许在参考域中使用球谐展开，同时保持对物理行星复杂形状的精确描述。
数值实现采用伪谱/谱元离散化方法。球谐变换用于处理角度部分，而径向部分则采用谱元方法离散。这种方法结合了谱方法的高精度和有限元方法的几何灵活性。
为了避免大型稠密矩阵的计算，文章采用了矩阵自由方法和预条件共轭梯度求解器。这种方法显著降低了计算成本，使得在单核处理器上也能高效求解。
文章推导了引力势对参考密度和映射的扰动响应。通过线性化处理，可以将扰动问题表述为与原问题相同形式的泊松方程，仅力项发生变化。这使得可以利用已有的求解器框架高效计算灵敏度。
对于任意势能泛函，文章建立了伴随方程来计算灵敏度核。这些核描述了目标泛函对密度、映射和变形梯度的敏感性。伴随方程同样具有泊松方程的形式，可以利用前向求解器进行计算。
文章还导出了一个重要的参考灵敏度核恒等式，该恒等式必须在参考公式中使用不同于物理域的计算域时得到满足。
以火卫一为例，展示了该方法在均匀和非均匀密度模型下的应用。计算结果表明，即使对于形状高度不规则的天体，该方法也能提供精确的引力势计算。灵敏度核的计算揭示了引力势与内部密度分布之间的复杂关系。
这一系列研究表明，参考域方法为复杂形状天体的引力场计算提供了有效途径，未来可进一步应用于行星内部结构反演等问题。