
2
ing characterized the bath of two nearby electronic spins
in diamond, we are able to probe, with nanoscale spatial
resolution, the dominant source of noise common to both
qubits arising from the quasistatic many-body electronic
spin bath. The noise model reveals the local spin density
and timescale of spin bath dynamics with nanoscale vari-
ations, information which is inaccessible by conventional
nuclear magnetic resonance (NMR) or ensemble-sensor
techniques.
II. QUANTUM NOISE SPECTROSCOPY
Several protocols for noise spectroscopy have been de-
veloped thus far, ranging from simple sequences [12–
14] to more complex continuous [15–17] and pulsed [18–
21] control. They have successfully elucidated noise
sources (from local fluctuators [18, 22–25] to spin en-
vironments [12–14, 19, 20, 26]), and their accuracy to
reproduce a given classical noise has been evaluated [27].
However, much less attention has been paid to analyze
their predictive power especially when the reconstructed
noise spectrum is only an approximation to the real noise,
i.e., whether because it arises from a quantum system [28]
or a complex classical source [29–31]—or more simply due
to experimental limitations. Here, to achieve a predictive
noise model, we propose to build a self-consistent noise
spectrum by combining complementary approaches.
The simplest approach, which we call R-E-noise spec-
troscopy, utilizes only decoherence under the free evolu-
tion [Ramsey, (R)] and spin echo (E) experiments. The
knowledge of their decay functionals and decay times
T∗
2(R) and T2(E) may be sufficient to fully character-
ize a noise model S(ω|~p) with unknown model param-
eters ~p [32]. While minimal in experimental cost, this
method requires a noise model that is already known
and sufficiently simple to uniquely identify ~p [12–14].
Furthermore, it can only investigate low-frequency noise
(ω < T −1
2).
A more general approach based on dynamical-
decoupling sequences with equidistant πpulses [Carr-
Purcell-Meiboom-Gill (CPMG) pulse sequences] can in
principle reconstruct the full noise spectrum. Under
the filter-function formalism, each CPMG experiment of
inter-pulse length 2τmforms a filter |˜
fT(ω)|2that approx-
imates a delta function δ(ω−ωm), ωm= (2π)(4τm)−1.
This allows direct measurement of S(ωm) from the
simple-exponential decay χm(T) under CPMG pulse se-
quences, where
χm(T) = 1
2ZS(ω)|˜
fT(ω)|2dω
2π≈4
π2S(ωm)T. (1)
While this method can characterize arbitrary, unknown
noise spectra with high-resolution, it comes at increased
experimental cost, as one CPMG experiment is needed
per frequency. Furthermore, the bandwidth, while much
broader, is still bounded by the coherence time T2and
Rabi frequency Ω0,T−1
2< ωmΩ0[20]. In particular,
low frequencies are harder to reach in the presence of
strong noise.
III. SELF-CONSISTENT NOISE
CHARACTERIZATION
Combining these techniques, we demonstrate how to
obtain a self-consistent classical model. We start with a
minimal noise model, consistent with initial experimental
data, and incrementally refine it as necessary to be con-
sistent with additional experiments. While other strate-
gies are possible, this minimizes the experimental cost.
We first demonstrate the protocol in the concrete case of
an NV center in diamond (Fig. 1).
A. NV electronic spin qubit
The first step is to measure the NV Ramsey dynamics.
We used the ms={0,−1}states of the NV electronic
spin (electronic spin S= 1) in an external static magnetic
field of strength B0≈350 G aligned approximately along
the NV axis. The control was achieved with a single-tone,
resonant microwave of ΩNV
0≈6.9 MHz amplitude to drive
both 15NV hyperfine transitions (Azz ≈3.2 MHz).
Observing a Gaussian decay under Ramsey control
[Fig. 2(b)], we assume as our minimal model an Ornstein-
Uhlenbeck (OU) process
S(ω|b, τc) = b2(2τc)
1+(ωτc)2,(2)
characterized by two parameters (b, τc). Indeed, a qua-
sistatic or “slow” OU noise, (bsτs)1, predicts a Gaus-
sian decay, χR(T) = (bsT)2/2≡(T /T ∗
2)2. More gen-
erally, the slow-OU noise has successfully modeled noise
from a slowly fluctuating spin bath [13, 14, 26] and is
expected [6] to be the dominant noise in our system [33].
Then, fitting for T∗
2we identify one of two unknown pa-
rameters, bs= 0.56(2) MHz.
Given a working model S0=Ssconsistent with Ram-
sey dynamics, we can ask whether it is already pre-
dictive of echo dynamics. Unfortunately, we find that
it is not, as while S0predicts a stretched-exponential
χE(T)≈(b2
sT3)/(12τs)≡(T/T2)3, the NV echo is domi-
nantly simple exponential [Fig. 2(c)]. Note that similarly
we could have started with the knowledge of NV echo
decay to first search for a minimal (single-termed) noise
model consistent with echo dynamics and test whether
it is predictive of Ramsey dynamics. In such a case, we
would arrive at either a fast-OU noise Sf(τfT) or
white-noise Sw, which both yield an exponential decay.
However, neither are consistent with NV Ramsey dynam-
ics.
This suggests that the environment around the NV is
sufficiently complex so as not to be reduced to a sin-
gle independent noise process. We thus introduce min-
imal complexity to the working model by considering