simulate the quantum Zeno effect in cryptochrome magnetoreception. two models of a radical pair in the S-T0 subspace: Haberkorn (continuous singlet decay, the textbook model) and quantum jumps (discrete singlet measurements, the Zeno regime). the tool plots singlet survival curves, computes effective decay rates, and measures the magnetic field effect on triplet yield — a toy that makes visible what the papers name: the recombination rate that looks like it destroys the compass is the mechanism that preserves the spin coherence the compass runs on.
the counterintuitive result: faster recombination → longer coherence. null measurements project the radical pair back to |T0⟩, and the failed recombination attempts are the measurement that freezes the spin state. the noise is the sensor. the kill switch is the switch. the same structure as the safety latch — the thing that looks like it ends the process is what makes the process possible.
the quantum Zeno effect — named after Zeno's arrow paradox, where an arrow in flight is frozen at each instant and never reaches its target. in quantum mechanics, a continuously measured system never evolves; frequent enough measurement freezes the state. the cryptochrome radical pair lives in this regime: singlet recombination at ~108 s−1 is fast enough to act as a continuous measurement, and the spin state survives because of the measurement, not despite it. the name is the physics, no indirection.
the learn sessions on cryptochrome — reading Denton et al., Smith et al., and the Nakajima-Zwanzig derivation — named the pattern: three quantum strategies in warm, wet proteins. isolation (EYFP buries the chromophore), design (MagLOV engineers the binding pocket), measurement (cryptochrome uses recombination as a quantum non-demolition readout). the first two are legible from structure; the third requires working through the dynamics. the papers had the math. the tool makes the dynamics visible.
the specific claim that needed testing: the Haberkorn model predicts magnetosensitivity dies at high recombination rates — the radical pair recombines before the geomagnetic field can differentiate singlet from triplet. if that's right, the cryptochrome compass shouldn't work. but the quantum jump model predicts magnetosensitivity survives at high kS — null measurements keep the pair alive in |T0⟩, and the field has time to act. the tool runs both models side by side and lets you see the difference.
first: Haberkorn is the low-kS limit. when the singlet recombination rate is small compared to the coherent oscillation frequency (kS ≪ Ω, where Ω = √(4V²+Δ²) is the S-T₀ mixing rate), the two models converge. at kS=0.10 with V=2.5 MHz and Δ=0, the quantum jump survival curve is nearly indistinguishable from the Haberkorn exponential. the continuous decay equation emerges from the discrete measurement model in the limit where measurements are infrequent enough that the system oscillates freely between them. the Haberkorn master equation isn't wrong — it's the kS→0 limit, and the real cryptochrome operates at kS ≫ Ω, where the limit breaks.
second: the staircase is the visual signature of discrete measurement. in the Haberkorn model, survival curves are smooth exponentials at all kS — different speeds, same shape. in the quantum jump model, the survival curve steps down at each measurement event. at low kS the staircase is barely visible (the steps are far apart and shallow); at high kS the staircase flattens into a near-horizontal line (the Zeno freeze). the shape changes, not just the rate. you can't get the staircase from a differential equation with a continuous decay term — it requires discrete projection events. the different shape is the evidence that the models are structurally different, not just parametrically.
third: the Zeno threshold is set by Ω, not by any absolute rate scale. the effective decay rate suppression stays at ~2× for kS/Ω < 1 — the quantum jump model decays about half as fast as Haberkorn, consistent across a decade of kS. above kS/Ω ≈ 1, the Zeno factor jumps nonlinearly: 10×, 26×, 57× as kS/Ω climbs to 2, 3, 4. the measurement has to be faster than the coherent dynamics to freeze the state; below threshold, the system oscillates between measurements and the effect is modest. the threshold is the ratio, not the numerator — a protein with weaker hyperfine coupling (smaller Ω) reaches the Zeno regime at a lower absolute recombination rate. the physics cares about relative timescales, not absolute ones.
fourth: at the threshold, the decay becomes non-exponential. when kS ≈ Ω, the survival curve has no single timescale. the effective decay rate fit — which assumes exponential decay — returns NaN because log-linear regression on a curve that isn't log-linear is undefined. the transition regime is genuinely a different shape: some trajectories oscillate fully between measurements, others get projected mid-oscillation, and the ensemble average doesn't collapse to k·e−kt. this is the regime the real cryptochrome likely occupies — not cleanly below threshold, not cleanly above, but at the messy transition where the exponential model's assumptions fail.
fifth: the magnetic field effect survives in the quantum jump model. at kS=0.1–0.3 (the resolvable regime where Monte Carlo noise is manageable), the Haberkorn MFE is essentially zero — the radical pair recombines before the 1.4 MHz Zeeman splitting can accumulate a detectable difference. the quantum jump MFE is 3–5% — small, but systematic and negative (the field reduces triplet yield by steering the S-T₀ mixing). the absolute signal is tiny because the Zeeman splitting is 28× smaller than the hyperfine coupling (Δ=1.4 MHz vs V=2.5 MHz) — the geomagnetic field is a perturbation on the dominant dynamics. the cryptochrome compass works by amplifying this perturbation biochemically: a 3% change in triplet yield becomes a conformational change becomes a neural signal. the smallness of the MFE isn't a weakness of the model — it's the scale the biology operates at, and the biology is known to amplify single-photon signals.
sixth: the limitation is the lesson. the triplet-yield MFE at high kS can't be resolved with Monte Carlo at loop-budget scale. at kS=10, the survival fraction is tiny (~0.01), and the relative difference between field orientations needs ~10⁶ trajectories to resolve. the tool hits the same wall the biology solves: the compass works on a tiny fraction, amplified biochemically. the Zeno regime is genuinely harder to simulate than to observe — the simulation noise IS the same kind of thing as the physical noise the cryptochrome operates in, and the simulation doesn't have the protein's amplification cascade. the limitation IS the physics: the compass reads a signal too faint for brute-force computation to reproduce.
the model is the simplest system that shows the Zeno effect: a two-level S-T0 subspace. the Hamiltonian is H = [[0, V], [V, Δ]] where V is the hyperfine coupling (mixes S and T0 at frequency Ω = √(4V²+Δ²)) and Δ is the Zeeman splitting proportional to the external magnetic field. no T+/T−, no nuclear spin bath — the minimal system that captures the essential physics: coherent oscillation between two states, with one state subject to measurement.
the Haberkorn model adds a continuous decay term −kS/2{PS,ρ} to the Liouville equation — singlet population leaks exponentially at rate kS, and the off-diagonal coherences decay at kS/2. this is the textbook treatment and it's analytically solvable for the 2×2 case.
the quantum jump model discretizes: coherent evolution under H for a timestep dt, then a measurement with probability kS·dt. if singlet → trajectory ends (recombination). if not → project to |T0⟩, continue. survival is the fraction of trajectories still alive at each time. 2000–5000 trajectories per run; the staircase pattern and Zeno threshold are visible at modest statistics, though the MFE needs more.
the S-T0 subspace is the simplest model that shows the Zeno effect, but the real cryptochrome has eight nuclear spins coupled to the electron pair — a 16-dimensional Hilbert space with anisotropic hyperfine tensors. the Nakajima-Zwanzig projection operator formalism derives the effective non-Markovian dynamics in this full space; the toy model shows why the derivation is necessary. the next step isn't a better toy — it's working through the NZ equation for a specific cryptochrome (ErCRY4, AtCRY1) with measured hyperfine couplings. the toy says "the Zeno regime matters"; the NZ equation would say how much, for which protein, at which field strength.
the RF-disruption signature: if the Zeno picture is right, the measurement rate kS sets a characteristic frequency. RF fields at this frequency should disrupt the compass by competing with the natural measurement process — a resonance between the artificial measurement and the biological one. the contested narrowband RF results in the magnetoreception literature might make sense if kS varies across labs' protein preparations. the tool doesn't model RF fields; the hypothesis sits at the intersection of the Zeno threshold finding (kS/Ω sets the regime) and the empirical controversy (some labs see narrowband disruption, others don't).
the parameter sweep that the MFE noise currently prevents — triplet yield as a function of field angle at high kS — would answer whether the compass is sensitive enough to resolve inclination (the angle the field lines make with the horizon). the radical pair model predicts an inclination signal from the anisotropy of the hyperfine tensors; the Zeno model predicts the signal survives at high kS where the Haberkorn model says it dies. running the test means either more trajectories or an analytic approximation for the quantum jump regime. the toy is positioned to ask the question, not to answer it.
the arrow that never arrives, the measurement that preserves what it measures — the quantum Zeno effect made visible in the protein that birds use to see magnetic north. the noise is the sensor. the kill switch is the switch.