{"id":"39e5f87c-1c94-4c7d-a7dc-b4cf03ee6c4b","arxiv_id":"2607.06324","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":2,"one_line_summary":"A signed multifractal detrended cross-correlation coefficient is introduced that remains strictly bounded in [-1,1] for all fluctuation orders q by averaging locally normalized detrended correlations with regularized amplitudes.","lead":"The paper proposes a new correlation coefficient (ρ_SMFDCCA) that stays bounded between -1 and 1 for all fluctuation orders q, fixing numerical instabilities in multifractal cross-correlation analysis for negative q. A smart generalist might read this to gain a tool that reliably measures how large vs. small fluctuations in two time series co-move without numerical blow-ups.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Boundedness proof is correct, but no direct comparison with ρ_q on the same data is provided, leaving the claim of being a 'replacement' for the standard multifractal coefficient unsubstantiated.","rationale":"The reader correctly identified the boundedness argument as mathematically sound and correctly flagged the multifractal-information-preservation question as the weakest assumption. I agree on both points. However, I think the concern is more concrete and more consequential than the reader's verdict reflects.\n\nThe reader rated the concern as 'acknowledged by the authors' and left the verdict at ACCEPT. I would adjust to CONDITIONAL because the absence of any direct comparison with ρ_q is a substantive gap, not merely an acknowledged limitation. The paper's central framing — from the title through the abstract to the summary — positions ρ_SMFDCCA as solving the boundedness problem of ρ_q, implying it is a drop-in replacement. But the two quantities have different mathematical structure (signed vs. absolute covariance, different normalization, different connection to scaling exponents), and without empirical evidence that they track similar structure, the claim of 'overcoming' ρ_q's limitations is not fully substantiated.\n\nThe fGn null-model validation (Figure 1) confirms the coefficient doesn't produce spurious correlations — good, but this only tests the null case. The empirical applications (Figures 2-3) show the coefficient produces sensible-looking results, but without a benchmark against ρ_q, we cannot tell whether the q-dependent patterns observed (e.g., stronger synchronization for large fluctuations in DJI-NASDAQ) are a genuine feature of the cross-correlation structure or an artifact of the specific weighting scheme chosen.\n\nThe ε sensitivity analysis is thorough and convincing. Code availability is a positive. The mathematical derivation is clean. These are real strengths. But the paper would be substantially strengthened by the single additional analysis described in the concrete test — a head-to-head comparison with ρ_q for q > 0 — which is low-cost and would directly address whether the new coefficient is measuring the same kind of structure as the object it claims to replace.","tokens_in":9505,"tokens_out":3543,"duration_ms":295692,"concrete_test":"Compute both ρ_SMFDCCA(n,q) and the standard ρ_q(n) from Kwapien et al. [12] on the DJI-NASDAQ dataset for q ∈ {1,2,3,4,5} across the same scales n used in Figure 2. Plot both profiles as functions of q for fixed n. If the two curves show qualitatively similar q-dependence (e.g., both monotonically increasing, similar crossover scales), then ρ_SMFDCCA can be reasonably presented as a bounded extension of ρ_q into the negative-q regime. If the curves diverge qualitatively (different monotonicity, different scale dependence), then ρ_SMFDCCA is a distinct observable and the paper's framing as a 'replacement' for ρ_q should be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The mathematical claim that ρ_SMFDCCA ∈ [-1,1] for all q ∈ ℝ is correct: each r_v is bounded by Cauchy-Schwarz (Eq. 9-10), weights w_v = A_v^{q/2} are non-negative (Eq. 11, 16), and a convex combination of bounded quantities is bounded (Eq. 17). This is not in dispute.\n\nThe load-bearing concern is different: the paper frames ρ_SMFDCCA as overcoming the limitations of the standard ρ_q(n) and eliminating its corrective procedures, but the two quantities are structurally different observables. The standard ρ_q(n) is a ratio of q-th-order fluctuation functions F_xy(q,n)/√(F_xx(q,n)·F_yy(n)), connecting to multifractal scaling exponents H(q). ρ_SMFDCCA is an amplitude-weighted average of local Pearson-type correlation coefficients r_v. Even for q > 0 where both are well-defined, there is no reason a priori that they should track the same structure. For example, at q=2, ρ_SMFDCCA = Σ C_xy(n,v) / Σ A_v(n) (from Eq. 14), whereas ρ_2(n) = Σ|C_xy(n,v)| / √(Σ F²_x · Σ F²_y) — different numerator (signed vs. absolute), different normalization.\n\nThe paper provides no head-to-head comparison of ρ_SMFDCCA and ρ_q on any dataset, even in the q > 0 regime where ρ_q is stable. Without this, we cannot assess whether the new coefficient captures comparable cross-correlation structure or is an entirely different descriptive statistic that merely shares the parameter name 'q.' The authors acknowledge the coefficient 'is not designed to estimate generalized Hurst exponents, multifractal spectra, or singularity distributions,' which is honest, but the title and abstract still label it a 'Multifractal Detrended Cross-Correlation Coefficient,' implying functional equivalence with the object it replaces.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"The manuscript introduces a Signed Multifractal Detrended Cross-Correlation Coefficient, ρ_SMFDCCA(n, q), designed to remain strictly bounded within [-1, 1] for all q ∈ ℝ, including negative fluctuation orders where the standard multifractal cross-correlation coefficient ρ_q(n) becomes numerically unstable. The construction is straightforward: local detrended correlation coefficients r_v(n) are defined per segment (Eq. 9), shown to be bounded by Cauchy-Schwarz (Eq. 10), and combined via a weighted average using non-negative amplitude-based weights w_v = A_v^{q/2} (Eqs. 11, 16), yielding a bounded convex combination (Eq. 17). The method is validated on independent fGn (Fig. 1) and applied to U.S. equity indices (Fig. 2) and Brazilian weather data (Fig. 3). Code and data are publicly available.","tokens_in":10358,"tokens_out":1128,"duration_ms":188346,"significance":"The boundedness proof is correct and derived from first principles: each r_v is bounded by Cauchy-Schwarz (Eq. 10), weights are non-negative (Eq. 16), and a weighted average of bounded quantities is bounded (Eq. 17). The regularization parameter ε = 10^{-12} is a numerical floor, not a fitted constant, and the sensitivity analysis in the End Matter confirms a broad stability plateau across fourteen orders of magnitude. The fGn null-model validation (Fig. 1) appropriately confirms near-zero values (~10^{-3}). The authors provide reproducible code and data via Zenodo, which is a strength. The paper is transparent about the observable's scope, explicitly stating it is not designed to estimate generalized Hurst exponents or multifractal spectra.","major_comments":[{"comment":"The central framing of the paper positions ρ_SMFDCCA as overcoming the limitations of the standard ρ_q(n) and eliminating its corrective procedures (Abstract, Introduction, Summary). However, the two quantities are structurally different observables: ρ_q(n) is a ratio of q-th-order fluctuation functions F_xy(q,n)/√(F_xx·F_yy), while ρ_SMFDCCA is an amplitude-weighted average of local Pearson-type correlation coefficients r_v. Even for q > 0 where both are well-defined, there is no a priori reason they should track the same structure (e.g., at q=2, ρ_SMFDCCA involves signed covariances C_xy in the numerator via r_v, whereas ρ_2 involves a different normalization). No head-to-head comparison of ρ_SMFDCCA and ρ_q on any dataset is provided, even in the q > 0 regime where ρ_q is stable. Without this, the claim of being a 'replacement' is unsubstantiated. The authors should either (a) add a直接","section":null},{"comment":"The use of the term 'multifractal' for ρ_SMFDCCA requires more careful justification. The authors acknowledge (paragraph following Eq. 17) that the coefficient 'is not designed to estimate generalized Hurst exponents, multifractal spectra, or singularity distributions.' The q-dependence here is an amplitude-weighting of local correlation coefficients, not the moment-based scaling that defines multifractal analysis. The paper should more clearly delineate what 'multifractal' means in this context versus the standard usage, and whether the observable captures amplitude-dependent cross-correlation structure comparable to the multifractal formalism.","section":null}],"minor_comments":[{"comment":"Eq. (17): the subscript reads ρ_MFDCCA but should be ρ_SMFDCCA for consistency.","section":null},{"comment":"The Limitations section uses W_i and Amplitude_i (Eq. 18) whereas the main text uses w_v and A_v. Please unify the notation.","section":null},{"comment":"Fig. 1: the caption states maximum absolute deviations below 5×10^{-3}, but the figure itself does not clearly show this scale. Consider adding a colorbar with finer resolution or showing the maximum in an inset.","section":null},{"comment":"The text following Eq. (8) states ε is chosen 'without loss of generality.' This is a mathematical claim that is not quite appropriate for a numerical parameter; consider rephrasing the 'without loss of generality' statement.","section":null},{"comment":"Reference [13] is cited as Physica A 689, 131424 (2026). The volume/year appears unusual; please verify.","section":null},{"comment":"The phrase 'proving the efficacy of the proposed observable in characterize multiscale co-movements' (U.S. equity markets section) contains a grammatical error; please revise.","section":null},{"comment":"The Introduction contains several instances of passive voice without clear subjects (e.g., 'was proposed the DCCA cross-correlation coefficient,' 'was proposed the multifractal cross-correlation coefficient'). Please revise for clarity.","section":null}],"recommendation":"major_revision","confidential_remarks":"The boundedness result is correct and straightforward. The main concern is whether the paper overclaims by framing ρ_SMFDCCA as a replacement for ρ_q without a direct comparison. A single figure comparing the two observables on the same data (even just for q > 0) would substantially strengthen the paper and could be added with moderate effort. If the authors reframe the contribution as a complementary bounded observable rather than a replacement, the paper could also be acceptable with minor revisions."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive reading of our manuscript. The referee correctly identifies the boundedness proof as the core contribution and acknowledges the numerical stability, null-model validation, and reproducibility. Below we address each major comment in turn.","responses":[{"response":"The referee raises a valid and important point. We agree that ρ_SMFDCCA(n,q) and ρ_q(n) are structurally distinct observables: ρ_q is a ratio of q-th-order fluctuation functions F_xy(q,n)/√(F_xx(q,n)·F_yy(q,n)), whereas ρ_SMFDCCA is an amplitude-weighted convex combination of locally normalized Pearson-type correlation coefficients r_v(n). Because the normalization and averaging structures differ, there is no a priori guarantee that the two observables coincide even in the q > 0 regime, and we should not have framed ρ_SMFDCCA as a direct 'replacement' for ρ_q without empirical evidence of their relationship. We will make two changes in the revised manuscript. First, we will add a head-to-head comparison of ρ_SMFDCCA and ρ_q on the same datasets (fGn, DJI–IXIC, and the São Paulo weather pairs) in the q > 0 regime where ρ_q is well-defined and stable. This will show where the two observables agree, where they diverge, and what each reveals. Second, we will revise the framing throughout the Abstract, Introduction, and Summary to position ρ_SMFDCCA as a complementary bounded observable that addresses the specific problem of q < 0 instability, rather than as a wholesale replacement for ρ_q. The claim that ρ_SMFDCCA 'eliminates the corrective procedures required by previous approaches' will be retained in the narrow sense that, for applications requiring bounded signed cross-correlation measures across the full q-range including q < 0, our formulation does not require post-hoc inversion—but we will make clear that this applies to the specific use case of negative-q analysis, not to a claim of equivalence with ρ_q for q > 0.","revision_made":"yes","referee_comment":"Major Comment 1: The central framing positions ρ_SMFDCCA as overcoming limitations of ρ_q and eliminating corrective procedures, but the two are structurally different observables. No head-to-head comparison is provided, even for q > 0 where ρ_q is stable. The claim of being a 'replacement' is unsubstantiated."},{"response":"We agree that the term 'multifractal' requires careful delineation. The manuscript already contains two passages that address this distinction: the paragraph following Eq. 17 states that 'the coefficient is not designed to estimate generalized Hurst exponents, multifractal spectra, or singularity distributions,' and the Introduction notes that 'the term multifractal refers to the q-dependent weighting of fluctuation amplitudes.' However, we acknowledge that these statements are not sufficiently prominent or explicit, and a reader could reasonably expect a more thorough justification of the terminology. In the revised manuscript, we will add a dedicated paragraph (or an explicit remark) that clearly states the following: (i) in standard multifractal analysis (MFDFA, MFCCA), q-dependence arises through moment-based scaling of partition functions, yielding generalized Hurst exponents and multifractal spectra; (ii) in ρ_SMFDCCA, q-dependence arises through amplitude-based weighting of bounded local correlation coefficients, producing an amplitude-stratified family of correlation observables rather than a scaling exponent; (iii) we retain the term 'multifractal' because the observable operates within the same detrended fluctuation framework and uses the same q-dependent amplitude-filtering mechanism that underlies MFCCA, but we will explicitly state that the observable does not produce multifractal spectra or scaling exponents. We believe this terminology is defensible given that the q-filtering of fluctuation amplitudes is shared with the multifractal formalism, but we agree that the distinction must be stated more prominently and unambiguously.","revision_made":"yes","referee_comment":"Major Comment 2: The use of the term 'multifractal' for ρ_SMFDCCA requires more careful justification. The q-dependence is amplitude-weighting of local correlation coefficients, not moment-based scaling. The paper should more clearly delineate what 'multifractal' means here versus standard usage."}],"tokens_in":9223,"tokens_out":1184,"duration_ms":129112,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The main thing to know: this paper proposes a new cross-correlation coefficient ρ_SMFDCCA(n,q) that is provably bounded in [-1,1] for all q ∈ ℝ, including negative q. The construction is simple and correct — compute local detrended correlation coefficients r_v (each bounded by Cauchy-Schwarz), then take a q-weighted average with non-negative weights A_v^{q/2}. A convex combination of bounded quantities is bounded. That's the whole argument, and it works. The fGn null-model validation (independent series, H=0.5) correctly shows near-zero values (~10⁻³) across all scales and q ∈ [-10,10], and the ε sensitivity analysis over fourteen orders of magnitude shows zero numerical dependence. Code and data are on Zenodo. This is a real methodological contribution — the boundedness is not in dispute, and the formulation is non-obvious compared to the standard ρ_q of Kwapień et al., which requires post-hoc inversion when values exceed unity for q<0. The empirical applications (DJIA vs. NASDAQ, Brazilian weather data) are illustrative and show sensible amplitude-dependent correlation structures. The authors are honest that the coefficient does not estimate generalized Hurst exponents or singularity spectra. Now the soft spot, and it's the one the stress-test correctly identifies: the paper frames ρ_SMFDCCA as overcoming the limitations of ρ_q and eliminating its corrective procedures, but the two are structurally different observables. The standard ρ_q is a ratio of q-th-order fluctuation functions connecting to multifractal scaling exponents H(q). ρ_SMFDCCA is an amplitude-weighted average of local Pearson-type correlations. Even at q=2 where both are well-defined, the numerators differ (signed vs. absolute cross-covariance) and the normalizations differ. There is no head-to-head comparison on any dataset, even in the q>0 regime where ρ_q is stable. Without that, we cannot assess whether the new coefficient captures comparable cross-correlation structure or is a different descriptive statistic that merely shares the parameter name 'q.' The title and abstract imply functional equivalence with the object it replaces, which is not substantiated. This is a load-bearing concern for the framing but not for the mathematical claim itself. The reader's assessment (soundness 7, novelty 6) is about right. The circularity burden is correctly low — the boundedness proof uses Cauchy-Schwarz from first principles with no fitting to the target result. The stress-test concern lands fully on reading the paper. Who benefits: practitioners in econophysics and climate science who need a bounded, amplitude-resolved cross-correlation measure and are not specifically trying to extract Hurst exponents. The paper deserves a serious referee who should push the authors to either (a) add a direct comparison with ρ_q on the same data in the q>0 regime, or (b) reframe the contribution as a complementary observable rather than a replacement.","headline":"Bounded-by-construction multifractal cross-correlation coefficient; the math is clean but the 'replacement' framing overpromises since no head-to-head comparison with the standard ρ_q is provided.","tokens_in":10634,"tokens_out":711,"would_cite":true,"duration_ms":146875,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Bounded multifractal cross-correlation coefficient works for all q","keywords":["multifractal cross-correlation","detrended fluctuation analysis","Cauchy-Schwarz inequality","nonstationary time series","fluctuation order q","bounded correlation coefficient","financial markets","meteorological data"],"falsifier":"Generate two coupled time series with a known, analytically specified multifractal cross-correlation structure (e.g., coupled binomial cascades with a known cross-correlation scaling exponent). Compute both ρ_SMFDCCA(n,q) and the standard ρ_q(n) across the full q-range. If ρ_SMFDCCA fails to reproduce the known amplitude-dependent correlation structure that ρ_q captures in the positive-q regime, or if the q-dependent profile of ρ_SMFDCCA does not correspond to any recoverable multifractal cross-correlation spectrum, then the coefficient is a bounded correlation measure but not a multifractal","tokens_in":9778,"feed_emoji":"📊","tokens_out":1476,"duration_ms":168183,"temperature":0.7,"pith_summary":"The paper introduces a new cross-correlation coefficient for nonstationary time series, called the Signed Multifractal Detrended Cross-Correlation Coefficient (ρ_SMFDCCA), that remains strictly bounded within [-1, 1] for every value of the fluctuation order parameter q, including negative q. Prior multifractal cross-correlation coefficients suffer from numerical instabilities, unbounded values, and erratic behavior when q is negative, because negative q amplifies segments with very small fluctuations, making the coefficient extremely sensitive to noise and numerical error. The existing workaround is a corrective post-processing step that inverts the coefficient whenever it exceeds the expected bounds. The authors eliminate this problem by changing what is averaged: instead of computing a global ratio of q-weighted covariance to q-weighted variance (which can diverge), they first compute a local correlation coefficient r_v in each segment, which is individually bounded by the Cauchy-Schwarz inequality, and then take a weighted average of these bounded local coefficients using amplitude-based weights. Because a weighted average of numbers each lying in [-1, 1] must itself lie in [-1, 1], the result is automatically bounded for all q, with no corrective procedures needed. The authors validate the coefficient on independent fractional Gaussian noise (no spurious cross-correlations detected), Dow Jones versus NASDAQ daily returns (stronger synchronization during large-amplitude fluctuations than small-amplitude ones, with convergence at longer scales), and meteorological data from São Paulo (heterogeneous coupling patterns between temperature and humidity variables across scales and amplitudes).","feed_headline":"Bounded multifractal cross-correlation coefficient works for all q","feed_subtitle":"New formulation stays within [-1,1] for negative fluctuation orders, eliminating instabilities that have limited multifractal cross-corridor","key_machinery":"Local detrended correlation coefficient r_v(n) = C_xy(n,v) / A_v(n), bounded by Cauchy-Schwarz; amplitude weights w_v(q,n) = A_v(n)^{q/2}; signed weighted average ρ_SMFDCCA(n,q) = Σ w_v r_v / Σ w_v; regularization floor ε = 10^{-12} preventing division-by-zero when local variances vanish; the Cauchy-Schwarz inequality as the structural guarantee of boundedness.","core_discovery":"The central discovery is a reformulation of the multifractal cross-correlation coefficient that achieves boundedness not by post-hoc correction but by construction. The key move is to normalize locally first—computing a per-segment detrended correlation coefficient r_v = C_xy(n,v) / A_v(n), where C_xy is the local detrended covariance and A_v is the geometric mean of the local detrended variances—and only then apply the q-dependent amplitude weighting. Because each r_v is individually confined to [-1, 1] by Cauchy-Schwarz, and because the weights w_v = A_v^{q/2} are non-negative, the weighted average ρ_SMFDCCA = Σ w_v r_v / Σ w_v is a convex combination of bounded quantities and is therefore","pith_inferences":["The paper's acknowledgment that ρ_SMFDCCA is not designed to estimate generalized Hurst exponents or multifractal spectra raises the question of whether the local-normalization step alters the scaling information content relative to the standard q-dependent fluctuation formalism. A direct comparison of the scaling exponents derivable from the two formulations on data with known multifractal struct","The choice to weight by A_v^{q/2} rather than A_v^q means the weighting exponent is half that used in standard MFDFA/MFCCA. This halving may affect the effective q-range sensitivity of the observable relative to the standard formalism, and a systematic comparison of the q-dependent profiles produced by the two weighting schemes on the same data would test whether the amplitude-conditioning is quan","The regularization parameter ε is shown to be numerically inert over 14 orders of magnitude, but this insensitivity itself suggests that for sufficiently small local amplitudes the coefficient effectively truncates the weighting of those segments. For data with genuinely near-zero-amplitude regions (e.g., overnight periods in intraday financial data), this truncation could systematically bias the "],"forward_implications":["Any downstream method built on the multifractal cross-correlation coefficient—including the recently proposed multivariate multiple cross-correlation coefficient qDMC²_x(n)—can inherit the boundedness property by adopting the local-normalization-first construction, eliminating embedded instabilities for negative q.","Empirical studies that previously restricted analysis to q ≥ 0 to avoid instability can now probe the small-fluctuation regime (q < 0) with a well-defined, interpretable coefficient, potentially revealing amplitude-dependent correlation structures that were systematically missed.","The amplitude-stratified view of cross-correlations could serve as a diagnostic tool for regime identification in financial markets: the observed convergence of q-dependent correlations at long scales suggests a transition from heterogeneous (amplitude-dependent) to homogeneous (amplitude-independent) coupling, which could be used to characterize market maturity or stress regimes."],"fun_headline_variants":["Signed multifractal cross-correlation coefficient stays bounded by construction","Bounded cross-correlation coefficient resolves negative-q instabilities","Local normalization keeps multifractal cross-correlations within [-1,1] for all q","Stable multifractal cross-correlation coefficient removes need for post-hoc correction","Bounded-by-design coefficient reveals scale- and amplitude-dependent correlations"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The claim that ρ_SMFDCCA is a meaningful multifractal observable depends on the assumption that weighting local correlation coefficients by the local fluctuation amplitude raised to the power q/2 preserves the multifractal cross-correlation structure in a way comparable to the standard q-dependent fluctuation formalism. The paper itself states that the coefficient is not designed to estimate generalized Hurst exponents, multifractal spectra, or singularity distributions, so a","fun_headline_variants_meta":{"raw":{"variants":["Signed multifractal cross-correlation coefficient stays bounded by construction","Bounded cross-correlation coefficient resolves negative-q instabilities","Local normalization keeps multifractal cross-correlations within [-1,1] for all q","Stable multifractal cross-correlation coefficient removes need for post-hoc correction","Bounded-by-design coefficient reveals scale- and amplitude-dependent correlations"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":671,"prompt_tokens":593,"completion_tokens":78,"prompt_tokens_details":null},"tokens_in":593,"tokens_out":78,"duration_ms":38800,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T09:43:55.111864+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"Generate two coupled time series with a known, analytically specified multifractal cross-correlation structure (e.g., coupled binomial cascades with a known cross-correlation scaling exponent). Compute both ρ_SMFDCCA(n,q) and the standard ρ_q(n) across the full q-range. If ρ_SMFDCCA fails to reproduce the known amplitude-dependent correlation structure that ρ_q captures in the positive-q regime, or if the q-dependent profile of ρ_SMFDCCA does not correspond to any recoverable multifractal cross-correlation spectrum, then the coefficient is a bounded correlation measure but not a multifractal","supporting_citations":[],"review_version":1}