docs/research/quantum-sensing/11-quantum-level-sensors.md
Date: 2026-03-08 Domain: Quantum Sensing × RF Topology × Graph-Based Detection Status: Research Survey
Classical RF sensing using ESP32 WiFi mesh nodes operates at milliwatt power levels with sensitivity limited by thermal noise floors (~-90 dBm). Quantum sensors offer fundamentally different detection mechanisms that can surpass classical limits by orders of magnitude, potentially transforming RF topological sensing from room-scale detection to single-photon field measurement.
This document surveys quantum sensing technologies relevant to RF topological sensing, evaluates their integration potential with the existing RuVector/mincut architecture, and identifies near-term and long-term opportunities.
NV centers are point defects in diamond crystal lattice where a nitrogen atom replaces a carbon atom adjacent to a vacancy. Key properties:
Diamond Crystal with NV Center:
C---C---C---C
| | | |
C---N V---C N = Nitrogen atom
| | | V = Vacancy
C---C---C---C C = Carbon atoms
| | | |
C---C---C---C
ODMR Protocol:
Green Laser → NV → Red Fluorescence
↕
Microwave Drive
Resonance frequency shifts with local B-field
ΔfNV = γNV × B_local
γNV = 28 GHz/T
SQUID Loop:
┌──────[JJ1]──────┐
│ │ JJ = Josephson Junction
│ Φ_ext → │ Φ = Magnetic flux
│ (flux) │
│ │ V = Φ₀/(2π) × dφ/dt
└──────[JJ2]──────┘ Φ₀ = 2.07 × 10⁻¹⁵ Wb
Critical current: Ic = 2I₀|cos(πΦ_ext/Φ₀)|
Voltage oscillates with period Φ₀
Atoms excited to high principal quantum number (n > 30) become extraordinarily sensitive to electric fields:
Rydberg EIT Level Scheme:
|r⟩ -------- Rydberg state (n~50) ← RF field couples |r⟩↔|r'⟩
↕ Ωc (coupling laser)
|e⟩ -------- Excited state
↕ Ωp (probe laser)
|g⟩ -------- Ground state
Without RF: EIT window → transparent to probe
With RF: Autler-Townes splitting → absorption changes
Splitting: Ω_RF = μ_rr' × E_RF / ℏ
where μ_rr' = n² × e × a₀ (scales as n²!)
Spin-exchange relaxation-free (SERF) magnetometers using alkali vapor:
| Sensor Type | Sensitivity | Temp | Bandwidth | Size | Cost Est. |
|---|---|---|---|---|---|
| NV Diamond | ~1 pT/√Hz | 300K | DC-10 GHz | cm | $1K-10K |
| SQUID | ~1 fT/√Hz | 4-77K | DC-1 GHz | cm | $10K-100K |
| Rydberg | ~1 µV/m/√Hz | 300K | DC-THz | 10 cm | $5K-50K |
| SERF | ~0.16 fT/√Hz | 420K | DC-1 kHz | cm | $5K-50K |
| ESP32 (classical) | ~-90 dBm | 300K | 2.4/5 GHz | cm | $5 |
Classical RF detection is limited by thermal (Johnson-Nyquist) noise:
Classical thermal noise floor:
P_noise = k_B × T × B
At T = 300K, B = 20 MHz (WiFi channel):
P_noise = 1.38e-23 × 300 × 20e6 = 8.3 × 10⁻¹⁴ W
P_noise = -101 dBm
Shot noise limit (coherent state):
ΔE = √(ℏω/(2ε₀V)) per photon
SNR_shot ∝ √N_photons
Heisenberg limit (entangled state):
SNR_Heisenberg ∝ N_photons
Quantum advantage: √N improvement over shot noise
For N = 10⁶ photons → 1000× SNR improvement
The quantum advantage for RF sensing depends on the signal regime:
| Regime | Classical | Quantum | Advantage |
|---|---|---|---|
| Strong signal (>-60 dBm) | Adequate | Unnecessary | None |
| Medium (-60 to -90 dBm) | Noisy | Cleaner | 10-100× SNR |
| Weak (<-90 dBm) | Undetectable | Detectable | Enabling |
| Single-photon | Impossible | Feasible | Infinite |
For RF topological sensing, the quantum advantage is most relevant for:
Squeezed States: Reduce noise in one quadrature at expense of other:
ΔX₁ × ΔX₂ ≥ ℏ/2
Squeeze X₁: ΔX₁ = e⁻ʳ × √(ℏ/2) (reduced)
ΔX₂ = e⁺ʳ × √(ℏ/2) (increased)
For r = 2 (17.4 dB squeezing):
Noise reduction in amplitude: 7.4×
Demonstrated: 15 dB squeezing (LIGO)
Quantum Error Correction: Protect quantum states from decoherence:
Rydberg atoms provide the most promising near-term quantum RF sensor for topological sensing because:
Rydberg Sensor Architecture:
┌─────────────────────────────┐
│ Cesium Vapor Cell │
│ │
│ Probe (852nm) ───────→ │──→ Photodetector
│ Coupling (509nm) ───→ │
│ │
│ ↕ RF field enters │
└─────────────────────────────┘
Frequency tuning:
n=30: ~300 GHz transitions
n=50: ~50 GHz transitions
n=70: ~10 GHz transitions (WiFi band!)
n=100: ~1 GHz transitions
For 2.4 GHz detection using Rydberg states near n=70:
Transition dipole moment:
μ = n² × e × a₀ ≈ 70² × 1.6e-19 × 5.3e-11
μ ≈ 4.1 × 10⁻²⁶ C·m
Minimum detectable field:
E_min = ℏ × Γ / (2μ)
where Γ = EIT linewidth ≈ 1 MHz
E_min ≈ 1.05e-34 × 2π × 1e6 / (2 × 4.1e-26)
E_min ≈ 8 µV/m
Compare to ESP32 sensitivity: ~1 mV/m
Quantum advantage: ~125× in field sensitivity
Key milestones in Rydberg RF sensing:
Hybrid Rydberg-ESP32 Architecture:
Classical Layer (ESP32 mesh):
┌────┐ ┌────┐ ┌────┐
│ESP1│────│ESP2│────│ESP3│ 120 classical edges
└────┘ └────┘ └────┘ CSI coherence weights
│ │ │
│ ┌────┴────┐ │
└────│Rydberg │────┘ Quantum sensor node
│ Sensor │ High-sensitivity edges
└─────────┘
The Rydberg sensor provides:
1. Ultra-sensitive reference measurements
2. Ground truth calibration for classical edges
3. Detection of sub-threshold perturbations
4. Phase reference for coherence estimation
Quantum illumination uses entangled photon pairs to detect objects in noisy environments:
Protocol:
1. Generate entangled signal-idler pair: |Ψ⟩ = Σ cₙ|n⟩_S|n⟩_I
2. Send signal photon toward target, keep idler
3. Collect reflected signal (buried in thermal noise)
4. Joint measurement on returned signal + stored idler
Classical detection: SNR = N_S / N_B
Quantum detection: SNR = N_S × (N_B + 1) / N_B
Advantage: 6 dB in error exponent (factor of 4)
Critical: Advantage persists even when entanglement is destroyed
by the noisy channel (unlike most quantum protocols)
For RF topological sensing at 2.4 GHz:
Microwave entangled source:
Josephson Parametric Amplifier (JPA)
→ Generates entangled microwave-microwave pairs
→ Or microwave-optical pairs (for optical idler storage)
Challenge: thermal photon number at 2.4 GHz, 300K:
n_th = 1/(exp(hf/kT) - 1) = 1/(exp(4.8e-5) - 1) ≈ 2600
Background: ~2600 thermal photons per mode
→ Classical detection hopeless for single-photon signals
→ Quantum illumination still provides 6 dB advantage
Quantum illumination could enhance RF topological sensing by:
Quantum walks are the quantum analog of random walks, with superposition and interference:
Continuous-time quantum walk on graph G:
|ψ(t)⟩ = e^{-iHt} |ψ(0)⟩
where H = adjacency matrix A or Laplacian L
Key property: Quantum walk spreads quadratically faster
Classical: ⟨x²⟩ ~ t (diffusive)
Quantum: ⟨x²⟩ ~ t² (ballistic)
For graph topology detection:
- Walk dynamics encode graph structure
- Interference patterns reveal symmetries
- Hitting times indicate connectivity
Grover-accelerated graph search:
Classical min-cut (Stoer-Wagner): O(VE + V² log V)
For V=16, E=120: ~4,000 operations
Quantum search for min-cut:
Use Grover's algorithm to search over cuts
Number of possible cuts: 2^V = 2^16 = 65,536
Classical brute force: O(2^V) = 65,536 evaluations
Quantum (Grover): O(√(2^V)) = 256 evaluations
Quadratic speedup for brute-force approach
However: For V=16, Stoer-Wagner (4,000 ops) beats Grover (256 oracle calls)
because each oracle call has overhead
Quantum advantage threshold: V > ~100 nodes
Quantum spectral analysis:
Quantum Phase Estimation (QPE) for graph Laplacian:
Input: L = D - A (graph Laplacian)
Output: eigenvalues λ₁ ≤ λ₂ ≤ ... ≤ λ_V
Fiedler value λ₂ → algebraic connectivity
Cheeger inequality: λ₂/2 ≤ h(G) ≤ √(2λ₂)
where h(G) = min-cut / min-volume (Cheeger constant)
QPE complexity: O(poly(log V)) per eigenvalue
Classical: O(V³) for full eigendecomposition
Quantum advantage for spectral analysis: exponential
for V >> 100
Variational Quantum Eigensolver (VQE) for normalized cut:
Minimize: NCut = cut(A,B) × (1/vol(A) + 1/vol(B))
Encode as QUBO:
min x^T Q x where x ∈ {0,1}^V
Q_ij = -w_ij + d_i × δ_ij × balance_penalty
Map to Ising Hamiltonian:
H = Σ_ij J_ij σ_i^z σ_j^z + Σ_i h_i σ_i^z
Solve with:
- VQE (gate-based): variational ansatz circuit
- QAOA: alternating cost/mixer unitaries
- Quantum annealing (D-Wave): native QUBO solver
Not every edge in the RF sensing graph benefits from quantum sensing. The advantage is concentrated in specific scenarios:
| Scenario | Classical | Quantum | Benefit |
|---|---|---|---|
| Strong LOS links | Adequate | Overkill | None |
| Weak NLOS links | Noisy/lost | Detectable | Enables new edges |
| Sub-threshold perturbations | Invisible | Detectable | Breathing, heartbeat |
| Phase coherence measurement | Clock-limited | Fundamental | Better edge weights |
| Multi-target disambiguation | Ambiguous | Resolvable | More accurate cuts |
Three-Tier Hybrid Sensing:
Tier 1: ESP32 Classical Mesh (16 nodes, $80 total)
┌─────────────────────────────────────┐
│ Standard CSI extraction │
│ 120 TX-RX edges │
│ ~30-60 cm resolution │
│ Person-scale detection │
└──────────────┬──────────────────────┘
│
Tier 2: NV Diamond Enhancement (4 nodes, ~$20K)
┌──────────────┴──────────────────────┐
│ pT-level magnetic field sensing │
│ Room-temperature operation │
│ Complements RF with B-field edges │
│ Breathing/heartbeat detection │
└──────────────┬──────────────────────┘
│
Tier 3: Rydberg Reference (1 node, ~$50K)
┌──────────────┴──────────────────────┐
│ µV/m electric field sensitivity │
│ Self-calibrated SI-traceable │
│ Ground truth for classical edges │
│ Sub-threshold perturbation detect │
└─────────────────────────────────────┘
Graph construction:
G_hybrid = G_classical ∪ G_magnetic ∪ G_quantum
Edge weight fusion:
w_ij = α × w_classical + β × w_magnetic + γ × w_quantum
where α + β + γ = 1, learned per-edge
Classical edge weight (ESP32):
w_ij = coherence(CSI_i→j)
Noise floor: ~-90 dBm
Phase noise: ~5° RMS (clock drift limited)
Quantum-enhanced edge weight:
w_ij = f(CSI_ij, B_field_ij, E_field_ij)
NV contribution:
- Local magnetic field map at pT resolution
- Detects metallic object perturbations
- Measures eddy current signatures
Rydberg contribution:
- Electric field at µV/m resolution
- Phase-accurate reference measurement
- Calibrates classical CSI phase errors
Quantum sensors naturally measure their environment through decoherence:
NV Center Decoherence:
T₁ (spin-lattice relaxation): ~6 ms at 300K
T₂ (spin-spin dephasing): ~1 ms at 300K
T₂* (inhomogeneous): ~1 µs
Environmental perturbation → T₂* change
Sensitivity:
ΔB_min = (1/γ) × 1/(T₂* × √(η × T_meas))
where η = photon collection efficiency
T_meas = measurement time
At η=0.1, T_meas=1s:
ΔB_min ≈ 1 pT
The key insight: decoherence signatures encode environmental structure. Different objects and materials produce different decoherence profiles:
| Object | Decoherence Mechanism | Signature |
|---|---|---|
| Metal | Eddy currents, Johnson noise | T₂* reduction, broadband |
| Human body | Ionic currents, diamagnetism | T₁ modulation, low-freq |
| Water | Diamagnetic susceptibility | Subtle T₂ shift |
| Electronics | EM emission | Discrete frequency peaks |
Quantum Fisher Information (QFI):
F_Q(θ) = 4(⟨∂_θψ|∂_θψ⟩ - |⟨ψ|∂_θψ⟩|²)
Quantum Cramér-Rao Bound:
Var(θ̂) ≥ 1/(N × F_Q(θ))
For sensor placement optimization:
- Compute F_Q at each candidate position
- Place quantum sensors where F_Q is maximized
- Typically: room center, doorways, narrow passages
Optimal placement for V=16 classical + 4 quantum:
┌─────────────────────────┐
│ E E E E E E │ E = ESP32 (perimeter)
│ │
│ E Q Q E │ Q = Quantum sensor
│ │ (high-FI positions)
│ E Q Q E │
│ │
│ E E E E E E │
└─────────────────────────┘
Quantum Graph Neural Network:
Input: Edge weights w_ij from RF sensing graph
Encoding: Amplitude encoding of adjacency matrix
|ψ_G⟩ = Σ_ij w_ij |i⟩|j⟩ / ||w||
Variational circuit:
U(θ) = Π_l [U_entangle × U_rotation(θ_l)]
U_rotation: R_y(θ₁) ⊗ R_y(θ₂) ⊗ ... ⊗ R_y(θ_V)
U_entangle: CNOT cascade matching graph topology
Measurement: ⟨Z₁⟩ → occupancy classification
Training: Minimize L = Σ (y - ⟨Z₁⟩)² via parameter-shift rule
For V=16: Requires 16 qubits + ~100 variational parameters
→ Within reach of current NISQ devices (IBM Eagle: 127 qubits)
Quantum kernel for CSI feature space:
Encode CSI vector x into quantum state: |φ(x)⟩ = U(x)|0⟩
Kernel: K(x, x') = |⟨φ(x)|φ(x')⟩|²
Properties:
- Maps to exponentially large Hilbert space
- Can capture correlations classical kernels miss
- Computed on quantum hardware, used in classical SVM/GP
For edge classification (stable/unstable/transitioning):
- Encode temporal CSI window as quantum state
- Quantum kernel captures phase correlations
- Classical SVM classifies using quantum kernel values
Quantum Reservoir for Temporal RF Patterns:
RF Signal → Quantum System → Measurement → Classical Readout
Reservoir: N coupled qubits with natural dynamics
H_res = Σ_i h_i σ_i^z + Σ_ij J_ij σ_i^z σ_j^z + Σ_i Ω_i σ_i^x
Input: CSI values modulate h_i (local fields)
Dynamics: ρ(t+1) = U × ρ(t) × U† + noise
Output: Measure ⟨σ_i^z⟩ for all qubits → feature vector
Advantages for temporal RF sensing:
- Natural temporal memory (quantum coherence)
- No training of reservoir (only readout layer)
- Captures non-linear temporal correlations
- Matches temporal graph evolution naturally
Min-cut as QUBO on D-Wave:
Variables: x_i ∈ {0,1} (node partition assignment)
Objective: minimize Σ_ij w_ij × x_i × (1-x_j)
QUBO matrix:
Q_ij = -w_ij (off-diagonal)
Q_ii = Σ_j w_ij (diagonal)
D-Wave Advantage2: 7,000+ qubits
→ Can handle graphs up to ~3,500 nodes
→ Our V=16 graph trivially fits
Practical consideration:
- Cloud API access: ~$2K/month
- Annealing time: ~20 µs per sample
- 1000 samples for statistics: ~20 ms
- Compatible with 20 Hz update rate
Multi-cut extension (k-way):
Use k binary variables per node
→ 16 × k = 48 qubits for 3-person detection
Variational Quantum Eigensolver for Laplacian spectrum:
Goal: Find smallest eigenvalues of L = D - A
Ansatz: |ψ(θ)⟩ = U(θ)|0⟩^⊗n
Cost: E(θ) = ⟨ψ(θ)|L|ψ(θ)⟩
Optimization: θ* = argmin E(θ) via classical optimizer
For Fiedler value (λ₂):
1. Find ground state |v₁⟩ (constant vector, known)
2. Constrain ⟨v₁|ψ⟩ = 0
3. Minimize in orthogonal subspace → λ₂
Application: Track λ₂ over time
- λ₂ large → graph well-connected → no obstruction
- λ₂ drops → graph nearly disconnected → boundary detected
- Rate of λ₂ change → speed of perturbation
Quantum Approximate Optimization Algorithm:
Cost Hamiltonian: H_C = Σ_ij w_ij (1 - Z_i Z_j) / 2
Mixer Hamiltonian: H_M = Σ_i X_i
p-layer circuit:
|ψ(γ,β)⟩ = Π_l [e^{-iβ_l H_M} × e^{-iγ_l H_C}] |+⟩^⊗n
For p=1: Guaranteed approximation ratio r ≥ 0.6924 for MaxCut
For p=3-5: Near-optimal for small graphs
Our V=16 graph: 16 qubits, p=3 → 96 parameters
→ Trainable on current hardware
→ Could provide better-than-classical cuts in some cases
Integration Pipeline:
ESP32 Mesh Quantum Sensors
┌──────────┐ ┌──────────┐
│ CSI Data │ │ QSensor │
│ 120 edges│ │ 4 nodes │
│ 20 Hz │ │ 100 Hz │
└────┬─────┘ └────┬─────┘
│ │
▼ ▼
┌──────────────────────────────┐
│ Edge Weight Fusion │
│ │
│ w_ij = fuse( │
│ classical_coherence, │
│ magnetic_perturbation, │
│ quantum_phase_ref │
│ ) │
└──────────────┬───────────────┘
│
▼
┌──────────────────────────────┐
│ RfGraph Construction │
│ G = (V_classical ∪ V_quantum, E_fused)
└──────────────┬───────────────┘
│
▼
┌──────────────────────────────┐
│ Hybrid Mincut │
│ - Classical: Stoer-Wagner │
│ - Or quantum: D-Wave QUBO │
│ - Select based on graph size│
└──────────────┬───────────────┘
│
▼
┌──────────────────────────────┐
│ RuVector Temporal Store │
│ - Graph evolution history │
│ - Quantum measurement log │
│ - Attention-weighted fusion │
└──────────────────────────────┘
/// Quantum sensor integration for RF topological sensing
pub trait QuantumSensor: Send + Sync {
/// Get current measurement with uncertainty
fn measure(&self) -> QuantumMeasurement;
/// Sensor sensitivity in appropriate units
fn sensitivity(&self) -> f64;
/// Decoherence time (characterizes environment)
fn coherence_time(&self) -> Duration;
}
pub struct QuantumMeasurement {
pub value: f64,
pub uncertainty: f64, // Quantum uncertainty
pub fisher_information: f64, // QFI for this measurement
pub timestamp: Instant,
pub sensor_type: QuantumSensorType,
}
pub enum QuantumSensorType {
NVDiamond { t2_star: Duration },
Rydberg { principal_n: u32, transition_freq: f64 },
SQUID { flux_quantum: f64 },
SERF { vapor_temp: f64 },
}
/// Fuse classical and quantum edge weights
pub trait HybridEdgeWeightFusion {
fn fuse(
&self,
classical: &ClassicalEdgeWeight,
quantum: Option<&QuantumMeasurement>,
) -> FusedEdgeWeight;
}
pub struct FusedEdgeWeight {
pub weight: f64,
pub confidence: f64, // Higher with quantum data
pub classical_contribution: f64,
pub quantum_contribution: f64,
pub fisher_bound: f64, // QCRB on precision
}
| Technology | Current TRL | Field-Ready | Clinical | Notes |
|---|---|---|---|---|
| NV Diamond magnetometer | TRL 5-6 | 2026-2028 | 2030+ | Room temp, most practical |
| Chip-scale NV | TRL 3-4 | 2028-2030 | 2032+ | Integration with CMOS |
| Rydberg RF receiver | TRL 4-5 | 2027-2029 | N/A | Military interest high |
| Miniature SQUID | TRL 7-8 | Available | Available | Requires cryogenics |
| SERF magnetometer | TRL 5-6 | 2026-2028 | 2029+ | Needs shielding |
| Quantum annealer (D-Wave) | TRL 8-9 | Available | N/A | Cloud access now |
| NISQ processor (IBM/Google) | TRL 6-7 | 2026+ | N/A | 1000+ qubits by 2026 |
Current vs Projected SWaP:
NV Diamond Sensor (2025):
Size: 15 × 10 × 10 cm
Weight: 2 kg
Power: 5 W (laser + electronics)
NV Diamond Sensor (2028 projected):
Size: 5 × 3 × 3 cm
Weight: 200 g
Power: 1 W
Rydberg Vapor Cell (2025):
Size: 20 × 15 × 15 cm
Weight: 3 kg
Power: 10 W (two lasers + control)
Chip-Scale Rydberg (2030 projected):
Size: 3 × 3 × 1 cm
Weight: 50 g
Power: 0.5 W
Compare ESP32:
Size: 5 × 3 × 0.5 cm
Weight: 10 g
Power: 0.44 W
Phase 1 (2026): Classical-only RF topology
- 16 ESP32 nodes
- Stoer-Wagner mincut
- Proof of concept
Phase 2 (2027-2028): Quantum-enhanced
- 16 ESP32 + 2-4 NV diamond nodes
- Hybrid edge weights
- Sub-threshold detection (breathing)
Phase 3 (2029-2030): Full quantum integration
- 16 ESP32 + 4 NV + 1 Rydberg
- Quantum-classical graph fusion
- D-Wave cloud for multi-cut optimization
Phase 4 (2031+): Quantum-native
- Chip-scale quantum sensors at every node
- On-device quantum processing
- Room-scale coherence imaging
Quantum advantage threshold: At what graph size does quantum mincut outperform classical? Preliminary analysis suggests V > 100, but constant factors matter.
Decoherence as feature: Can quantum decoherence rates serve as edge weights directly, bypassing classical CSI entirely?
Entanglement distribution: Can entangled sensor pairs provide correlated edge weights with fundamentally lower uncertainty?
Quantum memory for temporal graphs: Can quantum memory store graph evolution states more efficiently than classical RuVector?
Noise budget: In a real room with WiFi, Bluetooth, and power line interference, what is the practical quantum advantage?
Calibration: How often do quantum sensors need recalibration in field deployment?
Cost trajectory: When will quantum sensor nodes reach $100/unit for mass deployment?
Hybrid optimization: What is the optimal ratio of classical to quantum nodes for a given room size and detection requirement?
Resolution limits: Does quantum sensing fundamentally change the 30-60 cm resolution bound, or only improve SNR within the same Fresnel-limited resolution?
Multi-room scaling: Can quantum entanglement between rooms provide correlated sensing that classical links cannot?
Adversarial robustness: Are quantum-enhanced edge weights more robust against deliberate spoofing or jamming?
Quantum sensing represents a paradigm shift for RF topological sensing. While the classical ESP32 mesh provides adequate sensitivity for person-scale detection, quantum sensors enable:
The most practical near-term integration path uses NV diamond sensors (room temperature, pT sensitivity) as enhancement nodes within the classical ESP32 mesh, with Rydberg sensors providing calibration references. Quantum computing (D-Wave, NISQ) offers immediate value for graph cut optimization at scale.
The long-term vision is a quantum-native sensing mesh where every node performs quantum measurements, edge weights encode quantum coherence between nodes, and graph algorithms run on quantum hardware — a true quantum radio nervous system.