{"id":"f10cd8cd-a558-4510-950b-6cafaff880d3","arxiv_id":"2507.07347","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A weighted least-squares filter with de-weighted padding, solved by conjugate gradient, reduces boundary artifacts in bandlimited spectroscopy filtering compared to symmetric padding.","lead":"Researchers describe a signal-filtering method, called a rubber band filter, that reduces edge distortions when smoothing spectroscopy data. The trick is to let the signal outside the measured window move freely during the fit, so the filtered trace matches the true signal more closely than standard symmetric padding.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported edge-artifact reduction is measured against a full-signal oracle, and the padding weight is tuned to that same oracle; without a blind test, the central claim of optimal artifact-free filtering is not established.","rationale":"The reader's conditional verdict and the weakest-assumption analysis converge on the same load-bearing issue: the validation protocol uses a full-signal ground truth that is unavailable in the intended use case, while the free weight parameter is tuned against that same ground truth. This is not a mathematical inconsistency in Eqs. (1)-(3); the weighted least-squares problem is convex and the CG implementation is plausible. The concern is empirical and evaluative: the paper's quantitative support for 'optimal padding without edge artifacts' depends on oracle access to the signal that the filter is supposed to reconstruct. In the absence of code or data, a blind synthetic test with a fixed, pre-specified weight is the minimal check that would separate a genuine filtering improvement from a tuned fit to the reference. I agree with the reader's assessment and do not see a reason to move the verdict. The manuscript itself flags the weight as 'somewhat arbitrarily' chosen and notes that optimal parameters vary with dataset, which reinforces the concern. The lack of comparison to established bandlimited-extrapolation methods is a separate limitation, but it is secondary to the oracle-tuning issue. The appropriate verdict remains conditional: the method is plausible and internally consistent, but the central claim needs stronger, non-oracle evidence before full acceptance.","tokens_in":5616,"tokens_out":4851,"duration_ms":63106,"concrete_test":"Construct a synthetic bandlimited signal with a known spectrum and a steep edge, truncate it at the steep-slope point, and apply both symmetric padding and the rubber-band filter using the a priori weight w_pad = 0.03 * w_Tukey(alpha = 0.2) exactly as in the main THz example, without any oracle tuning. Compute the RMS boundary error against the known full-signal filtered result over 100 random phase realizations; then repeat with w_pad = 0.06, the value Fig. 5a identifies as optimal. If the median boundary error is not substantially lower than symmetric padding in both cases, or if the ranking flips between the two weight values, the claimed advantage is an artifact of oracle-tuned regularization rather than a robust property of the filter.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central evaluation step in Figs. 2, 3, and 4 defines the true filtered signal as the filtered version of the complete, untruncated measurement, while the rubber-band filter operates only on a symmetric-padded truncated window. Equations (2) and (3) show that the method minimizes a weighted least-squares cost: the result is a regularized bandlimited extrapolation determined by the chosen weight w(t), not by the unavailable true boundary values. The examples therefore demonstrate that, for a weight chosen with knowledge of the oracle, the rubber-band extrapolation matches the full-signal filtered result better than symmetric padding does. In practice the true filtered signal is not available; a real filtering application must choose w a priori. The paper uses w_pad = 0.03 with a Tukey window in the main examples, but Fig. 5a reports that the optimal weight for the THz dataset is 0.06. Because this tuning target is exactly the oracle reference, the reported roughly tenfold error reduction may reflect oracle-assisted weight selection rather than a reproducible property of the algorithm. Thus the manuscript's headline claims of 'optimal padding' and 'artifact-free' filtering are not yet supported by evidence that would be available when the filter is actually used.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a 'rubber band filter' for bandlimited filtering of truncated time-domain signals. The method symmetrically pads the signal, assigns weight 1 to the valid region and a small weight to the padded region, and solves a weighted least-squares problem with a bandlimited constraint via conjugate gradient, using FFTs for the required convolutions. The authors demonstrate the approach on terahertz time-domain spectroscopy, FTIR interferogram, and dual-comb spectroscopy data, reporting boundary errors roughly three to ten times lower than conventional symmetric padding. The central claim is that the method achieves 'optimal padding' and artifact-free bandlimited filtering with only a few extra FFTs over conventional approaches.","tokens_in":5895,"tokens_out":3491,"duration_ms":37089,"significance":"If the performance claims hold, the method would provide a practical and efficient alternative to symmetric padding for spectroscopic filtering, particularly where preserving sharp edge features is important. The derivation of the weighted least-squares formulation and the CG/FFT solution strategy is clear and technically sound, and the paper usefully identifies a real problem: symmetric padding introduces artifacts when the boundary slope is steep. The manuscript also makes a welcome effort to demonstrate the method on three distinct spectroscopic modalities. However, the evaluation is limited: errors are measured against an oracle reference (the full unpadded signal), the weight parameter is tuned per dataset, no error bars or statistical replicates are provided, and no comparison is made to existing bandlimited extrapolation or alternative padding methods. These limitations leave the 'optimal padding' claim stronger than the current evidence supports.","major_comments":[{"comment":"The main examples in Figs. 2–4 use a padded-region weight of 0.03 times the Tukey window (α=0.2), but Fig. 5a reports that the optimal constant for the THz dataset is 0.06. The paper does not explain this discrepancy, leaving it unclear whether the reported 'approximately ten times lower' error is obtained at the chosen weight or at the tuned optimum, and whether the main examples are suboptimal.","section":"§2, Fig. 5a"},{"comment":"The 'true filtered signal' in Figs. 2–4 is defined as the full, unpadded measurement filtered with the same passband, an oracle that is not available in practical filtering. The text states that the padding weight is 'initially chosen somewhat arbitrarily' and that it acts as a regularization constant, and Fig. 5a tunes that weight against the same oracle. Without a blind test or a principled a priori selection rule, the reported error reductions may reflect oracle-assisted parameter selection rather than a reproducible property of the algorithm, so the claim of 'optimal padding' is not yet supported by evidence available at the point of use.","section":"§2, Figs. 2–4 and Eq. (2)"},{"comment":"No comparison is made to established bandlimited extrapolation or alternative padding schemes, such as Papoulis–Gerchberg iteration, antireflective padding, or linear-prediction-based extension. The claim of 'overcoming longstanding challenges' and the assertion of superiority rest on a single baseline (symmetric padding), which is insufficient to establish the method's advantage relative to existing techniques that also aim to reduce boundary artifacts.","section":"§2; §3"},{"comment":"The convergence behavior is reported only as the residual against the direct solution A−1b of the CG iterate, with no information on the number of iterations or wall-clock time for the actual datasets, nor on how convergence depends on the weight parameter. Consequently, the claim that only 'a handful of iterations' and 'a few extra Fourier transforms' are needed cannot be quantitatively assessed from the presented results.","section":"§3, Fig. 5b"}],"minor_comments":[{"comment":"The phrase 'effectively removing the redundancy information' is ungrammatical; it should be 'effectively removing redundant information' or 'removing the redundant signal components.'","section":"Introduction"},{"comment":"The text says Fig. 5 highlights three aspects, including a 'padded region comparison between our method and symmetric padding method,' but the figure as described contains only panels (a) and (b) for weight dependence and convergence. Please either add the third panel or correct the caption and referencing text.","section":"Fig. 5"},{"comment":"The definition wpad(n) = 0.03 wTukey(n, α=0.2) is ambiguous about whether n denotes the number of padded samples or an index over the padded region; please clarify the notation.","section":"§2, weight definition"},{"comment":"The reference list omits seminal work on bandlimited extrapolation, such as Papoulis (IEEE Trans. Circuits Syst., 1975) and Gerchberg (Optica Acta, 1974), which are directly relevant to the method's relation to prior art; adding these would help readers position the contribution.","section":"References"},{"comment":"The data availability statement says data 'will be published prior to publication'; for a reproducibility-oriented methods letter, consider making the datasets and example code available at the time of submission.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a technically sound weighted least-squares bandlimited extrapolation implemented with CG and FFTs, and the three spectroscopic examples are appropriate. The main weaknesses are in the evaluation and framing: the 'optimal' claim is undercut by the oracle-based tuning and the 0.03 vs. 0.06 discrepancy, and the absence of comparison to existing bandlimited extrapolation methods makes the novelty claim hard to judge. The required revisions—clarifying the weight choice, providing a blind or sensitivity analysis, and adding a comparison to at least one established method—should be feasible within the scope of a letter, possibly with a supplementary document."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a useful and clearly-written thing to know about if you filter spectroscopic time traces. The rubber band filter is basically a weighted least-squares projection onto a bandlimited subspace, with the weight deliberately made tiny in the symmetrically padded region. That is not a brand-new concept—it sits in the bandlimited extrapolation family (Papoulis-Gerchberg and relatives)—but the authors apply it cleanly to a practical spectroscopy problem, solve the normal equations with CG/FFTs, and demonstrate on THz, FTIR, and dual-comb phase data that it reduces edge error by roughly 3–10x versus symmetric padding. The derivation is correct, the convergence plot (Fig. 5b) is convincing, and the computational overhead is modest.\n\nWhere I part ways: the headline claims \"optimal padding\" and \"artifact-free\" are stronger than the evidence. The ground truth in Figs. 2–4 is the filtered version of the complete, untruncated measurement. That is a sensible oracle for a filter-evaluation paper, but it is not available in a real filtering job, and the paper tunes the weight constant (0.03 in the main examples; 0.06 \"optimal\" in Fig. 5a) against that same oracle. So the reported tenfold improvement could be partially an artifact of oracle-assisted weight selection. The stress-test note is right. What is missing is a blind test: pick the weight by a fixed rule or cross-validation, then compare against symmetric padding. Also missing: any comparison to established bandlimited extrapolation baselines, error bars or multiple realizations, and the promised code/data in the arXiv version. That last point is a real reproducibility gap.\n\nMinor: \"artifact-free\" is contradicted by their own error plots, which still show elevated error at the edges. \"Optimal padding\" is also a misnomer—the padding itself is still symmetric; only the weight is optimized.\n\nVerdict: the central method holds up as a useful engineering contribution, but the paper needs a softer title, a blind evaluation, and a baseline comparison before the strong quantitative claim is credible. It is still worth a serious referee—this kind of preprocessing is widely used in spectroscopy, and the method is clearly explained and likely to be adopted once its robustness is demonstrated. I would send it to review, but I would push for major revision.","headline":"A clean, useful least-squares bandlimited filtering method for spectroscopy, but the 'optimal' and 'artifact-free' claims outrun the oracle-assisted validation.","tokens_in":6351,"tokens_out":2226,"would_cite":true,"duration_ms":26020,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that bandlimited filtering can be made free of boundary artifacts by recasting it as a weighted least-squares fit over an optimally de-weighted padded region, and demonstrates boundary errors roughly ten times lower than…","keywords":["rubber band filter","bandlimited filtering","edge artifacts","symmetric padding","weighted least squares","conjugate gradient","spectroscopy","Fourier transform"],"falsifier":"If, on a synthetic signal with a known bandlimit, no single small weight (or simple weight rule) keeps the boundary error consistently below symmetric padding across noise and cut-off choices, then the claim that the method recovers the bandlimited signal without edge artifacts fails.","tokens_in":5427,"feed_emoji":"📉","tokens_out":6434,"duration_ms":66519,"temperature":0.7,"pith_summary":"The paper introduces a \"rubber band filter\" for bandpass filtering of spectroscopic signals, aimed at removing the edge artifacts that conventional symmetric padding leaves at signal boundaries. It frames filtering as a weighted least-squares fit: the measured trace is symmetrically padded, the padded region is strongly de-weighted, and the best bandlimited reconstruction is found by solving a convolution system in the frequency domain. The authors show on terahertz time-domain, FTIR, and dual-comb spectroscopy data that the method preserves sharp edges with substantially lower error, around ten times lower in the THz example and about three times lower for dual-comb phase, while needing only a handful of extra Fourier transforms. If the method holds up, it gives spectroscopists a drop-in filtering step that does not distort the very edges where important signal features live.","feed_headline":"Rubber band filter cuts edge artifacts tenfold in spectroscopy","feed_subtitle":"Weighted padding plus a few Fourier transforms preserves sharp edges in THz, FTIR, and dual-comb measurements.","key_machinery":"The central object is the weighted least-squares cost $J = \\sum_n w(t_n) \\left| s_f(t_n) - s(t_n) \\right|^2$, with the filtered signal restricted to passband frequencies $\\Omega$. Its stationary condition is the frequency-domain convolution equation $W * F = W * S$ restricted to $\\Omega$, solved by conjugate gradient where each matrix-vector product $A F = W * F$ is evaluated with FFTs; the result is the optimal set of passband Fourier amplitudes $F(\\omega_k)$. The work it does is letting the padded region absorb the boundary discontinuity smoothly while keeping the reconstruction strictly bandlimited, at $O(N \\log N)$ cost.","core_discovery":"On the paper's own terms, the discovery is that edge-free bandlimited filtering can be posed as a convex least-squares problem over the passband Fourier amplitudes, with the padded region included in the fit but heavily down-weighted. The stationary condition reduces to the convolution equation $W * F = W * S$, where $S$ is the Fourier transform of the padded signal and $W$ is the transform of the weight function; solving only for frequencies in the passband $\\Omega$ and setting the rest to zero yields a bandlimited trace that is allowed to relax in the padding instead of being forced into the discontinuity that symmetric padding creates. The paper validates this on three experimental records by comparing against a true filtered signal obtained from the full unpadded measurement, reporting consistently smaller boundary error, roughly ten times lower in the terahertz example and about three times lower for dual-comb phase, with a few FFT-based conjugate-gradient iterations.","pith_inferences":["A natural next test is to use the same rubber-band formulation for interpolation or denoising where the weights encode measurement confidence rather than padding, since the convex least-squares core only needs a weight function and a passband.","The method's dependence on a hand-tuned weight (0.03 vs 0.06 in the demonstrations) suggests a practical autotuning rule: choose the weight by cross-validating the residual outside the padding, which the paper does not propose.","For imaging data, the same weighted-least-squares padding could suppress ringing at image borders, although the paper only demonstrates one-dimensional spectroscopy traces.","Because convergence is fast only when the solution resembles the naively filtered spectrum, pathological passbands or extremely steep signal edges may require more iterations; a fixed iteration budget could be risky in real-time instruments."],"forward_implications":["For stored spectroscopic traces that are already windowed, rubber band filtering can replace symmetric padding as the default preprocessing step, preserving sharp features at both ends instead of rounding them off.","Because the output is itself a padded, bandlimited signal whose padding has been reshaped by the weights, downstream operations like phase extraction or derivative estimation see a smooth extension rather than a jump.","The same weighted-fit idea extends directly to any bandpass shape, not just ideal filters, by changing which Fourier amplitudes $\\Omega$ are free parameters.","The error scaling in the three examples suggests that the advantage grows when the signal has steep slopes at the window edges, which is exactly where symmetric padding is most distorted."],"supporting_citations":[{"why":"Supplies the symmetric-padding baseline method whose boundary artifacts the rubber band filter is designed to beat.","marker":"[11–14]"},{"why":"Supplies the conjugate-gradient solution method used to solve the frequency-domain system $AF=b$ with FFT evaluations.","marker":"[15]"},{"why":"Identifies dual-comb spectroscopy phase signals, where the nonzero phase slope makes apodization inappropriate and symmetric padding artifacts matter.","marker":"[9, 10]"},{"why":"Provides the terahertz time-domain, FTIR interferogram, and dual-comb data used in the three demonstration comparisons.","marker":"[16–18]"}],"fun_headline_variants":["Rubber band filter eliminates edge artifacts in spectroscopy","Rubber band filter: artifact-free filtering in few FFTs","Optimal padding ends boundary distortion in spectroscopy","Rubber band filter cuts THz edge errors tenfold","Edge-free bandpass filtering via weighted padding"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the weighted least-squares fit, with the padded region's weight chosen by hand (roughly 0.03 to 0.06 in the demonstrations), reproduces the true bandlimited signal; since the validation reference is the full measurement that practical filtering never has, the method's success depends on that weight being close to optimal for unseen data.","fun_headline_variants_meta":{"raw":{"variants":["Rubber band filter eliminates edge artifacts in spectroscopy","Rubber band filter: artifact-free filtering in few FFTs","Optimal padding ends boundary distortion in spectroscopy","Rubber band filter cuts THz edge errors tenfold","Edge-free bandpass filtering via weighted padding"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000692,"raw_usage":{"total_tokens":3066,"prompt_tokens":815,"completion_tokens":2251,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":431,"completion_tokens_details":{"reasoning_tokens":2176}},"tokens_in":431,"tokens_out":2251,"duration_ms":17495,"temperature":1.0,"reasoning_tokens":2176,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:42:43.762896+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If, on a synthetic signal with a known bandlimit, no single small weight (or simple weight rule) keeps the boundary error consistently below symmetric padding across noise and cut-off choices, then the claim that the method recovers the bandlimited signal without edge artifacts fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the conjugate-gradient solution method used to solve the frequency-domain system $AF=b$ with FFT evaluations."}],"review_version":1}