计算与复现代码

01可下载的演示计算

本目录汇集已经存在的脚本、复现模板和研究试运行。先确定要核对的物理关系，再选择相应材料；这些项目使用的输入和完成的计算各不相同。

LVE 演示脚本 ASD_LVE_demo.py只使用 Python 标准库。它计算 $\beta=v/c$、$\gamma=(1-\beta^2)^{-1/2}$、横向尺度比 $L'_\perp/L_\perp=\gamma$，以及同速平行、横向分离构型中的力比 $F/F_0=1-\beta^2$ 和磁作用差额 $F_m/F_0=\beta^2$。长度关系来自 LVE；力比还使用核外真值作用距离和逆平方本构输入，见库仑作用与磁力推导。脚本中的 force_ratio 对应这里的真值物理力 $F_T/F_0$，不能直接当作参考坐标广义力或探针读数。速度域为 $0\leq\beta<1$；这个速度范围不把给定构型推广为任意源与探针。

下载脚本后运行：

python ASD_LVE_demo.py --beta 0.6 --csv lve-model.csv
python ASD_LVE_demo.py --beta 0

第一条命令的终端结果为 $\gamma=1.25$、$F/F_0=0.64$、$F_m/F_0=0.36$；第二条给出静止极限 $1,1,0$。--beta 选择终端显示的单点，CSV 则固定扫描 0.00 至 0.95、步长 0.01，共 96 行，字段为 beta,gamma,force_ratio,magnetic_fraction。这些数值由给定公式生成，不含测量记录、误差估计或一般动力学求解。

02方向钟试运行：几何与可识别性

方向钟试运行对应 ASD-RUN-CLOCK-001。复算代码已经包含在 ASD_AI_Research_Pack.zip中：tasks/clock-pilot.py、装置记录 tasks/clock-pilot.json 和冻结结果 tasks/clock-pilot-expected.json。这里的 JSON 是输入与未知量的记录；脚本在代码内固定下面的几何例子，并不读取 JSON 来设置装置参数。

设纬度 $40^\circ$，用户选定的背景轴赤纬 $30^\circ$，A、B 轴分别水平指北、指东，一个名义恒星日内采样 96 次。脚本使用未由 ASD 原子响应推导的几何模板 $Q=(\mathbf n\cdot\mathbf b)^2-1/3$，计算候选读数中三个系数对应的设计矩阵 $X=(Q_A,-Q_B,1)$。解压资料包后，在包的根目录运行：

python tasks/clock-pilot.py --verify tasks/clock-pilot-expected.json --output clock-pilot-output.json

输出 JSON 给出默认构型秩 3、赤纬 $0^\circ$ 时秩 2、两轴相同时秩 2；同时检查轴反向、一个周期的重复和同轴等系数的抵消。比较的背景轴和采样均为假设输入，秩只说明给定模板中参数能否分辨。代码目录中的方向钟结果只给出几何可识别性，尚未计算 ASD 频移振幅，也未使用观测数据。所选 $^{87}\mathrm{Sr}$ 跃迁、参考钟及传输链的响应系数仍是钟响应任务的下一份输入；不能用满秩代替这些物理计算或检测功效。

03实验比较基线与复现模板

ASD_Experimental_Reproducibility_Kit.zip包含 README.md、equations.md、models.py、run.py、verify.py、notebook.ipynb、input/、schemas/、expected/ 和 checksums.json。它用 Python 3.10 及以上的标准库运行比较模型，输入文件与单位见数据目录及包内说明。解压后运行：

python run.py --output results
python verify.py

results/summary.json 汇总结果，runtime.json记录实际 Python 版本；clock_geometry.csv、muon_bins.csv 和 charge_response.csv 分别保存几何因子、合成 cohort 的参考存活计数、标准点电荷场与假定高斯响应。这一包的钟几何有 97 个半小时采样点；它与上面的 96 点双钟秩试运行是不同计算，不能混用结果文件。expected/ 提供冻结数值，verify.py 核对分发文件的 SHA-256、数值和静止/周期/对称等极限。

μ 子部分使用恒能量、平直大气的标准存活公式和三个明确标为 SYNTHETIC 的 bin；高 γ 部分使用匀速点电荷的标准归一化横向场和单位面积高斯卷积，不含实际束团、晶体与电光传递函数。四份 schemas/*.csv 只有表头，是接收未来输入的格式模板；S301 没有轨道积分或拟合。详细前提、输出含义和缺项见复现资源与数据目录。包内 asd_output、predicted_difference 和 fit 均保留 null。

包内已有诊断图。需要重绘时，另行安装说明中指定的 matplotlib==3.10.8，运行 python figure.py --output results；数值计算不依赖它。Notebook 调用相同函数，需要读者自己的 Jupyter 环境。本包不包含真实观测比较或 Ota 实验复现。

04验证哪些结果

上述命令分别核对公式演示、几何秩与比较基线。2026-10-03 使用 Python 3.12.14 复算了实际分发材料：LVE 的 96 行 CSV、双钟秩和三个比较模型的 JSON/CSV 与所附冻结数值相符。运行结果应结合输入、代码、冻结输出和校验清单一起保存；数值一致只说明这份计算可核对，不会补出尚未求得的 ASD 装置响应。

05完整复现仓库的要求

本目录沿用现有下载包和登记，当前没有经确认的公共计算仓库链接。进一步计算应明确问题、方程、边界条件、真实输入来源、环境、输出和核对方法，并把数值误差与测量不确定度分开。方向钟先求响应系数，μ 子先闭合路径到衰变概率的映射，高 γ 场先闭合源与探针的联合响应，S301 先核实对象、历元、参考系与协方差；相应入口保留在预测登记与开放问题。这些明确缺项完成后，再登记真实数据比较。