Concepts Intermediate Free 42/53 in series 16 min read

Quantum Sensing: How Quantum Systems Detect the Undetectable

Explore how quantum superposition and entanglement enable sensors that exceed classical noise limits, from NV center magnetometers to atomic clocks, with a Qiskit simulation of the Ramsey spectroscopy protocol.

What you'll learn

  • quantum sensing
  • quantum metrology
  • Heisenberg limit
  • NV centers
  • atomic clocks

Prerequisites

  • Python proficiency
  • Beginner quantum computing concepts (superposition, entanglement)
  • Linear algebra basics

Sensing at the Edge of What Physics Allows

Every measurement is ultimately limited by noise. Classical sensors hit the shot noise limit, the unavoidable statistical uncertainty from using N independent probes. Quantum sensing exploits superposition and entanglement to push beyond this boundary, reaching the Heisenberg limit and enabling measurements of unprecedented precision.

This is not speculative. Quantum sensors are already deployed today: atomic clocks in GPS satellites, gravitational wave detectors at LIGO, and medical imaging systems that use quantum interference to measure brain activity. Understanding the physics behind these instruments reveals both how they work and why classical engineering cannot replicate them.

The Shot Noise Limit and the Heisenberg Limit

Imagine you have N particles and want to estimate an unknown phase phi. If you measure them classically and independently, the statistical uncertainty averages down as the square root of N:

Delta_phi (shot noise) = 1 / sqrt(N)

This is the standard quantum limit (SQL) or shot noise limit. It is the best any classical strategy can achieve.

Now entangle those N particles. With the right entangled state (a GHZ state, for example), the phase imprints on all N particles coherently. The measurement uncertainty shrinks to:

Delta_phi (Heisenberg) = 1 / N

This is the Heisenberg limit. For N = 1000 particles, the Heisenberg limit gives 31x better precision than the shot noise limit. The gap grows as sqrt(N), so it becomes dramatic at large N.

The catch: entangled states are fragile. Any decoherence destroys the advantage. Current quantum sensors operate close to the shot noise limit with some quantum enhancement, and the research frontier is about extending coherence times to capture the full Heisenberg scaling.

Nitrogen-Vacancy Centers in Diamond

The nitrogen-vacancy (NV) center is a point defect in diamond where two adjacent carbon atoms are replaced by a nitrogen atom and a vacant lattice site. Its electron spin behaves like a qubit at room temperature, a property almost no other solid-state system can claim.

The NV center’s ground state has spin S=1 with three sublevels: m_s = 0, +1, -1. Applying a microwave pulse at 2.87 GHz drives transitions between these states. The key trick: the m_s = 0 state fluoresces more brightly under green laser illumination than the m_s = +1 or -1 states. This means you can read out the spin state optically, no dilution refrigerator needed.

A magnetic field splits the m_s = +1 and m_s = -1 energy levels by an amount proportional to the field strength (the Zeeman effect). By measuring the resonance frequency shift, you can infer the magnetic field with femtotesla sensitivity. For context, a femtotesla is 10^-15 Tesla; the Earth’s magnetic field is about 50 microtesla, so NV magnetometers can detect signals nine orders of magnitude weaker.

Applications: Detecting the magnetic fields produced by neural currents in the brain (magnetoencephalography), imaging magnetic domains in materials science, measuring fields inside individual living cells.

Atomic Clocks and GPS

The best atomic clocks use cesium or ytterbium atoms cooled to microkelvin temperatures. The clock transition (the hyperfine transition in cesium at 9.192631770 GHz) is the SI definition of the second.

Modern optical lattice clocks trap neutral atoms in standing-wave laser fields and probe a narrow optical transition. The stability of these clocks reaches 10^-18 fractional frequency uncertainty, meaning they would neither gain nor lose a second over the age of the universe.

GPS depends on atomic clocks: each satellite carries multiple cesium and rubidium clocks. A timing error of 1 nanosecond translates to a position error of 30 centimeters. The more precisely the clocks tick, the more accurately your phone knows where you are.

Entanglement-enhanced clocks using spin-squeezed states of many atoms have been demonstrated in research labs, pushing below the shot noise limit while maintaining long coherence times.

Quantum Gravimeters and Accelerometers

Atom interferometry uses matter waves to measure acceleration and gravity. A cloud of cold atoms is split into two paths by laser pulses, the two paths travel different routes through a gravitational field, then recombine. The interference pattern depends on the gravitational phase accumulated on each path.

This gives absolute measurements of g (gravitational acceleration) to better than 10^-9 g precision, roughly 1000x better than commercial gravimeters. Applications include:

  • Underground navigation: Detecting tunnels, voids, and density variations in rock, enabling navigation without GPS in underground environments
  • Earthquake detection: Measuring slow ground tilts and precursory strain fields before seismic events
  • Oil and gas exploration: Mapping subsurface density to locate reservoirs without drilling
  • Fundamental physics: Testing the equivalence principle and searching for new gravitational-scale interactions

The Ramsey Spectroscopy Protocol

Ramsey spectroscopy is the key measurement protocol underlying both atomic clocks and NV magnetometry. It converts a phase accumulated during free evolution into a measurable population difference.

The protocol has four steps:

  1. Initialize: Prepare the qubit in the |0> state (ground state)
  2. First pi/2 pulse: Rotate the qubit into an equal superposition (|0> + |1>) / sqrt(2); put the state on the equator of the Bloch sphere
  3. Free evolution: Let the qubit precess for time T under the field being measured; it accumulates a phase phi = omega * T
  4. Second pi/2 pulse: Apply another pi/2 rotation to convert the phase into a population difference
  5. Measurement: Read out the probability of the |1> state

The final |1> probability is sin^2(phi/2). By scanning the precession time T or the reference frequency, you map out the resonance and extract the unknown field.

Simulating Ramsey Spectroscopy with Qiskit

import numpy as np
from qiskit import QuantumCircuit
from qiskit_aer import AerSimulator
import matplotlib.pyplot as plt

def ramsey_circuit(phase: float) -> QuantumCircuit:
    """
    Ramsey sequence: H - phase rotation - H - measure.
    The phase rotation simulates precession under an external field.
    """
    qc = QuantumCircuit(1, 1)
    qc.h(0)           # First pi/2 pulse (Hadamard)
    qc.rz(phase, 0)   # Free precession: accumulate phase phi
    qc.h(0)           # Second pi/2 pulse
    qc.measure(0, 0)
    return qc

simulator = AerSimulator()
phases = np.linspace(0, 2 * np.pi, 50)
probabilities = []

for phi in phases:
    qc = ramsey_circuit(phi)
    job = simulator.run(qc, shots=2000)
    counts = job.result().get_counts()
    p1 = counts.get("1", 0) / 2000
    probabilities.append(p1)

# The result is a cosine fringe -- the signature of quantum interference
print("Phase -> P(|1>) [first 5 points]:")
for phi, p in zip(phases[:5], probabilities[:5]):
    print(f"  phi={phi:.3f} rad  ->  P(1)={p:.3f}  (theory: {0.5 * (1 - np.cos(phi)):.3f})")

The measured probability P(|1>) = (1 - cos(phi)) / 2 traces a cosine fringe as a function of phase. An unknown field shifts this fringe. By fitting the fringe, you extract the unknown frequency or field strength.

Phase Kickback as a Sensing Primitive

A related quantum sensing technique uses phase kickback, the same mechanism underlying the quantum phase estimation algorithm. Here a target qubit in an eigenstate imprints its eigenvalue phase onto a control qubit, which can then be read out.

from qiskit import QuantumCircuit

def phase_kickback_sensor(unknown_phase: float) -> QuantumCircuit:
    """
    Uses phase kickback to sense an unknown phase.
    Control qubit in superposition kicks back the phase from a U gate.
    """
    qc = QuantumCircuit(2, 1)

    # Prepare control in superposition
    qc.h(0)

    # Prepare target in eigenstate of U (here, |1> for a Z-type gate)
    qc.x(1)

    # Apply controlled-U (simulates controlled interaction with field)
    # The eigenvalue of U on |1> is exp(i * unknown_phase)
    qc.cp(unknown_phase, 0, 1)

    # Interfere and read out control
    qc.h(0)
    qc.measure(0, 0)

    return qc

# Test at several phases
test_phases = [0, np.pi/4, np.pi/2, np.pi]
simulator = AerSimulator()

for phi in test_phases:
    qc = phase_kickback_sensor(phi)
    job = simulator.run(qc, shots=4000)
    counts = job.result().get_counts()
    p1 = counts.get("1", 0) / 4000
    print(f"phi={phi:.4f} rad  ->  P(|1>)={p1:.3f}  theory={(1 - np.cos(phi))/2:.3f}")

This single-shot phase estimation circuit is the simplest instance of quantum metrology: one qubit, one entangling gate, one measurement, yet it extracts the phase with shot-noise-limited precision.

Real-World Deployment Status

Quantum sensing is the most commercially mature branch of quantum technology:

  • Atomic clocks are in GPS satellites and are the primary frequency standard worldwide
  • SQUID magnetometers (superconducting quantum interference devices) are used in hospital MEG (magnetoencephalography) systems for brain imaging
  • NV center diamonds have been used to image magnetic fields in individual neurons and to detect nuclear magnetic resonance signals from nanoscale samples
  • Atom interferometer gravimeters have been tested in underground mapping missions and aboard aircraft for terrain mapping

The transition from laboratory demonstrations to fielded instruments is well underway. The remaining engineering challenges (miniaturization, ruggedization, reducing laser and electronics complexity) are primarily classical engineering problems rather than quantum physics ones.

Was this tutorial helpful?

Continue learning with these guides