{"id":"26c17941-a0a1-47a3-996d-949b3ea7c01d","arxiv_id":"2411.14891","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"New self-similar solutions for supersonic non-equilibrium gray radiation diffusion are derived for power-law density profiles, enabling benchmarks for non-homogeneous media with beta not equal to 4 material energies.","lead":"Physicists derived a new family of analytic similarity solutions for non-equilibrium Marshak waves, the radiation heat fronts driven into cold material in high-energy-density experiments, now extended to media with non-uniform density and general material energy laws. The solutions provide ready-made benchmarks for verifying radiation-transport and diffusion codes in regimes where radiation and matter are not in equilibrium.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Special-case formulas for μ=1 (Eqs. 38-39) contradict the general expressions (17)-(18) and (24); substitution gives ω=2/(2+λ+λ′) and δ=(2+λ+λ′)/(2+2λ′), not the printed values.","rationale":"The reader's weakest assumption concerned the circularity of the transport benchmark: the bath temperature is derived from the diffusion solution itself, weakening the claim of independent transport validation. That is a valid limitation of the verification strategy, but it does not affect the mathematical existence of the similarity solutions; the gray diffusion simulations still provide an independent check of the ODE reduction. The concern I identify is more load-bearing because it is a direct internal inconsistency in the paper's own formulas. For μ=1, the paper states exponents that do not follow from the general equations, which undermines the correctness of a specific claimed consequence of the central result and could mislead a reader implementing the μ=1 case. The general derivation for μ<1 appears sound, and the main claim of new similarity solutions for β≠4 and ω≠0 stands. Therefore the appropriate verdict remains conditional: the paper needs a correction to the μ=1 special-case formulas, and the reader's concerns about quantitative error metrics and data availability remain. I disagree with the reader that the transport boundary circularity is the weakest point; the algebraic error is more concrete and more directly tied to the paper's analytical claims.","tokens_in":34485,"tokens_out":24096,"duration_ms":206004,"concrete_test":"Symbolically substitute μ=1 into Eqs. (17)-(18) and (24) and simplify. The result should be τ=0, ω=2/(2+λ+λ′), δ=(2+λ+λ′)/(2+2λ′). If the printed Eqs. (38)-(39) do not match this output—and cannot be recovered under an additional constraint such as λ=λ′—then the special-case formulas in the manuscript are wrong and must be corrected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The manuscript's stated consequence for a density-independent specific energy (μ=1) is internally inconsistent with the general derivation. Setting μ=1 in Eq. (18) makes the term (α+α′)(1−μ) vanish, yielding ω = 2(β−4)/[(β−4)(2+λ+λ′)] = 2/(2+λ+λ′), not ω=1/(1+λ) as printed in Eq. (38). Likewise, Eq. (24) with τ=0 gives δ = 1/[2−ω(1+λ)] = (2+λ+λ′)/[2(1+λ′)], not δ=(1+λ)/(1+λ′) as printed in Eq. (39). The printed δ does not follow from the printed ω either. This is a concrete algebraic error in a stated consequence of the central derivation. It matters because the paper's range-of-validity discussion and parameter classification explicitly include the μ=1 case; a user constructing a benchmark for a density-independent material energy would apply the wrong density profile exponent and similarity exponent. The general formulas for μ<1 are unaffected, but the special-case claims require correction.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a family of self-similar solutions for the non-equilibrium, supersonic Marshak wave problem in the gray diffusion limit, for a planar non-homogeneous medium with density profile ρ(x)=ρ0 x^{-ω}, surface temperature drive T_s(t)=T0 t^τ, and power-law material opacities and specific energy. Using dimensional analysis (Appendix A), the authors show that self-similarity fixes τ and ω through the material exponents, reducing the partial differential equations (12)-(13) to the ODE system (27)-(28). They classify solutions by their near-origin behavior, derive the equivalent Marshak bath-temperature boundary condition (48), and construct six benchmark problems. Numerical solutions of the similarity ODEs are tabulated, and gray diffusion, implicit Monte-Carlo, and discrete-ordinates transport simulations are compared with the profiles in optically thick cases.","tokens_in":34728,"tokens_out":13197,"duration_ms":120153,"significance":"If correct, this is a valuable extension of the homogeneous-medium nonlinear non-equilibrium Marshak solutions, allowing self-similar benchmarks for materials with β≠4 through the inhomogeneous density profile. The construction is first-principles dimensional analysis with no fitted parameters in the similarity profiles, and the tabulated numerical profiles plus approximate analytic fits are a practical resource for code verification. The transport benchmarks are less independent than a fully external test, because the bath drive is taken from the diffusion solution (Eq. 48) and the scattering opacity k_s=k_t-k_a is an effective, nonphysical material law; however, these points are acknowledged in the text and do not affect the mathematical derivation of the diffusion-limit solutions.","major_comments":[{"comment":"The μ=1 special-case exponents are algebraically inconsistent with the general formulas. Setting μ=1 and β≠4 in Eqs. (17)-(18) and (24) gives τ=0, ω=2/(2+λ+λ′), and δ=(2+λ+λ′)/(2(1+λ′)); the printed values ω=1/(1+λ) and δ=(1+λ)/(1+λ′) are recovered only under the additional assumption λ=λ′. Because the text presents these as the general β-independent exponents for a density-independent material energy density, the formulas should be corrected or the λ=λ′ restriction should be stated explicitly. Without this fix, a user constructing a benchmark for μ=1 would use incorrect density-profile and front-propagation exponents.","section":"Sec. III A, Eqs. (38)-(39)"},{"comment":"The verification claim rests on visual agreement: the comparisons with gray diffusion, IMC, and SN transport simulations are shown only as overlaid curves and described as \"great agreement,\" without numerical error measures. For a paper whose central deliverable is a set of code-verification benchmarks, quantitative agreement should be reported, for example the maximum and mean relative deviation of T_r and T from the analytic profiles for each test and time. This would also make the benchmarks reproducible and would allow other codes to be checked against the same tolerances.","section":"Sec. IV B, Figs. 13-24 and Sec. V"}],"minor_comments":[{"comment":"The caption quotes the material exponents as α=3.9 and α′=1.9, while Figs. 4 and 7, which are cited as the same case, use α=3.5 and α′=1.9; one of these is a typo.","section":"Fig. 8 caption"},{"comment":"The stated value β′_c=2.69565 does not follow from α=3, λ=0.2, μ=0.4; Eq. (35) gives β′_c=2.5. The mistake does not affect the computed solution but should be corrected.","section":"Sec. IV B, Test 2"},{"comment":"The sentence \"there is no solution with τ=0\" is stated under the preceding assumption μ<1, but the transition to the μ=1 case in Eqs. (37)-(39) is abrupt; the domain of each statement should be made explicit to avoid confusion.","section":"Sec. III A"},{"comment":"The effective scattering opacity k_s=k_t-k_a is nonphysical and can become negative; although the paper acknowledges this, the benchmark section should state more prominently that the transport tests are consistency checks of the transport solver against the diffusion-derived boundary condition for a synthetic material, not independent physics validation.","section":"Sec. IV A, Eq. (54)"}],"recommendation":"major_revision","confidential_remarks":"The main obstacle is the incorrect μ=1 special-case formulas in Eqs. (38)-(39); the general derivation appears sound and the benchmarks are potentially useful. If the authors correct the algebraic error and add quantitative simulation-vs-analytic errors, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth a look if you work on verification of radiation-hydrodynamics codes. It generalizes the non-equilibrium Marshak wave similarity solutions from homogeneous media (beta=4) to non-homogeneous density profiles with power-law density rho~x^{-omega}, which allows self-similar solutions for arbitrary beta in the material energy density u~T^beta. The key result is the pair of exponent relations (17)-(18), which show that beta!=4 requires omega!=0. That is a clean and non-trivial extension, and it opens up benchmarks for realistic materials with beta!=4. The derivation via dimensional analysis is careful, the solution classification (omega<0, omega=0, omega>0) is physically sensible, and the six test cases with tabulated profiles and fitted approximations are a genuine contribution.\n\nThe main flaw I see is in the special-case discussion for mu=1. Equations (37)-(39) state that for mu=1, tau=0, omega=1/(1+lambda), delta=(1+lambda)/(1+lambda'). Substituting mu=1 into the general formulas (17)-(18) and (24) gives tau=0, omega=2/(2+lambda+lambda'), delta=(2+lambda+lambda')/[2(1+lambda')]. The printed omega equals 2/(2+lambda+lambda') only if lambda=lambda', and the printed delta is inconsistent with both. This is a concrete algebraic error in a stated consequence of the central derivation, and it sits in the range-of-validity section where a user would look for guidance on a density-independent material energy. The general formulas are unaffected, and this does not undermine the benchmarks, but it should be corrected.\n\nTwo smaller issues. The transport benchmarks use a scattering opacity ks=kt-ka which can be negative, and the bath temperature is taken from the diffusion solution itself. The paper states this explicitly, so it is not a hidden assumption, but it means the comparison is a consistency check between transport and diffusion in the optically thick limit, not an independent validation of the diffusion solution. Some quantitative error measures (e.g., L1 or L_infinity errors at the final time) would make the claimed 'great agreement' more credible than the visual plots.\n\nWho is this for? People building or maintaining gray transport/diffusion codes who need non-trivial verification problems with non-equilibrium and non-homogeneous density. It deserves a proper review, and I would recommend acceptance after the mu=1 error is fixed and the comparisons are quantified. The core derivation and the benchmark suite are solid.","headline":"A solid benchmark paper with a genuinely new generalization; the mu=1 special-case formulas in Eqs. (38)-(39) contain an algebra error that needs fixing.","tokens_in":35315,"tokens_out":3992,"would_cite":true,"duration_ms":37071,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives self-similar solutions of the non-equilibrium supersonic Marshak wave problem in a planar medium with a power-law density profile, showing that such solutions exist exactly when the surface-temperature exponent and the…","keywords":["Marshak wave","non-equilibrium radiative transfer","self-similar solutions","gray diffusion approximation","radiation hydrodynamics","code verification benchmark","power-law material model","non-homogeneous media"],"falsifier":"One concrete test: rerun any of the six benchmarks with the same total and absorption opacities but with the scattering opacity replaced by a different positive function that still sums to the same total opacity; if the transport results no longer match the self-similar profiles, the agreement is an artifact of the $k_s = k_t - k_a$ construction. Another: for a material with $\\beta \\neq 4$, set $\\omega = 0$ and attempt the shooting solution of Eqs. (27)-(28); the paper predicts no self-similar solution, so a successful find would refute the claimed if-and-only-if condition.","tokens_in":34300,"feed_emoji":"🔥","tokens_out":11837,"duration_ms":108210,"temperature":0.7,"pith_summary":"This paper derives a family of self-similar solutions for a non-equilibrium Marshak wave—a radiation heat wave that travels through a cold material—in a planar medium whose density follows a spatial power law and whose surface is driven by a radiation temperature that follows a temporal power law. The solutions exist only when the surface-temperature exponent $\\tau$ and the density exponent $\\omega$ are locked to the material's opacity and heat-capacity exponents by two algebraic relations, Eqs. (17)-(18). The key gain over prior work is that the material energy density can now have a general power-law temperature dependence; for homogeneous media self-similarity required the material energy to scale as $T^4$. The paper uses the solutions to define six optically thick verification benchmarks and reports agreement with gray-diffusion, implicit Monte-Carlo, and discrete-ordinates transport simulations.","feed_headline":"Self-similar Marshak waves now cover non-uniform media","feed_subtitle":"Radiation heat waves in power-law density plasmas now have analytic benchmarks for verifying simulation codes.","key_machinery":"The engine of the argument is the self-similar ansatz of Eqs. (21)-(23), combined with Buckingham-Pi dimensional analysis. The ansatz collapses $E$ and $U$ into dimensionless profiles $f(\\xi)$ and $g(\\xi)$ with similarity coordinate $\\xi = x / [t^{\\delta}(K E_0^{\\alpha/4})^{1/(2-\\omega(1+\\lambda))}]$. Imposing that the two dimensionless constants $A$ and $B$ (Eqs. (19)-(20)) are independent of space and time yields the consistency conditions (17)-(18) that determine $\\tau$ and $\\omega$. Once those hold, the problem reduces to the ODE system (27)-(28), whose front coordinate $\\xi_0$ is fixed by shooting until $f(0)=1$. The same machinery yields the Marshak boundary condition, with a bath temperature involving the dimensionless surface flux $S(0)$, and the front law $x_F(t)\\propto t^{\\delta}$.","core_discovery":"The central claim is that the two-temperature gray-diffusion equations (12)-(13), with total and absorption opacities $k_t \\propto T^{-\\alpha}\\rho^{1+\\lambda}$, $k_a \\propto T^{-\\alpha'}\\rho^{1+\\lambda'}$, material energy $u \\propto T^{\\beta}\\rho^{1-\\mu}$, density $\\rho=\\rho_0 x^{-\\omega}$ and surface drive $T_r(0,t)=T_0 t^{\\tau}$, admit self-similar solutions $E=E_0 t^{4\\tau} f(\\xi)$, $U=E_0 t^{4\\tau} g(\\xi)$ if and only if $\\tau$ and $\\omega$ take the values prescribed by Eqs. (17)-(18). When these hold, all dimensional factors cancel and the PDE system reduces to the two ordinary differential equations (27)-(28) for the similarity profiles, with $f(0)=1$ and a finite heat front at $\\xi=\\xi_0$. Self-similarity for $\\beta\\neq 4$ is possible only because $\\omega\\neq 0$: a homogeneous medium forces $\\beta=4$. The paper also derives the equivalent Marshak boundary condition, giving a time-dependent bath temperature $T_{bath}(t)$, and shows that the resulting profiles, integrated numerically and optionally fitted to closed-form approximations, are reproduced by gray-diffusion simulations and by stochastic and deterministic transport codes in the optically thick limit.","pith_inferences":["We infer that the derived bath temperature $T_{bath}(t)$ is not a pure power law whenever $\\delta\\neq 1$ (Eq. (48)), so driving a Marshak wave with a heat bath instead of a surface temperature requires a non-power-law temporal modulation to stay self-similar; this could be tested in a simulation that prescribes $T_{bath}$ rather than $T_s$.","We infer that because the transport agreement uses the constructed scattering opacity $k_t-k_a$, extending the benchmarks to physical, positive scattering opacities, or to multi-group transport where the gray-limit equivalence breaks, would test whether the similarity solutions survive outside the exact gray-diffusion setting.","We infer that for real materials with $\\beta<4$, the self-similarity condition $\\omega\\neq 0$ suggests a design rule for experiments: choose a density profile with the $\\omega$ that satisfies Eq. (18) to realize a self-similar non-equilibrium wave in the laboratory, rather than a homogeneous target.","We infer that the family may also serve as a test bed for inverse problems: because $A$, $B$, the front position, and the profile shape all depend on material exponents, matching a measured $x_F(t)$ to the self-similar form could in principle constrain an unknown power-law opacity model—an inversion the paper does not take up."],"forward_implications":["Wherever the material exponents satisfy the validity conditions, the heat front advances as a power law $x_F(t) = \\xi_0 (K E_0^{\\alpha/4})^{1/(2-\\omega(1+\\lambda))} t^{\\delta}$, with $\\delta$ from Eq. (24), so the wave accelerates, decelerates, or travels at constant speed depending on the material.","The solution family exhibits three distinct near-origin behaviors: material temperature locked to zero when $\\omega<0$, finite but below the radiation temperature when $\\omega=0$, and equal to the radiation temperature when $\\omega>0$; these are directly testable signatures in simulations.","The six benchmarks, with tabulated profiles and fitted analytic forms accurate to about 0.5%, give code developers a non-trivial target for verifying non-equilibrium radiation diffusion and transport codes in the optically thick limit.","In the large-$B$ (strong coupling) limit the radiation and material temperatures approach equality, $g\\approx f$, so the family contains its own equilibrium-limit consistency check, while test 6 shows equilibrium can fail in a low-density inner region even when $B$ is large."],"supporting_citations":[{"why":"The homogeneous-media ($\\omega=0$, $\\beta=4$) solution family that this paper generalizes to inhomogeneous media.","marker":"[27]"},{"why":"Defines the non-equilibrium Marshak wave diffusion problem whose linear solution this work extends to nonlinear power-law opacities.","marker":"[47]"},{"why":"Supplies the standard non-equilibrium Marshak benchmark format and the linear benchmark results that the new benchmarks mirror.","marker":"[48]"},{"why":"Introduces the Marshak wave phenomenon and the propagating heat-front boundary-value problem.","marker":"[30]"},{"why":"Provides the unified self-similar supersonic Marshak wave theory whose exponents and validity analysis are extended here.","marker":"[28]"},{"why":"Demonstrates dimensional-analysis self-similar solutions for nonlinear radiation diffusion in non-homogeneous media, the technical template for the new derivation.","marker":"[23]"},{"why":"The Buckingham Pi theorem used in Appendix A to reduce the dimensional problem to similarity variables.","marker":"[71]"},{"why":"Source of the power-law material exponents (Table I) used to show that the listed real materials do not satisfy the self-similarity condition.","marker":"[5]"}],"fun_headline_variants":["Marshak wave solutions extend to non-uniform media","Similarity solutions for Marshak waves in non-homogeneous matter","Benchmarking code: Marshak waves in power-law plasmas","Non-equilibrium Marshak waves solved for inhomogeneous media","Marshak wave benchmarks for non-uniform density plasmas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the transport verification is a genuine independent check, even though its scattering opacity is defined as the difference between the model's total and absorption opacities, a quantity that need not be a physical non-negative opacity, and its boundary drive is computed from the very diffusion solution being tested.","fun_headline_variants_meta":{"raw":{"variants":["Marshak wave solutions extend to non-uniform media","Similarity solutions for Marshak waves in non-homogeneous matter","Benchmarking code: Marshak waves in power-law plasmas","Non-equilibrium Marshak waves solved for inhomogeneous media","Marshak wave benchmarks for non-uniform density plasmas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000233,"raw_usage":{"total_tokens":1555,"prompt_tokens":1071,"completion_tokens":484,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":687,"completion_tokens_details":{"reasoning_tokens":400}},"tokens_in":687,"tokens_out":484,"duration_ms":19405,"temperature":1.0,"reasoning_tokens":400,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:45:09.986684+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One concrete test: rerun any of the six benchmarks with the same total and absorption opacities but with the scattering opacity replaced by a different positive function that still sums to the same total opacity; if the transport results no longer match the self-similar profiles, the agreement is an artifact of the $k_s = k_t - k_a$ construction. Another: for a material with $\\beta \\neq 4$, set $\\omega = 0$ and attempt the shooting solution of Eqs. (27)-(28); the paper predicts no self-similar solution, so a successful find would refute the claimed if-and-only-if condition.","supporting_citations":[{"cited_title":"& McClarren, R","cited_arxiv_id":null,"evidence_quote":"The homogeneous-media ($\\omega=0$, $\\beta=4$) solution family that this paper generalizes to inhomogeneous media."},{"cited_title":"The non-equilibrium marshak wave prob- lem","cited_arxiv_id":null,"evidence_quote":"Defines the non-equilibrium Marshak wave diffusion problem whose linear solution this work extends to nonlinear power-law opacities."},{"cited_title":"& Olson, G","cited_arxiv_id":null,"evidence_quote":"Supplies the standard non-equilibrium Marshak benchmark format and the linear benchmark results that the new benchmarks mirror."},{"cited_title":"Effect of radiation on shock wave behavior","cited_arxiv_id":null,"evidence_quote":"Introduces the Marshak wave phenomenon and the propagating heat-front boundary-value problem."},{"cited_title":"& McClarren, R","cited_arxiv_id":null,"evidence_quote":"Provides the unified self-similar supersonic Marshak wave theory whose exponents and validity analysis are extended here."},{"cited_title":"Analytic solutions of the nonlinear radiation diffusion equation with an instantaneous point source in non-homogeneous media","cited_arxiv_id":null,"evidence_quote":"Demonstrates dimensional-analysis self-similar solutions for nonlinear radiation diffusion in non-homogeneous media, the technical template for the new derivation."},{"cited_title":"On physically similar systems; illustra- tions of the use of dimensional equations.Physical review 4, 345 (1914)","cited_arxiv_id":null,"evidence_quote":"The Buckingham Pi theorem used in Appendix A to reduce the dimensional problem to similarity variables."}],"review_version":1}