{"id":"cde730ca-8cc2-4f7d-ad85-c17e64d1b9c3","arxiv_id":"2412.20813","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Inviscid isothermal Couette flow has a discrete tower of compressible modes that undergo repeated level-crossing instabilities, creating alternating stable and unstable bands.","lead":"Inviscid isothermal compressible Couette flow is shown to support a tower of compressible eigenmodes that become unstable in zebra-like stripes on the Mach number versus wavenumber plane. The instabilities come from level crossings between neighboring modes, a mechanism previously seen in adiabatic shear flows.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The instability pattern rests on an unvalidated Fredholm search; the infinite zebra-stripe extrapolation is not independently confirmed.","rationale":"I read the paper in good faith. The analytical derivations are internally consistent: the reduction to the second-order ODEs in Section 4.1 and the stability theorems in Section 4.2 (critical-layer condition, semicircle bound, spectral gap for M<1) are rigorous and check out. The asymptotic limits of Section 5 provide genuine support for an infinite tower of neutrally stable modes in the small-k, small-M, and high-frequency regimes. The symmetries of Section 3.3 are exact and correctly constrain the eigenvalue spectrum. What is least secure is the numerical evidence for the level-crossing instabilities: the Fredholm-alternative search of Appendix B is described only schematically, with no discretization details, no grid-convergence tests, no estimate of the numerical error in the eigenvalue positions, and no external validation. Because the operator is non-self-adjoint and the interesting events happen precisely where the leading coefficient can vanish (critical layers), the search is not protected by standard spectral approximation theorems. The infinite-sequence and zebra-stripe claims are extrapolations from the first few arches visible in Figures 3-5; the paper itself explicitly defers a proof to future work. This does not make the paper's thesis wrong, but it makes the strongest claim conditional on numerical reliability and on the extrapolation continuing. The proposed test—an independent solver based on the confluent hypergeometric formulation of Appendix C—would settle whether the numerics are sound and, by extension, whether the first few arches and their apparent recurrence are genuine. If that test agrees, the reader's conditional verdict should stand, with the infinite-sequence extrapolation still requiring either an independent computation at higher modes or an analytical proof.","tokens_in":31023,"tokens_out":10299,"duration_ms":110426,"concrete_test":"Independently verify the numerical foundation by solving the transcendental eigenvalue condition (108) from Appendix C (Whittaker/confluent hypergeometric functions) with a high-precision complex root finder, initialized from the Section 5 asymptotics. Recompute the first three (0*,0) and (0,1) instability arches for the n=0 mode for both the M-scan at k=1 and the k-scan at M=4.75, and also the first arch for n=1, comparing the gamma_r and gamma_i curves and the transition points M_c and k_c to Figs. 3-4. Agreement to about 1% would validate the Fredholm search; disagreement would reveal spurious or missing modes. Then extend the same solver to n=2 through 10 and to (k,M) values a factor of two beyond the displayed range; if the number of arches per mode stops growing or the alternating stripe pattern breaks, the infinite-sequence extrapolation would be refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—an infinite sequence of level-crossing instability windows forming zebra stripes in the (k,M) plane—rests on the numerical spectra of Appendix B, but that search has no convergence study, no error bounds, no published code or data, and no comparison against an independent method. The operator T_gamma in (104) is non-self-adjoint and, at the neutrally stable eigenvalues that participate in each merger, the leading coefficient k^2+M^2Gamma^2 can vanish inside the channel (critical-layer condition |gamma_i|<k), so the Fredholm alternative is invoked in a regime where its hypotheses are not checked. For nonnormal operators, discretizations can produce spectral pollution (spurious eigenvalues), and the method's 'norm tends to unbounded' criterion is not shown to be robust under grid refinement. Moreover, Fig. 5's stripes are sketched from four parameter slices, and the endlessly repeating pattern for each mode is an extrapolation from the first few arches (Figs. 3-4). The paper itself says 'we hope to prove the key features of our numerical spectra elsewhere' (Section 2), acknowledging that no proof of the infinite sequence is offered.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the linear stability of inviscid, isothermal, compressible Couette flow in a finite channel. It derives a compressible analogue of Rayleigh's equation, reduces it to second-order ODEs for the vorticity and velocity perturbations, and establishes symmetry properties of the spectrum. For uniform-vorticity Couette flow it proves several stability criteria: a critical-layer condition for instability, a semicircle-type bound, lower bounds on |Im(gamma)| for neutral modes, and a condition for neutral modes to possess critical layers. It then constructs analytic infinite towers of neutral modes in the small-wavenumber, small-Mach-number, and high-frequency limits. Beyond these regimes, a numerical Fredholm-alternative search is used to find complex growth rates, showing repeated level crossings between adjacent modes that open and close instability arches and produce a claimed zebra-stripe pattern of instability regions in the (k,M) plane. A canonical square-root power law near crossings is proposed and fit to the numerics. A continuous spectrum of non-smooth neutral modes is also described.","tokens_in":31324,"tokens_out":4979,"duration_ms":59503,"significance":"If the numerical results are correct, the paper would provide a striking new picture of noninflectional compressional instabilities in inviscid Couette flow, complementing earlier work by Glatzel and Renardy and extending the classical Rayleigh/Howard framework to isothermal compressible flows. The analytic parts are the strongest contribution: the integral identities in Section 4.2 and the asymptotic towers in Section 5 are clean, internally consistent, and match the numerics in overlapping regimes. The paper also gives a useful semicircle theorem and explicit neutral-stability bounds. However, the central quantitative claims about infinite level-crossing sequences and zebra stripes currently rest on an unvalidated numerical search and on extrapolation from a handful of low-lying arches; no code, data, convergence study, or independent verification is provided. The canonical square-root law is an ansatz fitted to the same numerical spectra it aims to explain. With additional validation, the paper would be a solid contribution; in its present form the main instability-pattern claims remain insufficiently supported.","major_comments":[{"comment":"The entire level-crossing scenario of Section 6 rests on the numerical Fredholm-alternative search, but the manuscript gives no convergence study, no grid-refinement test, and no independent check of the criterion that an eigenvalue is found when the norm of the solution 'tends to become unbounded.' The operator T_gamma is non-self-adjoint, and precisely in the level-crossing regime |Im(gamma)| < k the leading coefficient k^2 + M^2 Gamma^2 in (104) can have real zeros inside the channel, so the Fredholm hypotheses are not verified where the method is used. Spectral pollution is a known hazard for such problems. I request a convergence study, error estimates, a comparison with an independent method such as a spectral collocation scheme or the hypergeometric transcendental equation (108), and release of the code and data used to produce Figures 3-7.","section":"Appendix B, Eq. (104)"},{"comment":"The claim of an 'infinite sequence' of alternating instability stripes is extrapolated from the first few arches of the n = 0 mode, computed on only four parameter slices: k = 1 and k = 1.6 for variable M, and M = 4.75 and M = 7.6 for variable k. Figure 5 is itself a sketch using solid segments from those slices, with no computations shown for higher arches, higher modes, or larger parameter values. The paper explicitly states in Section 2 that the authors 'hope to prove the key features of our numerical spectra elsewhere,' which confirms that the unbounded repetition is currently a conjecture. This point is load-bearing for the abstract's claims of an infinite sequence of stability transitions and a zebra-stripe pattern; it should either be proven or clearly presented as a conjecture.","section":"Section 6.3, Fig. 5"},{"comment":"The canonical square-root law is proposed as an ansatz rather than derived from the eigenvalue problem. The coefficients A and B in (86)-(89) and the interpolating arch parameters in (91) are fitted to the same numerical spectra that the ansatz is supposed to explain, so the agreement in Figures 6 and 7 is a consistency check rather than independent confirmation. In particular, the assumption that only the 2x2 block of the two crossing modes contributes to the characteristic equation near a crossing is introduced without justification, and the biquadratic form (90) is not derived from the operator Q. The manuscript should either derive the 2x2 reduction from the resolvent structure or demonstrate that the fitted A-values are parameter-free predictions across multiple crossings and parameter values.","section":"Section 6.4, Eqs. (86)-(92)"}],"minor_comments":[{"comment":"The heading contains a typo: 'F redholm alternative' should read 'Fredholm alternative.'","section":"Appendix B title"},{"comment":"The regularity class 'C2' in the text and the remark should be written as C^2, and the sentence 'which is a regularity class' should be clarified: 'C^2' is standard notation, but 'C2' is ambiguous.","section":"Section 7"},{"comment":"The caption refers to solid and dashed boundary curves, but the figure as printed does not clearly distinguish them; a legend or explicit labels for the lines M = 1/2 and the curve M^2 = 1/4 + 1/(8k^2) would improve readability.","section":"Fig. 5 caption"},{"comment":"The notation Mc1, Mc2, and gamma_mid in Eq. (91) is introduced only in the surrounding prose; defining all symbols in the displayed equation itself would make the interpolation formula easier to follow.","section":"Section 6.4"}],"recommendation":"major_revision","confidential_remarks":"The analytic stability theorems and asymptotic mode towers are valuable and appear sound. The main risk is the numerical component: without a convergence study or reproducibility materials, the referee cannot independently assess whether the level-crossing spectra and the extrapolated zebra-stripe pattern are genuine or numerical artifacts. A revision that adds validation and clearly separates proven results from conjectured ones would make the paper suitable for publication. The manuscript fits the scope of the journal well."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here is my take. The analytical core of this paper is solid: the stability theorems in Section 4.2 are derived from integral identities and look correct, and the asymptotic formulas in Section 5 (small k, small M, constant coefficient) are explicit and checkable. The isothermal inviscid setting is genuinely new relative to Glatzel and Renardy. The paper is also honest about the limits of the numerical evidence. The soft spot is exactly where the reader put it: the infinite sequence of level-crossing instabilities and the zebra-stripe pattern rest on the Fredholm-alternative search in Appendix B, which has no convergence study, no error bounds, and no published code or data. The operator is non-self-adjoint, and at the neutral eigenvalues participating in mergers the leading coefficient can vanish inside the channel, so the hypotheses of the Fredholm alternative are not checked. The extrapolation from the first few arches to an infinite sequence is plausible but unproven. The square-root law is a fitted ansatz, not an independent prediction. Still, the numerics match the analytic limits in overlapping regimes, which gives some confidence that the low-lying spectrum is real. The stress-test note goes too far if it implies the whole paper is doubtful. The central analytical results stand. Who is this for? Hydrodynamic stability researchers working on compressible shear flows. It provides a clean benchmark and a map of instabilities for later nonlinear and viscous studies. I would send it to peer review, but a serious referee should demand numerical validation: convergence under grid refinement, comparison with an independent method, and ideally code and data. The infinite-sequence claim should be either proved for low modes or softened to a conjecture. With that revision, this is a useful contribution.","headline":"A solid analytical core with clean stability theorems and explicit mode formulas, but the headline infinite zebra-stripe pattern rests on a numerical search that lacks convergence checks and is extrapolated from low modes.","tokens_in":765,"tokens_out":1968,"would_cite":true,"duration_ms":36470,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76E05","76N10"],"pacs":["47.20.Ft","47.40.-x"],"model":"deepseek-v4-flash","headline":"Repeated level crossings make compressible Couette flow unstable in stripes.","keywords":["Couette flow","linear stability","compressible isothermal flow","level crossing","zebra stripes","continuous spectrum","noninflectional modes","growth rate"],"falsifier":"Calculate the full spectrum at larger $M$ and $k$ with an independent high-order method, for example by solving the transcendental equation from the confluent hypergeometric form, and check whether the first instability arch for $\\tilde{k}=1$ opens at $M\\approx 4.203$ and whether the imaginary and real parts of $\\tilde{\\gamma}$ follow the square-root law with fitted slope near $0.5$; a missing arch, a spurious arch, or a mode with $|\\Im\\tilde{\\gamma}|<\\tilde{k}$ that never crosses would show that the infinite sequence of stripes is not real.","tokens_in":30852,"feed_emoji":"🦓","tokens_out":6333,"duration_ms":61800,"temperature":0.7,"pith_summary":"The paper studies the linear stability of an inviscid, isothermal, compressible shear flow with a linear velocity profile, inviscid Couette flow, in a finite channel. In the incompressible limit this flow has no unstable modes, but once compressible perturbations are allowed, an infinite tower of noninflectional discrete modes appears. The paper argues that as the Mach number $M$ or wavenumber $k$ grows, neighboring purely oscillatory modes repeatedly collide on the imaginary axis, split into complex pairs, and thereby open finite windows of instability. These windows alternate with stable regions and form zebra-like stripes in the $k$-$M$ plane, with growth rates growing like the square root of distance from each transition. The picture matters because it shows that compressibility alone, without viscosity or inflection points, can produce a rich, organized pattern of linear instabilities.","feed_headline":"Level crossings stripe the stability map of compressible Couette flow","feed_subtitle":"Colliding sound-driven modes open zebra-striped instability windows as Mach number or wavenumber grows.","key_machinery":"The central object is the vortex-perturbation equation $\\tilde{w}'' - \\frac{2i\\tilde{k}}{\\Gamma}\\tilde{w}' - (\\tilde{k}^2 + M^2\\Gamma^2)\\tilde{w} = 0$, where $\\Gamma = \\tilde{\\gamma} + i\\tilde{k}\\tilde{y}$ is the Doppler-shifted growth rate and $\\tilde{w}$ is the vorticity perturbation; equivalently, a second-order equation for the normal velocity with Dirichlet boundary conditions at the channel walls. These equations are transformable into a confluent hypergeometric (Whittaker) equation, which expresses the same spectrum but is not solved analytically. The argument is carried by the merger-demerger of adjacent imaginary eigenvalues as $\\tilde{k}$ or $M$ varies, captured near each transition by the biquadratic characteristic equation $(\\tilde{\\gamma}^2 - \\tilde{\\gamma}_c^2 + A^2(M_c-M))^2 + 4\\tilde{\\gamma}_c^2 A^2(M_c-M) = 0$, whose two branches produce the square-root growth law. A $\\mathbb{Z}_2\\times\\mathbb{Z}_2$ symmetry constrains eigenvalues to four-tuples and fixes the allowed eigenfunction forms.","core_discovery":"For inviscid isothermal compressible Couette flow, the linearized eigenvalue problem reduces to a single second-order ODE for vorticity perturbations with the growth rate appearing nonlinearly. In the small-wavenumber, small-Mach, and high-frequency limits, the equation has an infinite tower of neutrally stable modes labeled by an integer $n$; each has a conjugate partner, and modes come in four-tuples. Beyond those limits, the purely imaginary eigenvalues of neighboring modes meet on the imaginary axis, leave the axis as a complex-conjugate pair, and later merge again, producing an arch-shaped interval of instability. The paper identifies this repeated level-crossing mechanism as the source of an infinite sequence of stability transitions, derives analytic stability theorems and bounds that delineate the stable regions, and fits the numerically computed eigenvalues near each crossing to a canonical biquadratic characteristic equation that gives square-root power-law behavior. A separate continuous spectrum of neutrally stable eigenmodes exists along the imaginary axis, with eigenfunctions that are only piecewise smooth across critical layers.","pith_inferences":["If the infinite-tower extrapolation is correct, a similar zebra-stripe pattern should appear for every mode, shifted toward larger $k$ and $M$; locating the first unstable band for a high mode $n$ would test the extrapolation.","Because the discrete operator is nonnormal, the paper's exponential-mode analysis does not rule out finite-time transient growth that could blur the sharp linear stability boundaries in numerical or experimental settings.","The level-crossing mechanism is not tied to isothermal thermodynamics or two-dimensional perturbations: the same merger-demerger route should appear in adiabatic and three-dimensional perturbations, so the stripes may be a general organizing feature of compressible shear-flow stability.","A direct probe of the $k$-$M$ plane by measuring the response to monochromatic perturbations could look for the predicted alternating stable and unstable bands in an experiment or simulation."],"forward_implications":["For small $k$ or small $M$ the flow is guaranteed neutrally stable; the paper proves $M<1/2$ is sufficient for stability at all wavenumbers.","Each mode $n$ undergoes an infinite sequence of stable-unstable transitions as $M$ or $k$ grows, so the number of participating modes increases while the maximum growth rate in each arch decreases.","The first instability always involves the ground mode colliding with its conjugate, and for the ground mode the resulting windows alternate between $(0^*,0)$ and $(0,1)$ mergers.","Near every transition, the real and imaginary parts of the growth rate scale as the square root of distance in $M$ or $k$, with a coefficient fixed by the biquadratic ansatz.","The continuous spectrum fills the imaginary interval $-\\tilde{k}<\\Im\\tilde{\\gamma}<\\tilde{k}$ and overlaps the discrete spectrum in the level-crossing regime, but its eigenfunctions are not smooth across critical layers."],"supporting_citations":[{"why":"Provides the standard incompressible hydrodynamic-stability background and the statement that inviscid Couette flow has no unstable modes.","marker":"[1]"},{"why":"Gives the incompressible Rayleigh equation and inflection-point criterion that the compressible problem reduces to in the $M\\to0$ limit.","marker":"[4]"},{"why":"Contains a qualitative phase-plane analysis of the same vorticity-perturbation equation, distinguishing subsonic and supersonic disturbances.","marker":"[16]"},{"why":"Supplies the continuous-spectrum construction for inviscid Couette flow that the paper adapts to the compressible case.","marker":"[17]"},{"why":"Found a similar Mach-number pattern of instabilities in an adiabatic linear shear layer, giving independent evidence for the level-crossing mechanism.","marker":"[21]"},{"why":"Provides a proof of instability arising from eigenvalue crossing in transonic shear flows, supporting the interpretation of level crossings as a general instability route.","marker":"[23]"}],"fun_headline_variants":["Zebra-striped instability windows from level crossings in gas Couette flow","Mode collisions paint striped stability maps in compressible Couette flow","Level crossing cascades open infinite instability windows in Couette flow","Compressible Couette flow: colliding modes stripe the stability map","Infinite tower of modes collides to stripe Couette stability plane"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that instabilities repeat indefinitely rests on a numerical eigenvalue search that is assumed to track the true spectrum of a non-self-adjoint operator without missing or spurious modes, extrapolating from the first few arches to all modes and all parameter ranges.","fun_headline_variants_meta":{"raw":{"variants":["Zebra-striped instability windows from level crossings in gas Couette flow","Mode collisions paint striped stability maps in compressible Couette flow","Level crossing cascades open infinite instability windows in Couette flow","Compressible Couette flow: colliding modes stripe the stability map","Infinite tower of modes collides to stripe Couette stability plane"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00045,"raw_usage":{"total_tokens":2317,"prompt_tokens":1040,"completion_tokens":1277,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":1185}},"tokens_in":656,"tokens_out":1277,"duration_ms":11721,"temperature":1.0,"reasoning_tokens":1185,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:09:30.912232+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Calculate the full spectrum at larger $M$ and $k$ with an independent high-order method, for example by solving the transcendental equation from the confluent hypergeometric form, and check whether the first instability arch for $\\tilde{k}=1$ opens at $M\\approx 4.203$ and whether the imaginary and real parts of $\\tilde{\\gamma}$ follow the square-root law with fitted slope near $0.5$; a missing arch, a spurious arch, or a mode with $|\\Im\\tilde{\\gamma}|<\\tilde{k}$ that never crosses would show that the infinite sequence of stripes is not real.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard incompressible hydrodynamic-stability background and the statement that inviscid Couette flow has no unstable modes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the incompressible Rayleigh equation and inflection-point criterion that the compressible problem reduces to in the $M\\to0$ limit."},{"cited_title":"2, 5, 13","cited_arxiv_id":null,"evidence_quote":"Contains a qualitative phase-plane analysis of the same vorticity-perturbation equation, distinguishing subsonic and supersonic disturbances."},{"cited_title":"Fluids, 8, 143 (1960)","cited_arxiv_id":null,"evidence_quote":"Supplies the continuous-spectrum construction for inviscid Couette flow that the paper adapts to the compressible case."},{"cited_title":"3, 8, 13, 31","cited_arxiv_id":null,"evidence_quote":"Found a similar Mach-number pattern of instabilities in an adiabatic linear shear layer, giving independent evidence for the level-crossing mechanism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides a proof of instability arising from eigenvalue crossing in transonic shear flows, supporting the interpretation of level crossings as a general instability route."}],"review_version":1}