{"id":"73e351bf-2e3e-4552-b07f-e51d9004f517","arxiv_id":"2607.27440","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Linearizing Boltzmann–BGK about kinetic shock base flows shows non-equilibrium VDFs shift 1D normal-shock spectra toward less stable branches than Maxwellian closures, enabled by a large-scale SLEPc/PETSc solver.","lead":"The authors build a kinetic linear-stability method that linearizes the Boltzmann–BGK equation on non-equilibrium shock profiles, keeping bi-modal velocity distributions inside the shock. At Mach 3–4 those kinetic effects push the spectrum toward less damping than Maxwellian-based analyses, and a parallel eigensolver reaches systems with ~280k unknowns.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"EQ vs NONEQ comparison may not isolate pure VDF non-equilibrium because the linearized collision coefficients still depend on how base g_c, h_c enter the operator.","rationale":"The reader correctly flags BGK/Pr=1 and the lack of an independent high-Mach benchmark, and the CONDITIONAL verdict is appropriate: the numerical framework, Couette verification, mesh studies, and HPC path are solid, and no unstable modes are claimed. The stronger load-bearing issue for the physics headline, however, is internal to the EQ/NONEQ design. The paper’s own linearization (Eqs. 17–21) makes the residual (g^e_c−g_c) terms part of the stability operator; nulling them when constructing EQ is not the same as “same macros, Maxwellian VDF only.” A hybrid spectrum would settle whether bi-modality in the base VDF, as opposed to the change in linearized collision coefficients, drives the shift that continuum closures are said to miss. Until that check (or an equivalent ES-BGK / full Boltzmann comparison) is done, the strongest claim should remain conditional on a clearer attribution, not only on Prandtl-number fidelity. This does not overturn the methods contribution or require rejection; it tightens the same CONDITIONAL posture the reader already reached.","tokens_in":40167,"tokens_out":868,"duration_ms":15369,"concrete_test":"Recompute the M∞=3, β=16 spectrum (Fig. 10c) with a hybrid operator: keep non-equilibrium g_c, h_c only in the streaming/ξ·∇ terms and in the moment recovery of macro perturbations, but evaluate all G_*, H_* and A_g1/A_h1 from the local Maxwellian built on the same macros (i.e., force the (g^e_c−g_c) residual blocks to zero). If the least-stable continuous branch moves back to within ~5–10% of the pure EQ branch (in Im(ω)), the published NONEQ shift is largely a collision-linearization artifact rather than bi-modal streaming; if it stays near the published NONEQ branch, the isolation holds.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that non-equilibrium base VDFs systematically shift shock eigenspectra toward less stable regions relative to Maxwellian-enforced spectra built from the same macroscopic base state (§V, Figs. 10–13; Abstract). That isolation rests on the EQ construction: force g_c = g^e_c and h_c = h^e_c inside the kLST operator while keeping macroscopic ρ_c, u_c, T_c from the BE-BGK base flow. In the linearized BGK system (Eqs. 17–21 and the assembled blocks in Eq. 23), the coefficients G_ρ, G_u, … and A_g1 / A_h1 contain both equilibrium derivatives and explicit (g^e_c − g_c) residual terms weighted by pressure/viscosity. Setting g_c = g^e_c therefore simultaneously (i) removes bi-modality from the base state that multiplies the streaming/collision terms and (ii) nulls those residual linearization pieces. The reported spectral shift therefore conflates genuine bi-modal streaming effects with a change in the linearized collision operator. The paper never shows a third spectrum that retains the non-equilibrium base VDF only in the streaming operator while restoring Maxwellian collision linearization (or the converse). Without that split, the claim that “kinetic effects” / “translational non-equilibrium” are what continuum closures miss is only partially supported; part of the shift could be an artifact of how EQ is defined inside BGK linearization. BGK/Pr=1 (the reader’s weakest assumption) is a real limitation but secondary: even inside pure BGK the EQ/NONEQ contrast is not cleanly attributed.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.5","summary":"The manuscript formulates a kinetic linear stability theory (kLST) by linearizing the Boltzmann–BGK equation in reduced distribution functions about BE-BGK base flows for isolated 1D argon normal shocks, so the stability operator acts on the VDF rather than a closed continuum system. After verifying the discrete operator against compressible Couette eigenvalues of Zou et al., the authors compute spectra at M∞=1.2, 3.0, and 4.0, comparing constructions that enforce local Maxwellian base VDFs (EQ) against those that retain the non-equilibrium BE-BGK VDFs (NONEQ). At higher Mach they report a systematic shift of continuous branches toward less stable regions when non-equilibrium VDFs are retained, and they support the large discrete systems (up to ~2.8×10^5 unknowns) with a parallel SLEPc/PETSc stack (shift-and-invert Arnoldi+MUMPS and Jacobi–Davidson with block-Jacobi ILU), including fill-in and scaling benchmarks on Frontera.","tokens_in":40549,"tokens_out":1405,"duration_ms":45037,"significance":"If the reported EQ/NONEQ spectral shifts hold under the stated operator definitions, the work is a genuine first step beyond NS/moment-based shock LST and beyond prior linearized-BGK studies confined to channel flows: it keeps bi-modal translational non-equilibrium inside both the base state and the perturbation problem for finite-thickness shocks, and it pushes kinetic modal analysis to M∞=4 with a documented HPC path. Strengths that should be credited include external Couette eigenvalue checks (Appendix A, ~1% on least-stable modes), documented base-flow and spectral mesh studies (Appendix B), an explicit shared-grid EQ vs NONEQ protocol, and quantitative solver/fill-in evidence rather than a methods sketch alone. The main scientific payload is the high-Mach spectral comparison; the solver infrastructure is enabling and carefully reported.","major_comments":[{"comment":"§II (Eqs. 17–21, 23) and §V (EQ vs NONEQ): The central claim attributes the high-M spectral shift to translational non-equilibrium / bi-modal base VDFs. In the BGK linearization, setting g_c=g^e_c and h_c=h^e_c simultaneously (i) replaces the base VDF in streaming and moment couplings and (ii) nulls the explicit residual pieces proportional to (g^e_c−g_c) and (h^e_c−h_c) that appear in Eq. 17 and feed the G_ρ, G_u, … / A_g1 blocks. The manuscript never states which of these channels dominates the shift in Figs. 10–13, nor that the EQ macroscopic field with forced Maxwellian VDFs is not itself a steady BE-BGK shock solution. Please add a short operator-level accounting of what EQ changes, and temper language that equates the EQ construction with “continuum” or pure “Maxwellian streaming” unless a controlled split (or an argument why a split is inconsistent) is given.","section":"§II, Eqs. 17–21; §V"},{"comment":"Abstract and §V: The abstract concludes that “continuum predictions can miss important changes even when the macroscopic profiles appear well resolved.” The high-M evidence is EQ vs NONEQ kinetic spectra (Figs. 10–13), not a side-by-side LNSE or G&P-based cLST spectrum at M∞=3 and 4 on the same Δ-based nondimensionalization. Duck & Balakumar are cited for low-M continuous branches; that does not by itself underwrite the continuum-miss claim at M=3–4. Either add a continuum LST comparison at least at one high-M case, or revise the abstract/conclusion wording to what is actually shown: within BGK-kLST, retaining non-equilibrium base VDFs shifts spectra relative to a Maxwellian-enforced kinetic operator built from the same macros.","section":"Abstract; §V"}],"minor_comments":[{"comment":"§VI already flags BGK/Pr=1; please mirror that limitation once in the abstract or early §I so readers do not over-read the argon spectra as quantitative real-gas predictions.","section":"Abstract / §I / §VI"},{"comment":"Notation: ϖ vs ω, A_g1 vs A1, and β as both spanwise wavenumber and G&P viscosity constant appear in close proximity; a single consistency pass would help.","section":"Nomenclature; §II"},{"comment":"Fig. 5 compares an incompressible NS extract (n=387) to a kLST extract for qualitative sparsity only; state explicitly in the caption that dimensions and physics differ so the panel is not read as a like-for-like cost comparison.","section":"§IV.C, Fig. 5"},{"comment":"Several encoding artifacts remain in the text (e.g., “shockâ=C™s”, “solverâ=C™s”); clean these in production.","section":"§I; §V"},{"comment":"Table 14 / length scale: using G&P Δ while base flows are BE-BGK is practical but should be restated when quoting Re_Δ and Kn_Δ next to BGK profiles so nondimensional spectra are unambiguous.","section":"§V, Table 14"},{"comment":"Appendix C solver tables are valuable; consider pointing to them from the M=4 JD run in §V so the 281,088-unknown result is reproducible from the main text alone.","section":"§V; Appendix C"}],"recommendation":"minor_revision","confidential_remarks":"I view this as a solid methods-plus-first-results paper for a fluids/HPC audience. The skeptic point on EQ conflating residual collision terms with bi-modal streaming is real and should be answered in revision, but it does not sink the literal NONEQ-vs-EQ eigenvalue claim; demanding an inconsistent hybrid operator is not necessary if the authors clarify the linearization. I would not reject on BGK/Pr=1 given that it is already disclosed. Fit is appropriate for a serious physics.flu-dyn / computational fluid dynamics venue."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new piece is real: they linearize BE-BGK about kinetic shock base flows, keep the bi-modal reduced VDFs in the stability operator, and produce the first reported kLST spectra for isolated finite-thickness normal shocks, including a 281k-unknown M=4 solve. That is not just Couette-style rarefied LST plus a shock profile.\n\nWhat they do well is concrete. Couette checks match Zou et al. to ~1% on the least-stable modes. Base-flow and spectral mesh studies are in the appendices. The HPC story is useful rather than decorative: multi-band velocity–space sparsity, LU fill-in, when Arnoldi+MUMPS dies, and when JD+block-Jacobi ILU is the only path. At M=1.2 the EQ and NONEQ spectra collapse as they should; at M=3 and 4 the NONEQ continuous branches sit systematically less damped, and the eigenvectors pick up extra upstream support that tracks the bimodal VDF. All reported modes stay stable—the claim is a relative spectral shift, not discovery of unstable shocks.\n\nSoft spots, in proportion. BGK with Pr=1 is the main physics limit; they say so and set G&P to Pr=1 for a fair comparison. The stress-test point has some bite: forcing g_c=g^e_c (and h) nulls both bi-modality in the base state and the residual (g^e−g) pieces inside the linearized collision coefficients, so the EQ/NONEQ split is not a pure “streaming with bi-modal VDF” experiment. They never show the crossed operator (non-eq base only in streaming, Maxwellian collision linearization, or the reverse). I would not call the central claim broken—NONEQ is still the consistent kinetic operator about the true base VDF—but the abstract’s “kinetic effects / translational non-equilibrium” language over-cleanly attributes the shift. No independent high-Mach stability benchmark exists, and code/data are not shipped; settings are documented enough to re-implement with effort.\n\nThis is for people who care about hypersonic LST inside shock layers and about large kinetic eigenproblems. Worth a serious referee. I would engage, cite the method and the M=3/4 shift with the BGK and EQ-definition caveats attached, and not treat continuum shock LST as automatically safe at these Mach numbers.","headline":"First kinetic LST of finite-thickness shocks with a real EQ/NONEQ spectral shift and a working O(10^5) HPC path; the isolation of “VDF non-equilibrium” is a bit messier than the abstract claims, and BGK/Pr=1 is the main physics caveat.","tokens_in":41249,"tokens_out":616,"would_cite":true,"duration_ms":11793,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Inside finite-thickness shocks, non-equilibrium velocity distributions systematically push linear stability spectra toward less stable regions than Maxwellian-based analyses predict.","keywords":["kinetic linear stability","normal shocks","Boltzmann-BGK","translational non-equilibrium","velocity distribution function","high-Mach compressible flow","parallel eigensolvers","Jacobi-Davidson"],"falsifier":"Repeat the Mach-3 and Mach-4 calculations with a collision operator that recovers the physical Prandtl number (or with a fully kinetic DSMC-linearized operator) and check whether the systematic shift of continuous branches toward weaker damping disappears or reverses.","tokens_in":41003,"feed_emoji":"💨","tokens_out":862,"duration_ms":23484,"temperature":0.7,"pith_summary":"This paper builds the first kinetic linear stability theory for isolated one-dimensional normal shocks. Instead of closing the problem on continuum variables, it linearizes the Boltzmann–BGK equation about kinetic base flows so both the mean shock and the disturbances live at the level of the velocity distribution function. At low Mach number the spectra recover the familiar stable continuous branches. At Mach 3 and 4, keeping the true non-Maxwellian distributions inside the shock shifts those branches toward weaker damping relative to spectra built from the same macroscopic profiles but forced to local Maxwellians. The authors also supply a parallel eigensolver stack so the resulting matrices—hundreds of thousands of unknowns—can actually be solved. A sympathetic reader cares because high-speed vehicle design often trusts continuum or moment-based stability tools; the work says those tools can miss a kinetic imprint even when the mean density and velocity look well resolved.","feed_headline":"Shock kinetics push high-Mach layers toward less stable spectra","feed_subtitle":"First kinetic LST of finite-thickness shocks shows non-Maxwellian VDFs weaken damping.","key_machinery":"Kinetic linear stability theory (kLST): the Boltzmann–BGK equation is linearized about a BE-BGK base flow, with disturbances posed in reduced distribution functions and macroscopic fields recovered only by velocity-space moments, so the stability operator acts on the VDF itself rather than on a closed continuum system.","core_discovery":"For argon normal shocks at Mach 3 and 4, eigenspectra constructed from non-equilibrium BE-BGK base velocity distribution functions are systematically less stable than spectra that use the same macroscopic base state but enforce local Maxwellian distributions. Continuum-style closures can therefore miss important spectral changes even when macroscopic shock profiles appear adequate.","pith_inferences":["If the less-stable kinetic shift survives better collision models, transition-prediction tools for hypersonic inlets and control surfaces may need a kinetic correction inside strong shocks.","The continuous, densely packed shock spectrum is what makes cheap ILU-based Jacobi–Davidson viable; discrete-mode problems such as Couette may still require exact shift-invert.","Comparing EQ versus NONEQ spectra is analogous to comparing a single-temperature Boltzmann emission model against a full non-Boltzmann state-to-state model: the macro state alone hides the stability-relevant microphysics."],"forward_implications":["Stability conclusions drawn from Navier–Stokes or Maxwellian-reconstructed bases for high-Mach shocks can be optimistic relative to a kinetic treatment.","Both the base flow and the perturbation problem must retain micro-velocity structure; using only macroscopic kinetic profiles is not enough.","Memory-efficient iterative eigensolvers (Jacobi–Davidson with local ILU) become necessary once velocity grids push unknowns past roughly 10^5.","The same operator framework can be extended to multidimensional and reacting shock layers once collision models improve."],"fun_headline_variants":["Non-Maxwellian VDFs shift high-Mach shock spectra toward less stable regions","Kinetic LST shows BE-BGK base states weaken damping at Mach 3–4","Continuum closures miss stability shifts inside finite-thickness shocks","kLST spectra less stable when base flows keep bi-modal VDFs","Argon shocks at Mach 4: kinetic base VDFs erode spectral damping"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The BGK collision model with Prandtl number fixed at one is assumed to capture the collision physics that set the linearized spectrum inside the shock.","fun_headline_variants_meta":{"raw":{"variants":["Non-Maxwellian VDFs shift high-Mach shock spectra toward less stable regions","Kinetic LST shows BE-BGK base states weaken damping at Mach 3–4","Continuum closures miss stability shifts inside finite-thickness shocks","kLST spectra less stable when base flows keep bi-modal VDFs","Argon shocks at Mach 4: kinetic base VDFs erode spectral damping"]},"model":"grok-4.5","effort":"low","cost_usd":0.003798,"raw_usage":{"total_tokens":1278,"prompt_tokens":913,"num_sources_used":0,"completion_tokens":84,"cost_in_usd_ticks":37984000,"prompt_tokens_details":{"text_tokens":913,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":281,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":913,"tokens_out":84,"duration_ms":5192,"temperature":1.0,"reasoning_tokens":281,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T00:49:16.648635+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Repeat the Mach-3 and Mach-4 calculations with a collision operator that recovers the physical Prandtl number (or with a fully kinetic DSMC-linearized operator) and check whether the systematic shift of continuous branches toward weaker damping disappears or reverses.","supporting_citations":[],"review_version":1}