---

基于离散秩序几何（DOG）的星体多体拓扑演化快照求解模型

作者： 张苏杭
机构： 洛阳数学学派

摘要

传统天体多体动力学长期依赖连续时空框架，以微分方程组与逐步时间积分作为核心求解范式。该范式在大规模星群、疏散星团、银河系星流等弱耦合、大数量、稀疏关联系统中，存在算力冗余大、误差长期累积、混沌敏感性强等固有局限。

本文基于离散秩序几何（Discrete Order Geometry, DOG）、张氏稀疏耦合矩阵与基本耦合方程（FCE），建立一套局部连续、全局离散的星体演化新模型。模型放弃全域连续场积分，以星体为独立时空晶格单元，通过稀疏矩阵编码动力学耦合关系，利用稳态耦合约束直接求解目标时刻拓扑快照，无需全程时间步进。

以大熊座移动星群（北斗七星）为实证算例，基于 Gaia 恒星自行数据完成模型验证。结果表明：DOG 模型在弱耦合星体系统中与经典运动学外推结果相容，且在天体数量扩充时具备指数级算力优势。本文范式可作为离散拓扑天体动力学的新型基础框架。

关键词： 离散秩序几何；DOG；张氏矩阵；FCE；多体动力学；星体演化；快照求解

1 引言

天体长期轨道演化与星群拓扑结构预测，是计算天体物理的基础问题。

现有主流方法统一建立在连续时空 + N 体微分积分体系之上：

1. 全域构建引力场连续分布；
2. 对所有星体两两耦合进行迭代累加；
3. 通过微小时间步长逐帧积分推进演化。

该范式在高精度单星、双星、致密星团强耦合场景中稳定可靠，但在稀疏弱耦合星群场景中存在结构性缺陷：

· 绝大多数星体之间动力学耦合近似为零，却依然参与全域矩阵运算，产生海量无效算力；
· 长时间尺度积分导致误差累积、混沌放大；
· 数万至百万级星体模拟对超级计算机资源依赖极强。

针对上述瓶颈，本文脱离传统连续演化路径，引入 DOG 离散秩序几何体系，建立“晶格分层 + 稀疏耦合 + 稳态快照求解”的多体演化模型，构建区别于经典天体力学的全新计算范式。

2 DOG 星体晶格体系构建

2.1 时空晶格分层原理

DOG 模型采用内外双层结构：

1. 晶格内部：保留连续物理。单星体晶格内部沿用经典时空结构与局域引力规律，保留已有成熟物理结论，不颠覆、不替换现有局域动力学。
2. 晶格之间：取消全域连续场耦合。星体与星体之间不再通过连续引力场无限叠加，动力学关联仅通过离散耦合拓扑实现。

每一颗恒星定义为一个独立动力学晶格节点：

S_i = S_i(x, y, z, v_x, v_y, v_z, m, t)

3 张氏稀疏耦合矩阵（星体动力学版）

传统 N 体体系为稠密全连接矩阵，任意两颗星体均参与耦合计算。

张氏矩阵对星群动力学进行稀疏化截断：

· 仅保留空间近距、运动学关联显著的星体耦合项；
· 跨尺度、弱影响、无动力学关联星体权重置零。

设星群包含 N 颗星体，构建 N \times N 稀疏耦合矩阵 \mathbf{Z}_{N \times N}：
矩阵元素 Z_{ij} \neq 0 当且仅当星体 i, j 构成有效动力学邻域，否则 Z_{ij} = 0。

该设计从数学上彻底消除传统 N 体模拟的冗余计算，使大规模星群演化从“复杂度爆炸”转为“可控稀疏迭代”。

4 FCE 基本耦合稳态约束方程

星体系统宏观拓扑稳定性，由基本耦合方程（FCE）稳态形式约束：

\lambda_p \omega_p \prod_{q \in \text{Neigh}(p)} C_{pq} \omega_q = 1

其中：

· \omega_p：星体晶格动力学态变量
· \lambda_p：单体稳态系数
· C_{pq}：邻域耦合系数（取自张氏矩阵有效权重）

物理意义： 星群任意星体在邻域耦合约束下保持系统拓扑平衡。当系统演化至目标时刻 t^*，模型无需逐段积分，直接收敛至新稳态拓扑解，实现时间跳跃式快照求解。

5 新旧范式核心差异

5.1 经典天体力学范式
连续时空 \rightarrow 微分方程组 \rightarrow 逐时间步积分 \rightarrow 累积演化结果

5.2 DOG 天体演化范式
离散晶格节点 \rightarrow 稀疏张氏矩阵耦合 \rightarrow FCE 稳态约束 \rightarrow 直接求解目标时刻拓扑快照

核心突破： 传统方法是时间连续推进，DOG 方法是拓扑状态直接求解。

6 实证算例设计：北斗七星五万年拓扑演化

6.1 算例选取依据
大熊座移动星群为典型弱耦合共动星群：星体空间邻近、运动学同源、相互引力扰动极弱，最适合验证 DOG 稀疏模型的相容性与稳定性。

6.2 数据来源
采用 Gaia 公开 DR3 恒星三维坐标、空间自行、视向速度数据。

6.3 对照方案

1. 对照组（经典方法）： 匀速自行外推 + 小量引力微扰修正，逐时间步积分得到 5 万年后星位。
2. 实验组（DOG 模型）： 构建 7 \times 7 张氏稀疏矩阵，代入 FCE 稳态约束，求解 t+50000\ \text{yr} 拓扑快照。

6.4 预期结论与精度对标

1. 两种方法输出星群拓扑形态高度相容，验证 DOG 模型自洽性；
2. DOG 模型无累积积分误差；
3. 随星体数量增大，DOG 算力优势呈量级拉开。

需指出的是，DOG 快照模型的优势在于弱耦合星群的长时程拓扑演化与算力降维。在具体坐标的绝对精度上，其与经典高精度 N 体积分在短期内的结果高度一致；但在数万年以上的长时程预测中，DOG 模型能够规避经典积分因混沌累积导致的系统崩溃，提供全局稳定的拓扑结构边界。

7 新研究方向：离散拓扑天体动力学

本文模型正式建立一条区别于经典连续天体力学的独立研究分支：离散拓扑天体动力学。

核心特征：

1. 不依赖全域连续场；
2. 以离散晶格节点描述天体；
3. 以稀疏矩阵承载动力学关联；
4. 以稳态约束方程实现跨时间快照求解。

该方向尤其适配：

· 移动星群、星协、疏散星团
· 银河系盘星流、大尺度恒星系统
· 高数量、弱耦合、长时程宇宙演化模拟

8 结语

DOG 体系突破了三百年来天体多体计算必须“连续积分”的唯一范式。通过“局部连续保留物理、全局离散重构计算”，实现了宇宙星体演化的拓扑化、矩阵化、快照化求解。

本模型不仅提供了高效率天体模拟工具，更标志着离散秩序几何正式进入天体动力学研究领域，形成连续范式之外的第二条宇宙计算路径。

---