--- { "path": "/en/research/code", "title": "Demonstration calculations", "category": "Research resources", "summary": "Existing LVE demonstration, clock geometry pilot and reproducibility templates, with inputs, commands, outputs and unfinished work.", "related": [ "/en/theory/lve", "/en/research/data", "/en/collaborate/how-to-contribute" ], "id": "ASD-DOC-RESEARCH-CODE-EN", "template": "document", "version": "1.1", "updated": "2026-10-03", "canonical": false, "source_tier": "website-guide", "depends_on": [], "open_problems": [], "equations": [], "sources": [ { "label": "ASD 建站蓝图与页面模板（用户提供）", "kind": "supplied-blueprint" } ], "theory_status": "FOUNDATION", "evidence_status": "DOCUMENT", "title_en": "Demonstration calculations", "locale": "en", "translation_path": "/research/code", "translation_of": "/research/code", "translation_version": "1.1", "canonical_url": "https://asd-physics.org/en/research/code", "section_anchors": [ "s0", "clock-pilot", "comparison-kit", "s1", "s2" ] } --- ## Downloadable demonstration This catalog collects existing scripts, reproducibility templates and research pilots. Choose the material for the physical relation you want to check; their inputs and completed calculations differ. The [LVE demonstration script ASD_LVE_demo.py](/downloads/ASD_LVE_demo.py) uses only the Python standard library. It calculates $\beta=v/c$, $\gamma=(1-\beta^2)^{-1/2}$, the transverse scale ratio $L'_\perp/L_\perp=\gamma$, and $F/F_0=1-\beta^2$ with magnetic deficit $F_m/F_0=\beta^2$ for equal-speed parallel motion with transverse separation. The length relation comes from [LVE](/en/theory/lve); the force ratio additionally uses the exterior true interaction distance and inverse-square constitutive input, described in [Coulomb interaction](/en/theory/electromagnetism/coulomb-interaction) and the [magnetic-force derivation](/en/theory/electromagnetism/magnetic-force). The script's `force_ratio` corresponds to the true physical force $F_T/F_0$, not directly to a reference-coordinate generalized force or a detector reading. Its speed domain is $0\leq\beta<1$; this does not extend the specified configuration to arbitrary sources and probes. After downloading the script, run: ```sh python ASD_LVE_demo.py --beta 0.6 --csv lve-model.csv python ASD_LVE_demo.py --beta 0 ``` The first command prints $\gamma=1.25$, $F/F_0=0.64$ and $F_m/F_0=0.36$; the second gives the rest limit $1,1,0$. `--beta` selects the single point printed to the terminal. The CSV always scans 0.00 through 0.95 in steps of 0.01, giving 96 rows with columns `beta,gamma,force_ratio,magnetic_fraction`. These numbers are generated from the stated formulas; no measurements, uncertainty estimates or general dynamics solver are included. ## Directional-clock pilot: geometry and identifiability The [clock pilot](/en/research/clock-response-pilot) is `ASD-RUN-CLOCK-001`. Its runnable materials are already included in [ASD_AI_Research_Pack.zip](/downloads/ASD_AI_Research_Pack.zip): `tasks/clock-pilot.py`, the apparatus record `tasks/clock-pilot.json` and frozen results `tasks/clock-pilot-expected.json`. The JSON records inputs and unknowns; the script fixes the following geometric examples in its code rather than reading apparatus settings from that file. The example uses latitude $40^\circ$, a user-selected background-axis declination of $30^\circ$, horizontal north/east axes for clocks A/B and 96 samples over one nominal sidereal day. The geometric template $Q=(\mathbf n\cdot\mathbf b)^2-1/3$ has not been derived from ASD atomic response. The script constructs the design matrix $X=(Q_A,-Q_B,1)$ for three coefficients in the candidate readout model. After extracting the research pack, run from its root: ```sh python tasks/clock-pilot.py --verify tasks/clock-pilot-expected.json --output clock-pilot-output.json ``` The JSON reports rank 3 for the default configuration, rank 2 at background declination $0^\circ$ and rank 2 for identical axes. It also checks axis reversal, periodicity and cancellation for identical axes with equal coefficients. The background and sampling are hypothetical inputs; rank establishes distinguishability only within this template. The catalog's clock result establishes geometric identifiability only: no ASD frequency amplitude has been calculated and no observations have been used. The selected $^{87}\mathrm{Sr}$ transition, reference clock and transfer-link response coefficients remain inputs to the [clock-response requirements](/en/tests/predictions#ASD-PRED-001); full rank cannot replace those physical calculations or establish detection power. ## Experimental comparison baselines and templates [ASD_Experimental_Reproducibility_Kit.zip](/downloads/ASD_Experimental_Reproducibility_Kit.zip) contains `README.md`, `equations.md`, `models.py`, `run.py`, `verify.py`, `notebook.ipynb`, `input/`, `schemas/`, `expected/` and `checksums.json`. Its comparison calculations require Python 3.10 or later and only the standard library. Inputs and units are documented in the [data catalog](/en/research/data) and the bundled instructions. After extraction, run: ```sh python run.py --output results python verify.py ``` `results/summary.json` summarizes the calculation, and `runtime.json` records the actual Python version. `clock_geometry.csv`, `muon_bins.csv` and `charge_response.csv` contain geometric factors, reference survival counts for synthetic cohorts, and the standard point-charge field with an assumed Gaussian response. This kit's clock geometry has 97 half-hour samples; it is a different calculation from the 96-sample two-clock rank pilot above, so their result files are not interchangeable. `expected/` contains frozen values. `verify.py` checks distributed-file SHA-256, numerical values and relevant limits, periodicity and symmetry. The muon calculation uses a constant-energy, flat-atmosphere reference survival formula and three bins explicitly labelled `SYNTHETIC`. The high-γ calculation uses a standard normalized transverse field of a uniformly moving point charge and convolution with a unit-area Gaussian, without an actual bunch, crystal or electro-optic transfer function. Four `schemas/*.csv` files contain headers only, as templates for future inputs; S301 has no orbital integration or fit. Assumptions, output meanings and missing work are detailed in [reproducibility resources](/en/research/reproducibility) and the [data catalog](/en/research/data). `asd_output`, `predicted_difference` and `fit` remain `null`. Diagnostic figures are bundled. To regenerate them, separately install the documented `matplotlib==3.10.8` and run `python figure.py --output results`; numerical calculations do not require it. The notebook calls the same functions and requires your own Jupyter environment. This kit contains no observational comparison or reproduction of the Ota experiment. ## What can be checked These commands check a formula demonstration, geometric ranks and comparison baselines respectively. On 2026-10-03, the actual distributed materials were rerun with Python 3.12.14: all 96 LVE CSV rows, the two-clock ranks and the three comparison models' JSON/CSV agreed with their supplied frozen values. Retain inputs, code, frozen outputs and checksum lists together. Numerical agreement makes this calculation checkable; it does not establish a missing ASD apparatus response. ## Requirements for full reproducibility This catalog uses the existing downloads and records; no confirmed public calculation-repository link is currently provided. Further calculations need a defined question, equations, boundary conditions, real-input provenance, environment, outputs and verification method, separating numerical error from measurement uncertainty. The clock needs response coefficients, muon decay needs a path-to-probability mapping, high-γ fields need joint source/probe response, and S301 needs verified object identity, epochs, frame and covariance. Their entries remain in the [prediction registry](/en/tests/predictions) and [open problems](/en/open-problems). Register real observational comparisons after those specific missing inputs are supplied.