{"id":"725df352-ac6b-423f-9e60-8ae895ea4e83","arxiv_id":"2607.17066","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Mapping solution snapshots to mass-coordinate (CDT) space before POD compresses transport-dominated 1D conservative PDEs into far fewer modes, with proven Kolmogorov-width bounds.","lead":"This paper develops a reduced-order modeling method that transforms each snapshot of a 1D conservative PDE into mass coordinates before applying POD, then transforms back. The theoretical bounds and numerical experiments show this can represent transport-dominated solution families with far fewer modes than standard POD.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pre-shock condition limits the sharper hyperbolic n-width claim to times before generic shock formation; post-shock only O(1/n) is proven.","rationale":"The reader's weakest_assumption identified two issues: the pre-shock hypothesis for the sharper hyperbolic rates, and the monotonicity of the truncated CDT map. I focus on the pre-shock assumption because it directly controls the strongest theoretical claim (Theorem 4.5), whereas monotonicity affects only the practical inverse step. The derivations in the paper appear internally consistent and the robust O(1/n) bound (Theorem 4.2) is independent support for entropy solutions; the concern is about the scope and advertised strength of the sharper result, not its correctness for the stated regime. Since the paper openly acknowledges the missing post-shock higher-rate bound, the appropriate verdict remains CONDITIONAL: the central theoretical claim is partially proven, with its sharpest form restricted to a generically finite pre-shock window. The proposed numerical test can indicate whether this is a practically essential restriction or merely a missing theoretical result.","tokens_in":45269,"tokens_out":16330,"duration_ms":167216,"concrete_test":"Take inviscid Burgers with a smooth compactly supported initial datum and known shock time t_* (using the provided code). For T > t_*, form CDT snapshots and compute the discrete Kolmogorov width / singular-value decay of the transformed snapshot matrix for n = 1,...,20; fit the decay exponent. If the exponent is ≥ 2, the pre-shock restriction is not practically binding; if it is ≈ 1, the O(1/n) bound is the correct post-shock ceiling and the stronger advertised claim applies only before shocks.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The advertised O(n^-2) hyperbolic bound (Theorem 4.5) rests on the hypothesis 1 + t f''(u0(ξ))u0'(ξ) ≥ α > 0 for all t ∈ [0,T]. For any C^1 initial datum with min u0' < 0 and f''(u0) ≠ 0, this quantity crosses zero at the first shock time t_*; for T > t_* the hypothesis fails identically. Thus the sharper estimate covers only the pre-shock portion of the trajectory. For post-shock entropy solutions, the only proven estimate is Theorem 4.2's O(1/n), which is an upper bound in the Wasserstein/L2(dr) metric, not in native L2, and is weaker than the claimed 'sharper compressibility.' The paper is explicit about this (Section 6: 'obtaining a higher-regularity post-shock counterpart ... remains open'), so this is a limitation rather than an internal inconsistency. But it is the most load-bearing point of the central claim: without the pre-shock condition, the theoretical advantage of CDT for the very shock-dominated regime motivating the method reduces to O(1/n), and any observed faster decay after shocks is unexplained by the theory.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a Cumulative Distribution Transform (CDT) based reduced-order modeling pipeline for one-dimensional scalar conservative PDEs with nonnegative, equal-mass states. After recalling the CDT and its isometry between the 1D Wasserstein metric and weighted L2(dr), the authors show analytically that linear transport becomes an affine segment in CDT space with zero 2-width, and they prove three families of estimates: (i) for entropy solutions of scalar conservation laws, a robust O(1/n) Kolmogorov-width bound in L2(dr) via the Wasserstein dynamic formulation (Theorem 4.2); (ii) a sharper O(1/n^2) pre-shock bound under the characteristic non-crossing condition 1+t f''(u0)u0' >= alpha > 0 (Theorem 4.5); and (iii) for advection-diffusion, d2 <= sqrt(2DT) (Theorem 5.1) plus Taylor-remainder based O(D^2T^2) bounds (Theorem 5.6). The proposed CDT-POD method is tested on Burgers, traffic flow, Buckley-Leverett, and advection-diffusion, with code provided. The paper is explicit that the sharper hyperbolic estimate is pre-shock and that a post-shock counterpart remains open.","tokens_in":45545,"tokens_out":23045,"duration_ms":208309,"significance":"If correct, the results provide a rigorous explanation for the empirical success of transport-based coordinates in ROM: the exact low-dimensional structure for linear transport, the robust O(1/n) bound for entropy solutions, and the quantification of diffusion's departure from the rank-two transport plane are clean and original contributions. The proofs are largely self-contained and use standard tools correctly; Theorem 4.2 is a particularly nice application of the Benamou-Brenier dynamic formulation with v=f(u)/u, and Theorem 5.1 is a simple Brownian-coupling argument. The paper is honest about its limitations: the O(n^{-2}) hyperbolic bound is conditional on a pre-shock hypothesis, the post-shock sharp estimate is explicitly left open (Section 6), and the inverse CDT requires monotonicity, which linear truncation may destroy. Reproducible code and explicit transform/inverse-transform error floors strengthen the numerical claims. The main caveat is that the theoretical width bounds are in the Wasserstein/L2(dr) metric rather than native L2, so the theoretical comparison with native POD is not a direct L2 statement; the authors are careful about this in the main text.","major_comments":[{"comment":"The advertised O(n^{-2}) hyperbolic width is conditional on 1+t f''(u0(ξ))u0'(ξ) >= alpha > 0 on [0,T]. For generic C^1 data with min u0' < 0 and f'' not identically zero on the range of u0, this quantity reaches zero at the first shock time, so the estimate covers only pre-shock times; after shocks, the only proven bound is the O(1/n) one in Theorem 4.2. Section 6 acknowledges this gap. Because the numerical examples in §4.2 (e.g., Burgers in Fig. 7) explicitly include post-shock times, the statement that the experiments 'support the theoretical picture' should be qualified: the observed faster decay after shocks is not covered by Theorem 4.5. I recommend making the pre-shock condition and the open post-shock question prominent in the abstract and contributions.","section":"§4.1, Theorem 4.5"},{"comment":"The phrase 'sharper O(D^2T^2) estimates under additional regularity or away from initial layers' in the abstract can be misread as an asymptotic statement in D. For t0>0, the constant C(u(·,t0)) in (60) is bounded in Remark 5.5 by a quantity proportional to (D t0)^{-3}, so for fixed t0 the right-hand side of (60) grows as D→0; only the t0=0 smooth-data case (Remark 5.7) yields D→0 recovery. Please state clearly in Theorem 5.6 that the D^2 factor comes with a D-dependent constant for t0>0, and that the zero-width recovery as D→0 is restricted to the t0=0 case.","section":"§5.1, Theorem 5.6 and Remark 5.5"}],"minor_comments":[{"comment":"The numerical experiments do not specify which monotone projection or rearrangement is applied when the truncated CDT map is not monotone. Please document this, and state the reference density r used in each experiment, for reproducibility of the reported reconstruction errors.","section":"§2.5, Algorithm 1"},{"comment":"The display bounding d2 by ||pu(t)-pu(t0)-...|| should include a supremum over t in [t0,T]; as written it appears to be missing the sup.","section":"§5.1, proof of Theorem 5.6"},{"comment":"For the Burgers initial condition u0 = A g(x)+0.1, please state explicitly how the total mass is normalized and how the CDT reference mass is matched; this affects the interpretation of the transformed snapshots.","section":"§4.2, Eq. (50)"},{"comment":"The labels 'formal' and the statement that the native-space estimates are not exact widths are welcome. In Section 6, avoid wording that suggests the native-space comparison is a theorem; the current text is mostly careful but could be tightened.","section":"§4.1, formal native-space comparison"}],"recommendation":"minor_revision","confidential_remarks":"The paper is heavily self-referential (CDT foundation papers, PyTransKit), but the Kolmogorov-width estimates are new and the central derivations appear correct. The scope is modest — 1D, nonnegative densities — but the results are solid and clearly presented. The main risk is that readers over-interpret the pre-shock O(n^{-2}) bound as applying to shock-dominated regimes; the revision should make the limitations impossible to miss. No concerns about novelty or authorship integrity beyond the usual expectation that the CDT-related citations are appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper gives real new estimates: O(1/n) Kolmogorov width in Wasserstein/L2(dr) for entropy solutions, O(n^-2) for pre-shock smooth regimes, sqrt(DT) and D^2T^2 bounds for advection-diffusion, plus a transform-reduce-invert CDT-POD pipeline with code. Second, the sharpest hyperbolic rate is limited by an explicit pre-shock hypothesis; for generic data with decreasing slope, characteristics cross in finite time and only O(1/n) is proved afterwards. That is a real scope limit, not an internal inconsistency.\n\nWhat is new: the zero 2-width for linear transport follows from known translation equivariance, but the entropy-solution width bound and the advection-diffusion deviation estimates are not in the cited literature. The proofs are mostly clean: Theorem 4.2 uses the Wasserstein dynamic formulation and Kruzhkov theory correctly; Theorem 5.1 is a simple Brownian coupling; Theorem 5.6 follows from Taylor remainder and Jensen-type convolution lemma. The paper does not hide the shock limitation; Section 6 explicitly says a higher-regularity post-shock counterpart is open. It also flags the monotonicity issue in linear CDT truncation, which matters for the inverse transform. Self-citations are heavy, but the estimates do not assume those cited results in a circular way; the cited CDT isometry and composition properties are standard and needed.\n\nSoft spots: the numerical section is broad but I could not independently reproduce the figures from the text alone; the code link helps, but constant tracking in Corollary 5.4 is loose and the C(u0) condition is strong. The bigger concern is proportionality: the abstract sells sharper O(n^-2) bounds without the pre-shock qualifier in the first summary, so a reader could overestimate the post-shock theory. That should be fixed by making the limitation more prominent. It does not invalidate the main theorem.\n\nWho should read: people doing reduced-order modeling for transport-dominated 1D problems, and anyone wanting a concrete example of transport coordinates helping with the Kolmogorov barrier. It deserves a serious referee; I would send it out, with a request to tighten the presentation of the pre-shock scope and to provide a reproducibility README or script for the numerics. My verdict would be conditional rather than unconditional, but it is a useful paper.","headline":"Solid transport-coordinate ROM paper; the advertised O(n^-2) shrink applies only before shock formation, and the paper says so, but that should not stop a serious referee.","tokens_in":46047,"tokens_out":1492,"would_cite":true,"duration_ms":15010,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["41A46","49Q22","65F55","35L65","35Q49"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that representing 1D conservative-PDE solutions by mass coordinates rather than spatial coordinates collapses transport-dominated solution manifolds into small linear subspaces, with rigorous Kolmogorov-width bounds and a","keywords":["cumulative distribution transform","optimal transport","Wasserstein distance","Kolmogorov width","reduced order modeling","proper orthogonal decomposition","scalar conservation laws","advection-diffusion"],"falsifier":"Take inviscid Burgers with a smooth initial profile that has negative slope somewhere, choose T before the shock time, compute the CDT snapshots of the exact solution, and measure the singular values of the transformed snapshot matrix; if they do not decay at least like (n-1)^{-2}, Theorem 4.5 is violated, and if after shock formation they decay exactly like n^{-1}, that confirms the general entropy-solution bound is the true post-shock rate.","tokens_in":45187,"feed_emoji":"🌊","tokens_out":7147,"duration_ms":66486,"temperature":0.7,"pith_summary":"This paper aims to show that one-dimensional conservative PDEs, whose solutions are moving, nonnegative, mass-preserving densities, become far more compressible if you change coordinates before applying linear model reduction. The change of coordinates is the cumulative distribution transform (CDT), which labels each solution by the position of accumulated mass rather than by the spatial coordinate; in these coordinates a translation becomes an affine shift, and the Wasserstein distance becomes an ordinary weighted L2 distance. The authors prove that for linear transport the transformed solution manifold lies in a two-dimensional plane, that for general hyperbolic conservation laws its Kolmogorov n-width decays at worst like O(1/n) even after shocks and like O(1/n^2) in smooth pre-shock regimes, and that for advection-diffusion the transformed trajectory stays within O(sqrt(DT)) of that plane. On top of this, they build a CDT-POD scheme that maps snapshots to transform space, performs POD there, and maps back, and they show numerically that it needs substantially fewer modes than Eulerian POD for transport-dominated examples. A sympathetic reader would care because the failure of linear reduced models on transported fronts is a known bottleneck, and this paper gives both a mechanism for that failure and a concrete fix with quantitative guarantees.","feed_headline":"Tracking mass, not position, shrinks moving-wave data to few modes","feed_subtitle":"Moving fronts in conservative PDEs become low-dimensional in mass coordinates, with proven error bounds for fewer modes.","key_machinery":"The central object is the cumulative distribution transform (CDT): the monotone optimal transport map hat-u from a reference density r to a target density u, defined implicitly by the mass-balance identity integral_{-inf}^{hat-u(xi)} u(x) dx = integral_{-inf}^{xi} r(s) ds, i.e., hat-u = F_u^{-1} o F_r. Three properties carry the argument: the CDT isometrically embeds the 1D Wasserstein space into the weighted Hilbert space L^2(dr); it converts translations into additive shifts (hat-{T_tau u} = hat-u + tau), so transported families become affine lines; and its mass labels are globally defined even after shocks, unlike classical characteristics. The width estimates then follow from a Hilbert-v","core_discovery":"The central claim is that the CDT is the right coordinate system for linear compression of conservative-transport data. In the CDT, a density u is represented by the monotone map hat-u that pushes a fixed reference density onto u; the transform is an isometry from the 2-Wasserstein space into L^2(dr). For linear advection u_t + A u_x = 0, hat-u(t) = hat-u_0 + A t, so the solution manifold lies in span{hat-u_0, 1} and d_n = 0 for all n >= 2. For a scalar conservation law u_t + f(u)_x = 0, entropy solutions satisfy W2(u(t),u(s)) <= B_M |t-s|, giving d_n(hat-M; L^2(dr)) <= B_M T/(2n); under the pre-shock condition 1 + t f''(u0) u0' >= alpha > 0, the sharper bound d_n <= C T^2 / (8 sqrt(alpha) (","pith_inferences":["If these bounds hold for time evolution, the same transform-reduce-invert recipe should extend to parametric families of conservative densities, where one varies the initial condition or coefficients and expects similar compression for any snapshot set dominated by mass displacement.","Because the linear truncation can leave the monotone image of the CDT, an editor-level extension is that projecting the reduced transport maps onto the cone of nondecreasing functions should convert the proven transform-space n-width bounds into native-space reconstruction guarantees; this is a testable algorithmic modification the paper leaves open.","The robust O(1/n) bound after shocks suggests that post-shock entropy-solution manifolds are still much more compressible in Wasserstein coordinates than in native L2, where translated jumps have n-width decaying only like n^{-1/2}; a direct numerical comparison of native versus CDT widths at fixed n after shock would quantify the practical gain.","The construction is essentially one-dimensional because it relies on monotone transport maps; an inference is that higher-dimensional analogues would need slicing (for example, a Radon-CDT) or regularized transport maps, and the isometric width bounds should not be expected to transfer directly."],"forward_implications":["For linear advection, a rank-2 POD in CDT space is exact, so purely transported data can be represented by exactly two modes.","For entropy solutions of scalar conservation laws, n modes in CDT space guarantee a worst-case Wasserstein/weighted-L2 error of at most B_M T/(2n), even after shock formation.","In smooth pre-shock regimes, the guarantee improves to O(n^{-2}) with a constant that degrades only like alpha^{-1/2} as the minimum characteristic Jacobian tends to zero.","For conservative advection-diffusion, n=2 modes in CDT space stay within sqrt(2DT) of the exact trajectory, recovering the pure-transport d_2 = 0 as D approaches 0.","Numerically, CDT-POD gives faster singular-value decay and lower reconstruction error than Eulerian POD for transport-dominated dynamics, provided the reduced transport maps remain monotone so the inverse CDT is valid."],"fun_headline_variants":["Mass coords make moving-wave data low-rank","CDT compresses transport PDEs to few modes","Track mass, not waves, to shrink basis size","Transform shows why transport is low-dimensional","Fewer modes with proven bounds via mass transform"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The sharp O(n^-2) hyperbolic bound assumes no shock forms by time T (1 + t f''(u0) u0' >= alpha > 0); the paper itself flags that a higher-regularity post-shock analogue remains open, and the numerical pipeline further assumes the truncated CDT maps stay monotone so the inverse transform is valid.","fun_headline_variants_meta":{"raw":{"variants":["Mass coords make moving-wave data low-rank","CDT compresses transport PDEs to few modes","Track mass, not waves, to shrink basis size","Transform shows why transport is low-dimensional","Fewer modes with proven bounds via mass transform"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000355,"raw_usage":{"total_tokens":1861,"prompt_tokens":938,"completion_tokens":923,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":682,"completion_tokens_details":{"reasoning_tokens":852}},"tokens_in":682,"tokens_out":923,"duration_ms":9009,"temperature":1.0,"reasoning_tokens":852,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T19:08:12.386974+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take inviscid Burgers with a smooth initial profile that has negative slope somewhere, choose T before the shock time, compute the CDT snapshots of the exact solution, and measure the singular values of the transformed snapshot matrix; if they do not decay at least like (n-1)^{-2}, Theorem 4.5 is violated, and if after shock formation they decay exactly like n^{-1}, that confirms the general entropy-solution bound is the true post-shock rate.","supporting_citations":[],"review_version":1}