Pith. sign in

REVIEW 3 major objections 4 minor 26 references

Compressed sensing enabled high-bandwidth and large dynamic range magnetic sensing

T0 review · 3 major / 4 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Compressed sensing reconstructs the eight-line NV spectrum from ~100 random samples in a 650 MHz window, reaching ~0.5 MHz accuracy.

desk verdict A genuine and useful application of compressed sensing to NV ESR, with real experiments and matched simulations, but the central factor-of-2/3 comparison needs one more calibration before I'd trust the headline numbers. read the letter →

arxiv 2502.06070 v1 pith:LQMDNFZ5 submitted 2025-02-09 quant-ph

Galya Haim , Chris Mullarkey , John Howell , Nir Bar-Gill This is my paper · ORCID
classification quant-ph
keywords compressedsensingnitrogen-vacancycenterselectronspinresonancemagneticfieldsparsereconstructionLorentziandictionaryquantumhighdynamicrange
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that compressed sensing can replace dense raster scanning in nitrogen-vacancy (NV) diamond electron spin resonance (ESR) measurements used for magnetic field sensing. The idea is that an ESR spectrum is sparse—eight Lorentzian peaks in a broad frequency window—so instead of measuring every microwave frequency and averaging many times, one can sample a few random frequencies and reconstruct the peak locations with an overcomplete Lorentzian model. The paper reports that in a 650 MHz window at low signal-to-noise ratio (SNR≈3), about 100 random samples give a normalized peak-location error near 0.5 MHz, matching a full 650-point raster sweep and beating subsampled raster scans by a factor of roughly 2 to 3. If this holds, magnetic sensing in high-dynamic-range and low-SNR settings becomes faster and more sensitive, and the same random-sampling strategy could apply to other resonance-based sensors.

What carries the argument

The load-bearing object is a highly overcomplete, non-orthogonal dictionary of Lorentzian lineshapes: the signal is modeled as $A(\nu_j)=\sum_k L_{jk} a_k$, where each column of $L$ is a Lorentzian of unknown center frequency and width, and $a$ is a sparse vector of amplitudes. The reconstruction minimizes the total variation $\sum_i \|D_i a\|_1$ subject to $a \ge 0$ and consistency with the random frequency samples, then iterates until eight peaks have converged. The sparse structure of $a$ is what lets roughly 100 random measurements determine eight peak locations in a 650 MHz window; the overcomplete basis lets the peaks sit at arbitrary frequencies rather than only on the measurement grid.

What would settle it

Repeat the high-field experiment (100 G, 650 MHz, SNR≈3) and analyze subsampled raster data with the same eight-Lorentzian fitting code, width prior, and convergence criterion used in the CS loop; if the raster normalized error at 100 samples drops to about 0.5 MHz rather than staying several times larger, the claimed CS advantage is a product of unequal analysis.

Watch

Extended reading notes

Core claim

The paper's central claim is that, in the low-SNR, large-frequency-window regime where NV ESR is normally slow and error-prone, compressed sensing reconstructs the eight resonance lines more accurately than raster scanning does at the same data budget, and reaches the accuracy of a full raster sweep with only about 15% of the measurements. Experimentally, with a 100-G bias field, a 650 MHz window, SNR≈3, and 15 MHz Lorentzian widths, CS achieved a normalized error below 0.5 MHz with about 100 measurements, while the raster baseline reached that error only near the full 650-point sweep; the advantage persisted at up to a factor of five in error at 650 points. At lower field (50 G, 400 MHz window, SNR≈13) the two methods were comparable, though CS still kept the error below 1 MHz with only 10% of the data. Simulations with randomly varied field angles and Lorentzian widths reproduced the trend, and showed that simultaneous excitation of several frequencies improves CS efficiency further—for 150 measurements, four simultaneous frequencies cut the error by a factor of four relative to two.

Load-bearing premise

The comparison assumes the raster-scan baseline is a fair, optimally processed reference; if the subsampled raster data were analyzed with the same eight-Lorentzian model and fitting details used by the CS loop, the reported factor-of-2 to 3 improvement might shrink or disappear.

Editorial extensions

If this is right

  • High-field NV magnetometry with unknown, large bias fields can achieve the same resonance-position accuracy with roughly 15% as many frequency samples, directly improving measurement bandwidth and time resolution.
  • Since the method only assumes that the spectrum consists of a small number of Lorentzian peaks on a flat background, it transfers to other ESR and magnetic-resonance systems, not just diamond NV centers.
  • Simultaneous multi-frequency excitation, which reduces SNR in raster scans, becomes a useful lever under CS: the paper finds that more simultaneous frequencies lower the error at fixed total sample count.
  • The normalized-error gain maps to sensitivity, because sensitivity is proportional to the frequency error divided by the square root of total measurement time, so a lower error at fixed time means better magnetic-field sensitivity in the large-signal regime.
  • In high-SNR regimes the CS advantage shrinks to near parity, so the practical gain is specific to low-SNR, wide-window measurements.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, a natural stress test is to run the same subsampled raster data through the identical eight-Lorentzian fitting code used in the CS loop; if the raster error then matches CS, the reported gain is analysis-driven rather than sampling-driven.
  • The paper fixes the number of peaks at eight and uses known characteristic widths; allowing the number of peaks to be unknown would extend the method to more general field inhomogeneities, but would likely require a stronger sparsity prior.
  • The per-spectrum convergence rule suggests a direct extension to imaging: feeding each pixel's reconstructed peaks as a prior to neighboring pixels, as the paper mentions for video-like analysis, could reduce the required samples per pixel further.
  • An experimental sweep of Lorentzian width and SNR would map where the advantage vanishes; the paper's simulations indicate the benefit disappears for linewidths well above 15 MHz at SNR≈7.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes compressed sensing (CS) as a replacement for dense raster scanning in NV-center ESR magnetometry, targeting the high-dynamic-range, low-SNR regime where standard frequency sweeps are slow. The CS protocol randomly samples frequencies, reconstructs the spectrum in an overcomplete Lorentzian dictionary using total-variation minimization, and stops once eight peaks have converged. The authors define a normalized error δν/√P, where δν is the mean absolute peak-location error against a full-raster ground truth and P is the success probability. They support the method with simulations and with experiments at bias fields of roughly 100 G and 50 G (Figs. 2–5), claiming a factor of at least 2–3 improvement over raster subsampling in low-SNR conditions and equal accuracy with about 15% of the full-sweep data points.

Significance. If the central comparison is sound, the result is practically valuable: it offers a concrete route to faster, higher-bandwidth NV ESR in regimes where the resonance locations are unknown and SNR is low. The paper has genuine strengths: experimental validation on two different setups, independent ground truth from full raster scans, and simulations that reproduce the experimental trends. However, the load-bearing comparison with raster scanning is not fully specified: the raster baseline is subsampled in post processing, the peak-extraction procedure for that baseline is not stated, and the success-probability metric is defined in a way that is tightly coupled to the CS algorithm's own 8-Lorentzian model. These omissions make the claimed factor-of-2–3 improvement difficult to verify from the manuscript as written.

major comments (3)
  1. [Section II and Fig. 2] The raster-scan baseline is not sufficiently specified to support the central comparison. The text states that complete raster scans were 'subsampled in post processing' and plotted as triangles in Figs. 2(a) and 2(b), but it never states how peaks are extracted from the subsampled raster data, whether an 8-Lorentzian fit is used, how the success probability P is computed for raster, or how the mean absolute error δν is accumulated across runs. Since the claim of a factor-of-2–3 improvement is derived from the normalized error, the raster analysis must be described with the same level of detail as the CS algorithm, or the advantage could reflect a weaker raster peak-picking routine rather than an intrinsic property of CS.
  2. [Methods IV.A and IV.B] The success and error metrics are asymmetric between CS and raster. The paper defines a successful measurement as one that 'retrieves 8 Lorentzians' and, for CS, additionally requires the Lorentzian widths to match diamond-characteristic values within error. No equivalent criterion is stated for raster scans. The normalized error δν/√P therefore depends on a CS-specific success definition, and it is unclear whether δν is the mean error over successful runs only, over all runs, or over all identified peaks. The authors should give explicit formulas for P and δν as functions of the number of measurement points, state how failed runs are treated, and apply exactly the same peak-count and width criteria to both methods.
  3. [Section II, Fig. 5, and Discussion] The counting of 'measurements' is ambiguous and affects the headline data-reduction claim. All results are said to use 3 simultaneous frequencies, but the x-axis of Figs. 2 and 5 is not defined as either the number of projection rounds or the total number of individual frequency samples. An abstract claim of 'same accuracy with only 15% of the data points' for 100 measurements in a 650 MHz window is consistent only if a 'measurement' is a single frequency sample; if it is a projection containing 3 simultaneous frequencies, the sampled-data fraction is roughly 46%, not 15% (and analogously for the 40-measurement/10% statement in the low-field case). Please define the unit of the horizontal axis explicitly and recompute the percentage claims accordingly.
minor comments (4)
  1. [Methods IV.B] The reconstruction uses total-variation minimization on the coefficient vector a of an overcomplete Lorentzian dictionary rather than a standard L1 minimization of a itself. The relationship of this objective to the paper's 'compressed sensing' framing and to the sparsity assumption should be clarified, since the dictionary is coherent and the usual RIP-based guarantees do not directly apply.
  2. [Fig. 3 and Fig. 4] The simulation curves in Figs. 3 and 4 would benefit from error bars or confidence intervals and from a precise statement of what is averaged (number of samples, random field angles, and whether the plotted quantity is the normalized error or a different metric).
  3. [Throughout] The text contains several typographical and grammatical errors, including a missing period before 'Additionally' in the Introduction, 'T he' in the Methods heading, and 'fluerescence' for 'fluorescence'. The figure captions should also state the number of independent runs used for each experimental and simulated point.
  4. [Data availability] No data or code availability statement is included. Given that the central comparison depends on the raster post-processing procedure, releasing the analysis code and processed data would substantially strengthen reproducibility and would allow reviewers and readers to verify the baseline.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the CS-vs-raster comparison is anchored to independent ground truth, and the reported improvements are empirical benchmarking results rather than constructions from the algorithm's own assumptions.

full rationale

The central claim of the paper is an empirical performance comparison between compressed sensing and conventional raster scanning for NV ESR spectra, and it is anchored to an external ground truth. The paper states that 'The error is derived from "ground truth" data obtained with a complete raster scan,' and the normalized error metric δν/√P is defined from measured peak-location deviations, not from the CS model itself. The CS algorithm does use the prior that exactly eight Lorentzian peaks are present and defines successful convergence as retrieving eight peaks with characteristic widths (Methods IV.A: 'A measurement is considered successful when it retrieves 8 Lorentzians'; Algorithm 1: 'if 8 peaks have converged then update basis matrix L at the converged locations'). This is a physically motivated modeling assumption about the NV ensemble spectrum, not a definition that forces the reported error values: the recovered peak positions must still match the ground-truth raster positions. No equation in the paper reduces to its own input by construction, no fitted parameter is relabeled as a prediction, and no load-bearing claim is justified solely by a self-citation. The only self-citations ([15], [16]) are incidental application references and do not support the central derivation. The main caveat—that the raster-scan subsampling and peak-fitting procedure is not specified in Methods IV.A/B—is a reproducibility and benchmark-fairness concern, not a circularity, because the ground truth remains independent of the CS reconstruction. Therefore, the paper contains no significant circular step.

Assumptions & free parameters 8 free parameters · 4 assumptions · 0 invented entities

The central empirical claim is an application result. It rests more on protocol choices (8 peaks, Lorentzian dictionary, stopping criteria, simultaneous frequencies) than on fitted physical constants. The main free parameters are measurement-protocol choices rather than physical parameters. No new physical entities, forces, or dimensions are introduced.

free parameters (8)
  • Number of simultaneous frequencies = 3
    Chosen as a balance between performance and practical implementation; simulations (Fig. 5) show more simultaneous frequencies improve CS error, so this choice affects the reported gains.
  • Peak convergence tolerance = 2 MHz (Alg. 1) / 3 MHz (Methods IV.B)
    Stopping rule for the CS loop; the two stated values are inconsistent, which affects both results and reproducibility.
  • Number of consecutive successes before termination = 4
    User-set stopping condition in Alg. 1; affects the number of points sampled and the reported error.
  • Initial number of random projections = unspecified
    The algorithm starts with a 'user-defined (small) number' of projections; this value is not reported and impacts reconstruction quality.
  • Maximum number of measurements = unspecified
    Termination cap; data collection was extended to the maximum, but the maximum itself is not stated.
  • Number of Lorentzians = 8
    The algorithm assumes exactly 8 resonances, which is physically motivated for NV ensembles but is an input, not learned; success is defined as finding 8 peaks.
  • Lorentzian width prior = 'value characterizing the diamond'
    Used in loop termination; the exact width or acceptance window is not reported.
  • Dictionary grid spacing = unspecified
    The overcomplete set of possible resonance locations {nu_k} is not specified; grid density affects resolution and runtime.
assumptions (4)
  • domain assumption The ESR spectrum is a sparse non-negative combination of Lorentzians centered on a small number of frequencies.
    The reconstruction model in Methods IV.B assumes every signal component has the same Lorentzian lineshape family; this is physically motivated but limits generalization to non-Lorentzian or overlapping spectra.
  • standard math Random frequency sampling is sufficiently incoherent with the overcomplete Lorentzian dictionary for compressed-sensing recovery.
    The paper invokes standard CS guarantees via refs [4,12] but does not verify RIP or coherence for its specific dictionary.
  • domain assumption The simplified Zeeman-only NV Hamiltonian is sufficient to produce simulation data that represents real spectra.
    Simulations use only the Zeeman shift, while real data contain strain, hyperfine, and orientation-dependent effects; the paper acknowledges minor discrepancies.
  • domain assumption Ground truth raster scans are treated as error-free references.
    CS errors are measured against complete raster scans (Methods IV.A); if the reference raster itself has bias, the reported normalized errors inherit it.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Compressed sensing enabled high-bandwidth and large dynamic range magnetic sensing." pith.science (2026). https://pith.science/paper/LQMDNFZ5

@misc{pith2026250206070,
  author       = {Pith},
  title        = {Pith review of: Compressed sensing enabled high-bandwidth and large dynamic range magnetic sensing},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LQMDNFZ5}},
  note         = {Machine review of arXiv:2502.06070}
}
read the original abstract

Electron Spin Resonance (ESR) is a widely common method in the field of quantum sensing. Specifically with the Nitrogen-Vacancy (NV) center in diamond, used for sensing magnetic and electric fields, strain and temperature. However, ESR measurements are limited in temporal resolution, primarily due to the large number of data points required especially in high dynamic range regimes and the need for extensive averaging caused by low signal-to-noise ratio (SNR). This study introduces a novel application of compressed sensing (CS) for magnetic sensing using NV centers. By comparing CS with conventional raster scanning, we demonstrate the potential of CS to enhance sensing applications. Experimental results, supported by simulations, show an improvement of factor 3 in measurement accuracy in low SNR data, which also translates to achieving the same accuracy with only 15% of the data points. Moreover, the proposed approach is not confined to NV centers but can be extended to ESR measurements in other systems, broadening its applicability in quantum sensing.

Figures

Figures reproduced from arXiv: 2502.06070 by the authors.

Figure 1
Figure 1. FIG. 1. (a) Experimental setup schematic, NV’s are excited with green laser, manipulated with MW, and red [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Experimental data (blue) and simulation re [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. (a) Simulated CS (orange circles) and raster’s [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (2 more)
Figure 5
Figure 5. Figure 5: FIG. 5. Simulated average normalized error in raster [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 4
Figure 4. Figure 4: FIG. 4. Simulated averaged error as a function of [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 23 canonical work pages

  1. [1]

    R. G. Baraniuk, E. Candes, M. Elad, and Y. Ma, Applications of sparse representation and com- pressive sensing [scanning the issue], Proceedings of the IEEE 98, 906 (2010)

    work page 2010

  2. [2]

    Y. C. Eldar, Sampling Theory (Cambridge Univer- sity press, 2015)

    work page 2015

  3. [3]

    H. B. et al., Compressed Sensing and its Applica- tions (Birkhauser, 2019)

    work page 2019

  4. [4]

    E. J. Candes and M. B. Wakin, An introduction to compressive sampling, IEEE Signal Process. Mag. 25, 21 (2008)

    2008

  5. [5]

    M. F. Duarte and Y. C. Eldar, Structured compressed sensing: From theory to appli- cations, IEEE Signal Process. Mag. (2011), arXiv:1106.6224 [cs.IT]

    work page Pith review arXiv 2011

  6. [6]

    Craven, B

    D. Craven, B. McGinley, L. Kilmartin, M. Glavin, and E. Jones, Compressed sensing for bioelectric signals: a review, IEEE J Biomed Health Inform 19, 529 (2015)

    work page 2015

  7. [7]

    O. N. Jaspan, R. Fleysher, and M. L. Lipton, Com- pressed sensing MRI: a review of the clinical liter- ature, Br. J. Radiol. 88, 20150487 (2015)

    work page 2015

  8. [8]

    G. A. Howland, D. J. Lum, M. R. Ware, and J. C. Howell, Photon counting compressive depth map- ping, Opt. Express 21, 23822 (2013)

    work page 2013

Show all 26 references
  1. [9]

    Y. Sher, L. Cohen, D. Istrati, and H. S. Eisenberg, Low intensity LiDAR using depth aware compres- sive sensing and a photon number resolving detec- tor, in Imaging and Applied Optics 2019 (COSI, IS, MATH, pcAOP) (Optical Society of America,

    2019

  2. [10]

    Gross, Y.-K

    D. Gross, Y.-K. Liu, S. T. Flammia, S. Becker, and J. Eisert, Quantum state tomography via com- pressed sensing, Phys. Rev. Lett. 105, 150401 (2010)

    2010

  3. [11]

    Sadeghi, M

    M. Sadeghi, M. Babaie-Zadeh, and C. Jut- ten, Learning overcomplete dictionaries based on Atom-by-Atom updating, IEEE Trans. Signal Pro- cess. 62, 883 (2014)

    2014

  4. [12]

    E. J. Cand` es, Y. C. Eldar, D. Needell, and P. Ran- dall, Compressed sensing with coherent and redun- dant dictionaries, Appl. Comput. Harmon. Anal. 31, 59 (2011)

    2011

  5. [13]

    Real data was acquired on a wide field setup (see methods), where 100 consecutive CS measurements ran, and one full raster scan that was sub-sampled in post processing

    Simulations generated 57 data samples on which CS and raster scans ran. Real data was acquired on a wide field setup (see methods), where 100 consecutive CS measurements ran, and one full raster scan that was sub-sampled in post processing. ment of at least 2. The CS’s normali...

    2020

  6. [14]

    L. M. Pham, Magnetic Field Sensing with Nitrogen-Vacancy Color Centers in Diamond , Ph.D. thesis, Harvard (2013)

    2013

  7. [15]

    J. F. Barry, J. M. Schloss, E. Bauch, M. J. Turner, C. A. Hart, L. M. Pham, and R. L. Walsworth, Sensitivity optimization for nv-diamond magne- tometry, Rev. Mod. Phys. 92, 015004 (2020)

    2020

  8. [16]

    Meirzada, N

    I. Meirzada, N. Sukenik, G. Haim, S. Yochelis, L. T. Baczewski, Y. Paltiel, and N. Bar-Gill, Long- time-scale magnetization ordering induced by an adsorbed chiral monolayer on ferromagnets, ACS Nano 15, 5574 (2021)

    2021

  9. [17]

    Ninio, N

    Y. Ninio, N. Waiskopf, I. Meirzada, Y. Romach, G. Haim, S. Yochelis, U. Banin, and N. Bar-Gill, High-Sensitivity, High-Resolution detection of re- active oxygen species concentration using NV cen- ters, ACS Photonics 8, 1917 (2021)

    2021

  10. [18]

    E. S. Grant, L. T. Hall, L. C. L. Hollenberg, G. Mc- Coll, and D. A. Simpson, Nonmonotonic super- paramagnetic behavior of the ferritin iron core re- vealed via quantum spin relaxometry, ACS Nano 17, 372 (2023)

    2023

  11. [19]

    D. J. McCloskey, N. Dontschuk, A. Stacey, C. Pat- tinson, A. Nadarajah, L. T. Hall, L. C. L. Hollen- berg, S. Prawer, and D. A. Simpson, A diamond voltage imaging microscope, Nature Photonics 16, 730 (2022)

    2022

  12. [20]

    W. S. Huxter, M. L. Palm, M. L. Davis, P. Wel- ter, C.-H. Lambert, M. Trassin, and C. L. Degen, Scanning gradiometry with a single spin quantum magnetometer, Nature Communications 13, 3761 (2022)

    2022

  13. [21]

    R. Li, T. A. Vedelaar, A. Sigaeva, Y. Zhang, K. Wu, H. Wang, X. Wu, P. Olinga, M. K. Wlodarzyk-Biegun, and R. Schirhagl, Fluorescent nanodiamonds for tracking single polymer parti- cles in cells and tissues, Analytical Chemistry 95, 13046 (2023)

    2023

  14. [22]

    M. W. Doherty, N. B. Manson, P. Delaney, F. Jelezko, J. Wrachtrup, and L. C. Hollenberg, The nitrogen-vacancy colour centre in diamond, Physics Reports 528, 1 (2013), the nitrogen- vacancy colour centre in diamond

    2013

  15. [23]

    B. J. Maertz, A. P. Wijnheijmer, G. D. Fuchs, M. E. Nowakowski, and D. D. Awschalom, Vector magnetic field mi- croscopy using nitrogen vacancy centers in diamond, Applied Physics Letters 96, 092504 (2010), https://pubs.aip.org/aip/apl/article- pdf/doi/10.1063/1.3337096/1442641...

    2010 doi

  16. [24]

    Clevenson, M

    H. Clevenson, M. E. Trusheim, C. Teale, T. Schr¨ oder, D. Braje, and D. Englund, Broad- band magnetometry and temperature sensing with a light-trapping diamond waveguide, Nature Physics 11, 393 (2015)

    2015

  17. [25]

    A. H. Marblestone*, B. M. Zamft*, Y. G. Maguire, M. G. Shapiro, T. R. Cybulski, J. I. Glaser, D. Amodei, P. B. Stranges, R. Kalhor, D. A. Dalrymple, D. Seo, E. Alon, M. M. Maharbiz, J. M. Carmena, J. M. Rabaey, E. S. Boyden**, G. M. Church**, and K. P. Kording**, Physical prin...

    2013

  18. [26]

    M. S. Asif and J. Romberg, Sparse recovery of streaming signals using l1-homotopy (2013), arXiv:1306.3331

    2013 arXiv

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.