Quantum Mechanics in Biology

Quantum phenomena in living systems: from enzyme catalysis to photosynthesis

Introduction: Quantum Effects in Life

For decades, quantum mechanics was thought to be relevant only at atomic and subatomic scales, with biological systems operating purely according to classical physics. However, mounting evidence reveals that quantum phenomena play crucial roles in fundamental biological processes.

Quantum biology explores how quantum mechanical effects—superposition, tunneling, entanglement, and coherence— influence biological functions that classical models cannot fully explain.

Key Quantum Phenomena in Biology

Quantum Tunneling: Particles passing through energy barriers in enzyme reactions
Quantum Coherence: Long-lived quantum states in photosynthetic complexes
Superposition: Vibrational states in olfactory receptors
Entanglement: Potential role in avian magnetoreception and radical pair mechanisms

Advanced Modules

Dive deeper into quantum biology with these interactive, multi-page modules.

The Green Network: Stomatal Dynamics & Quantum Transport

Explore how guard-cell signalling, photosynthetic energy transfer, and quantum tunnelling converge in the stomatal complex. Includes interactive phase portraits, cellular automata, 3-D manifold viewers, and an in-browser Python runtime.

Open module →

📹 Video Lectures

Jim Al-Khalili: The Fundamentals of Quantum Biology

Physicist Jim Al-Khalili introduces the emerging field of quantum biology, explaining how quantum mechanics operates in living systems including photosynthesis, enzyme catalysis, and bird navigation.

Philip Ball: An Introduction to Quantum Biology

Science writer Philip Ball provides an accessible introduction to quantum biology, covering photosynthesis, enzyme reactions, bird navigation, and other quantum phenomena in living systems.

Quantum Biology: The Hidden Nature of Nature

Can the spooky world of quantum physics explain bird navigation, photosynthesis and even our delicate sense of smell? Clues are mounting that the rules governing the subatomic realm may play an unexpectedly pivotal role in the visible world. Join leading thinkers in the emerging field of quantum biology as they explore the hidden hand of quantum physics in everyday life and discuss how these insights may one day revolutionize thinking on everything from the energy crisis to quantum computers.

Topics Covered:

Bird navigation and magnetoreception
Quantum coherence in photosynthesis
Quantum vibration theory of olfaction
Applications to energy crisis and quantum computing
Future directions in quantum biology research

More Lectures & Talks

How Quantum Biology Might Explain Life's Biggest Questions

Jim Al-Khalili • TED Talk (16 min)

Quantum Life: How Physics Can Revolutionise Biology

Jim Al-Khalili • The Royal Institution (1h)

A Future with Quantum Biology

Alexandra Olaya-Castro • The Royal Institution (1h)

How Nature Harnesses Quantum Processes to Function Optimally

Clarice Aiello (UCLA) • Fields Institute (59 min)

Quantum Biology and Consciousness

Stuart Hameroff • TSC Conference (1h 22min)

Decoding the Universe: Quantum

NOVA PBS Official (54 min)

Quantum Coherence in Photosynthesis

Quantum Biology & Nanobiophysics: Photosynthesis

ICTP Condensed Matter & Statistical Physics (30 min)

Photosynthesis — Quantum Life

Keio University (1h 11min)

Quantum Coherence & Entanglement in Photosynthetic Light-Harvesting

P. Nalbach • ICAM (33 min)

The Quantum Design of Photosynthesis

Dr. Rienk Van Grondelle • PARC (1h)

How Quantum Coherence Assists Light Harvesting

American Chemical Society (3 min)

Are There Non-Trivial Quantum Effects in Biology?

Susana Huelga • IFISC (1h 17min)

Quantum Tunneling in Enzymes

Professor Judith Klinman on Quantum Tunneling in Enzymes

The Guy Foundation (1h 12min)

Quantum Biology: Enzymes, the Engines of Life

Corporis (12 min)

Essentials of Quantum Physics for Quantum Biology

Khwarizmi Science Society • Physics of Life (1h 6min)

Quantum Tunneling Explained

Perimeter Institute for Theoretical Physics (5 min)

Avian Magnetoreception & Radical Pair Mechanism

Radical Pair Mechanism of Magnetoreception

Peter Hore • FENS (1h 30min)

Biology of Magnetoreception in Night-Migratory Songbirds

Henrik Mouritsen • FENS (1h 20min)

Avian Magnetoreception — A Radical Sense of Direction

Prof. Peter Hore • Oxbridge BioSoc (1h 17min)

Quantum Mechanics of Magnetoreception

C. Rodgers • ICAM (41 min)

Magnetoreception in Birds (Part 2)

Khwarizmi Science Society • Physics of Life (56 min)

How Quantum Mechanics Help Birds Find Their Way

Nature Video (6 min)

1. Quantum Tunneling in Enzyme Catalysis

Enzymes are biological catalysts that accelerate chemical reactions by factors of $10^{10}$ to $10^{23}$. Classical transition state theory cannot fully explain these extraordinary rate enhancements. Quantum tunneling plays a critical role.

The Hydrogen Transfer Problem

Many enzyme reactions involve hydrogen (proton or hydride) transfer. Classically, the transferred particle must overcome an activation energy barrier. Quantum mechanically, light particles like hydrogen can tunnel throughthe barrier rather than going over it.

Derivation: Tunneling Through a Rectangular Barrier

We begin with the one-dimensional time-independent Schrodinger equation for a particle of mass $m$ encountering a rectangular potential barrier of height $V_0$ and width $a$:

$$-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} + V(x)\psi = E\psi$$

where the potential is defined as $V(x) = V_0$ for $0 \le x \le a$ and $V(x) = 0$ elsewhere.

Region I ($x < 0$): The particle is free, so the solution is a superposition of incident and reflected plane waves:

$$\psi_I = A e^{ikx} + B e^{-ikx}, \qquad k = \frac{\sqrt{2mE}}{\hbar}$$

Region II ($0 \le x \le a$): Inside the barrier where $E < V_0$, the kinetic energy is negative and the wave function decays exponentially:

$$\psi_{II} = C e^{-\kappa x} + D e^{\kappa x}, \qquad \kappa = \frac{\sqrt{2m(V_0 - E)}}{\hbar}$$

Region III ($x > a$): Beyond the barrier, only a transmitted wave propagates:

$$\psi_{III} = F e^{ikx}$$

We apply boundary conditions requiring continuity of $\psi$ and $d\psi/dx$ at $x = 0$ and $x = a$. This yields four equations for the five coefficients. Setting $A = 1$ (unit incident amplitude), we solve for the transmission amplitude $t = F/A$.

At $x = 0$:

$$A + B = C + D, \qquad ik(A - B) = -\kappa(C - D)$$

At $x = a$:

$$C e^{-\kappa a} + D e^{\kappa a} = F e^{ika}, \qquad -\kappa(C e^{-\kappa a} - D e^{\kappa a}) = ik F e^{ika}$$

Solving the system of equations via transfer matrices, the exact transmission coefficient is:

$$T = |t|^2 = \left[1 + \frac{V_0^2 \sinh^2(\kappa a)}{4E(V_0 - E)}\right]^{-1}$$

In the thick barrier (WKB) limit where $\kappa a \gg 1$, we have $\sinh(\kappa a) \approx \frac{1}{2}e^{\kappa a}$, so the prefactor becomes negligible compared to the exponential:

$$T \approx \frac{16E(V_0 - E)}{V_0^2} e^{-2\kappa a} \approx e^{-2\kappa a}$$

Substituting back $\kappa = \sqrt{2m(V_0 - E)}/\hbar$, we arrive at the standard WKB tunneling formula used in enzyme catalysis:

The tunneling probability depends on the barrier width $a$ and height $V_0$:

$$T \approx \exp\left(-\frac{2a}{\hbar}\sqrt{2m(V_0 - E)}\right)$$

Quantum tunneling transmission coefficient

Experimental Evidence

Kinetic isotope effects (KIE) > 100
Temperature-independent reaction rates
Non-Arrhenius behavior at low T
H/D/T isotope substitution studies

Example Enzymes

Alcohol dehydrogenase
Methylamine dehydrogenase
Aromatic amine dehydrogenase
Soybean lipoxygenase

2. Quantum Coherence in Photosynthesis

Photosynthesis achieves near-perfect quantum efficiency (~95%) in light harvesting. Excitation energy is transferred from light-harvesting antenna complexes to reaction centers with minimal loss. How?

The 2007 Breakthrough

Fleming and colleagues (2007) discovered long-lived quantum coherence in the Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria at physiological temperatures. Energy transfer occurs via quantum superposition states that "sample" multiple pathways simultaneously.

Derivation: Excitonic Hamiltonian and Energy Transfer

We start from the electronic structure of individual chromophores (pigment molecules). Each chromophore $n$ has a ground state $|g_n\rangle$ and an excited state $|e_n\rangle$ with excitation energy $\epsilon_n$. In the Frenkel exciton model, we assume excitons are localized on individual molecules (valid when intermolecular overlap is small).

We define the single-exciton basis state $|n\rangle$ as the state where chromophore $n$ is excited and all others are in the ground state:

$$|n\rangle = |g_1\rangle \otimes \cdots \otimes |e_n\rangle \otimes \cdots \otimes |g_N\rangle$$

The diagonal (site) energy of state $|n\rangle$ comes from the isolated chromophore Hamiltonian:

$$H_{\text{site}} = \sum_{n=1}^N \epsilon_n |n\rangle\langle n|$$

The off-diagonal couplings arise from electrostatic (Coulomb) interactions between transition dipole moments of different chromophores. In the point-dipole approximation, the coupling between chromophores $n$ and $m$ is:

$$J_{nm} = \frac{1}{4\pi\epsilon_0} \left[\frac{\boldsymbol{\mu}_n \cdot \boldsymbol{\mu}_m}{|\mathbf{r}_{nm}|^3} - \frac{3(\boldsymbol{\mu}_n \cdot \hat{\mathbf{r}}_{nm})(\boldsymbol{\mu}_m \cdot \hat{\mathbf{r}}_{nm})}{|\mathbf{r}_{nm}|^3}\right]$$

where $\boldsymbol{\mu}_n$ is the transition dipole moment of chromophore $n$ and $\mathbf{r}_{nm}$ is the vector connecting the two chromophores. These couplings enable excitation transfer between sites.

Combining diagonal and off-diagonal terms gives the full excitonic Hamiltonian in second-quantized (site) representation:

$$H = \sum_{n=1}^N \epsilon_n |n\rangle\langle n| + \sum_{n \neq m} J_{nm} |n\rangle\langle m|$$

In a real biological environment, the chromophores are embedded in a fluctuating protein scaffold. To account for dissipation and decoherence, we use the Lindblad master equation for the reduced density matrix $\rho$ of the excitonic system:

$$\frac{d\rho}{dt} = -\frac{i}{\hbar}[H, \rho] + \mathcal{L}[\rho]$$

where $\mathcal{L}[\rho] = \sum_k \left(L_k \rho L_k^\dagger - \frac{1}{2}\{L_k^\dagger L_k, \rho\}\right)$ describes the dissipative coupling to the environment via Lindblad operators $L_k$.

In the weak-coupling limit where $J_{nm} \ll k_BT$, the excitation transfer rate between sites reduces to the celebrated Forster resonance energy transfer (FRET) rate, which scales as $k_{\text{FRET}} \propto J_{nm}^2 \propto 1/R^6$. The remarkable finding in photosynthesis is that the system operates in an intermediate regime where quantum coherence persists, enabling environment-assisted quantum transport (ENAQT) that outperforms both purely coherent and purely incoherent transfer.

The excitonic Hamiltonian for N pigments:

$$H = \sum_{n=1}^N \epsilon_n |n\rangle\langle n| + \sum_{n \neq m} J_{nm} |n\rangle\langle m|$$

$\epsilon_n$ = site energies, $J_{nm}$ = excitonic couplings

Key Findings

Coherence times: 660 fs at 77 K, 400 fs at 277 K (remarkably long for warm, wet biology!)
Quantum beats: Oscillations in 2D electronic spectra revealing coherent superposition
Environment-assisted quantum transport (ENAQT): Noise helps maintain optimal energy transfer
Quantum walk: Exciton explores network topology via quantum superposition

3. Quantum Vibration Theory of Olfaction

How do we distinguish molecules with identical shapes but different smells? The controversial but compelling vibration theory of olfaction proposes that olfactory receptors detect molecular vibrations via inelastic electron tunneling.

Turin's Mechanism (1996)

Luca Turin proposed that odorant molecules act as "bridges" enabling electrons to tunnel across olfactory receptor proteins. The tunneling rate depends on vibrational modes of the odorant matching the energy gap.

Derivation: Inelastic Electron Tunneling Current

We begin with the Bardeen transfer Hamiltonian formalism. The system consists of a left electrode (donor, L), a right electrode (acceptor, R), and a molecular bridge. The total Hamiltonian is:

$$H = H_L + H_R + H_{\text{mol}} + H_T$$

where $H_T = \sum_{k,q} T_{kq} c_k^\dagger c_q + \text{h.c.}$ is the tunneling Hamiltonian with transfer matrix elements $T_{kq}$ coupling states $k$ in L to states $q$ in R.

Using Fermi golden rule (first-order perturbation theory in $H_T$), the elastic tunneling current between two metallic electrodes at bias voltage $V$ is:

$$I_{\text{el}} = \frac{2\pi e}{\hbar} \int_{-\infty}^{\infty} |T|^2 \rho_L(E) \rho_R(E - eV) \left[f_L(E) - f_R(E - eV)\right] dE$$

where $\rho_{L,R}(E)$ are the electronic densities of states and $f_{L,R}(E) = [1 + e^{(E-\mu_{L,R})/k_BT}]^{-1}$ are the Fermi-Dirac distributions with chemical potentials split by the bias: $\mu_L - \mu_R = eV$.

Now we introduce electron-phonon (vibrational) coupling. When a molecule with vibrational mode of frequency $\omega$ sits between the electrodes, the electron can lose energy $\hbar\omega$ to excite a vibration during tunneling. The coupling Hamiltonian is:

$$H_{e\text{-ph}} = \sum_{\alpha} \lambda_\alpha \hbar\omega_\alpha \, c^\dagger c \,(b_\alpha^\dagger + b_\alpha)$$

where $\lambda_\alpha$ is the dimensionless coupling constant and $b_\alpha^\dagger, b_\alpha$ are phonon creation/annihilation operators for mode $\alpha$.

The inelastic channel opens when the bias energy $eV$ matches a vibrational quantum $\hbar\omega_\alpha$. The total current becomes $I = I_{\text{el}} + I_{\text{inel}}$, where the inelastic contribution adds a step at each vibrational threshold. The key signature is in the second derivative:

$$\frac{d^2I}{dV^2} \propto \sum_\alpha \lambda_\alpha^2 \,\delta(eV - \hbar\omega_\alpha)$$

Each vibrational mode produces a peak in $d^2I/dV^2$ at the voltage corresponding to its energy. This is exactly the IETS spectrum. In Turin's olfaction model, the olfactory receptor acts as the two-electrode junction and the odorant molecule is the bridge whose vibrational spectrum determines the tunneling current and thus the perceived smell. The full current integral takes the form:

Inelastic electron tunneling spectroscopy (IETS) current:

$$I \propto \int_{-\infty}^{\infty} \rho_L(E) \rho_R(E-eV) \left[f_L(E) - f_R(E-eV)\right] dE$$

Enhanced tunneling when $eV = \hbar\omega$ (vibrational mode energy)

Supporting Evidence

H/D isotope discrimination in flies
Similar vibrational spectra → similar smell
Drosophila behavioral studies (2011)
Human psychophysical tests

Challenges

Conflicting experimental results
Shape theory still explains many odors
Molecular dynamics question tunneling
Debate remains active (2024)

4. Quantum Entanglement in Avian Navigation

Migratory birds detect Earth's magnetic field with extraordinary sensitivity. The leading theory involves radical pair mechanism with potentially entangled electron spins.

The Cryptochrome Hypothesis

Cryptochrome proteins in bird retinas undergo photochemical reactions creating radical pairs—molecules with unpaired electrons. External magnetic fields influence the singlet/triplet interconversion rates.

Derivation: Radical Pair Spin Hamiltonian

Consider a radical pair: two molecules each with one unpaired electron, with spin operators $\mathbf{S}_1$ and $\mathbf{S}_2$. The combined two-electron spin space is 4-dimensional, spanned by the singlet and triplet states:

$$|S\rangle = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle - |\downarrow\uparrow\rangle), \qquad S = 0$$

$$|T_+\rangle = |\uparrow\uparrow\rangle, \quad |T_0\rangle = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle + |\downarrow\uparrow\rangle), \quad |T_-\rangle = |\downarrow\downarrow\rangle, \qquad S = 1$$

Exchange interaction: When the two radicals are close, their electronic wave functions overlap, leading to an energy splitting between singlet and triplet states proportional to the exchange coupling $J$:

$$H_{\text{exchange}} = J\, \mathbf{S}_1 \cdot \mathbf{S}_2$$

This Heisenberg interaction gives eigenvalues $-3J/4$ for the singlet and $+J/4$ for each triplet state (using $\mathbf{S}_1 \cdot \mathbf{S}_2 = \frac{1}{2}[S(S+1) - s_1(s_1+1) - s_2(s_2+1)]$).

Zeeman interaction: An external magnetic field $\mathbf{B}$ couples to each electron spin via the magnetic moment:

$$H_{\text{Zeeman}} = \gamma_e \mathbf{B} \cdot \mathbf{S}_1 + \gamma_e \mathbf{B} \cdot \mathbf{S}_2$$

where $\gamma_e = g_e \mu_B / \hbar$ is the electron gyromagnetic ratio. The Zeeman interaction mixes $|S\rangle$ and $|T_0\rangle$ when the two electrons have different local fields (different g-factors or hyperfine environments), which is the essential mechanism for magnetic field sensitivity.

Hyperfine coupling: Each unpaired electron interacts with nearby nuclear spins $\mathbf{I}_i$ (e.g., nitrogen and hydrogen nuclei in the cryptochrome flavin radical). The hyperfine tensor $\mathbf{A}_i$ couples electron and nuclear spins:

$$H_{\text{HF}} = \sum_{i=1}^{2} \mathbf{S}_i \cdot \mathbf{A}_i \cdot \mathbf{I}_i$$

This is the critical term: hyperfine interactions create different effective magnetic fields at each radical, driving singlet-triplet interconversion whose rate depends on the external field direction.

Combining all three contributions, the full spin Hamiltonian is:

$$H = J\, \mathbf{S}_1 \cdot \mathbf{S}_2 + \sum_{i=1}^{2} \gamma_e \mathbf{B} \cdot \mathbf{S}_i + \sum_{i=1}^{2} \mathbf{S}_i \cdot \mathbf{A}_i \cdot \mathbf{I}_i$$

The experimentally observable quantity is the singlet yield, which determines the chemical outcome. If the singlet state reacts with rate constant $k_S$ and both states decay with rate $k$, the singlet yield is:

$$\Phi_S = \int_0^\infty P_S(t)\, k_S\, e^{-kt}\, dt$$

where $P_S(t) = \text{Tr}[\hat{P}_S \rho(t)]$ is the time-dependent singlet probability and $\hat{P}_S = |S\rangle\langle S|$ is the singlet projection operator. Since $P_S(t)$ depends on the direction of $\mathbf{B}$ through the Hamiltonian, so does $\Phi_S$, giving the bird directional magnetic field information. The resulting Hamiltonian takes the compact form:

Singlet-triplet mixing Hamiltonian:

$$H = J \mathbf{S}_1 \cdot \mathbf{S}_2 + \sum_{i=1}^{2} \gamma_e \mathbf{B} \cdot \mathbf{S}_i + \sum_{i=1}^{2} \mathbf{S}_i \cdot \mathbf{A}_i \cdot \mathbf{I}_i$$

J = exchange, B = magnetic field, A = hyperfine coupling

Quantum Features

Entangled spin states: $|S\rangle = \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle)$ vs $|T_0\rangle = \frac{1}{\sqrt{2}}(|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle)$
Coherence times: ~100 μs (sufficient for navigation)
Sensitivity: Detect ~50 μT (Earth's field) with directional information
Experimental support: Behavioral disruption by RF fields, cryptochrome localization

5. Quantum Tunneling in DNA Mutations

Spontaneous mutations are the raw material of evolution. While many arise from chemical damage or replication errors, a fundamental class of mutations originates from quantum mechanical proton tunneling along the hydrogen bonds that hold DNA base pairs together. This mechanism, first proposed by Per-Olov Löwdin in 1963, connects quantum physics directly to the molecular basis of heredity and evolution.

5.1 Watson-Crick Base Pairing and Hydrogen Bonds

DNA's double helix is stabilised by hydrogen bonds between complementary base pairs: adenine (A) pairs with thymine (T) via two hydrogen bonds, and guanine (G) pairs with cytosine (C) via three. Each hydrogen bond involves a proton shared between a donor atom (typically N-H) and an acceptor atom (N or O) on the complementary base.

The potential energy landscape along each hydrogen bond can be modelled as a double-well potential. The proton normally resides in the deeper well (the canonical tautomeric form), but quantum tunneling allows it to transfer to the shallower well (the rare tautomeric form).

The Double-Well Potential Model

The proton in a hydrogen bond N-H···N sits in a one-dimensional double-well potential $V(x)$. The left well (at position $x_L$) corresponds to the canonical tautomer and the right well (at $x_R$) to the rare tautomer. A convenient model is the quartic double well:

$$V(x) = V_0\left[\left(\frac{x}{a}\right)^2 - 1\right]^2 - \frac{\Delta\epsilon}{2}\frac{x}{a}$$

where $V_0$ is the barrier height, $2a$ is the distance between the two minima, and $\Delta\epsilon$ is the asymmetry energy (the energy difference between the canonical and rare forms). For the symmetric case ($\Delta\epsilon = 0$), the barrier height at $x = 0$ is exactly $V_0$.

The time-independent Schrödinger equation for the proton of mass $m_p$ in this potential is:

$$-\frac{\hbar^2}{2m_p}\frac{d^2\psi}{dx^2} + V(x)\psi(x) = E\psi(x)$$

For a proton ($m_p = 1.673 \times 10^{-27}$ kg) in a hydrogen bond with typical parameters ($V_0 \approx 0.3$-$0.7$ eV, $a \approx 0.3$-$0.5$ Å), the de Broglie wavelength $\lambda = h/\sqrt{2m_p E}$ is comparable to the barrier width, making tunneling non-negligible.

5.2 Löwdin's Mechanism (1963)

Per-Olov Löwdin proposed that protons in the hydrogen bonds of Watson-Crick base pairs can tunnel simultaneously from their normal positions to create rare tautomeric forms. In the G·C pair, which has three hydrogen bonds, the concerted transfer of two protons produces G*·C* (where * denotes the rare tautomer). Similarly, the A·T pair undergoes a double proton transfer to yield A*·T*.

Tautomeric Shifts in Base Pairs

G·C Base Pair

Normal G: keto form (C=O at position 6)
Rare G*: enol form (C-OH at position 6)
Normal C: amino form (-NH$_2$)
Rare C*: imino form (=NH)
G*·C* mispairs as A·T during replication

A·T Base Pair

Normal A: amino form (-NH$_2$)
Rare A*: imino form (=NH)
Normal T: keto form (C=O)
Rare T*: enol form (C-OH)
A*·T* mispairs as G·C during replication

The key insight is that if the rare tautomer exists at the moment DNA polymerase arrives, it will be read as a different base, causing a transition mutation:

$$\text{G}\cdot\text{C} \xrightarrow{\text{tunneling}} \text{G}^*\cdot\text{C}^* \xrightarrow{\text{replication}} \text{A}\cdot\text{T} \quad (\text{transition mutation})$$

$$\text{A}\cdot\text{T} \xrightarrow{\text{tunneling}} \text{A}^*\cdot\text{T}^* \xrightarrow{\text{replication}} \text{G}\cdot\text{C} \quad (\text{transition mutation})$$

5.3 Quantum Mechanical Derivation

Step 1: WKB Tunneling Probability

We start from Gamow's tunneling theory. Consider a proton of mass $m_p$ and energy $E$ encountering a potential barrier $V(x)$. The WKB approximation gives the tunneling probability through a barrier of arbitrary shape:

$$\mathcal{T} = \exp\left(-\frac{2}{\hbar}\int_{x_1}^{x_2} \sqrt{2m_p\left[V(x) - E\right]}\, dx\right)$$

where $x_1$ and $x_2$ are the classical turning points defined by $V(x_1) = V(x_2) = E$. This is the Gamow factor, valid when the barrier is slowly varying compared to the de Broglie wavelength.

Step 2: Attempt Frequency

The proton oscillates in its potential well (the hydrogen bond) with a characteristic attempt frequency $\nu_0$. This is the vibrational frequency of the N-H or O-H stretching mode. For a harmonic well of force constant $k_f$:

$$\nu_0 = \frac{1}{2\pi}\sqrt{\frac{k_f}{m_p}} \sim 10^{13}\ \text{Hz}$$

Each oscillation gives the proton one chance to tunnel, so the tunneling rate is:

$$k_{\text{tunnel}} = \nu_0 \times \mathcal{T} = \nu_0 \exp\left(-\frac{2}{\hbar}\int_{x_1}^{x_2} \sqrt{2m_p[V(x) - E]}\, dx\right)$$

Step 3: Rectangular Barrier Approximation

For a rectangular barrier of effective height $\Delta E = V_0 - E$ and width $a$, the integral simplifies since $V(x) - E = \Delta E$ is constant:

$$\int_{x_1}^{x_2} \sqrt{2m_p \Delta E}\, dx = a\sqrt{2m_p \Delta E}$$

Substituting back, the tunneling rate for a rectangular barrier becomes:

$$\boxed{k_{\text{tunnel}} = \nu_0 \exp\left(-\frac{2a}{\hbar}\sqrt{2m_p \Delta E}\right)}$$

$\nu_0$ = attempt frequency, $a$ = barrier width, $m_p$ = proton mass, $\Delta E$ = barrier height above proton energy

Step 4: Temperature-Dependent Rate (Marcus-like Theory)

At finite temperature, the proton can access higher vibrational states, which see a thinner effective barrier. Combining WKB tunneling with transition state theory, the full rate is:

$$k(T) = \frac{\nu_0}{2\pi} \exp\left(-\frac{2a}{\hbar}\sqrt{2m_p \Delta E}\right) \exp\left(-\frac{E_a}{k_BT}\right)$$

where the first exponential is the quantum tunneling factor and $\exp(-E_a/k_BT)$ is the Boltzmann factor for thermal activation to the tunneling energy $E_a$. Three temperature regimes emerge:

High T: Classical over-barrier (Arrhenius): $k \propto \exp(-E_a/k_BT)$
Intermediate T: Thermally-assisted tunneling: both factors contribute
Low T: Pure quantum tunneling: $k$ becomes temperature-independent (Boltzmann saturates)

The crossover temperature $T_c$ below which tunneling dominates is:

$$T_c = \frac{\hbar\omega_b}{2\pi k_B}$$

where $\omega_b$ is the imaginary frequency at the barrier top (the curvature of the inverted barrier). For proton transfer in DNA, $T_c \approx 200$-$350$ K, meaning tunneling is significant at biological temperatures.

5.4 Concerted Double Proton Transfer

Löwdin's original proposal involved simultaneous tunneling of two protonsin opposite directions along the hydrogen bonds of a base pair. For concerted transfer, both protons must tunnel at the same time, which modifies the effective tunneling rate.

Two-Proton Tunneling Hamiltonian

For two protons at positions $x_1$ and $x_2$ along their respective hydrogen bonds, the Hamiltonian is:

$$H = -\frac{\hbar^2}{2m_p}\frac{\partial^2}{\partial x_1^2} - \frac{\hbar^2}{2m_p}\frac{\partial^2}{\partial x_2^2} + V_1(x_1) + V_2(x_2) + V_{12}(x_1, x_2)$$

where $V_1$ and $V_2$ are the individual double-well potentials and$V_{12}$ is the coupling between the two protons. The coupling term is crucial: it encodes the electrostatic interaction and the structural reorganisation of the base pair during transfer.

In the instanton (semiclassical path integral) formalism, the concerted tunneling rate for$n$ protons transferring through barriers of similar height is approximately:

$$k_{\text{concerted}} \approx \nu_0 \exp\left(-n \cdot \frac{2a}{\hbar}\sqrt{2m_p \Delta E}\right) = \nu_0\, \mathcal{T}^n$$

For the G·C pair ($n = 2$ protons transferring), the concerted rate is the square of the single-proton tunneling probability: $\mathcal{T}^2$. This makes concerted transfer significantly slower than single-proton transfer, but recent ab initio calculations show the coupling $V_{12}$can lower the effective barrier substantially, partially compensating for this suppression.

5.5 Numerical Estimates for DNA

Computing the Tunneling Rate

Using parameters from ab initio quantum chemistry (DFT/MP2) calculations for the G·C base pair:

Parameter	Symbol	Value
Barrier height	$\Delta E$	0.4 - 0.7 eV
Barrier width	$a$	0.3 - 0.5 Å
Proton mass	$m_p$	$1.673 \times 10^{-27}$ kg
Attempt frequency	$\nu_0$	$\sim 10^{13}$ Hz
Asymmetry energy (G·C)	$\Delta\epsilon$	0.02 - 0.05 eV
Tautomer lifetime	$\tau^*$	$10^{-13}$ - $10^{-10}$ s

The tunneling exponent evaluates to:

$$\frac{2a}{\hbar}\sqrt{2m_p \Delta E} = \frac{2 \times 0.5 \times 10^{-10}}{1.055 \times 10^{-34}}\sqrt{2 \times 1.673 \times 10^{-27} \times 0.4 \times 1.602 \times 10^{-19}}$$

$$= 9.49 \times 10^{23} \times 4.65 \times 10^{-24} \approx 4.4$$

Therefore:

$$k_{\text{tunnel}} = 10^{13} \times e^{-4.4} \approx 10^{13} \times 0.012 \approx 1.2 \times 10^{11}\ \text{s}^{-1}$$

This is the forward tunneling rate. However, the rare tautomer is thermodynamically unstable, so it tunnels back at a rate $k_{\text{back}} = k_{\text{tunnel}} \times \exp(\Delta\epsilon / k_BT)$. The equilibrium population of the rare tautomer at body temperature is:

$$P^* = \frac{k_{\text{forward}}}{k_{\text{forward}} + k_{\text{back}}} \approx \exp\left(-\frac{\Delta\epsilon}{k_BT}\right) \sim 10^{-4}\text{ to }10^{-5}$$

5.6 From Tunneling to Mutation Rate

The tunneling event alone does not guarantee a mutation. Several conditions must be satisfied simultaneously:

The probability that a tunneling event at a specific base pair leads to a mutation during one cell division is:

$$P_{\text{mut}} = P^* \times P_{\text{replication}} \times (1 - P_{\text{repair}})$$

where:

$P^*$ = equilibrium population of the rare tautomer ($\sim 10^{-4}$ to $10^{-5}$)
$P_{\text{replication}}$ = probability that the polymerase encounters the tautomer during the $\sim 10^{-3}$ s it spends at each nucleotide
$P_{\text{repair}}$ = probability of successful mismatch repair ($\sim 0.999$ for most organisms)

The per-base-pair, per-generation mutation rate from quantum tunneling is therefore:

$$\mu_{\text{quantum}} \approx 10^{-4} \times 1 \times 10^{-3} = 10^{-7}\ \text{per base pair per generation}$$

This is remarkably close to the observed spontaneous transition mutation rate of$\sim 10^{-8}$ to $10^{-9}$ per base pair per generation in many organisms, suggesting that quantum tunneling is a significant contributor to the spontaneous mutation rate.

5.7 Kinetic Isotope Effect as Experimental Test

The most direct experimental signature of quantum tunneling in DNA is the kinetic isotope effect (KIE). Replacing hydrogen (H) with deuterium (D) doubles the tunneling mass, which exponentially suppresses the tunneling rate:

$$\text{KIE} = \frac{k_H}{k_D} = \frac{\nu_0^H}{\nu_0^D} \exp\left(\frac{2a}{\hbar}\left[\sqrt{2m_D \Delta E} - \sqrt{2m_H \Delta E}\right]\right)$$

Since $m_D = 2m_H$, the ratio of attempt frequencies is $\nu_0^H / \nu_0^D = \sqrt{2}$, and the exponential factor gives:

$$\text{KIE} = \sqrt{2}\, \exp\left(\frac{2a}{\hbar}\sqrt{2m_H \Delta E}\left[\sqrt{2} - 1\right]\right)$$

For typical DNA parameters, this predicts KIE values of 3-15, well above the classical maximum of $\sim \sqrt{2} \approx 1.4$. An observed KIE significantly greater than the classical limit is a hallmark of quantum tunneling. Replacing H with tritium (T, mass $3m_H$) provides even stronger discrimination.

5.8 Modern Computational Studies

Supporting Evidence

Godbeer et al. (2015): Open quantum systems approach shows decoherence extends tautomer lifetimes
Slocombe et al. (2021): DFT + instanton theory gives tunneling rates consistent with observed mutation rates
Löwdin's prediction confirmed: Rare tautomers detected in crystal structures and NMR
Temperature independence: Mutation rates in thermophiles match tunneling predictions
Brovarets & Hovorun (2014): Ab initio MD shows concerted double proton transfer in G·C

Open Questions

Concerted vs. stepwise: Do protons tunnel simultaneously or sequentially?
Role of the environment: How do water molecules, ions, and the protein environment modify barriers?
Decoherence effects: Does the warm, wet cellular environment suppress or enhance tunneling?
Epigenetic modifications: Do methylation and other modifications alter tunneling rates?
Evolutionary implications: Does quantum tunneling create a non-random component to mutation?

5.9 Biological Consequences

Implications for Biology and Medicine

Evolution

Quantum tunneling provides a basal spontaneous mutation rate
Acts as a universal mutation mechanism across all organisms
May create a non-uniform mutational spectrum (GC→AT bias)
Sets a fundamental physical limit on genome stability

Cancer and Ageing

Accumulation of tunneling-induced mutations over a lifetime
Hotspot mutations at CpG sites may have tunneling contribution
DNA repair pathways evolved to counteract quantum-induced errors
Age-related increase in mutations partly from stochastic tunneling

The Decoherence Paradox

A surprising result from open quantum systems theory is that decoherence from the biological environment can actually stabilise the rare tautomer rather than destroying it. The Lindblad master equation for the proton density matrix $\rho$ in the double well is:

$$\frac{d\rho}{dt} = -\frac{i}{\hbar}[H, \rho] + \gamma\left(a\rho a^\dagger - \frac{1}{2}\{a^\dagger a, \rho\}\right)$$

where $\gamma$ is the decoherence rate and $a$ is the lowering operator. In the quantum Zeno regime ($\gamma \gg \Delta/\hbar$, where $\Delta$ is the tunnel splitting), frequent environmental measurements suppress coherent oscillations between wells, effectively trapping the proton in whichever well it happens to be in. This means that once a proton tunnels to the rare tautomer, environmental decoherence can prevent it from tunneling back, extending the tautomer lifetime by orders of magnitude.

Key References

Löwdin, P.-O. (1963). Proton tunneling in DNA and its biological implications. Rev. Mod. Phys. 35(3), 724-732.
Slocombe, L., Al-Khalili, J. & Mayall, M. (2021). An open quantum systems approach to proton tunnelling in DNA. Commun. Phys. 4, 196.
Godbeer, A.D., Al-Khalili, J. & Sherrill, D.C. (2015). Modelling proton tunnelling in the adenine-thymine base pair. Phys. Chem. Chem. Phys. 17, 13034.
Brovarets', O.O. & Hovorun, D.M. (2014). Why the tautomerization of the G·C Watson-Crick base pair via double proton transfer does not result in stable tautomers. J. Biomol. Struct. Dyn. 33, 925.
Tolosa, S. et al. (2020). Theoretical study of the proton transfer in the adenine-thymine base pair. J. Mol. Model. 26, 244.
Fang, W. et al. (2021). Inverse temperature dependence of nuclear quantum effects in DNA base pairs. J. Phys. Chem. Lett. 12, 10.1021/acs.jpclett.1c02013.

🌐 Research Platforms & Tools

Quantum-Proteins.ai

Computational platform for quantum simulations of biological systems, including enzyme tunneling, photosynthetic complexes, and DNA quantum mechanics.

Explore Platform →

Pattern.Quantum-Proteins.ai

Pattern recognition and analysis tools for identifying quantum signatures in biological data, including coherence detection and tunneling analysis.

Explore Platform →

💻 Computational Example

Let's simulate quantum tunneling through a biological energy barrier to see how quantum effects enable reactions that would be classically forbidden:

Quantum Tunneling in Enzyme Catalysis

Calculate tunneling probabilities for proton transfer in biological systems. Click Run to execute.

import numpy as np
import matplotlib.pyplot as plt

# Physical constants
hbar = 1.055e-34  # J·s
m_p = 1.673e-27   # proton mass (kg)
eV = 1.602e-19    # electron volt in Joules

# Barrier parameters (typical for enzyme catalysis)
V0 = 0.5 * eV     # barrier height (0.5 eV)
a = 1.0e-10       # barrier width (1 Angstrom)

def tunneling_probability(E, V0, a, m):
    """
    Calculate quantum tunneling probability using WKB approximation
    """
    if E >= V0:
        return 1.0
    kappa = np.sqrt(2 * m * (V0 - E)) / hbar
    T = np.exp(-2 * kappa * a)
    return T

# Energy range (from 0 to barrier height)
E_range = np.linspace(0.01 * eV, V0, 100)

# Calculate tunneling probabilities
T_proton = [tunneling_probability(E, V0, a, m_p) for E in E_range]
T_classical = [1.0 if E >= V0 else 0.0 for E in E_range]

print("Quantum Tunneling in Enzyme Catalysis")
print("=" * 55)
print(f"Barrier height: {V0/eV:.2f} eV")
print(f"Barrier width:  {a*1e10:.2f} Angstrom")
print(f"Particle: proton (mass = {m_p:.2e} kg)")
print()

test_energies = [0.1, 0.2, 0.3, 0.4, 0.45]
print("Tunneling Probabilities:")
print("-" * 55)
for e_frac in test_energies:
    E = e_frac * V0
    T = tunneling_probability(E, V0, a, m_p)
    print(f"E = {e_frac:.2f}*V0 = {E/eV:.3f} eV:  T = {T:.2e}")

print()
print("Temperature Dependence:")
print("-" * 55)
kB = 1.381e-23
for T_kelvin in [273, 300, 310, 350]:
    E_thermal = kB * T_kelvin
    if E_thermal < V0:
        T_tunnel = tunneling_probability(E_thermal, V0, a, m_p)
        print(f"T = {T_kelvin}K: E = {E_thermal/eV:.4f} eV, Tunneling = {T_tunnel:.2e}")

# Visualization
fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 5))

ax1.semilogy(E_range/eV, T_proton, 'b-', linewidth=2, label='Quantum (proton)')
ax1.plot(E_range/eV, T_classical, 'r--', linewidth=2, label='Classical')
ax1.axvline(x=V0/eV, color='gray', linestyle=':', alpha=0.5, label='Barrier height')
ax1.set_xlabel('Energy (eV)', color='white')
ax1.set_ylabel('Transmission Probability', color='white')
ax1.set_title('Quantum Tunneling vs Classical Barrier', color='white')
ax1.legend()
ax1.grid(True, alpha=0.3)
ax1.tick_params(colors='white')

x = np.linspace(-2*a, 3*a, 1000)
V = np.where((x >= 0) & (x <= a), V0, 0)
ax2.plot(x*1e10, V/eV, 'k-', linewidth=2, label='Potential barrier')
ax2.axhline(y=0.3, color='blue', linestyle='--', alpha=0.7, label='E = 0.3*V0')
ax2.fill_between(x*1e10, 0, V0/eV, where=((x>=0) & (x<=a)),
                  alpha=0.2, color='red', label='Classically forbidden')
ax2.set_xlabel('Position (Angstrom)', color='white')
ax2.set_ylabel('Energy (eV)', color='white')
ax2.set_title('Energy Barrier in Biological System', color='white')
ax2.legend()
ax2.grid(True, alpha=0.3)
ax2.set_ylim(-0.1, 0.7)
ax2.tick_params(colors='white')

plt.tight_layout()

Click Run to execute the Python code

First run will download Python environment (~15MB)

📚 Key References

McFadden & Al-Khalili (2014). Life on the Edge: The Coming of Age of Quantum Biology. Crown Publishers.
Engel et al. (2007). Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems. Nature 446, 782-786.
Ball (2011). Physics of life: The dawn of quantum biology. Nature 474, 272-274.
Hore & Mouritsen (2016). The Radical-Pair Mechanism of Magnetoreception. Annu. Rev. Biophys. 45, 299-344.
Cao et al. (2020). Quantum biology revisited. Science Advances 6(14), eaaz4888.