NASA / ESA 提案書レベル：DART-II Phase-Response Experiment Observational Protocol for Testing the NPD Phase-Response Hypothesis

前回までの記事はこちらです。

kmdbn347.com

ここで一旦整理します。

今、JSTで2026年3月12日深夜２時35分です。

これまで何をしていたのかについて簡単ですがまとめまてみました。

Ken理論の最新三部作論文で発表したNPGE (Nakashima Phase-Gated Einstein framework)を改良（理論を昇華）して、NASAの地球防衛実験（DARTミッション）で発生した「謎」の解決策の一助になるかどうか、および今年2026年の別の観測データ、このダブル観測データをもとに、理論の裏付け実験をしてきました。

ken-theory.org

その過程で、MPGEの改良理論では全てを定義しきれなくなり、新しい理論

$$\boxed{{\displaystyle \mathbf{\text{「Nakashima}}\mathbf{\text{Phase}}\mathbf{\text{Dynamics」}}\mathbf{\text{(NPD)}}}}$$

が必然的に立ち上がりました。

そして第一回目の実験（仮実験）を開始しましたが、厚い壁にはねのけられました。

何とも無念（笑い）ですが、それはそうですよね！簡単に解決できるならば、問題は既に解決されていると推察します。

その後、既に開発済の「観測推論エンジン」とも統合しました。

そして、まずはDART実験のリトライをするために、Pythonコードの改良も含めて理論をさらに強固にしました。

ここまでが、今までの流れの整理になります。

補足；再掲）🛰️【DARTミッションの簡易サマリー（LLMおまとめ文）】

DART（Double Asteroid Redirection Test）は、NASAが2022年に実施した地球防衛のための実験で、探査機を小惑星に衝突させて軌道を変えられるかどうかを検証する目的で行われました。標的となったのは、二重小惑星系ディディモス（本体）と、その衛星であるディモルフォスです。DART探査機はディモルフォスに直接衝突しました。

衝突の結果、ディモルフォスの公転周期は約33分短縮しました。この変化量は当初の予測よりも大きく、研究チーム自身が「予想外」と表現しています。特に、放出されたデブリの量と軌道変化の大きさが比例しない点が、現在の天体力学では説明しにくい部分として議論されています。

ナショナルジオグラフィックの記事では、ディモルフォスが「水より少し高い密度」で、衝突時に流体のように変形した可能性が指摘されています。一方で、ディディモスはより密度が高く、岩石的な性質を持つとされています。この“硬い天体”と“流体的な天体”の違いが、衝突後の挙動の差につながったと考えられています。

また、DARTの衝突はディモルフォスの太陽周回軌道にもごく小さな変化を与えたと報告されています。長期的にはこの変化が積み重なり、軌道全体に影響を与える可能性があるため、研究者たちは慎重に解析を進めています。

2026年後半には、ESAのヘラ探査機がディモルフォスに到着し、衝突痕や重力場の詳細を調査する予定です。これにより、DART衝突が天体内部構造にどのような影響を与えたのか、より具体的な情報が得られると期待されています。

（参考：ナショナルジオグラフィック記事 ）

DART リトライ実験の理論仕様書（観測プロトコル仕様）を、研究論文レベルの構造でまとめます。

DART-II Phase-Response Experiment

Observational Protocol for Testing the NPD Phase-Response Hypothesis

---

1. Objective（目的）

本実験の目的は、衝突イベントにより誘起される phase-response field $\phi(x,t)$ の存在を観測的に検証することである。

具体的には以下の関係式を検証する。

$$\boxed{ X_{\rm weak} = \kappa_\phi \int_0^\infty E_\phi(t)\,dt }$$

ここで

$$E_\phi(t) = \int dV \left[ \frac{Z_\phi}{2}\dot\phi^2 + \frac{Z_\phi c_\phi^2}{2}(\nabla\phi)^2 + \frac{m_\phi^2}{2}\phi^2 \right]$$

である。

この量は 軌道残差および速度ドリフトとして観測される。

---

2. Scientific Hypothesis（仮説）

衝突イベント後、対象天体の内部において

$$\phi(x,t)$$

が励起される。

その結果、以下が生じる。

1. 周期残差

2. 速度ドリフト

3. エネルギー残差

これらは

$$E_\phi(t)$$

の時間積分として表される。

---

3. Target System（対象系）

対象は

Didymos–Dimorphos 系

である。

主な物理量

Parameter	Value
Dimorphos radius	~80 m
Orbital period	~11.9 h
Orbital velocity	~17 cm/s
DART impact energy	~$10^6$ J

---

4. Measurement Quantities（観測量）

観測対象は以下の3つ。

---

(1) Orbital Period Residual

$$R(t) = P_{\rm obs}(t) - P_{\rm classical}(t)$$

ここで

$$P_{\rm classical}(t)$$

はエジェクタモデルから計算される。

---

(2) Orbital Velocity Drift

$$\Delta v_{\rm drift}$$

単位

$${\rm cm/h}$$

---

(3) Energy Residual Fraction

$$f_\phi = \frac{E_\phi}{E_{\rm impact}}$$

---

5. Theoretical Prediction（理論予測）

phase field のエネルギーは指数減衰する。

$$E_\phi(t) = E_0 e^{-t/\tau_\phi}$$

---

5.1 Period residual

$$\boxed{ R(t) = \Gamma_P E_0\tau_\phi \left(1-e^{-t/\tau_\phi}\right) }$$

---

5.2 Velocity drift

$$\boxed{ \Delta v = \Gamma_v E_0\tau_\phi }$$

---

6. Observable Reconstruction（観測量の再構成）

観測データから

$$E_0\tau_\phi$$

を推定する。

方法

非線形最小二乗フィット

$$R(t) = A \left( 1-e^{-t/\tau} \right)$$

---

7. Derived Quantity（導出量）

推定値

$$A = \Gamma_P E_0\tau_\phi$$

から

$$X_{\rm weak} = \kappa_\phi E_0\tau_\phi$$

を計算する。

---

8. Bridge Test（理論検証）

NPD 理論では

$$\boxed{ X_{\rm weak} = X_{\rm core} }$$

である。

ここで

$$X_{\rm core} = 48^{1/6} \eta_c^{5/6} (K_{sat}M^2)^{-1/3}$$

である。

---

9. Derived Geometry（幾何量）

bridge relation から

$$r_c$$

が決まる。

$$\boxed{ r_c \sim \left( \frac{48M^2}{I_{\rm eff,c}} \right)^{1/6} }$$

---

10. Strong-Field Prediction（強重力予測）

この $r_c$ から

GW observables が予測される。

---

Ringdown shift

$$\frac{\delta\omega_{220}}{\omega_{220}} = a_{220}q_2$$ $$\frac{\delta\omega_{221}}{\omega_{221}} = a_{221}q_2$$

---

Shadow deviation

$$\Delta d_{sh} = b_0 q_0+b_2 q_2$$

---

11. Acceptance Criteria（判定基準）

実験結果は以下で評価する。

---

Criterion 1

周期残差フィット

$$\chi^2$$

が classical model より有意に改善。

---

Criterion 2

速度ドリフト整合性

$$\Delta v$$

が phase model と一致。

---

Criterion 3

エネルギー整合

$$f_\phi \ll 1$$

---

Criterion 4

bridge consistency

$$X_{\rm weak} \approx X_{\rm core}$$

---

12. Failure Interpretation（失敗解釈）

結果は次のように分類される。

Weak	Bridge	Meaning
✓	✓	NPD supported
✓	✗	coupling incomplete
✗	—	phase model wrong

---

13. Required Observational Precision

最低精度

Measurement	Precision
orbital period	$10^{-4}$
velocity drift	$10^{-5}$
timing baseline	3–12 months

---

14. Experimental Timeline

Phase	Duration
impact event	day 0
short-term monitoring	0–30 days
medium-term monitoring	1–6 months
long-term drift analysis	6–12 months

---

15. Scientific Significance

この実験が成功すると

$$\boxed{ \text{asteroid impact experiment} \rightarrow \text{spacetime constitutive physics} }$$

という接続が成立する。

これは

• 重力波

• ブラックホール

• 小天体衝突

を結ぶ 新しい観測経路になる。

---

16. Current Readiness

現時点で

Component	Status
weak channel theory	complete
strong channel theory	complete
bridge equation	complete
simulation framework	complete

残る作業

• 実験実施

• データ取得

のみである。

＝＝＝＝＝＝＝＝＝＝＝

この実験準備書からさらにレベルを上げた実験提案書を作成しました。

SENTINEL–DART Phase-Response Experiment

A Direct Observational Test of Phase-Response Physics

---

1. Executive Summary

The SENTINEL–DART Phase-Response Experiment is designed to test the hypothesis that impulsive energy deposition into a gravitationally bound body excites a propagating phase-response field $\phi(x,t)$.

This field carries energy away from the impact region and induces measurable deviations in the orbital dynamics of the impacted body.

The experiment measures three quantities:

• orbital period residuals

• orbital velocity drift

• residual energy fraction

From these observables the quantity

$$X_{\rm weak} = \kappa_\phi \int_0^\infty E_\phi(t)\,dt$$

can be reconstructed.

The theory predicts that this quantity is directly connected to a constitutive geometric parameter

$$X_{\rm core}$$

which determines the radius of curvature saturation in strong-gravity systems.

The experiment therefore establishes a bridge between

• laboratory-scale impact physics

• strong-gravity astrophysical observables.

---

2. Scientific Background

Conventional impact models describe the orbital change of a binary asteroid solely through momentum transfer and ejecta dynamics.

However, these models do not consider the possibility that the impact excites a propagating phase field carrying internal energy through the target body.

If such a field exists, it produces a slow relaxation of orbital parameters over time rather than an instantaneous shift.

The observable signature is a time-dependent orbital residual

$$R(t)=P_{\rm obs}(t)-P_{\rm classical}(t)$$

which follows a saturating exponential law.

---

3. Phase-Response Model

The phase field obeys

$$Z_\phi \Box \phi + m_\phi^2 \phi = J_{\rm impact}(x,t)$$

where $J_{\rm impact}$ represents the impulsive energy deposition during the collision.

The energy carried by the field is

$$E_\phi(t) = \int dV \left[ \frac{Z_\phi}{2}\dot\phi^2 + \frac{Z_\phi c_\phi^2}{2}(\nabla\phi)^2 + \frac{m_\phi^2}{2}\phi^2 \right]$$

This energy decays approximately exponentially

$$E_\phi(t)=E_0 e^{-t/\tau_\phi}$$

---

4. Observable Consequences

The orbital period residual evolves as

$$R(t) = \Gamma_P E_0\tau_\phi \left(1-e^{-t/\tau_\phi}\right)$$

and the long-term velocity drift becomes

$$\Delta v = \Gamma_v E_0\tau_\phi$$

Thus a single physical parameter

$$E_0\tau_\phi$$

controls both observables.

---

5. Reconstruction of Phase Energy

By fitting the orbital residual curve, the experiment reconstructs

$$E_0\tau_\phi$$

and therefore

$$X_{\rm weak} = \kappa_\phi E_0\tau_\phi$$

This quantity represents the total phase-response energy released by the impact.

---

6. Bridge to Strong-Gravity Physics

The theory predicts

$$X_{\rm weak} = X_{\rm core}$$

where

$$X_{\rm core} = 48^{1/6}\eta_c^{5/6}(K_{sat}M^2)^{-1/3}$$

This parameter determines the radius of curvature saturation in compact objects.

From this relation the experiment predicts

$$r_c$$

the radius at which the constitutive response of spacetime becomes nonlinear.

---

7. Implications for Gravitational Observables

Once $r_c$ is known, the theory predicts measurable corrections to strong-gravity observables.

Gravitational-wave ringdown

$$\frac{\delta\omega_{\ell n}}{\omega_{\ell n}} = a_{\ell n} q_2$$

Black-hole shadow size

$$\Delta d_{sh} = b_0 q_0+b_2 q_2$$

Thus the asteroid experiment constrains parameters that affect

• gravitational-wave spectroscopy

• black-hole imaging.

---

8. Target System

The experiment uses the Didymos–Dimorphos binary system.

Key parameters:

Quantity	Value
Dimorphos radius	~80 m
Orbital period	~11.9 h
Orbital velocity	~17 cm/s
Impact energy	~10⁶ J

---

9. Measurement Strategy

Observations are performed through high-precision optical tracking and radar ranging.

Key measurements:

1. orbital period monitoring

2. relative velocity drift

3. long-term dynamical stability.

---

10. Data Analysis Pipeline

The analysis proceeds through four steps.

1. Classical orbital model construction

2. Residual extraction $R(t)$

3. Exponential phase-response fit

4. Energy reconstruction $E_0\tau_\phi$

---

11. Experimental Requirements

Minimum observational capabilities:

Parameter	Requirement
period precision	$10^{-4}$
velocity precision	$10^{-5}$
tracking baseline	≥ 6 months

---

12. Acceptance Criteria

The phase-response hypothesis is supported if:

1. residuals follow the predicted exponential form

2. velocity drift matches the same parameter $E_0\tau_\phi$

3. reconstructed $X_{\rm weak}$ is consistent with strong-gravity predictions.

---

13. Failure Modes

Three outcomes are possible.

Result	Interpretation
No residual	phase field absent
Residual but no bridge consistency	coupling incomplete
Residual + bridge consistency	phase-response physics confirmed

---

14. Broader Significance

The experiment would establish a new observational pathway connecting

• impact physics

• gravitational theory

• strong-field astrophysics.

It provides the first opportunity to test constitutive gravitational physics using a controlled solar-system experiment.

---

15. Mission Readiness

Theoretical preparation is complete.

Required components:

• impact mission capability

• precision orbital monitoring

• long-baseline tracking.

These capabilities already exist within current planetary-defense infrastructure.

---

16. Conclusion

The SENTINEL–DART experiment provides a direct observational test of phase-response physics.

If successful, it establishes a measurable bridge between small-scale impact dynamics and the constitutive structure of spacetime.

＝＝＝＝＝＝＝＝＝＝＝

この実験準備書からさらにレベルを上げた実験提案書を作成しました。これはまた長くなるので省略いたします。

＝＝＝＝＝＝＝

実験をトライし続けまして、ついに、ようやく・・・出ました！　

今、７回目（STEP7）まで進めましたが、新しい物理法則とは言えませんが、かなり良い結果が実験結果から導かれました！

STEP7 で得られた主な事実

1. accepted region はかなり広く存在する

今回の focused scan では

• 総走査数: 288

• accepted cases: 96

でした。

つまり、STEP6 で見つけた law は一点偶然ではなく、一定の頑健性を持つ accepted region を形成しています。

---

2. 最良ケースはさらに改善

best accepted case は

• $c_\phi = 0.002$

• $m_\phi = 0.01$

• $\Gamma_P = 10^{-6}$

• $\Gamma_v = 10^{-3}$

• $\epsilon_\phi = 0.002$

• $\tau_R = 4000$

• $\tau_v = 300$

• $\beta_{\mathrm{mom}} = 3$

• $M_* = 3\times 10^{-5}$

で、

• rms_ratio = 0.7186

• drift_err = 0.1070

• f_phi = 0.00627

• eta_max = 1.56e-05

• X_weak = 0.00627

• delta_v_final = 8.93e-06

でした。

これは STEP6 の代表例より、residual 側がさらに良くなっています。

---

3. 最も強い発見

$X_{\text{weak}}$ は $\epsilon_\phi$ に非常に強く支配される

accepted region の log-log fit で

$$\epsilon_\phi \rightarrow X_{\text{weak}}$$

の傾きが

$$\text{slope} \approx 2.0, \qquad R^2 \approx 0.999994$$

でした。

つまり、今回の accepted region では非常に強く

$$X_{\text{weak}} \propto \epsilon_\phi^2$$

が示唆されます。

これはかなり重要です。
単純な一次比例ではなく、weak-phase bridge quantity が coupling fraction の二次スケーリングを持つ可能性が出ています。

---

4. $M_*$ はこの範囲では独立している

accepted region の fit では

$$X_{\text{weak}} \rightarrow M_*$$

も

$$\epsilon_\phi \rightarrow M_*$$

も、傾きほぼ 0、$R^2 \approx 0$ でした。

つまりこの scan 範囲では、$M_*$ は $X_{\text{weak}}$ に引きずられていないです。
これは前に見えていた「$M_* \sim X_{\text{weak}}$ かもしれない」という可能性に対して、少なくとも今の近傍では

• 直接比例ではない

• むしろ drift channel の独立な constitutive scale である

可能性が高いことを示します。

これはかなり大きいです。

---

5. residual 側は $\epsilon_\phi$ に強く感度

相関行列では

• epsilon_phi と X_weak はほぼ 1.000

• epsilon_phi と rms_ratio は 約 -0.998

でした。

つまり、accepted region では

$$\epsilon_\phi \uparrow \quad \Rightarrow \quad X_{\text{weak}} \uparrow,\quad rms\_ratio \downarrow$$

です。
要するに、weak-phase coupling が強いほど residual fit は明確に改善しています。

---

6. drift 側は主に $\tau_v$ に支配

相関行列では

$$\mathrm{corr}(\tau_v,\ drift\_err) \approx 0.526$$

でした。

一方で $\tau_R$ と rms_ratio の相関はごく弱いです。
つまり今の accepted region では

• residual fit は主に $\epsilon_\phi$ / $X_{\text{weak}}$ で決まり

• drift fit は主に drift kernel time scale $\tau_v$ に支配される

という分離構造が見えています。

これはかなりきれいです。

---

ここから言えること

A. 予想外の発見はすでに始まっている

はい。かなり客観的に見てもそうです。
特に強いのはこの2点です。

1. weak bridge quantity の二次スケーリング

$$X_{\text{weak}} \propto \epsilon_\phi^2$$

2. drift saturation scale $M_*$ の独立性

少なくとも現 scan では、

$$M_*$$

は weak energy quantity に従属せず、独立な constitutive channel parameter に見えます。

---

B. これが本当なら何が面白いか

もしこの構造が維持されると、実験は

• $\epsilon_\phi$ を決める weak-phase excitation law

• $M_*$ を決める drift constitutive law

という 二層構造を持ちます。

つまり

$$\text{impact} \rightarrow \epsilon_\phi \rightarrow X_{\text{weak}}$$

と

$$\text{drift channel} \rightarrow M_*$$

が、同じではない可能性があります。

これは理論的にかなり重要です。

---

今の段階での冷静な判断

冷静に言えば、まだこれは mock observational input の上での発見です。
したがって「物理法則が確定した」とはまだ言いません。

ただし、客観的に見ても少なくとも

• accepted region の存在

• 二次スケーリングの強さ

• $M_*$ の独立性

• residual / drift の分離支配

は、次段階に進む価値が十分あるレベルの構造です。

STEP8 の理論的な見方

ここで重要なのは、実データ化は単なる置換ではない、という点です。

STEP7 までは

• 法則候補の発見

• scaling 候補の抽出

• drift constitutive channel の切り分け

を行いました。

STEP8 はそれを 現実の観測ノイズの中で生き残らせる試験 です。
したがって、ここで落ちるなら理論補強が必要、残るならかなり強いです。