# 空间运动与物质耦合模型阶段一至六计算研究报告

作者：UFT Labs  
日期：2026 年 10 月 22 日  
对象：[理论框架 V1.0](../空间运动与物质耦合理论_研究框架V1.0.md) 第 22 节的第 1—6 项计算任务。

## 研究范围与完成标准

本次研究采用 V1.0 已明确提出的经典有效作用量，完成该模型内可解析的六项问题，并为非线性与控制问题提供可复算的低维数值算例。“完成”指得到明确推导、可重复的计算及适用边界，不指完成未知参数的实验测定、任意三维装置设计、全部非线性背景的稳定性证明或宏观载人飞行理论。

计算中的长度、时间、能量和场幅均为无量纲示例。未由实验拟合任何参数，不能直接换算为装置的伏特、安培、赫兹或牛顿。

配套文件：

- [计算脚本](reproduce.py)
- [依赖版本](requirements.txt)
- [完整数值结果](data/results.json)
- [数值表格与收敛摘要](数值结果摘要.md)
- `data/*.csv`：原始计算序列。
- `figures/*.png` 与 `figures/*.svg`：计算图。

## 1 作用量与变分复核

### 1.1 约定

取号差 \((-+++)\)，自然单位 \(c=\hbar=1\)，并定义：

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

主体数值算例取 \(Y=0\)。所有粒子先取中性；这既减少普通电磁力混入，也限定了结果不能直接推广到电子、导电材料或完整核结构。

\[
S=\int\sqrt{-g}\,d^4x\left[
\frac{M_{\rm Pl}^2}{2}R-\frac Z4F^2
-\frac{f^2}{2}(\nabla\chi)^2-U-\frac Y4F\widetilde F
\right]+S_m+S_{\rm src}+S_k,
\]

\[
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].
\]

这里 \(f=f_\chi>0\)，\(e\) 是辅助变量。引力变分取紧支撑扰动，或加入与所用边界条件相容的引力边界项；不能将含二阶导数的曲率项当作完全没有边界问题。

### 1.2 状态场变分

\[
\delta S_{\chi}=\int\sqrt{-g}\,d^4x\left[
f^2\Box\chi-U'-\frac{Z'}4F^2-\frac{Y'}4F\widetilde F-J_\chi
\right]\delta\chi,
\]

其中：

\[
J_\chi=-\frac1{\sqrt{-g}}\frac{\delta(S_m+S_{\rm src})}{\delta\chi},
\qquad
J_{\chi,m}=\frac1{\sqrt{-g}}\sum_a\int e\,m_am'_a\delta^4(x-x_a)\,d\lambda.
\]

因此 V1.0 的状态场方程符号正确。对 \(m>0\) 取固有时参数，\(e=1/m\)，点粒子源化为 \(\int m'\delta^4d\tau\)，与对 \(-\int m(\chi)d\tau\) 直接变分一致。质量归零时则保留一阶形式，不能使用 \(e=1/m\)。

### 1.3 电磁变分

由 \(\delta F_{\mu\nu}=2\nabla_{[\mu}\delta A_{\nu]}\) 和反对称性：

\[
\delta S_{\rm EM}
=\int\sqrt{-g}\,\nabla_\mu(ZF^{\mu\nu}+Y\widetilde F^{\mu\nu})\delta A_\nu\,d^4x.
\]

沿用 V1.0 的源约定 \(j^\nu=-\delta S_{\rm src}/(\sqrt{-g}\delta A_\nu)\)，得到：

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

源电流的正负号与源作用量的定义配套，不能只更换一边。

### 1.4 物质变分与哈密顿量

分别对 \(p_\mu,e,x^\mu\) 变分：

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

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

局部平直背景、正能量支：

\[
H=\sqrt{p^2+m^2},\qquad
\dot{\mathbf x}=\frac{\mathbf p}{H},\qquad
\dot{\mathbf p}=-\frac mH\nabla m,\qquad
\dot H=\frac mH\partial_t m.
\]

后一式是哈密顿量沿解的全导数，右侧是其显式时间依赖；不是将质量沿轨迹的全导数代入。

### 1.5 应力能量与源交换

取 \(T_{\mu\nu}=-2\delta S/(\sqrt{-g}\delta g^{\mu\nu})\)。得到 V1.0 中的标量场和电磁场应力能量；点粒子项为：

\[
T^{\mu\nu}_m=\frac1{\sqrt{-g}}\sum_a\int e\,p^\mu p^\nu\delta^4(x-x_a)\,d\lambda.
\]

\(YF\widetilde F\) 在本写法中没有独立显式度规应力项，因为 \(\sqrt{-g}\varepsilon^{\mu\nu\alpha\beta}\) 是相应的符号密度。但它仍改变场的解和场间交换。

完整系统满足 \(\nabla_\mu T^{\mu\nu}_{\rm total}=0\)。在 \(Y=0\)、无电磁源的简化情况下：

\[
\nabla_\mu T^{\mu\nu}_{\chi}
=\left(\frac{Z'}4F^2+J_\chi\right)\nabla^\nu\chi,
\]

\[
\nabla_\mu T^{\mu\nu}_{\rm EM}
=-\frac{Z'}4F^2\nabla^\nu\chi,
\qquad
\nabla_\mu T^{\mu\nu}_m=-J_\chi\nabla^\nu\chi.
\]

三者相消。若把驱动电源当作固定外场而不计入源，则试验粒子子系统可以不守恒，但不能据此宣布整体违反守恒。

### 1.6 约束与量纲检查

- \(p^2+m^2=0\) 是重参数不变性的约束；固定坐标时间后选择正能量支。
- \(Z>0,f^2>0\) 保证所写经典动能符号正确，但不足以证明所有背景稳定。
- \([\chi]=[Z]=[Y]=0,[f]=[\Lambda]=[m]=1,[F]=2,[U]=4\)，每项拉格朗日密度均为能量四次方。
- 独立零流作用量 \(S_k=-\frac12\int\sqrt{-g}\rho(\nabla\vartheta)^2\) 给出零约束、连续性方程和零测地线。其梯度形式无旋，不是任意螺旋流。

脚本独立核验七项标量与哈密顿恒等式。完整张量变分是上述手工推导，不能把七项符号检查描述成对全部引力方程的自动证明。

## 2 真空、扰动与有限区域边界

### 2.1 均匀真空

\[
U'=2\Lambda^4\chi(1-\chi)(1-2\chi),\qquad
U''=2\Lambda^4(1-6\chi+6\chi^2).
\]

\(\chi=0,1\) 是简并局部极小值，\(\chi=1/2\) 是势垒顶。真空小扰动满足：

\[
(\Box-\mu^2)\delta\chi=0,\qquad
\mu^2=2\Lambda^4/f^2,
\qquad \omega^2=k^2+\mu^2.
\]

因此 \(\chi\) 不是处处以光速传播的无质量场。四项核心命题中的“空间光速运动”只能在本模型中指零轨迹扇区；若要求全部自由度都以光速传播，则与该双稳势模型不相容。

### 2.2 平面边界精确解

静态方程 \(f^2\chi''=U'\) 在 \(\chi(-\infty)=0,\chi(+\infty)=1\) 下有：

\[
\chi_w(z)=\frac1{1+e^{-\mu(z-z_0)}}.
\]

一阶关系与单位面积能量为：

\[
\frac{f^2}{2}(\chi'_w)^2=U(\chi_w),
\quad
\sigma=\int_{-\infty}^{\infty}\left[\frac{f^2}{2}(\chi'_w)^2+U\right]dz
=\frac{\sqrt2}{6}f\Lambda^2.
\]

数值边值求解采用 \([-12,12]\)、精确尾部 Dirichlet 值、独立初猜和 `solve_bvp`，不把解析解直接当作数值结果。误差与张力核对见数值摘要。

### 2.3 平面壁线性谱

对 \(\chi=\chi_w+\eta(z)e^{-i\omega t}\)：

\[
\left[-\partial_z^2+\mu^2\left(1-\frac32\operatorname{sech}^2\frac{\mu z}{2}\right)\right]\eta
=\omega^2\eta.
\]

离散谱包含平移零模 \(\omega^2=0\) 和形状模 \(\omega^2=3\mu^2/4\)，连续谱阈值为 \(\mu^2\)。加入沿壁波数后增加非负 \(k_\parallel^2\)。这证明该孤立平面壁在线性扰动层面无负模，不是全非线性、含物质和电磁背景的普遍稳定性证明。

非负性也可直接由因子化得到：令 \(Q=\partial_z+\mu\tanh(\mu z/2)\)，则扰动算符为 \(Q^\dagger Q\)。在平方可积、边界项消失的扰动上，\(\langle\eta,Q^\dagger Q\eta\rangle=\|Q\eta\|^2\geq0\)；零模正比于 \(\operatorname{sech}^2(\mu z/2)\)。

有限差分的最低本征值会出现趋于零的小负值；网格加密后按离散误差缩小，不能把它误判为物理不稳定。

### 2.4 三维孤立有限泡不能静态自持

对无外源、\(U\geq0\) 的三维静态有限能量构型，令 \(\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
\]

对非平凡构型无法满足驻值条件。这是本势能下的 Derrick 缩放障碍，适用于正则有限能量、无外部支持的静态标量区域。

薄壁半径 \(R\) 下：

\[
E(R)=4\pi\sigma R^2+\frac{4\pi}{3}\Delta u\,R^3.
\]

两个真空简并时 \(\Delta u=0\)，能量随半径增加，泡倾向收缩。若 \(\Delta u<0\)，临界半径 \(2\sigma/|\Delta u|\) 是势垒极大值而非稳定极小值。电磁压力、外源、荷守恒或时间依赖可能改变结论，但必须重新计算。

因此第二项任务的结果是：获得稳定平面壁及其谱，同时证明当前无源模型不存在静态自持的有限三维“零质量舱”。

## 3 静态电磁弱场响应

### 3.1 线性化与格林函数

在 \(\chi_0=0\)、物质源可忽略、外部电磁场较弱时，定义：

\[
S(\mathbf x)=\frac\beta2(B^2-E^2)-\alpha\mathbf E\cdot\mathbf B.
\]

状态方程给出：

\[
(\nabla^2-\mu^2)\delta\chi=\frac S{f^2}.
\]

三维无限空间解为：

\[
\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'.
\]

对 \(\beta>0,\alpha=0\)，电场主导的负源提高 \(\chi\)，磁场主导的正源降低 \(\chi\)。作用长度受 \(\mu^{-1}\) 控制。引力响应还须计算物质和度规，不能把 \(\delta\chi\) 自身称为引力加速度。

### 3.2 无限平板局域源

若仅依赖 \(z\)，且 \(S=S_0\) 于 \(|z|<a\)，其余为零，则：

\[
\delta\chi(z)=-\frac{S_0}{f^2\mu^2}
\begin{cases}
1-e^{-\mu a}\cosh\mu z,& |z|<a,\\
e^{-\mu|z|}\sinh\mu a,& |z|\geq a.
\end{cases}
\]

函数与一阶导数在边界连续。它说明场效应取决于源厚度及状态场传播尺度，而不只取决于场强。

### 3.3 数值验证

为避免阶跃源带来的边界离散混淆，数值选取 \(S/f^2=-0.02e^{-z^2/2}\)，在 \([-10,10]\) 上使用二阶中心差分。边界由独立无限线格林积分提供，再以 400、800、1600 个内点检验收敛。正反源对应相反线性响应，误差见摘要。

### 3.4 固定电荷与固定电压必须区分

\(Y=0\) 时，静电位移 \(\mathbf D=Z\mathbf E\)。固定自由电荷的局部电能为 \(D^2/(2Z)\)，其对 \(\chi\) 的导数为负。固定电压的电源维持场，需要使用包含电源做功的适当势，例如均匀情况下的 \(-ZE^2/2\)。两个系综的一阶源在同一瞬间相容，但二阶反馈和稳定性不同。

把一个维持固定电压的电源排除在系统能量之外，会产生虚假的能量增益。

### 3.5 单束波、背景位移与阈值

单束理想真空平面电磁波有 \(F^2=F\widetilde F=0\)，不经本模型所选的直接耦合给 \(\chi\) 提供非零源。近场、驻波和叠加场可不同。

\(\chi=1\) 下虽然 \(U'=0\)，但持续纯电场一般仍有非零源。因此电场并不自动把系统停在 \(\chi=1\)；可能推过该点并使质量重新增加。需要关闭、平衡或反馈调节驱动。

线性小扰动解只在 \(|\delta\chi|\ll1\) 有效，不能把它外推到 \(\chi=1\) 得出转换电压或功率。

## 4 外部控制单粒子数值积分

### 4.1 控制函数

设 \(q(s)=6s^5-15s^4+10s^3\)，在 \(0\leq s\leq1\)。取：

\[
\chi(t)=\begin{cases}
q(t/10),&0\leq t<10,\\
1,&10\leq t\leq15,\\
1-q((t-15)/10),&15<t\leq25.
\end{cases}
\]

该函数在连接处一、二阶导数连续。取 \(m_0=1,p_0=0.3\)，积分 \((x,p,W)\)：

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

独立解析基准为 \(p=0.3\) 与 \(H(t)-H(0)=W(t)\)。使用 DOP853，以三档容差检验。

### 4.2 物理结果

初始 \(H_0=\sqrt{1.09}\)，零质量平台 \(H=0.3,v=1\)。进入平台时粒子能量降低，差额转入控制源；恢复质量时同样能量需要返还。完整循环结束后恢复初始速度和能量，但粒子位置已经改变。

此算例证明一个指定的平滑背景中存在数学上正则的非零动量进入及退出过程。它不证明电磁装置已经能生成该控制曲线，也不证明宏观结构可承受转换。

### 4.3 空间质量阱

另取静态 \(\chi(x)=0.35e^{-x^2/(2\times1.5^2)}\)，初态 \(x=-6,p=0.6\)。积分完整相对论方程：

\[
\dot x=p/H,\qquad \dot p=-m\partial_xm/H.
\]

粒子进入低质量区加速，离开时减速，总能量保持。有限起终点处的质量不完全相同，所以不能把末动量与初动量的微小差异全部称为数值误差；守恒误差单独计算。

## 5 动态控制源与整机动量

### 5.1 一般结论

在平直隔离系统中：

\[
\frac{dP^i_{\rm total}}{dt}=-\oint T^{ji}n_j\,dS.
\]

没有边界通量时总动量不变。若能量密度局域且总能量有限，能量中心的速度也由守恒总动量固定。这是由模型对称性得到的全系统约束，不依赖某次数值结果。

试样受力不是整机推力。若存在地球外场、辐射、外部空间流或其他外部源，则可交换动量，但相应反作用与通量必须明确。

### 5.2 闭合数值模型

为把控制场的反作用真正纳入积分，采用平直、周期边界的一维横向电磁场 \(A=A_y(x,t)\)，\(E_y=-\dot A,B_z=A_x\)，以及动态状态场和两个中性平滑试样。

这不是三维线圈或飞行器仿真，而是完整纳入控制能量与动量的受限模型。取 \(L=2\pi,f=\Lambda=1,\beta=0.4,Y=0\)。

\[
H_{\rm field}=\int_0^L\left[
\frac{\Pi_A^2}{2Z}+\frac Z2 A_x^2+
\frac{\pi_\chi^2}{2}+\frac{\chi_x^2}{2}+U(\chi)
\right]dx.
\]

为避免网格上的点源奇异，定义归一化周期高斯核 \(K_\sigma\)，\(\sigma=0.35\)：

\[
\bar\chi_a=\int_0^L K_\sigma(x-X_a)\chi(x)dx,
\quad m_a=m_{a0}(1-\bar\chi_a)^2,
\]

\[
H=H_{\rm field}+\sum_a\sqrt{p_a^2+m_a^2},\quad
m_{10}=0.03,\quad m_{20}=0.05.
\]

平滑化是明确的有限尺度正则化；它在这一坐标系保留空间平移对称性，但不能声称在核尺度保留完整三维洛伦兹协变性。宽度固定，未作点粒子 \(\sigma\to0\) 极限。

### 5.3 自洽演化方程

\[
\dot A=\Pi_A/Z,\qquad
\dot\Pi_A=\partial_x(ZA_x),\qquad
\dot\chi=\pi_\chi,
\]

\[
\dot\pi_\chi=\chi_{xx}-U'
+\frac{\beta}{2}(\Pi_A^2/Z-ZA_x^2)-\sum_a q_aK_\sigma(x-X_a),
\quad q_a=\frac{m_am'_a}{\sqrt{p_a^2+m_a^2}},
\]

\[
\dot X_a=\frac{p_a}{\sqrt{p_a^2+m_a^2}},\qquad
\dot p_a=-q_a\int K_\sigma(x-X_a)\chi_x(x)dx.
\]

总动量为：

\[
P=\sum_a p_a-\int_0^L(\Pi_AA_x+\pi_\chi\chi_x)dx.
\]

此处场动量负号来自物理应力能量的 \(T^{0x}\) 与所用号差。试样增益与场损失应相消。

### 5.4 初态、收敛与结果

\[
\chi_0=0.015\cos(x-0.4),\quad\pi_{\chi0}=0,
\]

\[
A_0=0.15\sin x+0.05\sin(2x+0.2),\quad
\Pi_{A0}=-Z(\chi_0)\partial_xA_0,
\]

\[
(X_1,X_2)=(1,4),\qquad(p_1,p_2)=(0,0).
\]

初始电磁场具有非零总动量，系统不是从零总动量凭空启动。使用 Fourier 配点法及 DOP853，在 \(N=64,128,256\) 和递紧容差下积分到 \(t=10\)。周期边界无净外部通量。

结果见数值摘要：试样获得动量，电磁及状态场承担补偿；总能量、总动量误差随数值控制保持很小。连续方程的守恒是解析结果，数值只验证实现和具体交换过程。

这一模型不产生隔离整机的持续无反作用推力。它不能排除与外界交换动量的推进，但此类推进必须另外给出外界模型。

## 6 零质量端点、方向锁定与退出

### 6.1 正则零质量端点

令 \(\delta=1-\chi\)，则：

\[
m=m_0\delta^2,\quad
H=\sqrt{p^2+m_0^2\delta^4},\quad
\dot{\mathbf p}=\frac{2m_0^2\delta^3}{H}\nabla\chi.
\]

当 \(p\neq0\) 且背景光滑，\(H\geq|p|>0\)，方程在 \(\delta=0\) 正则；速度趋于光速，质量梯度力趋于零。不能使用 \(\nabla\ln m\) 在此处产生的形式发散。

若动量方向为 \(\hat{\mathbf p}\)：

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

在有限 \(p\) 和有限场梯度下，接近零质量状态的该类横向控制逐渐变弱，而不是无限增强。

### 6.2 双零极限不是唯一速度

若同时令 \(p=a\delta^2\)，则：

\[
v\longrightarrow\frac{a}{\sqrt{a^2+m_0^2}}.
\]

该值可以小于 1；这是趋向 \(p=m=H=0\) 的路径极限，端点本身不是非零能量的零质量传播粒子。故“零质量必光速”必须保留非零能量条件。

### 6.3 自动方向锁定未出现

V1.0 的 \(S_k\) 与物质只有通过共同度规的间接联系，没有强制 \(p^\mu\parallel k^\mu\) 的直接耦合。在平直背景，物质可沿其自身初始方向光速传播，而零流沿另一方向传播。

将二者画成平行不能替代作用量中的机制。“平行即正比”是成立时的关系，不是把任意两个独立矢量变平行的动力。

### 6.4 退出是可能的，但驱动尚未工程闭合

第四项算例给出明确反例：保持 \(p\neq0\)，先将 \(\chi\) 设为 1，再平滑移开，物质即可重新获得质量。零质量处 \(m'=0\) 只意味着该瞬间的源和力为零，并不强迫外部场永久停在 \(\chi=1\)。

因此必须修正 V1.0 可能引起的“零质量是不可退出状态”的读法。当前模型允许受控退出，尚未证明实际电磁装置能实现整段控制曲线及其边界。

## 阶段结论与后续目标

本次六项计算给出一组正面结果和一组约束：作用量的主要变分成立；双稳状态场有可解析平面壁；电磁弱场响应可计算；非零动量粒子可在指定背景中通过零质量平台并恢复；控制场纳入后总动量守恒；有限无源静态三维泡、自动方向锁定和宏观飞行没有被当前模型实现。

后续目标应集中在：实验参数限制、实际边界条件下的动态状态区、带反作用的方向锁定机制、多粒子束缚态，以及在这些限制内是否还存在可行的外部动量交换飞行方案。

## 参考资料

1. UFT Labs 项目研究稿，《空间运动与物质耦合理论研究框架 V1.0》，2026 年 10 月 22 日。本次计算的直接模型来源。
2. S. M. Carroll, *Lecture Notes on General Relativity*, 1997, [arXiv:gr-qc/9712019](https://arxiv.org/abs/gr-qc/9712019)。用于几何与协变守恒的标准背景；不是本候选模型的实验支持。
3. G. H. Derrick, *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)。本稿对所用非负双稳势自行写出缩放论证。

参考资料的书目信息经过检索核对；本阶段没有完成候选模型全部实验约束的系统综述。