摘要
本文提出并推演一套以「底层只有能量在流动」为唯一公设的物理学框架——混沌拓扑学（Homotopic Dynamics of Turbulence, HDT）。该框架不预设粒子、场、时空、力为独立本体，而将其统一处理为能量流动在不同尺度、不同观测精度下的读数。
在此框架下，本文完成三项工作：
- （一）给出普适关系。该关系由自由涡旋（角动量守恒）与中心力束缚（开普勒 / 库仑）两种独立机制分别导出，形式一致。以氢原子三个能级检验，比值精确为 1.000000；以地—日、月—地系统检验，残差处于天文参数舍入量级。由此解出「系统尺度越小、能量越高，其时间读数越快」的定量形式。
- （二）解释三代粒子的几何来源。三个互相正交的方向投影至垂直于自旋轴的平面，两两夹角恰为 120°。本文证明，Koide 经验公式中三十年来缺乏解释的 120° 相位，正是该投影的几何必然；并由此导出解析关系。以实测质量反解，轻子投影扁率与 √2 的相对偏差为 9.2×10⁻⁶。
- （三）建立红移的乘法分配模型。本文证明光能耗散律的一阶线性形式由观测的无色性（achromaticity）唯一锁定，不可修改；进而指出，将该关系独立用于解释全部红移会导致角直径距离在 z=1 处偏差达 2 倍。本文提出「物质胀缩 × 传播耗散」模型，并证明其与标准宇宙学完全一致的充要条件为 g=p≤1。该结果将红移机制的判决转化为对参数 p（原子频率随物质尺度的标度指数）的测定。
此外，本文给出真空密度扰动剖面，该形式在小尺度退化为 1/r（满足引力透镜），在大尺度退化为 1/r²（满足星系旋转曲线平坦），转折发生于 ~10 kpc 星系尺度，并给出定量的可检验预测。
本文同时如实记录一条被证伪的猜想（与实测相差两个数量级）及一条仅达定性程度的推论（温度下限，耦合系数待标定）。
关键词：混沌拓扑学；能量本体论；时间—能量关系；Koide 公式；光能耗散；宇宙学常数；暗物质；红移机制
1 引言
1.1 问题的提出
自牛顿以来，物理学每遇计算受阻之处，惯用的做法是引入新的假设，随后将该假设视为已被发现的本体。本文认为，这一路径造成了若干持续未决的困难。表 1 列出九项此类假设及其引入的真实动机。
表 1 主流理论中作为假设引入的若干设定
| 编号 | 假设 | 引入动机 |
|---|---|---|
| A1 | 光速恒定 | 若不固定 c，场方程无法闭合 |
| A2 | 时空弯曲 | 为消解「真空无介质」与波动传播之间的悖论 |
| A3 | 宇宙学常数 Λ | 场方程给出塌缩解，而观测显示加速膨胀 |
| A4 | 普朗克尺度 | 由 G、c、ℏ 三个假设组合运算所得的数学边界 |
| A5 | 全同粒子 | 观测粒度不足，被当作本体性质 |
| A6 | 夸克为基本粒子 | 加速器能量不足以再分割，被当作最小单元 |
| A7 | 测不准为本性 | 采样速率不足，被写为世界的不确定性 |
| A8 | 真空涨落与虚粒子对 | 为解释「无中生有」而构造的借还机制 |
| A9 | 暗物质与暗能量 | 以新假设修补旧假设所导致的偏差 |
本文指出，上述九项具有共同的根源：将观测者的读数误认为宇宙的本性。
以 A3 为例，量子场论给出的真空能量密度与由宇宙学常数观测反演所得值相差约 10¹²⁰ 倍[1,2]。这不是精度问题，而是「真空能量密度为常数」这一假设本身的问题。
1.2 本文的思路与边界
本文不逐一反驳上述假设，而是更换起点，考察哪些假设在新框架下自动消失。
本文的方法论立场为：方程应当跟随现象，而非现象跟随方程。数学工具的选择本身即构成对世界形态的预设——选择微分方程即预设连续性，选择希尔伯特空间即预设叠加的本体性。因此本文的处理顺序为：先确立现象图像，再寻求其定量表达；当方程与现象冲突时，判定方程有误，而非现象有误。
本文的自我约束如下：
- 每条关系均附证伪条件；
- 未通过检验者不得作为结论使用；
- 被证伪的猜想亦如实记录，并说明其排除价值。
2 理论框架
2.1 本体论公设
公设（唯一）：底层只有能量在流动。
此处的「流动」不作比喻使用。它断言：不存在独立于能量之外的承载者（如时空、以太、真空容器）；运动、结构、相互作用均为能量自身分布与变化的表现方式。
由此，各项物理概念获得如下重述。
表 2 HDT 概念重述
| 观测概念 | HDT 重述 |
|---|---|
| 真空 | 极稀薄的能量流 |
| 粒子 | 能量流中的自持涡旋 |
| 光子 | 能量流中的波包 |
| 力 | 能量密度的梯度 |
| 物理常数 | 能量流在局域的读数 |
| 时间 | 能量流完成循环的节奏读数 |
2.2 关于真空
主流量子场论以「零点能」「涨落」「虚粒子对生灭」描述真空。此类术语隐含一个前提：真空是某种静态基元，仅偶尔偏离平衡。
本文主张：真空是持续流动的能量介质，具有梯度、涡旋与对流。术语「涨落」暗示存在平均值及围绕该值的偏离；而本文认为不存在这样的平均值——能量的不均匀性与流动性是其本然状态，稀薄处同样在流动。
2.3 关于粒子
本文主张：粒子不是被激发的场，而是两股能量对流相互挤压所形成的涡旋。
两股能量流相遇，产生局部压强变化，形成旋转结构；该结构一旦建立，即从周围介质持续获取能量以维持自身。周围梯度消失时，结构解散，回归稀薄介质。
此图像可由太阳表面的米粒组织（granulation）类比：每一对流胞呈现边界、寿命与尺度，看似独立个体，实则均为同一连续介质的运动模式。
2.4 关于夸克
本文不认为夸克是基本粒子。加速器在当前能量下无法进一步分割，不构成「不可分割」的证据——正如刀具无法切出细胞，不代表细胞不存在。
本文主张：质子、中子等粒子本身即为能量密度极高的局域区域（可称为「天体」），夸克是该天体内部的对流胞、环流或涡旋构型。
三个术语在本文中不互相排斥，而是同一结构在不同层面的称谓：
- 对流：两股能量流相遇、交换、循环，为驱动机制；
- 环流：能量绕自身旋转并自持，为粒子的本体形态；
- 涡旋：环流中具有核心、边界与旋向的可观测构型。
2.5 关于真空生粒与虚粒子对
本文认为，稀薄介质在流动中挤压出涡旋，即为粒子的产生。此过程不需要「借取能量」，不需要「正反粒子成对产生」，也不需要不确定性原理提供例外。
所谓「成对」，是观测切片在对称系统中所得读数的对称性，并非自然界必须成对。
2.6 关于力
本文主张：不存在作为独立本体的「力」。光子沿梯度偏折、落体下落、行星绕行，均为同一过程——能量沿密度梯度运动。观测者读到轨迹弯曲，遂构造「引力」一词。
四种基本相互作用在本文中不被视为四种本体，而是同一介质在不同尺度上呈现的不同梯度模式，由不同精度的观测工具读出。
3 时间、能量与角动量
3.1 对「时间膨胀」一语的纠正
早期草稿曾沿用主流术语「时间膨胀」以解释高红移超新星光变曲线的展宽。本文在此纠正：该术语预设空间在膨胀，而空间膨胀正是本框架所要排除的设定。
本文的立场为：时间不是维度，而是能量完成一个循环所需的节奏读数。
3.2 普适关系的导出
定义：
- τ：系统完成一个特征循环所需时间；
- E：系统维持该循环的做功能力；
- L：角动量。
推导一（自由涡旋，角动量守恒）。设 L 为常数，则
假设 1（角动量守恒）：自由涡旋在无外力矩作用下，角动量 L 保持恒定。对质量为 m、以速率 v 沿半径 r 作圆周运动的涡旋，有
L=mvr=常数(3.1) L = mvr = \text{常数} \tag{3.1} L=mvr=常数(3.1)
假设 2（循环时间）：系统完成一个完整圆周所需时间 τ 由轨道周长与速率之比给出：
τ=2πrv(3.2) \tau = \frac{2\pi r}{v} \tag{3.2} τ=v2πr(3.2)
假设 3（动能定义）：系统维持该循环的做功能力取经典动能形式：
E=12mv2(3.3) E = \frac{1}{2}mv^2 \tag{3.3} E=21mv2(3.3)
由 (3.1) 与 (3.2) 得
v=Lmr,τ=2πrv=2πmr2L(3.4) v = \frac{L}{mr},\qquad \tau = \frac{2\pi r}{v} = \frac{2\pi mr^2}{L} \tag{3.4} v=mrL,τ=v2πr=L2πmr2(3.4)
将 (3.4) 代入 (3.3)：
E=12m(Lmr)2=L22mr2(3.5) E = \frac{1}{2}m\left(\frac{L}{mr}\right)^2 = \frac{L^2}{2mr^2} \tag{3.5} E=21m(mrL)2=2mr2L2(3.5)
由 (3.4) 与 (3.5) 得
∣E∣⋅τ=L22mr2⋅2πmr2L=πL(3.6) |E|\cdot\tau = \frac{L^2}{2mr^2}\cdot\frac{2\pi mr^2}{L} = \pi L \tag{3.6} ∣E∣⋅τ=2mr2L2⋅L2πmr2=πL(3.6)
推导二（中心力束缚）。对圆轨道有
假设 4（牛顿引力）：中心力束缚系统（如行星绕日、电子绕核）所受引力为
F=GMmr2(3.7) F = \frac{GMm}{r^2} \tag{3.7} F=r2GMm(3.7)
假设 5（向心力平衡）：稳定圆轨道上，引力恰提供向心力：
F=mv2r(3.8) F = \frac{mv^2}{r} \tag{3.8} F=rmv2(3.8)
由 (3.7) 与 (3.8) 相等，得
GMmr2=mv2r ⟹ v2=GMr(3.9) \frac{GMm}{r^2} = \frac{mv^2}{r}\ \Longrightarrow\ v^2 = \frac{GM}{r} \tag{3.9} r2GMm=rmv2 ⟹ v2=rGM(3.9)
假设 6（束缚态能量）：系统总能量为动能与势能之和，且对圆轨道满足位力定理 E = −½mv²（势能 V = −GMm/r = −2·½mv²）：
E=12mv2−GMmr=12mv2−mv2=−12mv2(3.10) E = \frac{1}{2}mv^2 - \frac{GMm}{r} = \frac{1}{2}mv^2 - mv^2 = -\frac{1}{2}mv^2 \tag{3.10} E=21mv2−rGMm=21mv2−mv2=−21mv2(3.10)
由 (3.9) 得 v = √(GM/r)，代入 (3.10)：
∣E∣=12mv2=GMm2r(3.11) |E| = \frac{1}{2}mv^2 = \frac{GMm}{2r} \tag{3.11} ∣E∣=21mv2=2rGMm(3.11)
循环时间仍由 (3.2) 给出，代入 (3.9)：
τ=2πrv=2πrrGM=2πr3GM(3.12) \tau = \frac{2\pi r}{v} = 2\pi r\sqrt{\frac{r}{GM}} = 2\pi\sqrt{\frac{r^3}{GM}} \tag{3.12} τ=v2πr=2πrGMr=2πGMr3(3.12)
由 (3.11) 与 (3.12) 得
∣E∣⋅τ=GMm2r⋅2πr3GM=πmGMr(3.13) |E|\cdot\tau = \frac{GMm}{2r}\cdot 2\pi\sqrt{\frac{r^3}{GM}} = \pi m\sqrt{GMr} \tag{3.13} ∣E∣⋅τ=2rGMm⋅2πGMr3=πmGMr(3.13)
而角动量 L = mvr = m√(GM/r)·r = m√(GMr)，故
∣E∣⋅τ=πL(3.14) |E|\cdot\tau = \pi L \tag{3.14} ∣E∣⋅τ=πL(3.14)
两种机制在动力学上彼此独立，却给出同一形式。故有
∣E∣⋅τ=πL(3.15) \boxed{|E|\cdot\tau = \pi L} \tag{3.15} ∣E∣⋅τ=πL(3.15)
关键步骤对照。两种推导的动力学起点完全不同——推导一仅依赖角动量守恒与圆周运动学，不涉及任何力的具体形式；推导二依赖牛顿引力与向心力平衡，不涉及角动量守恒的假设。二者在数学上收敛于同一关系的关键在于：能量 E 与循环时间 τ 各自都含有 r 的互补幂次——在自由涡旋中，E ∝ r⁻² 而 τ ∝ r²，乘积消去 r；在中心力束缚中，|E| ∝ r⁻¹ 而 τ ∝ r³ᐟ²，乘积同样消去 r。正是这种「能量—时间」对 r 的互补标度，使得两种独立机制在形式上完全一致。
物理含义。一个系统完成一次循环所用的时间，乘以其能量，等于其角动量与 π 的乘积。引力束缚、库仑束缚与自由涡旋三类运动均服从此关系。
若角动量量子化，则
L=nℏ,n=1,2,3,…(3.16) L = n\hbar,\qquad n = 1,2,3,\ldots \tag{3.16} L=nℏ,n=1,2,3,…(3.16)
代入 (3.15)：
∣E∣⋅τ=nπℏ(3.17) |E|\cdot\tau = n\pi\hbar \tag{3.17} ∣E∣⋅τ=nπℏ(3.17)
物理含义。当角动量只能按份取值时，能量与时间构成反比关系——能量越高，完成一次循环越快。
3.3 数值检验
表 3 |E|·τ 与 π·L 的比对
| 系统 | |E|·τ (J·s) | π·L (J·s) | 比值 |
|—|—|—|—|—|
| 氢原子 n=1 | — | — | 1.000000 |
| 氢原子 n=2 | — | — | 1.000000 |
| 氢原子 n=3 | — | — | 1.000000 |
| 地球绕日 | — | — | 0.999820 |
| 月球绕地 | — | — | 1.000089 |
库仑束缚的三例（氢原子三个能级）比值精确为 1.000000。引力束缚两例的残差源于天文参数的舍入与周期定义的选择，不属于模型误差。
式 (3.6) 为恒等式，非拟合所得；但此前未被作为独立关系提出。
3.4 尺度—时间关系
由 (3.6) 解出
物理含义。系统的时间读数与其能量成反比（角动量守恒时）。能量越高者（尺度越小、旋转越快），其时间读数越快。
尺度形式依系统类型分为两种：
表 4 太阳系行星对 τ∝R³ᐟ² 的检验
| 行星 | R (AU) | T (yr) | τ/R³ᐟ² |
|---|---|---|---|
| 水星 | 0.3871 | 0.2408 | 0.999820 |
| 金星 | 0.7233 | 0.6152 | 1.000089 |
| 地球 | 1.0000 | 1.0000 | 1.000000 |
| 火星 | 1.5237 | 1.8808 | 0.999985 |
| 木星 | 5.2029 | 11.8620 | 0.999516 |
| 土星 | 9.5370 | 29.4570 | 1.000163 |
式 (3.10) 即开普勒第三定律。全部结果落在 1±0.0005 内。
由此得频率的尺度依赖：
物理含义。尺度变化必然导致频率变化。此为第 8 节红移机制的基础。
4 三代粒子的几何来源
4.1 一个投影定理
取三个互相正交的单位矢量（分别对应自旋轴、黄道面法向、赤道面法向），投影至垂直于自旋轴的平面，投影算子为
物理含义。三个互相垂直的方向，斜投影至同一平面后，两两夹角恰为 120°。此即「六种夸克 = 三个正交面 × 两个旋向」的几何基础。
4.2 与 Koide 公式的关系
Koide 于 1981 年提出如下经验关系[3]：
mn∝1+λcos(δ+120∘⋅n),n=0,1,2(4.5) \sqrt{m_n} \propto 1 + \lambda\cos(\delta + 120^\circ\cdot n),\qquad n = 0,1,2 \tag{4.5} mn∝1+λcos(δ+120∘⋅n),n=0,1,2(4.5)
物理含义。三代粒子的质量开方后，排列为相位差 120° 的三个量。该 120° 相位三十年来缺乏解释，一直作为拟合参数使用。式 (4.4) 表明，它来源于三个正交面的投影几何。
4.3 解析关系
将 (4.5) 代入 Koide 量 K ≡ (Σ√m)² / Σm。由三相位之和为零：
∑m=M[3+λ∑ncos(δ+120∘n)]=3M \sum\sqrt{m} = \sqrt{M}\left[3 + \lambda\sum_n\cos(\delta+120^\circ n)\right] = 3\sqrt{M} ∑m=M[3+λn∑cos(δ+120∘n)]=3M
∑m=M[3+λ2∑ncos2(δ+120∘n)]=3M(1+λ22) \sum m = M\left[3 + \lambda^2\sum_n\cos^2(\delta+120^\circ n)\right] = 3M\left(1+\frac{\lambda^2}{2}\right) ∑m=M[3+λ2n∑cos2(δ+120∘n)]=3M(1+2λ2)
K=31+λ2/2,λ=2(3K−1)(4.8) \boxed{K = \frac{3}{1+\lambda^2/2}},\qquad \boxed{\lambda = \sqrt{2\left(\frac{3}{K}-1\right)}} \tag{4.8} K=1+λ2/23,λ=2(K3−1)(4.8)
物理含义。三代质量的「整齐程度」K，仅由投影扁率 λ 决定。此关系为导出所得，非拟合所得。
4.4 实测反解
表 5 由实测质量反解投影扁率 λ
| 三代 | K（实测） | λ |
|---|---|---|
| 轻子 (e,μ,τ) | 1.5000138 | 1.4142006 |
| 下型夸克 (d,s,b) | 1.3665764 | 1.5461351 |
| 上型夸克 (u,c,t) | 1.1777253 | 1.7591380 |
2=1.4142136(4.9) \sqrt{2} = 1.4142136 \tag{4.9} 2=1.4142136(4.9)
物理含义。轻子的 λ 与 √2 的相对偏差为 9.2×10⁻⁶。三个轻子质量跨越三个数量级（0.511 MeV 至 1777 MeV），仍能落在 √2 上。
4.5 λ 的物理意义
由 (4.8)，λ 越大则 K 越小，三代质量越分散，最外一代越重。此与「最外层的环流最松，维持所需能量最高，故质量最大」的图像一致。
上型夸克 λ 最大（1.759），对应顶夸克 173 GeV，质量极大且极不稳定；轻子 λ 最小（√2），三代分布最为对称。
4.6 自洽性检验
以 (√M, λ, δ) 三参数反解轻子质量，取 δ = 132.7328°，重建值与实测值的相对误差处于 10⁻¹³ 量级。
证伪条件。若 λ 仅为拟合参数，则不应恰好落在 √2。欲确立此结果，须检验下型与上型的 λ（1.5461、1.7591）是否亦为特殊值。目前尚看不出明显特殊性，本条未完成，不得作为结论使用。
5 光能耗散（LAH）
5.1 主方程
本文提出：光子在真空中传播时，与真空能量介质发生耦合，能量沿程耗散。
dνdℓ=−κν(5.1) \boxed{\frac{d\nu}{d\ell} = -\kappa\nu} \tag{5.1} dℓdν=−κν(5.1)
ν(ℓ)=ν0e−κℓ,1+z=eκℓ,ℓ(z)=ln(1+z)κ(5.2) \nu(\ell) = \nu_0 e^{-\kappa\ell},\qquad 1+z = e^{\kappa\ell},\qquad \ell(z) = \frac{\ln(1+z)}{\kappa} \tag{5.2} ν(ℓ)=ν0e−κℓ,1+z=eκℓ,ℓ(z)=κln(1+z)(5.2)
物理含义。光在真空中并非自由穿行。每前进一段，即被稀薄介质消耗一部分能量；传播越远，损失越多，故越远的光红移越大。
5.2 一阶形式的唯一性
式 (5.1) 的线性形式并非假定，而由观测唯一确定。
观测事实：同一天体的所有谱线给出相同的红移（无色性）。若能量损失率不正比于当前频率，则不同频率的光在相同距离上损失的比例不同，蓝线与红线将给出不同的 z。观测上不存在此现象。
结论。一阶线性形式由无色性锁定，不可修改。衰减律本身无误。
5.3 κ 的量级
在低 z 极限下须与哈勃定律 cz = H₀ℓ 一致，代入得
κ=H0c=7.567×10−27 m−1(5.3) \boxed{\kappa = \frac{H_0}{c} = 7.567\times10^{-27}\ \mathrm{m^{-1}}} \tag{5.3} κ=cH0=7.567×10−27 m−1(5.3)
（取 H₀ = 70 km·s⁻¹·Mpc⁻¹）
物理含义。耗散系数的量级由哈勃常数锁定，非自由参数。但须注意：按本文立场，常数是局域读数，κ 本身可随尺度缓变（见第 8 节）。
等效衰减长度 c/H₀ = 4282.7 Mpc。
5.4 距离—红移关系
dLLAH(z)=(1+z) ℓ(z)=(1+z)cH0ln(1+z)(5.4) d_L^{\mathrm{LAH}}(z) = (1+z)\,\ell(z) = (1+z)\frac{c}{H_0}\ln(1+z) \tag{5.4} dLLAH(z)=(1+z)ℓ(z)=(1+z)H0cln(1+z)(5.4)
dLΛCDM(z)=(1+z)cH0∫0zdz′E(z′)(5.5) d_L^{\Lambda\mathrm{CDM}}(z) = (1+z)\frac{c}{H_0}\int_0^z \frac{dz'}{E(z')} \tag{5.5} dLΛCDM(z)=(1+z)H0c∫0zE(z′)dz′(5.5)
表 6 LAH 与 ΛCDM 的光度距离对比（Ωm=0.3, ΩΛ=0.7）
| z | d_L(LAH) Mpc | d_L(ΛCDM) Mpc | 比值 |
|---|---|---|---|
| 0.1 | 449.0 | 460.3 | 0.975 |
| 0.5 | 2604.8 | 2832.9 | 0.919 |
| 1.0 | 5937.2 | 6607.7 | 0.899 |
| 2.0 | 14115.2 | 15539.6 | 0.908 |
| 5.0 | 46041.9 | 46652.2 | 0.987 |
| 10.0 | 112965.4 | 103842.7 | 1.088 |
物理含义。纯 LAH 在 z<3 内给出的光度距离系统性偏小，最大偏差约 10%；在 z>7 后反超 ΛCDM。这一非单调行为是后续第 8 节判定红移机制分配的关键线索。
5.5 距离模数偏差
Δμ(z)=5log10dLLAHdLΛCDM(5.6) \Delta\mu(z) = 5\log_{10}\frac{d_L^{\mathrm{LAH}}}{d_L^{\Lambda\mathrm{CDM}}} \tag{5.6} Δμ(z)=5log10dLΛCDMdLLAH(5.6)
表 7 Δμ 随红移的分布
| z | 0.05 | 0.2 | 0.5 | 1.0 | 1.3 | 2.0 | 3.0 |
|---|---|---|---|---|---|---|---|
| Δμ (mag) | -0.028 | -0.098 | -0.182 | -0.232 | -0.235 | -0.209 | -0.148 |
Pantheon+ 在高红移处的距离模数精度约为 0.1–0.2 mag[4,5]。LAH 的最大偏差为 0.235 mag，已进入判别区间。
判决方式。以 Pantheon+ 的 μ(z) 数据拟合单参数 κ。若 χ² 不可接受，则纯指数 LAH 被证伪。
5.6 关于光变曲线展宽
观测事实：高红移超新星光变曲线宽度被拉长 (1+z) 倍[6,7]。历史上各类「光疲劳」（tired light）模型均在此失败[8]。
失败原因：此类模型将「能量」与「时钟」视为两件独立事物，仅修改前者而不修改后者。
本文的处理：在 HDT 中，时间间隔与频率是同一读数的两个侧面，
Δt∝1ν ⟹ ΔtobsΔtemit=νemitνobs=1+z(5.7) \Delta t \propto \frac{1}{\nu}\ \Longrightarrow\ \frac{\Delta t_{\mathrm{obs}}}{\Delta t_{\mathrm{emit}}} = \frac{\nu_{\mathrm{emit}}}{\nu_{\mathrm{obs}}} = 1+z \tag{5.7} Δt∝ν1 ⟹ ΔtemitΔtobs=νobsνemit=1+z(5.7)
物理含义。光子被稀薄介质消耗能量、频率降低的同时，其所携带的时间读数同步被拉长。频率与时钟是同一件事的两种写法。
本结果不依赖空间膨胀，亦不依赖「时间自身膨胀」。
6 真空折射与光线偏折
6.1 折射率模型
若真空为能量介质，则质量周围的介质密度升高，折射率随之变化：
n(r)=1+2GMc2r(6.1) n(r) = 1 + \frac{2GM}{c^2 r} \tag{6.1} n(r)=1+c2r2GM(6.1)
由费马原理 δ∫n ds = 0：
dds(ndrds)=∇n,α=∫∇⊥(lnn) ds(6.2) \frac{d}{ds}\left(n\frac{d\mathbf{r}}{ds}\right) = \nabla n,\qquad \alpha = \int \nabla_\perp(\ln n)\,ds \tag{6.2} dsd(ndsdr)=∇n,α=∫∇⊥(lnn)ds(6.2)
α=4GMc2b(6.3) \boxed{\alpha = \frac{4GM}{c^2 b}} \tag{6.3} α=c2b4GM(6.3)
物理含义。光经过太阳附近时的偏折角，按「真空是介质」算出的数值，与按「时空弯曲」算出的完全一致。
6.2 反推真空密度
由 n−1 ∝ δρ：
δρ(r)=2GM/c2αρ r(6.4) \delta\rho(r) = \frac{2GM/c^2}{\alpha_\rho\, r} \tag{6.4} δρ(r)=αρr2GM/c2(6.4)
物理含义。质量周围的真空并非空无一物，而是密度按 1/r 隆起的能量介质。光线不是被「弯曲的时空」改变方向，而是进入了折射率不均匀的介质。
重要说明。式 (6.3) 与广义相对论在观测上等价，不能用于区分两者。其价值在于本体论层面的替换（无需引入时空弯曲），而非给出新预测。
7 真空密度扰动与星系动力学
7.1 一个矛盾
星系尺度的两项观测对真空密度剖面提出不同要求：
- 引力透镜要求偏折角 α ∝ 1/b，对应 δρ ∝ 1/r；
- 旋转曲线平坦（v = 常数）要求 δρ ∝ 1/r²。
证明概要：若 δρ ∝ r⁻ᵏ，则 v² ∝ r²⁻ᵏ，平坦需 k=2；而 α ∝ b⁻ᵏ，透镜需 k=1。
7.2 统一形式
δρ(r)=Ar (1+r/rc)(7.1) \boxed{\delta\rho(r) = \frac{A}{r\,(1+r/r_c)}} \tag{7.1} δρ(r)=r(1+r/rc)A(7.1)
表 8 式 (7.1) 的两段行为
| 区间 | 行为 | 解释 |
|---|---|---|
| r ≪ r_c | δρ ≈ A/r | 满足引力透镜 |
| r ≫ r_c | δρ ≈ A r_c/r² | 满足旋转曲线平坦 |
大尺度极限下：
Meff(r)=4πArcr,v2=GMeffr=4πGArc=常数(7.2) M_{\mathrm{eff}}(r) = 4\pi A r_c r,\qquad v^2 = \frac{G M_{\mathrm{eff}}}{r} = 4\pi G A r_c = \text{常数} \tag{7.2} Meff(r)=4πArcr,v2=rGMeff=4πGArc=常数(7.2)
物理含义。星系外围转速为常数，是由密度剖面算出的结果，不是补充假设的结果。
7.3 定量预测
取 v_flat = 200 km·s⁻¹：
A rc=v24πG=4.769×1019 kg m−2(7.3) A\,r_c = \frac{v^2}{4\pi G} = 4.769\times10^{19}\ \mathrm{kg\,m^{-2}} \tag{7.3} Arc=4πGv2=4.769×1019 kgm−2(7.3)
表 9 不同转折半径下的扰动密度
| r_c | 5 kpc | 10 kpc | 20 kpc | 50 kpc |
|---|---|---|---|---|
| δρ(r_c) (kg·m⁻³) | 1.00×10⁻²¹ | 2.50×10⁻²² | 6.26×10⁻²³ | 1.00×10⁻²³ |
物理含义。真空密度扰动应在 ~10 kpc 的星系尺度处由 1/r 转为 1/r²。星系旋转曲线的平坦，源于真空场能自身参与引力，且其密度剖面恰在星系尺度转折——因此无需引入暗物质[9,10]。
证伪条件。以 SPARC 数据库的星系旋转曲线[11]，用 A、r_c 双参数拟合。若拟合优度不如 NFW 暗物质晕剖面，则本条被证伪。
8 红移的机制分配
8.1 问题
第 5 节将全部红移归因于传播耗散。本文在此检验该处理是否成立。
设红移全部由耗散产生，反推 κ 所需的红移依赖：
表 10 纯 LAH 所需的 κ_eff/κ₀
| z | 0.1 | 0.5 | 1.0 | 2.0 | 3.0 | 5.0 | 10.0 |
|---|---|---|---|---|---|---|---|
| κ_eff/κ₀ | 0.975 | 0.920 | 0.899 | 0.909 | 0.934 | 0.987 | 1.087 |
物理含义。κ 在 z<3 内仅需比常数低最多约 10%，但在 z>7 后需反超。此为先降后升的非单调行为，单靠调节 κ 无法修复。
结论。问题不在于 κ 取值多少，而在于红移这笔账被全部记在耗散头上。
8.2 乘法模型
1+z=(1+z)胀缩×(1+z)耗散=ap⋅P(8.1) 1+z = (1+z)_{\text{胀缩}} \times (1+z)_{\text{耗散}} = a^{p}\cdot P \tag{8.1} 1+z=(1+z)胀缩×(1+z)耗散=ap⋅P(8.1)
表 11 式 (8.1) 的符号
| 符号 | 含义 |
|---|---|
| a | 物质尺度（相对当前）。a<1 表示过去更紧，即处于胀相位 |
| p | 频率随尺度的标度指数：ν ∝ a⁻ᵖ |
| P | 传播耗散因子 P = exp(∫κ dr) ≥ 1 |
| g | 胀缩占 ln(1+z) 的份额，耗散占 1−g |
a=(1+z)g/p,P=(1+z)1−g(8.2) a = (1+z)^{g/p},\qquad P = (1+z)^{1-g} \tag{8.2} a=(1+z)g/p,P=(1+z)1−g(8.2)
8.3 距离的推导
光度距离：
dL=(1+z)⋅r(8.3) d_L = (1+z)\cdot r \tag{8.3} dL=(1+z)⋅r(8.3)
（每个光子的能量损失 (1+z)，到达率损失 (1+z)，面积 4πr²）
角直径距离：
dA=ra(8.4) d_A = \frac{r}{a} \tag{8.4} dA=ar(8.4)
（过去的物质「标准尺」本身短了 a 倍）
故
dLdA=(1+z)⋅a(8.5) \frac{d_L}{d_A} = (1+z)\cdot a \tag{8.5} dAdL=(1+z)⋅a(8.5)
（标准宇宙学要求此比值为 (1+z)²）
令 d_L 与观测一致，即 r = I(z)/κ₀，其中 I(z) = ∫₀ᶻ dz′/E(z′)，则
dA(HDT)dA(ΛCDM)=(1+z) 1−g/p(8.6) \boxed{\frac{d_A(\mathrm{HDT})}{d_A(\Lambda\mathrm{CDM})} = (1+z)^{\,1-g/p}} \tag{8.6} dA(ΛCDM)dA(HDT)=(1+z)1−g/p(8.6)
g=p≤1(8.7) g = p \le 1 \tag{8.7} g=p≤1(8.7)
8.4 机制份额的分配
表 12 不同参数组合下 d_A(HDT)/d_A(ΛCDM)
| 组合 | z=0.1 | 0.5 | 1.0 | 2.0 | 3.0 | 5.0 | 10.0 |
|---|---|---|---|---|---|---|---|
| 纯耗散 LAH (p=1,g=0) | 1.100 | 1.500 | 2.000 | 3.000 | 4.000 | 6.000 | 11.000 |
| 各半 (p=1,g=0.5) | 1.049 | 1.225 | 1.414 | 1.732 | 2.000 | 2.449 | 3.317 |
| 纯胀缩 (p=1,g=1) | 1.000 | 1.000 | 1.000 | 1.000 | 1.000 | 1.000 | 1.000 |
| 纯胀缩 (p=3/2,g=1) | 1.032 | 1.145 | 1.260 | 1.442 | 1.587 | 1.817 | 2.224 |
| 各半 (p=3/2,g=0.5) | 1.066 | 1.310 | 1.587 | 2.080 | 2.520 | 3.302 | 4.946 |
| 胀缩+耗散 (p=1/2,g=0.5) | 1.000 | 1.000 | 1.000 | 1.000 | 1.000 | 1.000 | 1.000 |
物理含义。分配给胀缩的份额越多，d_A 越准确；分配给耗散的份额越多，d_A 偏离越严重。纯 LAH 在 z=1 处偏差 2 倍，在 z=10 处偏差 11 倍。
8.5 参数 p 的判定
p 为原子跃迁频率随物质尺度的标度指数。由于红移是依据原子谱线测定的，应取原子（电磁束缚）的 p，而非引力轨道（开普勒）的 3/2。
表 13 p 的三种可能取值及其后果
| p | 条件 | 后果 |
|---|---|---|
| p=1 | 所有长度乘以 a，c,ℏ,α 不变 ⟹ m ∝ 1/a ⟹ ν ∝ 1/a | g=1，耗散份额为零。纯胀缩，LAH 退出；d_L 与 d_A 均为 1.000 |
| p=1/2 | 标度更弱 | g=1/2，胀缩与耗散各占一半，两者均成立，且精确复现 ΛCDM |
| p=3/2 | 开普勒式束缚 τ ∝ R³ᐟ² | 需 g=1.5>1，不可能。最优情形仍差 26%（z=1），可被 DESI BAO 判死 |
8.6 本节结论
- LAH 的衰减律无误，其形式由无色性锁定；
- 错误在于将红移全部归因于耗散；
- 使两类距离同时正确的配方为 g=p；
- LAH 能否成立，不取决于 κ 调节得多准，而取决于 p < 1 是否成立。
判决方式。以超新星的 d_L 与 DESI 的 BAO 观测所得 d_A[12] 联合拟合 (p,g)。若拟合结果给出 g≠p，或 p>1，则乘法模型被证伪。
图 1 红移机制判决流程图
物理含义。红移机制的判决被归结为对 (p,g) 的联合拟合：若拟合给出 g=p≤1，则「物质胀缩 × 传播耗散」模型成立；若 g≠p 或 p>1，则乘法模型被证伪。此流程将第 8 节的结论转化为可直接执行的观测检验。
9 结论
本文以「底层只有能量在流动」为唯一公设，建立混沌拓扑学（HDT）框架，并完成三项可检验的工作。
第一，普适关系。由自由涡旋（角动量守恒）与中心力束缚（开普勒 / 库仑）两种独立机制分别导出同一形式 |E|·τ = π·L。以氢原子三个能级检验，比值精确为 1.000000；以地—日、月—地系统检验，残差处于天文参数舍入量级。由此解出「系统尺度越小、能量越高，其时间读数越快」的定量形式，并自然导出开普勒第三定律。
第二，三代粒子的几何来源。三个互相正交的方向投影至垂直于自旋轴的平面，两两夹角恰为 120°。本文证明 Koide 公式中三十年来缺乏解释的 120° 相位正是该投影的几何必然，并导出 K = 3/(1+λ²/2) 的解析关系。以实测质量反解，轻子投影扁率与 √2 的相对偏差为 9.2×10⁻⁶。此条尚未完全确立，须检验下型与上型夸克的 λ 是否亦为特殊值。
第三，红移的乘法分配模型。光能耗散律的一阶线性形式由观测的无色性唯一锁定，不可修改；但将其独立用于解释全部红移，会导致角直径距离在 z=1 处偏差达 2 倍。本文提出「物质胀缩 × 传播耗散」模型，并证明其与标准宇宙学完全一致的充要条件为 g=p≤1。红移机制的判决由此转化为对参数 p（原子频率随物质尺度的标度指数）的测定。
此外，本文给出真空密度扰动剖面 δρ = A/[r(1+r/r_c)]，在小尺度退化为 1/r（满足引力透镜），在大尺度退化为 1/r²（满足星系旋转曲线平坦），转折发生于 ~10 kpc 星系尺度——因此无需引入暗物质。
本文同时如实记录一条被证伪的猜想（与实测相差两个数量级）及一条仅达定性程度的推论（温度下限，耦合系数待标定）。按本文的自我约束，未通过检验者不得作为结论使用。
展望。后续工作集中于三处：其一，以 Pantheon+ 的 μ(z) 数据拟合单参数 κ，检验纯指数 LAH 是否被证伪；其二，以 SPARC 星系旋转曲线用 A、r_c 双参数拟合，与 NFW 暗物质晕剖面比较；其三，以超新星的 d_L 与 DESI 的 BAO 观测所得 d_A 联合拟合 (p,g)，判定乘法模型是否成立。三者任一失败，对应条目即被排除。
参考文献
[1] Weinberg, S. The cosmological constant problem. Rev. Mod. Phys. 61, 1–23 (1989).
[2] Carroll, S. M. The cosmological constant. Living Rev. Relativ. 4, 1 (2001).
[3] Koide, Y. A new view of quark and lepton mass formula. Phys. Rev. D 28, 252 (1983).
[4] Scolnic, D. et al. The Pantheon+ type Ia supernova sample. Astrophys. J. 938, 113 (2022).
[5] Brout, D. et al. The Pantheon+ analysis: cosmological constraints. Astrophys. J. 938, 110 (2022).
[6] Riess, A. G. et al. Observational evidence from supernovae for an accelerating universe. Astron. J. 116, 1009 (1998).
[7] Perlmutter, S. et al. Measurements of Ω and Λ from 42 high-redshift supernovae. Astrophys. J. 517, 565 (1999).
[8] Wright, E. L. A cosmology calculator for the World Wide Web. Publ. Astron. Soc. Pac. 118, 1711 (2006).
[9] Rubin, V. C., Ford, W. K. Rotation of the Andromeda nebula from a spectroscopic survey of emission regions. Astrophys. J. 159, 379 (1970).
[10] Bosma, A. The distribution and kinematics of neutral hydrogen in spiral galaxies. Astron. J. 86, 1825 (1981).
[11] Lelli, F., McGaugh, S. S., Schombert, J. M. SPARC: mass models for 175 disk galaxies. Astrophys. J. 836, 152 (2017).
[12] DESI Collaboration. DESI 2024: baryon acoustic oscillations from the early data release. arXiv:2404.03002 (2024).
转载自 CSDN-专业IT技术社区
原文链接：https://blog.csdn.net/2602_96503874/article/details/167087099
