{"id":"903364ea-b7a3-450e-80e2-16ef888f19ac","arxiv_id":"2509.10198","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"The log-likelihood-ratio detector for analytic Gaussian signals is expressed as correlations of cross-Bertrand time-frequency distributions, with a single simulation suggesting gains for phase-agnostic power-law chirp detection.","lead":"This paper rewrites the optimal detector for Gaussian signals as a correlation of Bertrand time-frequency distributions, a family suited to power-law chirps such as dispersed pulsar pulses. The author also reports a simulation in which this representation beats Wigner-Ville and spectrogram methods for phase-agnostic detection.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unproved cross-Bertrand identity (43) is internally inconsistent; since (46)-(47) depend on it, the central representation is unsupported.","rationale":"Reading the manuscript in good faith, the goal is clearly to give a Bertrand-class counterpart of Flandrin's Cohen-class representation of the optimal detector. For that, the argument must connect the double time-frequency integral of cross-Bertrand distributions to the product of ordinary time-domain inner products that appear in the detector. That connection is entirely made by equations (43)-(45). I find the reader's objection accurate: the printed identity has an unexplained weight mismatch and an inconsistent signal pairing. These are not minor typographical details, because the choice of pairing determines whether the right-hand side equals the squared projection |∫ r φ_i^*|^2 that constitutes the detector. I also note that the numerical simulation in Section VIII does not resolve the issue: it tests a discrete version of the same formulas, and its 'phase-agnostic' substitution of |·| for Re{·} is itself an unproved heuristic. Thus the central claim is not established as stated. A correct cross-Bertrand identity might exist for different parameters or a different pairing, but that would require rewriting the main equations. For these reasons, I agree with the reader's REJECT and recommend no change to that verdict. If the identity can be proved and corrected, the paper could be reconsidered, but as submitted the load-bearing step is missing.","tokens_in":16471,"tokens_out":11657,"duration_ms":111892,"concrete_test":"Re-derive (43) from the definition (42) for the Unterberger case k=-1, r=1/2, q=0, using the substitution a=f e^{u/2}, b=f e^{-u/2}. Verify whether the resulting left-hand side factorizes as the printed RHS, and which of the two pairings, (X1,X3)(X2,X4) or (X1,X2)(X3,X4), appears. If the correct factorization is not the printed one, or if the weight is f^0 rather than f^2, then the identity is false and the central representation is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, equations (46) and (47), is derived solely through the cross-Bertrand identity (43)-(45), introduced as 'a straightforward computation' with no derivation. As printed, the identity is internally inconsistent. With r=1/2 and q=0, equation (43) has f^{2q}=1 on the left but f^{2r+1}=f^2 inside each frequency integral on the right; equation (44) silently drops this factor and uses plain df. Independently, the Parseval statement (45) pairs time-domain integrals as (x1,x2) and (x3,x4), whereas the frequency-domain right-hand side of (43) pairs (X1,X3) and (X2,X4). These two pairings cannot both be correct. Since (46) and (47) are obtained by applying the time-domain pairing of (45) to the products in (38)-(39), any error here invalidates the claimed representation. No lemma, derivation, or reference supports (43), so the central bridge of the paper is unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives an optimal detector for analytic Gaussian signals using a generalized multivariate complex normal distribution and a likelihood-ratio argument. It then claims to represent this detector in the time-frequency domain as correlations of Bertrand-class distributions, specifically equations (46)-(47). The derivation relies on a new cross-Bertrand identity (43) which is said to follow by 'a straightforward computation.' The paper also gives two analytic examples (fully known signal and Rayleigh fading signal) and a numerical simulation for a power-law chirp with a phase-agnostic detector.","tokens_in":16782,"tokens_out":9297,"duration_ms":89566,"significance":"If the main representation were correct, it would extend Flandrin's Cohen-class formulation of optimal detection to Bertrand's class, which is better adapted to power-law chirps and hence relevant to pulsar and fast radio burst detection. The derivation of an optimal detector for analytic signals in Sections III-V is a useful contribution in itself. However, the central bridge between the detector and Bertrand's class, equations (43)-(45), is unproved and, as printed, inconsistent. The claimed representation (46)-(47) is therefore not established, and the numerical conclusions rest on an unsupported heuristic.","major_comments":[{"comment":"The cross-Bertrand identity (43) is asserted without proof and is internally inconsistent as printed. With r=1/2 and q=0, the factor f^{2r+1}=f^2 on the right-hand side of (43) is silently dropped in (44). Moreover, the time-domain pairing in (45) is (x1,x2),(x3,x4), whereas the frequency-domain right-hand side of (43) pairs (X1,X3),(X2,X4). These two pairings cannot both follow from the same identity. Since (46)-(47) are based on (44)-(45), this is a load-bearing error. A correct derivation of a cross-Bertrand Moyal-type formula, with explicit hypotheses, is required.","section":"Section VI, Eqs. (43)-(45)"},{"comment":"Even if a corrected version of (43) were supplied, the identity yields frequency integrals weighted by f^{2r+1}, not the unweighted inner products appearing in (38)-(39). For r=1/2 the weight is f^2, and no parameter choice in the printed formulas removes it. Thus the equality between (46) and (38) is not established. The detector would need to be reformulated with weighted inner products, or a different Bertrand-class parameterization would be needed.","section":"Section VI, Eq. (46)"},{"comment":"For a fully known signal (K=0), the Karhunen-Loève eigenvalues η_i^2 are zero, not E_g. The coefficient in (55) should therefore be 1/(2N0), not 1/(2N0+η_g^2). The statement that η_g^2 equals the energy of the expected signal confuses an eigenvalue of the covariance matrix with a signal energy, making the example incorrect.","section":"Section VII-A, Eq. (55)"},{"comment":"The phase-agnostic detectors are obtained by replacing Re{·} with |·| in (66)-(68). This heuristic is not derived from the likelihood ratio and is not claimed to be optimal. The numerical comparison uses a single signal (64) with one value of b and one noise model; the conclusion that Bertrand's representation is 'superior' for phase-agnostic detection is therefore not supported by the simulations as presented.","section":"Section VIII"}],"minor_comments":[{"comment":"Typo: 'more effective then' should be 'more effective than'.","section":"Abstract"},{"comment":"Typo: 'complex conjuration' should be 'complex conjugation'; 'power-low chirps' should be 'power-law chirps' (two occurrences).","section":"Section II"},{"comment":"In the enumeration of equations, '(66), (66)' should be '(66), (67), (68)'.","section":"Section VIII"},{"comment":"Typo: 'cross-Bernard's class' should be 'cross-Bertrand's class'.","section":"Section VII-A"}],"recommendation":"reject","confidential_remarks":"The central identity (43) is the linchpin of the paper. It is unproved, internally inconsistent, and the resulting representation (46)-(47) does not follow from the detector equations (38)-(39). The f^2 weighting and the pairing mismatch are not mere typos; they invalidate the claimed equivalence. A successful revision would require a correct derivation of a cross-Bertrand Moyal formula and a reworking of the representation, examples, and simulation. Given the scope of these changes, I do not see a path to acceptance within the manuscript's current framing."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: the new result the paper claims—representing the optimal detector as correlations of cross-Bertrand distributions—is not actually established. The identity that does all the work, eq. (43), is asserted without proof, and as printed it is internally inconsistent. The paper still has a decent section deriving the optimal detector for analytic signals, which is worth keeping.\n\nWhat's new and good: the derivation in Sections III–V is careful and self-contained. Starting from the generalized complex normal, it correctly handles the fact that analytic signals have correlated real and imaginary parts, does the Karhunen–Loève expansion, and lands on (38)–(39). That is a real contribution and I don't see a problem with it. The idea of looking for a Bertrand-class representation is also motivated reasonably, given that power-law chirps are the natural signals for this class.\n\nThe soft spots are load-bearing. Eq. (43) is introduced as a 'straightforward computation' with no derivation, no lemma, no reference. Then the specialization (44) drops the f^{2r+1}=f^2 weight that (43) puts in the frequency integrals. Even if that's a typo, the Parseval step (45) uses a different pairing—time-domain x1x2 times x3x4—whereas the right side of (43) pairs X1X3 with X2X4. These two cannot both be right. Since (46)–(47) are just (45) applied to the detector terms, the central representation collapses unless the identity is fixed.\n\nThe examples don't help. In VII-A, a fully known signal has K=0, so the KL eigenvalues should be zero, not E_g; the formula (55) uses 2N0+eta_g^2 as if there were signal energy. That's another inconsistency. The simulation is one favorable power-law chirp that the Unterberger distribution localizes perfectly; there is no code or data, and the phase-agnostic trick of replacing Re with |.| is asserted without justification.\n\nWho should read this: anyone working on time-frequency detection for pulsar/FRB searches will find the analytic-signal detector derivation useful, and the Bertrand representation idea is worth pursuing. But the paper as submitted is not ready. It deserves a serious referee—the core idea is plausible and the flaws look repairable—but the referee should ask for a complete proof of the identity, or the paper should be cut down to the detector derivation.\n\nMy recommendation: send it to peer review, but with a clear expectation of major revision. If the identity can't be repaired, the representation claim should be withdrawn.","headline":"The detector derivation is solid, but the cross-Bertrand identity that carries the main result is unproved and as printed inconsistent; the paper needs major repair before the representation claim can be trusted.","tokens_in":89,"tokens_out":3096,"would_cite":false,"duration_ms":70162,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["94A12","62M15"],"pacs":[],"model":"deepseek-v4-flash","headline":"A new derivation expresses the optimal detector for Gaussian analytic signals as a correlation of Bertrand-class time-frequency distributions, a form well-suited to power-law chirp detection.","keywords":["optimal detection","Bertrand's class","time-frequency distributions","power-law chirp","analytic signals","likelihood ratio","Gaussian signals","phase-insensitive detection"],"falsifier":"Evaluate numerically both sides of equation (43) for several simple analytic signals (e.g., two overlapping chirps with different delays) at r=1/2, q=0. If the integral of B^L_{X1X2} B^{A*}_{X3X4} dtdf is not equal to (∫ X1 X3* df)(∫ X2 X4* df)^*, then the representation (46)-(47) fails. A second check: compare the pairing structure of (45) with (43) — the printed version pairs x1x3 with x2x4 in one place and x1x2 with x3x4 in the other, so a direct derivation from the definitions of B^L and B^A is needed to identify the correct identity.","tokens_in":16318,"feed_emoji":"📡","tokens_out":5755,"duration_ms":54366,"temperature":0.7,"pith_summary":"The paper claims that the optimal detector—the likelihood-ratio test—for Gaussian analytic signals can be rewritten exactly as a correlation of Bertrand-class time-frequency distributions. This matters because Bertrand's class is built to localize power-law chirps, the type of dispersed signals seen in pulsar and fast-radio-burst observations, whereas the usual Wigner-Ville-based (Cohen-class) representation is tied to linear chirps. The author derives the optimal detector for analytic signals from first principles, then uses a cross-Bertrand Moyal-like identity to convert each pair of inner products into an integral over the time-frequency plane. If the identity holds, the detector becomes a sum of weighted time-frequency correlations, and simulations suggest it outperforms both the Wigner-Ville representation and spectrogram correlation when the signal phase is unknown.","feed_headline":"Optimal detection reduces to Bertrand time-frequency correlations","feed_subtitle":"Rewriting the likelihood-ratio test as a correlation of Bertrand-class distributions could sharpen phase-blind chirp detection.","key_machinery":"The load-bearing tool is the cross-Bertrand class (eq. 42), which generalizes Bertrand's time-frequency distribution to two different signals X and Y. The decisive step is a Moyal-like identity (eq. 43): the integral of B^L_{X1X2} (conjugated B^A_{X3X4}) times f^{2q} over the half-plane equals the product of two inner products, with the signal pairings permuted. By picking r=1/2, q=0 and applying Parseval's theorem, this identity turns the pair of inner products in the optimal detector into a single time-frequency correlation. The combination of the 'localized' weighting µ_L (which gives ideal power-law chirp localization) with the 'auxiliary' weighting µ_A (which supplies unitarity) is what","core_discovery":"The central result is equations (46) and (47). The random component of the log-likelihood ratio becomes a sum over the Karhunen-Loève modes: each term is a weight η_i^2/(2N0+η_i^2) times the time-frequency integral of B^L_{rr} with the conjugate of B^A_{φ_iφ_i}; the deterministic component is the analogous correlation between B^L_{rm} and B^A_{φ_iφ_i}. The derivation starts from the optimal detector for analytic signals under generalized complex Gaussian noise, expands inverse covariances via Karhunen-Loève, and then uses a cross-Bertrand Moyal-like identity to convert pairs of inner products into integrals over the time-frequency plane. For a fully known signal or a Rayleigh-fading signal,","pith_inferences":["Because the identity (43) is stated without proof and its printed specialization appears to contain a weight mismatch (f^2 vs plain df) and a signal-pairing mismatch, the entire representation hinges on a missing derivation; a direct numerical check of (43) on simple analytic signals would settle it.","If the identity is repaired to hold for general k, the same construction would yield optimal detectors for hyperbolic and other chirp models beyond the Unterberger case studied in the simulation.","The phase-agnostic improvement seen in simulation comes from replacing Re with |·|; that substitution is heuristic, so the reported advantage may depend on the specific signal, noise level, and threshold choice.","The single simulated chirp is perfectly localized by the Unterberger distribution; testing on non-localized, noisy, or multi-component chirps would clarify how general the performance gain is."],"forward_implications":["If equation (43) is correct, the optimal detector for Gaussian analytic signals can be computed as a weighted sum of time-frequency correlations, extending the known Cohen-class time-frequency formulation to power-law chirp signals.","For a fully known signal, the detector reduces to a correlation between B^L_{rg} and B^{A*}_{gg}, a generalization of earlier path-integration chirp detectors to arbitrary power-law chirps.","For a Rayleigh-fading signal, the detector is a single correlation between B^L_{rr} and B^{A*}_{gg}, which is naturally phase-insensitive.","The representation is not unique: µ_L and µ_A can be swapped, or replaced by the unitary µ_U, yielding equivalent detector forms.","Numerical simulations with a dispersed chirp signal indicate that the Bertrand representation matches the matched filter in phase-sensitive detection and outperforms both Wigner-Ville and spectrogram correlation when the phase is unknown."],"fun_headline_variants":["Optimal detector rewritten as Bertrand time-frequency correlation","Bertrand representation of optimal detector improves phase-blind chirp detection","Phase-blind chirp detection aided by Bertrand correlation","Log-likelihood ratio as sum of Bertrand correlations for chirp detection","Optimal Gaussian detection via Bertrand time-frequency inner products"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire representation depends on an unproved identity, introduced as a straightforward computation, that connects the integral of a product of cross-Bertrand distributions to a product of inner products; if that identity is wrong, the detector representation collapses.","fun_headline_variants_meta":{"raw":{"variants":["Optimal detector rewritten as Bertrand time-frequency correlation","Bertrand representation of optimal detector improves phase-blind chirp detection","Phase-blind chirp detection aided by Bertrand correlation","Log-likelihood ratio as sum of Bertrand correlations for chirp detection","Optimal Gaussian detection via Bertrand time-frequency inner products"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000697,"raw_usage":{"total_tokens":2928,"prompt_tokens":629,"completion_tokens":2299,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":373,"completion_tokens_details":{"reasoning_tokens":2219}},"tokens_in":373,"tokens_out":2299,"duration_ms":18675,"temperature":1.0,"reasoning_tokens":2219,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T18:07:08.469170+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate numerically both sides of equation (43) for several simple analytic signals (e.g., two overlapping chirps with different delays) at r=1/2, q=0. If the integral of B^L_{X1X2} B^{A*}_{X3X4} dtdf is not equal to (∫ X1 X3* df)(∫ X2 X4* df)^*, then the representation (46)-(47) fails. A second check: compare the pairing structure of (45) with (43) — the printed version pairs x1x3 with x2x4 in one place and x1x2 with x3x4 in the other, so a direct derivation from the definitions of B^L and B^A is needed to identify the correct identity.","supporting_citations":[],"review_version":1}