{"id":"dc253edb-14c4-4fa9-927a-ad7409d63ad1","arxiv_id":"2412.08254","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Padé approximants to the z-transform reveal separate complex-plane pole and zero patterns for signal and noise (Froissart doublets), a method the report tests on gravitational-wave templates and Virgo data.","lead":"A physics report argues that tracking the complex poles and zeros of a signal's z-transform, through noisy Padé approximants, can separate gravitational-wave signals from detector noise better than looking at Fourier spectrum peaks. It demonstrates the idea on toy models and on Virgo detector channels, compressing a seismic-to-mirror-suspension transfer function to a small rational function.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim depends on the faithful-pole versus Froissart-doublet separation that Section 17 itself leaves open; the applied proof of concept in Sections 13 and 16 lacks a quantitative, reproducible criterion, so the method is not yet demonstrated on noisy data.","rationale":"Section 17 explicitly designates the separation of faithful poles/zeros from Froissart doublets as the 'corner stone' of Padé Filtering and states the issue is 'still under discussion.' Section 16's central quantitative-looking result—the [5/5] skeleton for the Virgo seismic-to-suspension transfer function—is obtained by erasing doublets with no stated rule, and the 'check on other chunks' is not shown. Section 13's claim that local maxima of pole histograms must propagate between channels is asserted rather than derived; the relation for T_{O_{j+1}} includes an additive input term and possible pole-zero cancellations, so the claim is not automatic. These are not internal inconsistencies in the toy-model derivations, which are a legitimate strength of the report: the constant-signal, damped-cosine, ringdown, and AR(1)/AR(2) calculations are explicit and checkable. The gap is that none of these exercises validates the decision rule that is needed to apply the method to real data. A synthetic-chain test with known ground truth would settle whether the separation and the histogram-correspondence claims hold; if they do not, the applied claims in Sections 13 and 16 collapse. Since the concern is concrete but addressable by releasing data/code and adding quantitative validation, the reader's CONDITIONAL verdict stands.","tokens_in":26255,"tokens_out":5135,"duration_ms":55150,"concrete_test":"Run a synthetic two-stage chain with known rational transfer function H(z) (5 poles, 5 zeros), Gaussian inputs, and independent output noise, following the §13 relation exactly. Apply the paper's PF pipeline ([20/20] PAs, pole-density histograms, §16 doublet erasure) for ℳ=1300 draws over several chunk lengths and SNRs; compare recovered poles before/after erasure with true H poles, and repeat erasure with thresholds spanning the observed doublet-extension range. If recovered poles do not match true poles with high precision/recall, or if the skeleton varies with the erasure threshold, the central claim fails. For Virgo, estimate the [5/5] skeleton on one chunk and measure its prediction error against the empirical transfer function of a held-out chunk.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing premise is that, in practical noisy data, PA pole/zero statistics split into a 'faithful' class carrying signal information and a 'spurious' Froissart-doublet class that can be recognized and erased. The authors themselves call the clear separation 'the corner stone' and say it is 'still under discussion' (Section 17), citing robust Padé approximants and spurious-pole results. Two applications depend on this premise without supplying the missing criterion. (1) Section 13 asserts that singularities of T_{O_j} 'have to reappear' among those of T_{O_{j+1}} and hence that local maxima of pole histograms propagate between channels. This is not a consequence of the displayed z-transform relation: the second term T_{E_j} T_{L_{j+1}} T_{M_j} T_{I_j} contributes its own poles, and pole-zero cancellations can remove or shift maxima; no proof or numerical test on synthetic chains is given. (2) Section 16 derives a [5/5] transfer-function skeleton by 'erasing' doublets, but the erasure threshold is not specified and the promised check on other chunks of data is not reported. Without an operational rule, the [5/5] compression and the channel-to-channel pole correspondences are not reproducible, so the central claim that PF separates signal from noise in Virgo data is not demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript is an introductory report on 'Padé Filtering' (PF), the practice of using Padé Approximants of the z-Transform of a data set to locate complex poles and zeros and thereby separate signal from noise. Part A motivates the approach: a two-Lorentzian example (Section 3) and AR(1)/AR(2) computations (Section 8) show that real spectral peaks are ambiguous while complex-pole structures are stable, and the Kac and Froissart phenomena (Section 4) are presented as the probabilistic backdrop. Part B develops toy-model mathematics: the random z-Transform of white Gaussian data (Section 5), random polynomials (Section 7), and explicit noisy rational-interpolation and Padé models in which a 'Froissart polynomial' controls the statistics of pole-zero doublets (Sections 10-11). Part C applies PF to prototypical signals and noises (Section 14), to gravitational-wave templates (chirp, ringdown, pulsar; Section 15), and to Virgo data, where a [5/5] transfer-function 'skeleton' between a seismicity channel and a suspension channel is claimed (Sections 13 and 16), and the open state of the 'Froissart debate' is acknowledged (Section 17). The paper is explicitly a mixture of review and research report.","tokens_in":26703,"tokens_out":18622,"duration_ms":178243,"significance":"If it can be made operational, PF would give gravitational-wave analysts a genuinely different diagnostic: a data-compression and signal/noise-separation tool whose resolution is governed by complex-pole statistics rather than by spectral-peak separation. The paper's strengths are its explicit, checkable toy-model calculations; I verified the Lorentz two-line threshold S >= 2/sqrt(3) (Section 3), the AR(1)/AR(2) energy spectra and the AR(2) peak-transition curve r_2(phi) (Section 8), and the mean energy E[E_N] = 1 for white noise (Section 5). The Section 10 derivation of an explicit Cauchy doublet density with condensation probability 2/3 inside the interpolation interval is a concrete and useful result, and Section 17's candid statement of the open Froissart question is a model of scientific honesty. However, the two applied demonstrations -- the channel-to-channel pole reappearance (Section 13) and the Virgo [5/5] skeleton (Section 16) -- rest on exactly the faithful-versus-Froissart separation that Section 17 concedes is not settled, and neither demonstration is currently quantitative or reproducible.","major_comments":[{"comment":"Section 17, 'The Froissart debate': the manuscript itself states that the clear separation between faithful poles and spurious Froissart doublets is 'the corner stone of the Padé Filtering method' and that 'the situation is not clear cut and is still under discussion.' Both the channel-to-channel proof of concept (Section 13) and the Virgo [5/5] extraction (Section 16) interpret PA pole/zero statistics using precisely that separation, so the central applied claim of the report is conditional on a premise the paper leaves open. The report should either supply an operational, validated separation criterion for the noisy-data regime of Sections 13 and 16, or re-frame the Virgo and channel-cascade results as illustrations of an intended workflow whose validity awaits resolution of the Froissart question.","section":"Section 17"},{"comment":"Section 13, displayed z-transform relation: the sentence 'the singularities of T_{O_j}(z), for each random realization, have to reappear among the singularities of T_{O_{j+1}}(z)' does not follow from the displayed identity T_{O_{j+1}} = (T_{E_j} T_{L_{j+1}} / T_{L_j}) T_{O_j} + T_{E_j} T_{L_{j+1}} T_{M_j} T_{I_j}. A pole of T_{O_j} at z_0 is cancelled if T_{E_j} or T_{L_{j+1}} has a zero at z_0, and the second term contributes singularities of its own; moreover, for the finite data sets used here the z-transforms are polynomials, so the statement can only be meant at the level of PA poles, while PA's are nonlinear functions of the data. The histogram version ('the local maxima of the polar histogram ... have to reappear') is stronger still and is supported only by the two Virgo figures, whose content is qualitative. Please either prove the claim under stated genericity conditions (no pole-zero cancellations, known transfer functions) or test it on a synthetic cascade with known stages, and state what the histograms are expected to show in the presence of shared versus channel-specific poles.","section":"Section 13"},{"comment":"Section 16, [5/5] skeleton: the procedure is not reproducible as written. The doublet-erasure threshold is not specified ('we erase those doublets'), the number of erased doublets is not reported, the [20/20] PA is applied to a single realization, and the promised 'check on other chunks of data' is not presented. As a result, the claimed compression of the seismicity-to-suspension transfer function to a [5/5] rational fraction cannot be verified. Please specify the erasure criterion (distance in |z|, residue size, or histogram-based rule), report the counts and the variability across realizations, and show the held-out-chunk validation, or explicitly label the [5/5] as a single-realization illustration pending that validation.","section":"Section 16"},{"comment":"Section 16, Table 16: the correlation coefficients and Kullback-Leibler divergences are point estimates from single realizations with no uncertainty quantification and no stated null model. The within-channel KL entries are printed with an unclear exponent, and the histogram binning and number of draws are not given, so it is not possible to decide whether the contrasts (for example 6.14 versus 5.4) are statistically meaningful. For a statistical method whose output is a two-dimensional histogram, the proof of concept should include bootstrap or jackknife error bars (over draws or over disjoint time chunks) and a statement of the histogram resolution.","section":"Section 16, Table 16"}],"minor_comments":[{"comment":"Section 11: in the displayed formula for the [n-1/n] PA of the geometric series, the Froissart polynomial appears once as K_{n-1} and once as K_k; please use a single symbol and define its degree, and add an explicit inline citation to Gilewicz and Pindor (JCAM 105, 1999, already in the bibliography) at the point where the general rational case is asserted.","section":"Section 11"},{"comment":"Section 15, pulsar paragraph: the authors are explicit that the test uses b=4 as 'a model of the model' because the physical b=9 regime is not resolvable by the PA; the section heading 'testing on signals that are expected from the Gravitational Universe' should state up front that the pulsar test is a surrogate, so that readers do not mistake it for a test of the physical signal.","section":"Section 15"},{"comment":"Section 2: 'We prove the validity of the concept on actual Virgo data' overstates what Sections 13 and 16 establish given the caveats of Section 17; suggest 'provide a proof of concept subject to the open Froissart question' or an equivalent formulation.","section":"Section 2"},{"comment":"Table 16: print the diagonal KL entries in decimal form, state the histogram binning, the number of Monte Carlo draws, and the exact definition of the correlation factor between two pole distributions, so that the numbers can be reproduced.","section":"Table 16"},{"comment":"Bibliography and Sections 4/IX: the surname 'Kahn' should read 'Kahane' (J.-P. Kahane), and the year of the Porquerolles school is given as 2003 in Section 4 and in the bibliography but as 2005 in Section IX for E. Saff's lecture; please unify these references.","section":"Bibliography and Section 4"},{"comment":"Section 5a: the asymptotic statement that the z-Transform 'explodes like 2^{(N-1)/2}(1-rho)^{-1/2}' does not follow from the preceding variance expression E[T_N^2] = (1 - z^{2N})/(1 - z^2), whose square root in the regime rho^{2N} << 1 and rho roughly 1 scales as (1-rho)^{-1/2} with no exponential-in-N factor; please correct the prefactor or its typesetting.","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reads as lecture notes plus a research report; its archival novelty is concentrated in Sections 11-16, and the two Virgo-based sections are the weakest. The formal toy-model sections (3, 5, 8, 10) are sound and would be publishable on their own. I would ask the editor to weigh whether the journal can accept a paper whose central applied claim depends on a premise the authors themselves acknowledge in Section 17 to be unsettled; a major revision that reframes the applied sections as conditional illustrations and adds synthetic-cascade and held-out-chunk validation would resolve this. The reference list is dominated by the authors' own work and the small classical PF literature; this is understandable for a niche method, but it means the 'proof of concept' has not been independently replicated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take on 2412.08254, the Padé Filtering report by Fournier and Pichot du Mézeray. Read it as an introductory report, not a finished research paper, and you'll find it genuinely useful if you work on GW detector noise or rational approximation of noisy data. The core idea is coherent: stop reading real peaks of the energy spectrum and instead track the complex poles and zeros of the z-Transform via noisy Padé Approximants, treating Froissart doublets as the footprint of randomness. The paper is also honest — Section 17 admits the method's cornerstone, the clean separation of faithful poles from spurious doublets, is still under discussion.\n\nWhat's actually new and good. The AR(1) and AR(2) analysis in Section 8 is fully worked out and is the strongest part: it shows explicitly how the energy-spectrum peak structure changes (one peak, two peaks, flat transition) while the z-Transform poles sit at stable exterior locations. That is a real, checkable demonstration of the paper's motto. The noisy-interpolant toy model in Section 10 is exactly solvable — doublet position follows a Cauchy-Lorentz law, condensation inside the interval is computed. The ringdown template in Section 15 reduces the [49/50] PA to the exact [1/2] fraction with predicted pole positions. Citation pattern is proper, including the robust-Padé work that challenges the method.\n\nWhere it's soft. Section 13's proof of concept overstates things: the claim that singularities of one suspension channel 'have to reappear' in the next does not follow from the displayed z-Transform equation — the second term contributes its own poles, and the histogram claim is never tested on synthetic chains. Section 16's [5/5] transfer-function skeleton comes from erasing doublets with an unspecified threshold, and the promised check on other data chunks is not reported. Table 16's correlations and KL divergences have no error bars. None of this kills the report — it is explicitly pedagogical — but 'we prove the validity of the concept on actual Virgo data' is too strong as written.\n\nWho it is for: GW data analysts, and anyone wanting a real entry point into the Padé-filtering literature. It deserves a serious referee — the toy math is sound, and the right referee can push for the missing quantitative validation. Send to review, expect revisions.","headline":"Honest, useful introduction to Padé Filtering with solid toy-model mathematics and a candid admission of the open Froissart-separation problem, but the Virgo 'proof of concept' outruns what the data demonstration actually supports.","tokens_in":27172,"tokens_out":5793,"would_cite":true,"duration_ms":59366,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["41A21","30C15","62M10","94A12"],"pacs":["04.80.Nn","95.75.Wx"],"model":"deepseek-v4-flash","headline":"This report argues that the reliable way to separate a signal from noise in a data stream is to study the complex poles and zeros of the data's z-transform via Padé approximants, rather than the real peaks of its Fourier energy spectrum.","keywords":["Padé approximants","z-transform","Froissart doublets","gravitational wave data analysis","signal-noise separation","random polynomials","ARMA noise","Virgo interferometer"],"falsifier":"To test the central claim, synthesize a damped cosine with known resonance $z_0=e^{-\\lambda\\tau}e^{i\\omega\\tau}$, add white Gaussian noise, and run the paper's [10/10] Padé filtering with $N=20$ and many realizations: if the pole-density maximum does not converge to $z_0$ and remain separated from the Froissart doublets, the method fails. A second decisive check is the Virgo transfer property: the pole-histogram maxima of the seismic channel must reappear in the F0 suspension channel on independent non-overlapping data chunks; if they do not, the [5/5] skeleton is an artifact of the chosen chunk.","tokens_in":26082,"feed_emoji":"📡","tokens_out":9636,"duration_ms":91374,"temperature":0.7,"pith_summary":"Padé Filtering is a proposed way to separate signal from noise by reading the complex singularities of a data set's z-transform instead of the real peaks of its Fourier energy spectrum. The paper argues that Padé Approximants of the z-transform produce pole and zero patterns in the complex plane: the poles that reflect the true signal stay stable, while random noise shows up as nearby zero-pole pairs, the Froissart doublets. If this separation holds, gravitational-wave analysts could identify damped oscillations, binary chirp features, suspension resonances, and environmental couplings from pole statistics, and compress a long noisy relation between detector channels into a small rational transfer function—on Virgo data the seismic-to-suspension relation reduces to a [5/5] fraction. The report is pedagogical and includes a proof of concept on real Virgo channels, which makes the central claim testable.","feed_headline":"Look for complex poles, not spectral peaks","feed_subtitle":"Padé approximants map noisy gravitational-wave data into pole-zero pictures where noise leaves Froissart doublets.","key_machinery":"The load-bearing object is the Padé Approximant $[L/M](z)=P_L(z)/Q_M(z)$ of the z-transform $T_N(\\mathcal{D})(z)$, whose numerator zeros and denominator poles are tracked in the complex plane. Two probabilistic phenomena give the method its power: the Kac phenomenon (random polynomials have a positive density of real roots, so the real axis attracts PA poles in the noisy case) and the Froissart phenomenon (noise in the Taylor coefficients creates short zero-pole doublets, the Froissart doublets, located near the unit circle). The practical step that carries all the applications is separating these doublets from faithful poles and erasing them, which yields the compressed rational transfer functions used on Virgo data.","core_discovery":"The central claim is that the z-transform $T_N(\\mathcal{D})(z)=\\sum_{n=0}^{N-1} d_n z^n$ encodes the physics of a data set in its complex analytic structure, and that (possibly noisy) Padé Approximants $[L/M](z)=P_L(z)/Q_M(z)$ reveal that structure even for finite $N$. For the paper's model signals the faithful poles are explicit: a constant signal has a pole at $z=1$; a damped cosine has the complex-conjugate pair $1/a$ and $1/\\bar a$; AR(1) pink noise has a pole at $1/a$ outside the unit circle; AR(2) red noise has two poles at $r e^{\\pm i\\varphi}$. White Gaussian noise, by contrast, makes the unit circle a natural boundary of the infinite-$N$ z-transform, and the PA's respond by placing a dense set of poles there, broken near $\\pm 1$, together with Froissart doublets—zero-pole pairs with $O(\\varepsilon)$ separation that are the signature of randomness. The paper's motto is 'Better study the complex poles of the z-Transform of a process rather than the real peaks of its energy spectrum,' and it supports this by showing that Virgo pole-density maxima of one channel reappear in another and that erasing doublets leaves a [5/5] rational skeleton for the seismic-to-suspension transfer function.","pith_inferences":["If the Froissart doublet separation can be made quantitative, the density and angular spread of doublets could serve as a noise-level estimator, giving the method an error bar on the cleaned signal.","Comparing complex singularity patterns rather than time series might provide a new coincidence test between separate detectors (LIGO, Virgo, KAGRA), especially in the low signal-to-noise regime where spectral peaks are ambiguous.","The contrast between the real-axis pole attraction and the unit-circle condensation suggests a cheap diagnostic for non-Gaussian noise: real-axis PA pole density should be sensitive to non-Gaussian contaminants, a prediction that can be checked with simulated non-Gaussian noise."],"forward_implications":["Gravitational-wave searches could stop hunting for peaks in the energy spectrum and instead ask whether a particular pole of the z-transform appears consistently across realizations of the noisy data.","Noise types become identifiable by their complex-plane signatures: white noise fills the unit circle with poles except near $\\pm 1$, AR(1) noise concentrates poles near $1/a$, AR(2) noise creates a pair at $r e^{\\pm i\\varphi}$, and these signatures persist when a signal is added.","Erasing Froissart doublets gives a data-compression tool: the Virgo seismic-to-suspension coupling is represented as a [5/5] rational transfer function rather than 432,000-sample channels, and the same idea generalizes to other linear apparatus.","For modeled gravitational-wave signals the PA reproduces the expected analytic structures—the chirp's accumulating singularity, the ringdown's [1/2] fraction with one zero and complex-conjugate poles, and the pulsar's slowly drifting frequency as secondary peaks near the unit circle."],"supporting_citations":[{"why":"Supplies the Kac phenomenon—random polynomials have a positive density of real roots—which shapes the real-axis component of PA pole statistics in noisy data.","marker":"M. Kac (1943)"},{"why":"Introduced the zero-pole doublets near the unit circle that this report identifies as the footprint of noise in Padé Approximants.","marker":"M. Froissart (1969)"},{"why":"First demonstrated PA-based noise filtering on scattering data and the data-compression effect that motivates Padé Filtering.","marker":"D. Bessis (1996)"},{"why":"Proved that PA poles of the z-transform of complex Gaussian white noise condense uniformly on the unit circle, providing the white-noise baseline.","marker":"P. Barone (2005)"},{"why":"Worked out the Froissart doublet structure for noisy Padé Approximants of a geometric series, the simplest rational case.","marker":"J. Gilewicz and M. Pindor (1997)"},{"why":"Established the analogous Froissart phenomenon for rational interpolation from stochastic data, with doublets condensing inside the interpolation interval.","marker":"J.-D. Fournier and M. Pindor (2000)"},{"why":"Applied Padé Filtering to gravitational-wave burst identification using complexified two-detector data, extending the method to coincident searches.","marker":"L. Perotti, T. Regimbau, D. Vrinceanu, D. Bessis (2014)"},{"why":"Introduced robust Padé Approximants via SVD to suppress spurious poles, representing the alternative approach contested in the paper's Froissart debate.","marker":"P. Gonnet, S. Guttel, L. N. Trefethen (2013)"}],"fun_headline_variants":["Padé filtering: track poles, not spectral peaks","Complex poles beat spectral peaks in noise analysis","Padé approximants expose poles hidden in noise","Noise leaves Froissart doublets: read the poles","Padé filtering: z-transform poles reveal signal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole method rests on being able to tell the faithful poles from the spurious Froissart doublets in real noisy data; the paper explicitly notes that this separation is still under discussion and gives no precise rule for the threshold used when erasing doublets.","fun_headline_variants_meta":{"raw":{"variants":["Padé filtering: track poles, not spectral peaks","Complex poles beat spectral peaks in noise analysis","Padé approximants expose poles hidden in noise","Noise leaves Froissart doublets: read the poles","Padé filtering: z-transform poles reveal signal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000118,"raw_usage":{"total_tokens":1050,"prompt_tokens":878,"completion_tokens":172,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":98}},"tokens_in":494,"tokens_out":172,"duration_ms":2596,"temperature":1.0,"reasoning_tokens":98,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:01:10.820366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"To test the central claim, synthesize a damped cosine with known resonance $z_0=e^{-\\lambda\\tau}e^{i\\omega\\tau}$, add white Gaussian noise, and run the paper's [10/10] Padé filtering with $N=20$ and many realizations: if the pole-density maximum does not converge to $z_0$ and remain separated from the Froissart doublets, the method fails. A second decisive check is the Virgo transfer property: the pole-histogram maxima of the seismic channel must reappear in the F0 suspension channel on independent non-overlapping data chunks; if they do not, the [5/5] skeleton is an artifact of the chosen chunk.","supporting_citations":[],"review_version":1}