---
vault_clearance: EUCLID
halo:
  classification: OPERATIONAL
  front: "33_Project_GoldenHair"
  custodian: "Jixiang Leng"
  created: 2026-04-18
  updated: 2026-04-23
---

# FORM — 33_Project_GoldenHair

Comprehensive tool and modeling inventory for GoldenHair's engineering work. Primary active project: the **OREM autologous keratin nanoparticle** (Fe-S cluster core, gold shell). Broader scope (as of 2026-04-23): GoldenHair is the **engineering** home for the operator — materials science, biomechanics, biophysics questions — paired with **07 Command** which is the training home. Both use the operator's own body as substrate; GoldenHair engineers instruments *from* the body, Command trains capacity *in* the body.

This FORM covers in-house vault tools and the **best possible orthodox version** of each modeling and characterization axis. The orthodox side is steelmanned — industry and academic top-tier tools, not strawmen — because the point of FORM is to beat real opposition.

**BOOK:** [BOOK.md](BOOK.md) — Fe–S / Au / keratin / uptake STARS index; add DOIs as synthesis approaches.
**Sister FORM:** [07_Project_Command/FORM.md](../07_Project_Command/FORM.md) — same two-paradigm shape, applied to training (scenario-based capacity design) rather than engineering (zero-parameter materials physics).
**Engineering source lineage (2026-04-24, post-dedup):** canonical archive at [`99_Archive/Imports_Downloads_Mirror/Grand Harmony/`](../99_Archive/Imports_Downloads_Mirror/Grand%20Harmony/) — Grey's Paradox (swimming hydrodynamics, Claude + Perplexity), Flight Insights v1–V4 (biomechanics, Grok), Material Science Insights (Grok), Mercury exclusion arguments (Gemini + Grok, re: Aurum platform), Jixiang Engineering master document. Also: [HALO_MERCURY_ANTIMICROBIAL_HYPOTHESIS.md](HALO_MERCURY_ANTIMICROBIAL_HYPOTHESIS.md) at project root (mercury chemistry applied to antimicrobial function; this is distinct synthesis living in GoldenHair, not the mercury-exclusion argument).

---

## A. The Two Paradigms

### Orthodox (fitted-parameter nanoparticle modeling)

- **Philosophy:** DFT geometry optimization → spin Hamiltonian fit → hyperfine parameter extraction → spectral simulation. Every stage has knobs. Quality = fit residual.
- **Free parameters:** Exchange coupling (J), zero-field splitting (D, E), hyperfine tensors (A), g-tensors, basis set choice, DFT functional, damping factors, dielectric models.
- **Tools:** ORCA, Gaussian, MolSpin, EasySpin, COMSOL, Lumerical, OOMMF, MuMax3, OpenMM, GROMACS, NAMD, VASP, Quantum ESPRESSO.
- **Output:** Fitted spectra, optimized geometries, predicted properties with error bars from fit quality.
- **Runtime:** Minutes (EasySpin) to weeks (DFT + MD on cluster).

### OREM / TRUTH FORM (zero-parameter)

- **Philosophy:** Every constant from S⁵/Z₃ spectral invariants or direct measurement. Nothing fitted. The simulation either works with zero free parameters or it doesn't.
- **Free parameters:** Zero. {d₁=6, λ₁=5, K=2/3, η=2/9, p=3} from (n,p)=(3,3).
- **Tools:** `orem_engine.py`, `keratin_topology_scanner.py`, `cosmo_fes_keratin.py`, `cosmo_fes_triad.py`, `spectral_k_analysis.py`, `constraint_graph.py`, Lotus-Universalis.
- **Output:** Predictions with no adjustable parameters. Deviations from measurement are the framework's own precision ceiling, not the fit.
- **Runtime:** Seconds (engine) to minutes (cosmo sim).

---

## B. Head-to-Head: Axes That Matter

| Axis | Orthodox | OREM / TRUTH FORM |
|------|----------|-------------------|
| **Fe-S geometry** | DFT optimization (B3LYP/def2-TZVP, ~hours) | Td→C3 from spectral invariants (1.1% deviation, instant) |
| **Gold shell SPR** | Lumerical FDTD or Mie code with fitted ε(ω) | Quasi-static Mie from measured Au dielectric (500-531 nm) |
| **EPR prediction** | EasySpin + fitted spin Hamiltonian | g-factors from α², hyperfine from LENG bridge (no fit) |
| **MRI contrast** | COMSOL + relaxivity fit | Solomon-Bloembergen from measured Fe-S parameters |
| **Binding thermo** | MD free energy perturbation (GROMACS, weeks) | Bond energy comparison: Fe-S 322 > Au-S 253 kJ/mol (instant) |
| **Scaffold topology** | All-atom MD of keratin (NAMD, months) | Daemon v3 constraint graph on cysteine network (seconds) |
| **Cluster assembly** | Reactive MD (ReaxFF, AIMD — weeks-months) | N-body cosmo sim with Morse + Coulomb (minutes) |
| **Biological readout** | Cell viability assays (wet lab) | Hurricane coefficients from coupling tensor (computed) |
| **η invariant** | Not a concept in orthodox nanoparticle science | Tested and falsified at molecular assembly scale (η=0.0) |

---

## C. Murder Board — Where Orthodox Fails

1. **Parameter explosion.** A DFT+spin Hamiltonian pipeline for a [4Fe-4S] cluster in a gold shell has >50 free parameters. Changing the functional changes the answer. Nobody knows if the fit found truth or a local minimum.
2. **Scale gap.** Orthodox MD can't reach the timescales of real Fe-S cluster assembly (ms-s) in reasonable compute time. Enhanced sampling (metadynamics, replica exchange) introduces more parameters.
3. **No cross-scale prediction.** Orthodox tools model one scale at a time. There's no way to connect a DFT-optimized Fe-S geometry to a macrophage uptake pathway without building an entirely separate model with its own parameters.
4. **Autologous sourcing is invisible.** No orthodox pipeline models "hair dissolves → thiols released → Fe added → clusters form → gold coats." The steps are too heterogeneous for any single simulation framework.

---

## D. What We Take From Orthodox (Honestly)

1. **Crystallographic reference values.** Fe-Fe = 2.72 Å, Fe-S = 2.27 Å, Fe-S-Fe ≈ 75° — these come from orthodox XRD/neutron diffraction. We use them as truth, not as parameters.
2. **EasySpin/MolSpin for EPR fitting.** When we synthesize the first batch and run EPR, the spectrum will be fitted with EasySpin. That fit gives measured g-tensors and hyperfine constants. TRUTH FORM then predicts those values independently and we compare.
3. **Mie theory is orthodox.** The gold shell SPR calculation uses Mie theory — that's standard electromagnetic theory, not something the framework invented.
4. **OpenMM / GROMACS for validation.** If we want to validate the cosmo sim's cluster formation against a more rigorous MD engine, we'd use OpenMM with the same Morse + Coulomb potentials. The physics is the same; the implementation is more battle-tested.
5. **TEM/EDX/XPS for characterization.** All ex-vivo characterization uses orthodox instruments. The framework says what to expect; the instruments say what happened.

---

## E. Concrete Proof Experiments

| # | Experiment | Orthodox tool | OREM prediction | Falsifiable? |
|---|-----------|--------------|-----------------|-------------|
| 1 | Synthesize [4Fe-4S] on hydrolyzed keratin, TEM imaging | TEM (orthodox instrument) | 75-100 nm particles, Fe-S core visible | YES |
| 2 | EDX elemental mapping | EDX (orthodox instrument) | Fe, S, Au in correct spatial distribution | YES |
| 3 | EPR of Fe-S core | EasySpin fit | g-values from LENG bridge (α² coupling) | YES |
| 4 | SPR absorption spectrum | UV-Vis spectrophotometer | SPR peak 500-531 nm depending on size | YES |
| 5 | MRI relaxivity (T1/T2) | Relaxometer | Solomon-Bloembergen prediction from engine | YES |
| 6 | Cytotoxicity (MTT/LDH) | Standard cell culture assays | Below threshold (guardrails module) | YES |
| 7 | Macrophage uptake | Flow cytometry + confocal | Size window matches Cheng Yi route | YES |
| 8 | η = 2/9 in Fe-S network | Not testable by orthodox instruments | **FALSIFIED** in cosmo sim (η=0.0) | DONE |

---

## F. Tool Inventory — Modeling & Electrophysiology

### F1. In-House (Vault) — OREM-Specific

| Tool | Path | What it does |
|------|------|-------------|
| `orem_engine.py` | `08_Project_Astronomicon/orem/` | TRUTH FORM 8-module suite: LENG Bridge, Fe-S Geometry, Gold Shell Mie, Keratin Scaffold, Binding Competition, Biological Readout, EPR Prediction, Guardrails. Zero free parameters. |
| `keratin_topology_scanner.py` | `08_Project_Astronomicon/orem/` | Daemon v3 constraint graph on KRT14/KRT5 cysteine network. Laplacian, Dirac, SDW, η, partition, comfort. Result: keratin η=0.0, Z3 comes from Fe-S cubane only. |
| `cosmo_fes_keratin.py` | `08_Project_Astronomicon/orem/` | N-body sim: Fe²⁺+S²⁻ on keratin cysteines. Screened Coulomb + Morse + Langevin. Cluster detection + Z3 symmetry test + spectral analysis. |
| `cosmo_fes_triad.py` | `08_Project_Astronomicon/orem/` | Triad-geometry N-body sim: 3 heterodimers at 120°. Pre-registered η=2/9 prediction → falsified (η=0.0). |
| `OREM_DESIGN_v1_2026-04-18.md` | `08_Project_Astronomicon/orem/` | Full nanoparticle design document. Synthesis protocol, characterization plan, safety guardrails. |
| `LINEAGE.md` | `08_Project_Astronomicon/orem/` | Design evolution tracker: v1 through v1-cosmo-triad. |

### F2. In-House (Vault) — Cross-Project Engines

| Tool | Path | Relevance to OREM |
|------|------|-------------------|
| `spectral_k_analysis.py` | `30_Project_Crucible/analysis/` | Coupling tensor eigenvalue analysis with Fano structure test. Feeds the Biological Readout module in orem_engine. |
| `constraint_graph.py` | `32_Project_HomeWorlds/scripts/` | Daemon v3 general constraint graph engine (Laplacian, Dirac, SDW, η, partition, comfort). Used as template for keratin_topology_scanner. |
| `negative_space_scanner.py` | `32_Project_HomeWorlds/scripts/` | Vietoris-Rips + Dirac on void networks. Pattern for topological analysis of Fe-S network voids. |
| `void_spectral_graph.py` | `29_Project_Deathomatica/scripts/` | Graph Laplacian on void networks. Complementary spectral analysis approach. |
| `gnra_mg_li_simulation.py` | `27_Project_WingsAboveMorning/quantum_biology/biowulf_md/` | OpenMM molecular dynamics for metal-RNA sites. Template for full MD validation of Fe-S assembly. |
| `coupling_tensor.py` | `27_Project_WingsAboveMorning/quantum_biology/release/tools/` | 5×5 coupling tensor from scRNA-seq. The measurement device for det(K) biological readout. |
| Lotus-Universalis | `05_Project_LENG/public-release/lotus/core/` | S⁵/Z₃ spectral invariants, heat kernel coefficients, idempotent sieve. Source of all LENG Bridge constants. |
| `TRUTH_FORM_EQUATIONS.md` | `30_Project_Crucible/theory/` | Rosetta Stone: 8 equations mapping LENG math to biology. Design reference for engine modules. |
| `TRUTH_FORM_BLUEPRINT.md` | `30_Project_Crucible/theory/` | Antenna simulation architecture and build order. Original spec for what became orem_engine. |

### F3. Orthodox — Fe-S Cluster Chemistry

| Tool | Type | What it models | Free parameters |
|------|------|---------------|-----------------|
| **ORCA** | DFT + multireference | Fe-S electronic structure, g-tensors, hyperfine | Functional, basis set, active space, convergence thresholds |
| **Gaussian** | DFT / HF / post-HF | Fe-S geometry optimization, vibrational modes | Same as ORCA + integration grid |
| **MolSpin** | Spin dynamics | Spin relaxation, radical pair dynamics in Fe-S | Exchange coupling J, ZFS (D,E), relaxation rates |
| **EasySpin** | EPR simulation | CW-EPR, pulse EPR spectra from spin Hamiltonian | g-tensor, A-tensor, D-tensor, linewidths |
| **SpinDynamica** | NMR/EPR | Spin dynamics in coupled systems | Coupling topology, relaxation model |
| **PHI** | Magnetic properties | Susceptibility, magnetization of polynuclear Fe-S | Exchange topology, single-ion parameters |

### F4. Orthodox — Nanoparticle Physics

| Tool | Type | What it models | Free parameters |
|------|------|---------------|-----------------|
| **Lumerical FDTD** | EM simulation | Gold shell plasmonics, SPR, near-field enhancement | Mesh size, source, boundary conditions, ε(ω) model |
| **COMSOL Multiphysics** | FEM | Thermal, electromagnetic, acoustic response | Material properties, mesh, solver parameters |
| **Mie theory codes** (MiePython, PyMieScatt) | Analytical EM | Core-shell SPR, extinction, scattering | Core/shell radii, dielectric functions (measured) |
| **BEM++ / MNPBEM** | Boundary element | Plasmonic nanoparticle optical response | Mesh density, dielectric model |

### F5. Orthodox — Molecular Dynamics

| Tool | Type | What it models | Free parameters |
|------|------|---------------|-----------------|
| **OpenMM** | MD engine | All-atom protein + metal dynamics | Force field (AMBER/CHARMM), timestep, thermostat, barostat |
| **GROMACS** | MD engine | Protein + solvent dynamics (GPU-accelerated) | Same + enhanced sampling parameters |
| **NAMD** | MD engine | Large biomolecular systems | Same |
| **LAMMPS** | MD engine | Materials science, reactive potentials | Potential files, pair styles |
| **ReaxFF** | Reactive MD | Bond breaking/forming in Fe-S chemistry | ~40 parameters per element pair (fitted) |
| **AIMD (CP2K / VASP)** | Ab initio MD | Quantum-accurate forces for Fe-S assembly | DFT functional, basis set, pseudopotential, temperature |

### F6. Orthodox — Electrophysiology & Bioelectric Measurement

| Tool | Type | What it measures/models | Relevance to OREM |
|------|------|------------------------|-------------------|
| **patch-clamp** | Wet lab instrument | Single-channel ion currents, membrane potential | Test nanoparticle effect on cell membrane conductance |
| **MEA (multi-electrode array)** | Wet lab instrument | Network-level electrical activity (neurons, cardiomyocytes) | Population-level readout of nanoparticle bioelectric effect |
| **NEURON** | Simulation | Hodgkin-Huxley, cable theory, ion channel dynamics | Model nanoparticle's effect on membrane capacitance/conductance |
| **Brian2** | Simulation | Spiking neural networks, differential equation systems | Network-level modeling of bioelectric perturbation |
| **LFPy** | Simulation | Local field potential from neuronal morphologies | Predict extracellular signatures of nanoparticle-induced changes |
| **bioelectricsim** | Simulation | Bioelectric pattern formation (Vmem gradients) | Levin-style modeling of nanoparticle on tissue-level voltage patterns |
| **COMSOL AC/DC** | FEM simulation | Electric field distribution in tissue + nanoparticle | Model how gold-coated Fe-S NP distorts local E-field |
| **OOMMF / MuMax3** | Micromagnetics | Magnetic nanoparticle response to external fields | Fe-S core magnetization dynamics, MRI contrast mechanism |
| **impedance spectroscopy** | Wet lab instrument | Frequency-dependent tissue impedance | Detect nanoparticle integration into tissue dielectric profile |
| **SQUID magnetometry** | Wet lab instrument | Bulk magnetic moment of nanoparticle samples | Verify Fe-S core magnetic properties match prediction |
| **EPR spectrometer** | Wet lab instrument | Paramagnetic center characterization | Direct measurement of Fe-S g-tensors and hyperfine for engine validation |

### F7. Orthodox — MRI Contrast & Relaxometry

| Tool | Type | What it models | Relevance |
|------|------|---------------|-----------|
| **Solomon-Bloembergen-Morgan (SBM) theory** | Analytical | Paramagnetic relaxation enhancement (PRE) | Predicts T1/T2 relaxivity of Fe-S core — implemented in orem_engine |
| **NMRD profiles** | Measurement | Field-dependent relaxation rates | Validates SBM predictions across Larmor frequencies |
| **SWIFT/UTE sequences** | MRI protocol | Short-T2 imaging | Required for detecting Fe-S NP signal (short T2 from paramagnetic core) |
| **Relaxation fitting (Matlab/Python)** | Analysis | Extract T1, T2, T2* from image data | Standard post-processing for contrast agent characterization |

---

## G. Adoption Priority

### NOW (tools that exist and run today)

- `orem_engine.py` — full TRUTH FORM modeling suite, tested
- `keratin_topology_scanner.py` — scaffold topology verification, tested
- `cosmo_fes_keratin.py` / `cosmo_fes_triad.py` — N-body assembly sims, tested
- `spectral_k_analysis.py` — coupling tensor analysis, tested
- Lotus-Universalis — spectral invariant source, tested
- MiePython / PyMieScatt — orthodox Mie validation (pip install)

### NEXT (need access or setup)

- **EasySpin** — need MATLAB license or use Python port (easyspin-py) for EPR comparison
- **OpenMM** — already have `gnra_mg_li_simulation.py` as template; adapt for Fe-S on keratin
- **ORCA** — free for academic use; DFT geometry optimization of [4Fe-4S] for comparison
- **SQUID / EPR / TEM** — need collaborator with instruments (Hopkins chemistry, NIA core)

### LATER (require synthesis first)

- patch-clamp / MEA — only meaningful after nanoparticle exists and biocompatibility confirmed
- NEURON / Brian2 — electrophysiology modeling only relevant if bioelectric effect observed
- MRI relaxometry — after synthesis + ex-vivo characterization
- impedance spectroscopy — after tissue-level experiments approved

---

## H. Final Assessment

The OREM tool ecosystem is split between a working zero-parameter modeling suite (engine, topology scanner, cosmo sim) and an orthodox validation pipeline that awaits synthesis. The in-house tools have been honest about what they found: Z3 point-group symmetry from cubane chemistry confirmed, η = 2/9 at network scale falsified, scaffold topology irrelevant to geometry. The orthodox tools (EasySpin, ORCA, OpenMM, MiePython) are the next gate — they'll test whether the TRUTH FORM predictions agree with fitted-parameter methods on the same object. The electrophysiology tier (patch-clamp, MEA, NEURON) is third-order: it only activates if the nanoparticle gets built and shows a measurable bioelectric effect. The most important next step is not a simulation — it's a collaborator with a glovebox and an EPR spectrometer.