# 空间运动与可变质量耦合模型的动力学与传播约束

## 变分结构 状态边界 电磁响应和零质量转换的阶段性研究

**UFT Labs**

2026 年 10 月 22 日  
研究论文草稿 V1.0

## 摘要

本文研究一个由空间光速运动、零质量物质光速传播、垂直方向演化及平行比例关系四项命题启发的经典有效模型。模型采用状态场 \(\chi\)、质量函数 \(m=m_0(1-\chi)^2\)、非负双稳势 \(U=\Lambda^4\chi^2(1-\chi)^2\) 及电磁耦合 \(Z=e^{\beta\chi}\)。研究完成了六项基础计算：作用量变分复核、真空与边界解、电磁弱场响应、受控粒子运动、包含控制场反作用的动量分析，以及零质量端点和退出过程。结果表明，该模型具有可解析的平面状态壁与非负线性扰动谱，但无外源的三维有限静态状态泡受到缩放不稳定性的限制；静态电磁响应具有 Yukawa 型空间核；具有非零初始动量的粒子可以在指定的平滑背景中进入零质量传播平台并恢复质量，且能量交换闭合。在一维周期电磁场—状态场—双试样系统中，试样获得的动量由场动量补偿，没有出现封闭系统总动量增长。模型没有产生物质与独立空间流的自动方向锁定，也未建立宏观物体在零总不变质量状态下保持结构的机制。本文据此明确四项核心命题的数学实现范围、当前模型的可检验内容及下一阶段目标。

**关键词：**可变质量；状态场；零质量传播；电磁耦合；动量守恒；空间运动；有效场论

## 1 引言

“空间自身具有运动结构”“物质质量与空间状态相关”可以被提出为物理研究假设，但只有在这些陈述对应明确的变量、方程和测量时，才能讨论其与引力控制或飞行的关系。原始统一场论材料及人工场技术指南提出了这一方向，同时混用了几何关系、物理定义、既有定律与工程推测 [1,2]。其概念性启发不能替代动力学证明。

此前研究稿 [3] 提出四项核心命题：空间以光速运动；物质静止质量归零时随空间以光速运动；垂直关系生成运动方向的变化；平行关系表示比例耦合。本文保留这些命题，但区分其逻辑地位。零轨迹条件可以实现局部光速运动；垂直关系可以写成方向演化；平行关系可以写成比例表达式。然而，“任意物体自动降质量并沿指定空间流飞行”的强陈述仍需要独立耦合机制。

本文的目标是完成 [3] 所列计算安排的第 1—6 项，使讨论从概念组织推进到解析结果、数值结果和可复核的限制。本文采用标准洛伦兹几何与协变守恒语言 [4]，不主张已统一全部基本相互作用，不把标量场附加力直接等同于反引力，也不将低维数值解认定为实际飞行器方案。

## 2 四项核心命题的数学表述

### 2.1 空间光速轨迹

采用号差 \((-+++)\) 和自然单位 \(c=\hbar=1\)。空间运动方向由非零、未来指向的 \(k^\mu\) 表示：

\[
g_{\mu\nu}k^\mu k^\nu=0.\tag{1}
\]

相对于单位观察者四速度 \(u^\mu\)，有：

\[
k^\mu=\omega(u^\mu+n^\mu),\qquad
u\cdot n=0,\qquad n^2=1,\qquad\omega>0.\tag{2}
\]

因此该观察者测得轨迹速率为光速。式（1）本身也是既有光锥几何的结构，不能仅凭新的名称赋予它独立物理证据。

若采用独立零流实现：

\[
S_k=-\frac12\int\sqrt{-g}\,\rho\,(\nabla\vartheta)^2d^4x,
\qquad k_\mu=\partial_\mu\vartheta,\tag{3}
\]

则变分给出 \(k^2=0\)、\(\nabla_\mu(\rho k^\mu)=0\)，并有 \(k^\nu\nabla_\nu k_\mu=0\)。它是一种无旋零质量流，不自动生成任意螺旋轨迹，也没有在现有模型中自动携带物体。

### 2.2 垂直演化与平行比例

在局部惯性系写：

\[
\mathbf C=c\mathbf n,\qquad
\dot{\mathbf n}=\boldsymbol\Omega\times\mathbf n,\qquad
\mathbf n\cdot\dot{\mathbf n}=0.\tag{4}
\]

此式保留“垂直即运动”的方向演化内核。它保证方向变化与当前方向垂直，但不能代替 \(\boldsymbol\Omega\) 的动力学来源。在平直背景中，自由零质量粒子仍沿直线运动。

对同一点的非零平行矢量，可写：

\[
\mathbf A=\alpha(x,t)\mathbf B.\tag{5}
\]

平行不自动保证 \(\alpha\) 为常数，更不把原本不平行的自由度强制锁定。比例函数必须由构成关系或动力学决定。

当 \(\boldsymbol\Omega\) 的轴和大小恒定时，式（4）给出螺旋运动学，其横向和轴向速率满足：

\[
a^2\omega^2+v_\parallel^2=c^2.\tag{6}
\]

该结果不能被解释为“轴向与旋转分量各自均以光速运动”。

### 2.3 零质量传播命题

取局部质量壳与速度关系：

\[
\mathcal E^2=c^2p^2+m^2c^4,\qquad
\mathbf v=\frac{c^2\mathbf p}{\mathcal E}.\tag{7}
\]

在 \(m=0,\mathcal E>0\) 时，\(|\mathbf v|=c\)。若另有 \(\mathbf p=(\mathcal E/c)\mathbf n\)，便得到 \(\mathbf v=\mathbf C\)。因此“随空间运动”在这里是方向匹配状态的传播关系；自动匹配仍需证明。

## 3 有效作用量与变分结果

### 3.1 模型定义

引入无量纲状态场 \(\chi\)：

\[
m_a(\chi)=m_{a0}(1-\chi)^2,\quad
U(\chi)=\Lambda^4\chi^2(1-\chi)^2,\quad
Z(\chi)=e^{\beta\chi}.\tag{8}
\]

这三个函数是构成选择，不是从四项核心命题唯一推出的自然定律。它们分别提供非负质量、两个经典简并状态和正电磁动能系数。

场作用量为：

\[
S_f=\int\sqrt{-g}\,d^4x\left[
\frac{M_{\rm Pl}^2}{2}R-\frac Z4F_{\mu\nu}F^{\mu\nu}
-\frac{f^2}{2}(\nabla\chi)^2-U
-\frac Y4F_{\mu\nu}\widetilde F^{\mu\nu}
\right].\tag{9}
\]

主体数值研究设 \(Y=0\)；解析响应允许 \(Y=\alpha\chi\)。\(f>0\) 为状态场能标。中性粒子采用在零质量支仍可使用的一阶作用量：

\[
S_m=\sum_a\int d\lambda\left[
p_\mu\dot x^\mu-\frac e2(g^{\mu\nu}p_\mu p_\nu+m_a^2)
\right].\tag{10}
\]

完整系统还包括源装置和可选的式（3）。没有给出驱动装置自由度时，只能研究外部背景下的子系统，不能判断整机推进。

### 3.2 变分方程

采用 \(j^\nu=-\delta S_{\rm src}/(\sqrt{-g}\delta A_\nu)\)，场方程为：

\[
\nabla_\mu(ZF^{\mu\nu}+Y\widetilde F^{\mu\nu})=j^\nu,
\qquad \nabla_\mu\widetilde F^{\mu\nu}=0,\tag{11}
\]

\[
f^2\Box\chi-U'-\frac{Z'}4F^2-\frac{Y'}4F\widetilde F-J_\chi=0,\tag{12}
\]

\[
J_{\chi,m}=\frac1{\sqrt{-g}}\sum_a\int e\,m_am'_a\delta^4(x-x_a)d\lambda.\tag{13}
\]

粒子方程为：

\[
\dot x^\mu=e g^{\mu\nu}p_\nu,\qquad
p^2+m^2=0,\tag{14}
\]

\[
\dot p_\mu=-\frac e2\partial_\mu g^{\alpha\beta}p_\alpha p_\beta
-\frac e2\partial_\mu m^2.\tag{15}
\]

引力变分在相容边界条件下给出 \(M_{\rm Pl}^2G_{\mu\nu}=T^{\rm total}_{\mu\nu}\)。完整系统满足：

\[
\nabla_\mu T^{\mu\nu}_{\rm total}=0.\tag{16}
\]

上述方程经手工变分复核；另用符号计算独立核对七项势能、哈密顿量及边界解恒等式。该符号检查的范围有限，不作为全部张量方程或量子一致性的自动证明。

## 4 真空与状态边界

### 4.1 真空与小扰动

势能在 \(\chi=0,1\) 取极小值，在 \(\chi=1/2\) 取势垒极大值。极小值附近：

\[
(\Box-\mu^2)\delta\chi=0,\qquad
\mu^2=\frac{2\Lambda^4}{f^2},\qquad
\omega^2=k^2+\mu^2.\tag{17}
\]

因此状态场扰动一般不是光速群速度。这明确限定了核心光速命题适用于零轨迹扇区，而不是模型中所有自由度。

### 4.2 平面壁与线性谱

连接两真空的静态平面解为：

\[
\chi_w(z)=\frac1{1+e^{-\mu(z-z_0)}},\qquad
\sigma=\frac{\sqrt2}{6}f\Lambda^2.\tag{18}
\]

其扰动算符为：

\[
\mathcal L=-\partial_z^2+
\mu^2\left[1-\frac32\operatorname{sech}^2\left(\frac{\mu z}{2}\right)\right].\tag{19}
\]

离散本征值为 \(0\) 与 \(3\mu^2/4\)，连续谱阈值为 \(\mu^2\)。零模表示平移自由度；没有负本征值表示该平面解的线性稳定性。

在 \(f=\Lambda=1\) 时，独立边值求解与解析壁的最大剖面差为 \(4.69\times10^{-8}\)，张力为 \(0.235702260395\)。有限差分谱在 600、1200、2400 个内点下趋向解析值；最细网格的形状模为 \(1.49999129\)。最低模的小负离散值随网格加密趋零，是截断误差而非已发现的物理负模。

![平面状态壁及能量密度](Research_阶段1至6/figures/01_wall.png)

图 1　平面状态壁的独立边值解与解析解比较。图中量均为无量纲。平面边界稳定不代表有限三维舱体稳定。

### 4.3 有限三维区域的限制

对无外源、非负势能、正则有限能量的三维静态标量构型，缩放 \(\chi_\lambda(\mathbf x)=\chi(\mathbf x/\lambda)\) 给出：

\[
E(\lambda)=\lambda E_{\rm grad}+\lambda^3E_U,
\qquad E'(1)=E_{\rm grad}+3E_U>0.\tag{20}
\]

非平凡构型不可能驻定。这是本模型中的 Derrick 型障碍 [5]。薄壁近似下，简并真空泡的能量 \(4\pi\sigma R^2\) 随半径增加，表面张力导致收缩。

因此，当前模型不支持无外部条件维持的有限静态“零质量区域”。外场压力、物质、守恒荷或动态过程可能改变该结论，但必须进入新计算。

## 5 电磁弱场响应

### 5.1 线性响应核

取 \(F^2=2(B^2-E^2)\)、\(F\widetilde F=-4\mathbf E\cdot\mathbf B\)。在普通状态附近、忽略物质源并将电磁场作为弱背景：

\[
(\nabla^2-\mu^2)\delta\chi=\frac{S}{f^2},\qquad
S=\frac\beta2(B^2-E^2)-\alpha\mathbf E\cdot\mathbf B.\tag{21}
\]

无限空间解为：

\[
\delta\chi(\mathbf r)=-\frac1{4\pi f^2}\int
\frac{e^{-\mu|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|}
S(\mathbf r')d^3r'.\tag{22}
\]

在 \(\beta>0,\alpha=0\) 下，电场主导提高 \(\chi\)，磁场主导降低 \(\chi\)。这是一项模型特定预测，不是对所有电磁重力理论的结论。

一维高斯源 \(S/f^2=-0.02e^{-z^2/2}\) 的有限差分解与独立格林积分比较，采样最大误差由 \(4.36\times10^{-7}\) 降至 \(2.74\times10^{-8}\)，符合二阶空间离散预期。

![线性状态场响应](Research_阶段1至6/figures/02_weak_response.png)

图 2　相反局域源产生相反的线性响应。源是无量纲示例，未对应实际电压或线圈。

### 5.2 对控制构想的约束

第一，单束理想真空平面波有 \(F^2=F\widetilde F=0\)，不通过式（21）的直接耦合激发状态场。近场、驻波和叠加场须按各自不变量计算。

第二，纯 \(ZF^2\) 源对单独电场极性反转为偶性；若加入 \(YF\widetilde F\)，单独翻转一个场方向可改变交叉源符号。普通放电和材料效应也可能具有极性依赖，不能由极性效应单独确认新场。

第三，固定自由电荷与固定电压是不同控制条件。电源做功必须计入相应能量势。第四，持续纯电场通常在 \(\chi=1\) 仍提供非零源，不会自动将状态停在质量零点。第五，弱场线性解不能外推到 \(\delta\chi\sim1\) 计算转换阈值。

## 6 受控粒子转换与能量交换

### 6.1 计算方法

在平直、均匀但时间依赖的背景中：

\[
H=\sqrt{p^2+m^2},\quad
\dot x=p/H,\quad\dot p=0,\quad
\dot H=m\dot m/H.\tag{23}
\]

选择 \(q(s)=6s^5-15s^4+10s^3\)，使 \(\chi\) 在 \(t=0\) 至 10 从 0 平滑升至 1，在 10 至 15 保持 1，在 15 至 25 平滑恢复 0。取 \(m_0=1,p_0=0.3\)。采用 DOP853 积分，并同时累计控制场对粒子做功。

### 6.2 转换结果

| 量 | 初始状态 | 零质量平台 | 恢复状态 |
|---|---:|---:|---:|
| 质量 \(m\) | 1 | 0 | 1 |
| 动量 \(p\) | 0.3 | 0.3 | 0.3 |
| 能量 \(H\) | 1.04403065 | 0.3 | 1.04403065 |
| 速度 \(v/c\) | 0.28734789 | 1 | 0.28734789 |

进入平台时粒子释放无量纲能量 \(0.74403065\) 给控制源；恢复时需要返还。完整循环没有产生净能量。工作积分与能量差的诊断残差约为浮点舍入量级，在报告中不将其解释为物理预测达到十五位精度。

![质量转换与恢复循环](Research_阶段1至6/figures/03_conversion_cycle.png)

图 3　指定背景下的进入、零质量平台和恢复过程。绿色区域为零质量平台。“控制源增益”按粒子能量交换计算；此图没有模拟真实电源如何生成背景。

另一个静态空间质量阱算例取 \(\chi(x)=0.35e^{-x^2/(2\times1.5^2)}\)，初态 \(x=-6,p=0.6\)。粒子在质量较低区域加速，最大速度约 \(0.93204c\)，离开后减速，能量保持。此例核验的是质量梯度动力学，不是整机推进。

## 7 控制场反作用与全系统动量

### 7.1 解析约束

在平直封闭系统中：

\[
\frac{dP^i_{\rm total}}{dt}=-\oint T^{ji}n_jdS.\tag{24}
\]

若没有边界通量，总动量不变。内部质量变化和内部场交换不能绕开这一结果。与地球、外部场或辐射的交换可以产生子系统推力，但必须包含相应反作用。

### 7.2 动态闭合算例

为实际纳入控制场而不是只把它写在能量账目中，构造一维周期模型。设横向电磁势为 \(A_y(x,t)\)，并引入两个质量依赖于平滑局域状态的中性试样：

\[
H=\int_0^{2\pi}\left[
\frac{\Pi_A^2}{2Z}+\frac Z2A_x^2+
\frac{\pi_\chi^2}{2}+\frac{\chi_x^2}{2}+U
\right]dx+\sum_{a=1}^2\sqrt{p_a^2+m_a(\bar\chi_a)^2}.\tag{25}
\]

\(\bar\chi_a\) 为宽度 \(\sigma=0.35\) 的归一化周期高斯平均。该有限尺度正则化保持本算例的空间平移对称性，但不是三维洛伦兹协变微观理论。取 \(\beta=0.4,m_{10}=0.03,m_{20}=0.05\)。详细初值和方程在计算报告及脚本中给出。

所有控制场和试样共同演化，没有外加时间驱动。总动量为：

\[
P=p_1+p_2-\int_0^{2\pi}(\Pi_AA_x+\pi_\chi\chi_x)dx.\tag{26}
\]

初始电磁场携带非零动量，故结果不应被描述为系统从零动量自行启动。

### 7.3 数值结果与边界

使用 64、128、256 点 Fourier 配点法积分到 \(t=10\)。在最细网格，两个试样的末动量为：

\[
p_1=3.90018423\times10^{-4},\qquad
p_2=-5.44101564\times10^{-4}.\tag{27}
\]

试样总动量变化约 \(-1.54083141\times10^{-4}\)，由场动量增加补偿。总动量保持约 \(0.1023362610\)。主计算的守恒诊断达到浮点舍入量级，并补做不同最大时间步的独立检查；具体数值和误差定义见附件。

![闭合系统动量分配](Research_阶段1至6/figures/04_closed_momentum.png)

图 4　左：总动量及各部分动量。右：试样动量变化与场动量变化相互补偿。周期边界无净外部通量，控制能量和动量均属于系统。

此结果验证了所实现的受限模型中的反作用过程，不是全三维装置、强引力区或 \(\chi\approx1\) 闭合转换的模拟。在该弱激发算例中 \(\max|\chi|\approx0.057\)。完整进入零质量平台的闭合电磁控制过程仍是下一阶段问题。

## 8 零质量端点与方向控制

令 \(\delta=1-\chi\)。质量梯度力为：

\[
\dot{\mathbf p}=\frac{2m_0^2\delta^3}{\sqrt{p^2+m_0^2\delta^4}}\nabla\chi.\tag{28}
\]

在非零有限动量与光滑背景下，\(\delta\to0\) 时该力趋于零。相应方向变化：

\[
\frac{d\hat{\mathbf p}}{dt}
=\frac1{|p|}(I-\hat{\mathbf p}\hat{\mathbf p}^{\mathsf T})\dot{\mathbf p}
=O(\delta^3).\tag{29}
\]

因此低质量不意味着该质量梯度机制的转向能力无限增强。含 \(\nabla\ln m\) 的低速公式在零质量处失效，不能用于推导无限加速度。

若 \(p=a\delta^2\) 同时趋于零，则：

\[
v\to\frac{a}{\sqrt{a^2+m_0^2}}.\tag{30}
\]

这说明双零极限依赖路径，但端点的能量也为零，并不是一个非零能量的光速粒子。它不否定式（7）的零质量传播结论，而是要求保留非零能量条件。

![零质量端点的路径与力](Research_阶段1至6/figures/05_endpoint.png)

图 5　固定非零动量与动量同时趋零的不同极限；有限动量下质量梯度力按三次幂关断。

第 6 节还修正了一种可能的误读：零质量点的瞬时力与源为零，不意味着状态永久锁死。只要外部 \(\chi\) 后续改变，非零动量轨迹即可恢复质量。当前作用量却没有将物质方向自动锁定到独立 \(k^\mu\) 的项，所以“随某一指定空间流运动”仍未被动力学实现。

## 9 对引力控制与宏观飞行的含义

弱引力、低速且质量非零时，静态试样的近似加速度为：

\[
\mathbf a\simeq-\nabla\Phi-c^2\nabla\ln m.\tag{31}
\]

第二项是一种候选状态场附加作用，而不是已经证明的排斥引力。形式上抵消重量仍需证明目标梯度能由装置产生，且要计算发生器的反作用。

宏观系统的总不变质量满足：

\[
M^2c^4=E_{\rm total}^2-c^2P_{\rm total}^2.\tag{32}
\]

对自由正能量零质量粒子集合：

\[
M^2c^4=2\sum_{i<j}E_iE_j(1-\cos\theta_{ij}).\tag{33}
\]

组成粒子的质量归零并不保证系统总质量归零。内部热运动、非平行动量、结合能和控制场能量都必须计入。要求所有成分同向光速传播，与普通刚性结构及内部时钟的存在之间具有实质性困难。因此本文的点粒子转换不能被直接推广为载人飞行。

## 10 六项任务的阶段性结论

| 任务 | 本阶段完成的研究 | 仍不包含的内容 |
|---|---|---|
| 1 变分与约束复核 | 主要场方程、粒子方程、应力能量与交换关系；七项符号恒等式 | 全部背景的稳定性及量子完备性 |
| 2 真空与边界 | 真空、平面壁、线性谱、三维无源静态泡的缩放障碍 | 有外源有限区域的所有动态解 |
| 3 电磁弱场响应 | 格林核、平板解、数值收敛、控制条件区别 | 强场转换阈值及真实装置参数 |
| 4 单粒子控制 | 时间循环与空间质量阱；能量交换核对 | 宏观物质结构保持 |
| 5 全系统动量 | 一般守恒证明、含动态电磁控制场的闭合低维算例 | 三维线圈与整机工程仿真 |
| 6 零质量过程 | 正则端点、双零极限、力关断、可控退出及锁定缺失 | 自动锁定和闭合强场进入退出机制 |

这些结果表明，四项核心命题可以在明确限定下组织研究，但不能被当作飞行可实现性的共同证明。特别是，本模型支持条件性的零质量传播，不支持从均匀质量下降直接推出静止物体启动，也没有提供隔离系统无反作用推进。

## 11 下一阶段目标与判定标准

**目标一：将参数与证据连接。**系统检索等效原理、第五力、光谱、时钟、电磁耦合及天体约束，确定是否存在尚未被排除的参数区间。若没有，则当前函数选择被淘汰，而不是继续将其推广到工程设计。

**目标二：求解闭合的强控制状态区。**把电源和边界纳入动力学，在有限空间内尝试建立并撤销 \(\chi\approx1\) 区域。验收条件包括稳定时间、输入和回收能量、边界运动及所有能量动量流，而不只检查区域中心的 \(\chi\) 值。

**目标三：研究方向锁定。**若仍要求物质由特定空间流携带，应提出显式相互作用，并检查守恒、稳定性及可区分的观测结果。若只获得与既有零测地线相同的描述，则不将其宣传为新增机制。

**目标四：从点粒子扩展到束缚态。**先计算双粒子或简单场论束缚态的总不变质量，再分析尺寸、结合能和内部信息是否可保持。没有这一结果，不进入载人飞行推论。

**目标五：由独特预测指导实验。**优先选择源极性、材料组成、空间衰减长度和独立惯性响应等相互关联的观测量，要求同一组参数同时解释结果。普通静电、离子风、磁力、热漂移和机械耦合须与候选效应定量比较。

## 12 结论

本研究将空间运动与质量耦合的概念框架推进到六项可复核计算。主要正面结果是：作用量能够一致地产生所列经典方程；平面状态壁及其线性谱可解析求得；电磁弱场响应和受控单粒子转换可以计算；非零动量轨迹能平滑通过零质量平台并恢复质量。

同样重要的限制是：无源有限三维静态状态泡不能自持；控制场纳入后试样动量增益由其他部分补偿；独立空间流没有自动锁定物质方向；宏观结构在零总不变质量状态下的存在尚未建立。本文因此将理论定位为具有明确假设与淘汰条件的经典候选模型，后续进展须由受约束参数、闭合状态控制及多粒子结构计算决定。

## 数据与方法可得性

所有数值结果均由实际运行的计算脚本生成，而非实验数据。数据、参数、软件版本、收敛诊断及图件保存在同一项目中：

- [阶段一至六计算研究报告](Research_阶段1至6/阶段1至6_计算研究报告.md)
- [数值结果与收敛摘要](Research_阶段1至6/数值结果摘要.md)
- [计算脚本](Research_阶段1至6/reproduce.py)
- [机器可读结果](Research_阶段1至6/data/results.json)
- [复算说明](Research_阶段1至6/README.md)

图 1—5 均由配套程序生成。主算例为无量纲、经典、低维或试验粒子计算，没有进行实验标定，也没有进行完整三维广义相对论数值演化。误差诊断描述数值实现，不代表理论的物理可信度。

## 参考文献

[1] 张祥前．统一场论．本项目提供的原始文稿，版本及发表年份未在本研究中核定。[本地文件](Original/4统一场论.docx)。

[2] 人工场，飞碟技术制造指南，第九版．本项目提供的技术笔记，作者未在本研究中核定。[本地文件](Original/22人工场，飞碟技术制造指90度-第九版_.docx)。

[3] 空间运动与物质耦合理论研究框架 V1.0．项目工作稿，2026 年 10 月 22 日。[本地文件](空间运动与物质耦合理论_研究框架V1.0.md)。

[4] Carroll, S. M. Lecture Notes on General Relativity. 1997. [arXiv:gr-qc/9712019](https://arxiv.org/abs/gr-qc/9712019).

[5] Derrick, G. H. Comments on Nonlinear Wave Equations as Models for Elementary Particles. *Journal of Mathematical Physics*, 5, 1252–1254, 1964. [doi:10.1063/1.1704233](https://doi.org/10.1063/1.1704233).

文献 [1,2] 提供研究动机，不作为本模型的新物理证据；[4,5] 提供标准理论背景。本文没有声称这些参考资料支持人工反引力或光速飞行装置。
