{"id":"8310db0b-f950-4b06-8f50-24f76e4de980","arxiv_id":"2506.15805","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Quantum Invariant Filtering reverses the usual time-domain design of quantum control, mapping a chosen frequency filter directly into continuous drive fields, and is validated on a nitrogen-vacancy qubit.","lead":"Researchers built a quantum control scheme that starts from a desired frequency filter and computes the microwave waveform needed to realize it on a qubit. They demonstrate it on a single nitrogen-vacancy center, achieving tunable, dual-band, phase-sensitive filters that outperform standard pulse sequences.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The experimental linearization β(t)=−π/2+H(t) replaces the exact arcsin mapping; the realized filter kernel is sin H(t), so strong and multi-band filters develop unquantified spectral replicas that the 'arbitrary FIR' claim relies on.","rationale":"I read the central claim as the systematic mapping from a target frequency-domain filter to a time-dependent control Hamiltonian, with Eq. (S39) as the link that makes the measured σz response proportional to the convolution of the signal with the designed impulse response H(t). The load-bearing step is therefore the choice of β(t). The paper's exact construction is β(t)=−π/2+arcsin H(t), but the experimental implementation uses β(t)=−π/2+H(t), which changes the effective impulse response from H to sin H. This is an admitted simplification (Sec. S2 B) but its error is never quantified, and the central claim of arbitrary FIR filtering depends on it. The concern is concrete: for the dual-band kernel, sin H generates intermodulation products inside the scanned frequency range, not just a global amplitude rescaling. If the realized H amplitudes are small in the demonstrated data, the narrow-band demonstrations remain valid, but the 'arbitrary FIR' claim is not established without either using the exact arcsin mapping or bounding the error. I agree with the reader's weakest-assumption identification. The experimental demonstrations themselves—passband position, tunability, phase sensitivity, dual-band response, and amplitude robustness—are real evidence and are not called into question by this concern. The coherence comparison overclaim and missing error bars flagged by the reader are secondary; the linearization directly affects the core frequency-to-time mapping. The verdict should remain CONDITIONAL pending a quantitative analysis of the linearization error or a switch to the exact β relation.","tokens_in":21608,"tokens_out":8738,"duration_ms":93988,"concrete_test":"For the dual-band filter of Fig. 3e–f (tf=4 μs, f1=1.5 MHz, f2=2.5 MHz), compute the Fourier transform of the actually implemented kernel sin H(t) with H=(θ/2)[cos(2πf1t)+cos(2πf2t)] and compare its magnitude to the designed |H(f)|. If intermodulation peaks at 0.5 and 3.5 MHz exceed the measured noise floor, the linear choice β=−π/2+H is falsified for that demonstration; in addition, report max|H(t)| for every experimental filter and re-fit the measured spectra using the exact kernel sin H to see whether the passband widths and sidelobe levels shift by more than the error bars.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mapping is Eq. (S39)/(5): with the exact auxiliary field β(t)=−π/2+arcsin H(t), the second-order response is A(tf)=∫fin(s)H(s), so the designed frequency response H(ω) is enacted. The experiment instead uses β(t)=−π/2+H(t) (main text and Sec. S2 B). This is not a harmless reparameterization: the realized kernel is cosβ=sin H(t), so the actual response is |F[sin H]|², not |H(ω)|². The first correction is −H³/6. For a single bandpass, this creates a replica at 3f0. For the dual-band kernel H=(θ/2)[cos(2πf1t)+cos(2πf2t)] used in Fig. 3e–f, expansion of sin H produces intermodulation products at |2f1−f2| and |2f2−f1| (0.5 and 3.5 MHz for f1=1.5, f2=2.5 MHz), both inside the measured band. The paper provides no bound on max|H(t)| and no quantification of these replicas or of the validity range of the small-H linearization. Since the abstract claims arbitrary FIR profiles are mapped into control Hamiltonians, this unquantified approximation is load-bearing: the demonstrated passband positions may be correct while the general filter-shape claim is not yet supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces Quantum Invariant Filtering (QIF), a method that reverse-engineers continuous qubit control fields from a target frequency-domain filter profile. Within the dynamical-invariant framework and a second-order Dyson expansion, the authors derive that the final-time expectation value <σz(tf)> responds as a convolution of the input signal with an impulse response H(t) (Eq. S39). An exact construction sets the auxiliary field to β(t) = −π/2 + arcsin H(t); the experiment instead uses the linearized form β(t) = −π/2 + H(t), yielding the control field ε(t) = ∂t H(tf − t). The method is demonstrated on a single NV center in diamond, with experiments showing tunable bandpass filtering, phase-sensitive detection, amplitude detection, dual-band filtering, and robustness/coherence comparisons against CPMG sequences.","tokens_in":21912,"tokens_out":5560,"duration_ms":54486,"significance":"The theoretical idea is elegant and, if properly supported, would constitute a broadly applicable framework for frequency-domain-designed quantum control and sensing. Strengths include a clean analytic mapping from a specified filter profile to a control Hamiltonian, the numerical cross-check in Fig. S1, an explicit construction recipe, and a diverse set of experimental demonstrations. The novelty is clear and the level of ambition is appropriate. However, the central 'arbitrary FIR' claim currently rests on an unquantified linearization of the auxiliary field, and the experimental evidence lacks error bars, uses normalization and a fitted phase offset, and the coherence comparison may be confounded by the use of a different NV center. These issues are fixable, but they are load-bearing for the paper's headline claims.","major_comments":[{"comment":"The substitution β(t) = −π/2 + H(t) replaces the exact relation β(t) = −π/2 + arcsin H(t) from Eq. (S35). The realized second-order kernel is then cos β = sin H(t) ≈ H(t) − H(t)^3/6, so the implemented filter is |F[sin H]|^2 rather than |H(ω)|^2. For the dual-band kernel of Eq. (S56) with f1 = 1.5 MHz and f2 = 2.5 MHz, the leading third-order correction generates intermodulation products at |2f1 − f2| = 0.5 MHz and |2f2 − f1| = 3.5 MHz, both inside the measured band. The manuscript neither bounds max|H(t)| nor quantifies these replicas. Because the abstract's claim of arbitrary FIR mapping and the dual-band demonstration depend on the realized kernel, the authors should either implement the exact arcsin construction or provide a rigorous small-H error bound and demonstrate that the replicas are below the experimental sensitivity.","section":"Main text after Eq. (5); Supplementary S2 B, Eq. (S48)"},{"comment":"The supplement states that the T1 and T2 measurements were performed on a different NV center that was not polarized. If the QIF and CPMG decay curves in Fig. 4b were acquired on different NV centers, the reported factor-of-approximately-40 improvement in decay time cannot be attributed solely to the protocol. The authors should either measure all protocols on the same center or report per-center noise characterization (e.g., T2* or an independently measured noise spectral density) to make the coherence comparison valid.","section":"Fig. 4b; Supplementary S2 C, 'T1 and T2 Figure'"},{"comment":"Experimental points are presented without error bars, and the theoretical curves in Figs. 2b and 3f are explicitly normalized to the maximum experimental contrast. The normalization constant is a free parameter, so the agreement currently demonstrates qualitative line-shape correspondence rather than quantitative transfer-function fidelity. The authors should report statistical uncertainties (one million averages per point should allow tight error bars) and either show absolute, unnormalized data or state explicitly that only relative response is claimed.","section":"Figs. 2b, 2e, 3b, 3f, 4a"},{"comment":"The central response formula is second-order perturbative in the signal amplitude δ, and the numerical check in Fig. S1 shows that the perturbation deviates from exact dynamics at the response peaks even for δ = 0.5. The manuscript does not report the δ values used in the experimental scans of Figs. 2 and 3, nor does it give the domain of validity of Eq. (S39) as a function of δ·tf. Since Fig. 3d explicitly addresses signal amplitude extraction, the authors should state the experimental δ and verify that the quoted 'high fidelity' and amplitude calibration lie within the validated perturbative regime.","section":"Eqs. (S37)–(S39) and Fig. S1"}],"minor_comments":[{"comment":"There is a typo: 'of the the QIF protocol' should read 'of the QIF protocol'.","section":"Fig. 1 caption"},{"comment":"The symbol H is used both for the Hamiltonian (Eq. (1)) and for the impulse response (Eq. (5) and later). This is confusing; consider denoting the impulse response with a lowercase symbol such as h(t).","section":"Throughout main text"},{"comment":"The phrase 'arbitrary finite-impulse responses' is broader than the class actually implemented: the construction requires a smooth continuous H(t) with H(0)=H(tf)=0, |H(t)|<1, and time symmetry about tf/2. The authors should qualify this claim.","section":"Abstract and Discussion"},{"comment":"The phase offset Δφ_signal ≈ π/8 is introduced as a fitted adjustment and then corroborated by electronic measurement; the fit uncertainty and the uncertainty of the electronic corroboration should be reported.","section":"Fig. 3c"}],"recommendation":"major_revision","confidential_remarks":"I see no grounds for questioning novelty or author conduct. The theoretical core is sound under the stated assumptions, and the main risks are the unquantified experimental linearization and the cross-center coherence comparison; both are fixable within the scope of a revision. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What to know: QIF is a real advance. The idea of starting from an FIR filter specification and reverse-engineering a continuous control Hamiltonian via dynamical invariants is not in the cited literature, and the experiments back it up: the NV response tracks the designed passbands, the dual-band filter works, and the phase sensitivity is a nice lock-in-like feature. The theoretical derivation in Eq. (S39) is clean under the stated assumptions, and the numerical comparison in Fig. S1 gives me confidence the perturbative picture is right. I agree with the reader that circularity is low: the filter is the input, the measured response is the benchmark, and the only fitted adjustments are a normalization and a phase offset that is independently corroborated.\n\nThe soft spots are real but not fatal. The most important is the linearization. The theory uses beta = -pi/2 + arcsin(H), which makes cos(beta) = H exactly. The experiment uses beta = -pi/2 + H, so the realized kernel is sin(H), not H. For small H that's fine, but for the dual-band filter and any stronger filter, sin(H) generates third-order replicas and intermodulation products that are never quantified. The paper should either implement the exact relation (it's not that hard, just a sqrt in the denominator) or state a bound on max|H| and show the replicas are below the noise floor. As written, the claim of \"arbitrary FIR profiles\" is stronger than what is actually demonstrated.\n\nThe abstract also overstates the coherence gain: Fig. 4b shows QIF decay around 2.7 ms versus 0.07–0.14 ms for CPMG, which is roughly a factor of 20–40, not \"two orders of magnitude.\" And the coherence comparison uses different NV centers, with no error bars on any of the main figures. These are easy fixes and don't undermine the core result.\n\nWho gets value: anyone working on quantum control, NV-center sensing, or noise spectroscopy. The method gives a practical, analytic route to custom filter shapes that pulsed sequences can't produce. It deserves a serious referee. My recommendation: send it to review, but require the authors to quantify the linearization error, correct the coherence claim, and add error bars or acknowledge their absence. The central idea holds up; the presentation just needs to match the evidence.","headline":"A genuinely new way to go from frequency-domain filter specs to continuous control fields, with a solid experimental demo and a real but fixable linearization gap that the paper never quantifies.","tokens_in":22441,"tokens_out":2748,"would_cite":true,"duration_ms":28644,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Quantum Invariant Filtering turns a desired spectral filter into the continuous driving field that realizes it.","keywords":["quantum invariant filtering","dynamical invariants","quantum control","dynamical decoupling","frequency-domain filter design","nitrogen-vacancy center","quantum sensing","dual-bandpass filter"],"falsifier":"Scale a designed impulse response $H(t)$ to a maximum amplitude around $0.5$, measure $\\langle\\sigma_z(t_f)\\rangle$ versus probe frequency, and compare with two predictions: the linearized kernel gives $|H(\\omega)|^2$, while the exact $\\arcsin$ construction shows a third-order replica at three times the center frequency whose height grows as $\\max|H|^3/6$. A direct time-domain deconvolution of the realized kernel would also settle whether the linearization is acceptable.","tokens_in":21376,"feed_emoji":"⚛️","tokens_out":5258,"duration_ms":48953,"temperature":0.7,"pith_summary":"The paper proposes and demonstrates a method, Quantum Invariant Filtering, that reverses the usual order of quantum control design: instead of choosing a time-domain pulse sequence and then measuring its frequency response, the user starts with a desired spectral filter and derives the continuous Hamiltonian modulation that implements it. The claim is that any finite-impulse-response filter profile, whether single-band, multi-band, or phase-sensitive, can be mapped into an experimentally feasible drive on a qubit through the dynamical-invariant framework. On a single nitrogen-vacancy center in diamond, the method reproduces prescribed passbands, suppresses low-frequency noise, keeps coherence for milliseconds, and tolerates drive-amplitude errors far better than Carr-Purcell-Meiboom-Gill sequences. If correct, this gives quantum control and sensing a direct design path from the frequency domain to the waveform.","feed_headline":"Design the spectrum, get the quantum drive field","feed_subtitle":"NV-center experiment realizes multi-band and phase-sensitive filters that pulse sequences cannot.","key_machinery":"The central object is the dynamical invariant $I(t)$ of a driven qubit, parametrized by two auxiliary angles $\\alpha(t)$ and $\\beta(t)$. Under the constraint $\\alpha(t)=\\pi$ and $\\dot\\alpha(t)=0$, the relative Lewis-Riesenfeld phase accumulated between invariant eigenstates is exactly $\\beta(t)$, so the second-order response to a $z$-axis perturbation becomes a convolution with kernel $\\cos\\beta(t)=H(t)$. To make the target filter the response, one sets $\\beta(t) = -\\pi/2 + \\arcsin H(t)$, so that the Hamiltonian control $\\varepsilon(t)$ is obtained from the time derivative of the reversed impulse response. This identity, filter specification in frequency to inverse Fourier transform to auxiliary angle to drive field, is what carries the argument.","core_discovery":"The central discovery is that the final-time readout $\\langle\\sigma_z(t_f)\\rangle$ of a qubit driven by a reverse-engineered invariant control is, to second order in the signal amplitude, the squared convolution of the input signal with a kernel $H(t)$: $\\langle\\sigma_z(t_f)\\rangle \\simeq \\langle\\sigma_z(t_f)\\rangle_0 - \\frac{\\delta^2}{2}\\left(\\int f_{\\mathrm{in}}(s)H(t_f-s)\\,ds\\right)^2$. Since $H(t)$ is generated from the inverse Fourier transform of the target filter $H(\\omega)$, the measured contrast for a monochromatic probe directly reproduces $|H(\\omega)|^2$. The construction fixes the auxiliary invariant angle to $\\beta(t) = -\\pi/2 + \\arcsin H(t)$, and in the implemented linearized version $\\beta(t) = -\\pi/2 + H(t)$, the control field is $\\varepsilon(t) = \\partial_t H(t_f - t)$. Because the protocol starts from the spectral specification rather than from a pulse sequence, it can realize passband shapes and dual-band filters that pulsed dynamical decoupling cannot, and it does so with a smooth, continuous drive.","pith_inferences":["Editorial inference: the exact $\\arcsin$ mapping would produce kernel $\\sin(\\arcsin H(t))=H(t)$, but the linearized implementation uses $\\sin H(t)\\approx H(t)-H(t)^3/6$; for stronger filters the cubic term will create spectral replicas that the paper does not quantify.","Editorial inference: because the response is quadratic in the signal amplitude, very weak signals produce a quadratically suppressed readout, so practical sensitivity will require either larger $\\delta$ or additional amplification; a dedicated calibration curve would show where this becomes limiting.","Editorial inference: by symmetry, the same invariant construction should work for perturbations along other axes, turning QIF into a vector-field filter, though the paper only demonstrates $z$-axis noise and signal coupling."],"forward_implications":["A user can design quantum control for a qubit by specifying the spectral response directly, including multi-band passbands and phase offsets.","Noise concentrated outside the passband is suppressed, extending coherence much longer than Carr-Purcell-Meiboom-Gill sequences; the paper reports about two orders of magnitude.","The measurement of the final population yields both amplitude and phase of a signal near the filter frequency, enabling lock-in-style quantum sensing.","The same prescription applies to other qubit platforms such as superconducting qubits, trapped ions, and nuclear magnetic resonance.","Because the control is continuous, drive-amplitude miscalibration of $\\pm 50\\%$ still leaves $\\langle Z\\rangle$ above roughly $0.8$, unlike $\\pi$-pulse sequences."],"supporting_citations":[{"why":"Supplies the dynamical-invariant formalism whose eigenstate expansion gives the closed-form propagator used to derive the response.","marker":"[36]"},{"why":"Establishes reverse-engineering of controls from invariants in shortcuts to adiabaticity, which QIF adapts to filter design.","marker":"[37]"},{"why":"Gives the invariant-basis treatment of open-system response used for the second-order perturbation computation.","marker":"[12]"},{"why":"Provides the finite-impulse-response filter design methodology whose impulse response is mapped into the auxiliary angle.","marker":"[41]"},{"why":"Shows single-ion quantum lock-in phase sensitivity, the basis for the QIF phase-detection experiments.","marker":"[42]"},{"why":"Defines the Carr-Purcell-Meiboom-Gill sequence used as the baseline for the robustness and coherence comparisons.","marker":"[26]"},{"why":"Introduces dynamical decoupling, the pulsed paradigm against which QIF's continuous, frequency-first control is contrasted.","marker":"[24]"},{"why":"Characterizes the 1/f noise spectrum typical of solid-state qubit environments that the QIF is designed to suppress.","marker":"[38]"}],"fun_headline_variants":["Spectrum-first quantum control: filters pulses cannot make","Design the filter, get the drive: invariant-based qubit control","NV demo: control fields from spectral targets, robust to errors","Multi-band quantum filters from smooth drives, demonstrated"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The experimental filter uses the linearized auxiliary angle $\\beta(t) = -\\pi/2 + H(t)$ instead of the exact $\\beta(t) = -\\pi/2 + \\arcsin H(t)$, so the realized spectral kernel is $\\sin H(t)$ rather than $H(t)$, and the paper does not quantify how large $H(t)$ may be before the cubic error distorts the passband.","fun_headline_variants_meta":{"raw":{"variants":["Spectrum-first quantum control: filters pulses cannot make","Design the filter, get the drive: invariant-based qubit control","NV demo: control fields from spectral targets, robust to errors","Multi-band quantum filters from smooth drives, demonstrated"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000272,"raw_usage":{"total_tokens":1632,"prompt_tokens":944,"completion_tokens":688,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":621}},"tokens_in":560,"tokens_out":688,"duration_ms":7373,"temperature":1.0,"reasoning_tokens":621,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:31:28.858468+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Scale a designed impulse response $H(t)$ to a maximum amplitude around $0.5$, measure $\\langle\\sigma_z(t_f)\\rangle$ versus probe frequency, and compare with two predictions: the linearized kernel gives $|H(\\omega)|^2$, while the exact $\\arcsin$ construction shows a third-order replica at three times the center frequency whose height grows as $\\max|H|^3/6$. A direct time-domain deconvolution of the realized kernel would also settle whether the linearization is acceptable.","supporting_citations":[{"cited_title":"Lewis and W","cited_arxiv_id":null,"evidence_quote":"Supplies the dynamical-invariant formalism whose eigenstate expansion gives the closed-form propagator used to derive the response."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the invariant-basis treatment of open-system response used for the second-order perturbation computation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the finite-impulse-response filter design methodology whose impulse response is mapped into the auxiliary angle."},{"cited_title":"Viola and S","cited_arxiv_id":null,"evidence_quote":"Introduces dynamical decoupling, the pulsed paradigm against which QIF's continuous, frequency-first control is contrasted."},{"cited_title":"Paladino, Y","cited_arxiv_id":null,"evidence_quote":"Characterizes the 1/f noise spectrum typical of solid-state qubit environments that the QIF is designed to suppress."}],"review_version":2}