Project Quantum Shield · Reference Guide

Biophysical Verification & Bibliography

The published work the project leans on, the mathematics it uses, and the line between the two.

Compiled with computer assistance from peer-reviewed physical and biological publications. Version 5 · August 2026.

Abstract

Spontaneous mutagenesis is classically modeled as a thermodynamic, over-the-barrier process. Recent biophysical work argues that quantum tunneling is a leading contributor to point mutations instead. This reference guide couples that subatomic layer to the population genetics of cultural divergence.

An open-quantum-systems treatment of double proton transfer in a G–C base pair, run under environmental dephasing, settles at a steady-state mutagenic occupancy near 1.73 × 10⁻⁴ within about 1.5 picoseconds. A migration–selection model then asks what happens to those rare events when cultural barriers reduce effective migration toward zero: isolated pools, faster drift, and genetic differentiation climbing toward unity.

What is established here is the tunneling literature and the population-genetics algebra. What is proposed — and remains unproven — is the coupling between them at human timescales. The coupling is the claim to attack.

Conceptual mapping

Each macro-hypothesis in this citizen-science model, set against the peer-reviewed mechanism it depends on.

Macro-hypothesisEmpirical mechanismValidation
Subsurface point mutations accumulateDouble proton tunneling (DPT)Löwdin (1963); Slocombe, Sacchi & Al-Khalili (2022)
Mutations are captured during replicationStrand-separation “locking” of tautomersSlocombe et al. (2022, Communications Chemistry)
Active biological resistanceEnzymatic helicase suppressionWinokan, Slocombe et al. (2023, Scientific Reports)
Environmental feedback loopsAdaptive / directed mutation couplingMcFadden & Al-Khalili (1999)
Cultural silos drive divergenceWright’s island model and Fₛₜ driftWright (1931); parapatric speciation models

Annotated bibliography

  1. Löwdin, P.-O. (1963). Proton tunneling in DNA and its biological implications. Reviews of Modern Physics 35(3), 724–732. doi:10.1103/RevModPhys.35.724

    The founding paper of quantum biology. Because protons are quantum-mechanical particles, they can tunnel across the double-well potential of a Watson–Crick hydrogen bond, producing rare tautomeric base pairs (T* and G*). If those forms persist into replication, they become point mutations.

    Status · Established as a theoretical proposal; the mechanism it opened is still being quantified.

  2. Slocombe, L., Sacchi, M., & Al-Khalili, J. (2022). An open quantum systems approach to proton tunnelling in DNA. Communications Physics 5, 109. doi:10.1038/s42005-022-00881-8

    Density functional theory plus open-system dynamics put tunneling rates in G–C base pairs orders of magnitude above classical thermal hopping. Despite strong dephasing from the hydration shell, the system reaches thermal equilibrium inside 1.5 ps at a steady-state tautomeric occupancy near 1.73 × 10⁻⁴.

    Status · Published calculation. The number is model-dependent: it moves with the assumed barrier, coupling, and bath.

  3. Slocombe, L., Winokan, M., Al-Khalili, J., & Sacchi, M. (2022). Proton transfer during DNA strand separation as a source of mutagenic guanine–cytosine tautomers. Communications Chemistry 5, 144. doi:10.1038/s42004-022-00760-x

    How a transient state becomes a permanent one. Mechanical unzipping by the replication machinery alters the local potential landscape and can freeze a transferred proton in its tautomeric position, leaving a base mismatch behind.

    Status · Mechanism modeled. Direct in-vivo capture of a frozen tautomer has not been demonstrated.

  4. Winokan, M., Slocombe, L., Al-Khalili, J., & Sacchi, M. (2023). Multiscale simulations reveal the role of PcrA helicase in protecting against spontaneous point mutations in DNA. Scientific Reports 13, 21749. doi:10.1038/s41598-023-48472-z

    Evidence of a defence system. PcrA helicase appears stereochemically and electrostatically arranged to suppress the double-proton-transfer rate immediately before strand separation, rather than acting as a passive zipper.

    Status · Simulation-level result. A strong falsifier: disrupt the relevant geometry and find no change in mutation rate.

  5. McFadden, J., & Al-Khalili, J. (1999). A quantum mechanical model of adaptive mutation. Biosystems 50(3), 203–211. doi:10.1016/S0303-2647(99)00004-0

    The most speculative thread in the set: if environmental coupling shapes the quantum state of the base pair, mutation is not uniformly random with respect to environment. That is a bias in the mechanism, not intent or direction.

    Status · Early theoretical proposal. Hold loosely; treat as open.

Standard evolutionary foundations

Behavioural and cultural isolation. Shifts in mating preference, language, and social structure are recognised pre-zygotic isolating mechanisms. Sustained cultural polarization can function as a barrier to gene flow, producing demes without any geographic separation.

Parapatric speciation and F-statistics. Where gene flow is restricted behaviourally rather than geographically, isolated pools drift. Wright's F-statistics give the standard way to measure it, with differentiation between silos climbing as effective migration falls.

A caveat this guide keeps visible: humans remain one species, and measured differentiation between human populations is low. The model describes a direction of pressure, not an accomplished split.

Simulation results and speciation tipping points

New in version 5. A 1,000-generation, 100-locus stochastic run, Ne = 30, baseline migration m₀ = 0.05, mutation held at the tunneling-derived rate μ = 1.73 × 10⁻⁴. Thresholds: moderate differentiation at FST ≥ 0.50, lock-in at FST ≥ 0.75.

ScenariomeFST ≥ 0.50FST ≥ 0.75Plateau
Extreme ideological silo (C = 0.99)0.0005Generation 61Generation 169≈ 0.80
Severe ideological silo (C = 0.95)0.0025Generation 111Not reached in 1,000≈ 0.57
High ideological silo (C = 0.90)0.0050Generation 537Not reached in 1,000≈ 0.30
Integrated population (C = 0.00)0.0500NeverNever≈ 0.03

The picosecond layer and the millennial layer meet here. Tunneling supplies a steady baseline mutation rate; behavioural isolation decides whether those events stay local long enough to accumulate. In this parameterisation the transition is not gradual — below C ≈ 0.90 differentiation stays modest across the whole window, while at C ≥ 0.95 the moderate threshold falls inside about a hundred generations.

Read the numbers as properties of the model, not measurements of any living population. Ne = 30 is a small deme, the loci are neutral, and a cultural filter coefficient of 0.99 describes closure no human society has sustained. The falsifier stays the same: show that measured differentiation between strongly separated modern populations does not track effective migration this way, and the coupling fails.

The figure, parameter mapping, and source code are on the cultural isolation FST simulation page.

Mathematical appendix

1 · Two-state Hamiltonian of the base pair

The canonical Watson–Crick pairing is written as |L⟩ and the mutagenic tautomer as |R⟩:

H^S=(ELΔΔER)=ϵ2σzΔσx\hat{H}_S = \begin{pmatrix} E_L & -\Delta \\ -\Delta & E_R \end{pmatrix} = \frac{\epsilon}{2}\,\sigma_z - \Delta\,\sigma_x

with energy asymmetry ε = EL − ER ≈ 0.2246 eV and tunneling matrix element Δ ≈ 0.1 meV.

2 · Lindblad master equation

The base pair sits in a warm, wet cell, so it is treated as an open system:

dρ(t)dt=i[H^S,ρ(t)]+k(LkρLk12{LkLk,ρ})\frac{d\rho(t)}{dt} = -\frac{i}{\hbar}\left[\hat{H}_S, \rho(t)\right] + \sum_k \left( L_k \rho L_k^{\dagger} - \tfrac{1}{2}\left\{ L_k^{\dagger} L_k, \rho \right\} \right)
  • Pure dephasing from solvent and neighbouring nucleotides, L₁ = √γdephase σz, with γdephase ≈ 50 ps⁻¹.
  • Thermal relaxation |R⟩ → |L⟩, L₂ = √γdown σ₋, with γdown ≈ 10 ps⁻¹.
  • Thermal excitation |L⟩ → |R⟩, L₃ = √γup σ₊, fixed by detailed balance:
γup=γdownexp ⁣(ϵkBT)\gamma_{up} = \gamma_{down}\,\exp\!\left(-\frac{\epsilon}{k_B T}\right)

At body temperature, T = 310.15 K, this gives a steady-state tautomeric population ρRR(∞) ≈ 1.73 × 10⁻⁴, reached in under 1.5 ps.

These are the parameters the interactive model uses. You can move them and watch the steady state move with them: DNA Proton Tunneling Simulator.

3 · Population genetics of cultural isolation

Wright's island model, with effective migration scaled by a cultural filter coefficient C ∈ [0, 1] against a baseline rate m₀:

me=m0(1C)m_e = m_0\,(1 - C)

Under drift, steady-state differentiation between sub-populations of effective size Ne is:

FST11+4Neme=11+4Nem0(1C)F_{ST} \approx \frac{1}{1 + 4 N_e m_e} = \frac{1}{1 + 4 N_e m_0 (1 - C)}
limC1FST=1\lim_{C \to 1} F_{ST} = 1

The limit is the point of the algebra, not a prediction about the present: total cultural closure would mean total genetic isolation, and independent mutation lines would then accumulate separately. Any real value of C well short of 1 leaves differentiation small.

Source documents