{"id":"fe13fb67-3aeb-4a6c-b505-9ec0c6e19a19","arxiv_id":"2502.02258","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For small viscosity, exact Tollmien-Schlichting wave solutions exist for wavenumbers from ν^{1/8} to ν^{1/12}, with imaginary phase speed crossing zero at the lower and upper neutral branches.","lead":"This paper constructs stable, neutral, and unstable Tollmien-Schlichting waves for boundary layer flow over a flat plate, confirming the location of the neutral stability curve at the predicted Reynolds-number scalings. The result turns a classical asymptotic picture from fluid mechanics into a theorem under mild regularity assumptions on the shear profile.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lower-branch neutral crossing rests on an unproved figure-based property of the Tietjens function; without a rigorous quantitative statement, the existence of A_c with Im(c(A_cν^{1/8}))=0 is not established.","rationale":"The reader's verdict is CONDITIONAL with moderate confidence, and the reader's rationale explicitly mentions the 'textbook-figure property of the Tietjens function for the lower branch' as a point to address before acceptance. However, the reader's formal weakest_assumption field focuses on the modified Rayleigh operator approximation in §1.2.1. My stress-test identifies a different, more concrete hinge: the lower-branch neutral crossing A_c in Proposition 5.1 depends on a quantitative property of the Tietjens function that is only cited from a figure. This is load-bearing because Theorem 1.1's first bullet claims exactly a sign change of Im(c) across α = A_cν^{1/8}, and the proof of that sign change uses the figure-based property to produce a point with Im(c) < 0 at A0 ∼ 1. The modified Rayleigh construction is also central, but it is at least supported by detailed estimates and follows the framework of [27]; the Tietjens step has a clear gap between a figure and the required O(ν^{1/16}) error control. A numerical verification of the Tietjens property with a quantitative slope bound would settle whether the gap is merely a missing proof or a genuine obstruction. Since this concern reinforces the existing CONDITIONAL verdict rather than overturning it, the verdict should remain UNCHANGED.","tokens_in":64029,"tokens_out":19702,"duration_ms":160551,"concrete_test":"Compute the Tietjens function F(z) to high precision on a fine grid covering Re z ∈ [2, 2.5] using its standard representation from hydrodynamic stability theory (e.g., via the Airy/Hankel integral used in Fig 3.2 of [26]). Verify: (i) Im F(z) < 0 for Re z < z0, (ii) Im F(z0) = 0 for some z0 ∈ [2, 2.5], (iii) |d/dz Im F(z0)| > 1, and (iv) Re F(z0) ∈ [2, 3]. These quantitative bounds would justify the O(ν^{1/16}) absorbing step in Proposition 5.1. If the slope at z0 is smaller than the error threshold, the claimed A_c crossing is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Proposition 5.1, Case 1, the proof of the lower-branch sign change reduces to equation (5.9): F(−κη(0)) + O(ν^{1/16}) = u′(0)c/α. To obtain a point A0 with Im(c(A0ν^{1/8})) < 0, the paper invokes a property of the Tietjens function F taken from a figure: 'By the property of F(z) (see Fig 3.2 in [26]), we know that there exists z0 ∈ [2, 2.5] such that Im(F(z0)) = 0 and Im(F(z)) < 0 for any z ∈ [2, z0).' The proof then asserts existence of a small perturbation z1 with Re z1 < 0 and |Im z1| ≤ Cν^{1/16} such that Im(F(z0+z1)+O(ν^{1/16})) < −(1/8)ν^{1/16}. This requires a quantitative transversality statement: near z0, Im F must leave the real axis with a slope large enough to dominate the O(ν^{1/16}) error. The paper neither states this as a precise lemma nor proves it; a textbook figure confirms the qualitative shape but not the quantitative margin needed to absorb the error term. If the slope of Im F at z0 is small or zero, the O(ν^{1/16}) perturbation can flip the sign, and no A0 with Im c < 0 would follow. The upper-branch crossing B_c is obtained from explicit inequalities (5.10), so this concern is specific to the lower branch. Theorem 1.1's first bullet, which asserts a stable side below A_c and an unstable side above A_c, is exactly the point that depends on this Tietjens property. This is a genuine gap in the rigor of the central claim, even though the underlying physical picture may be correct.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies linear stability of boundary-layer shear flows over a flat plate. After rescaling, Tollmien-Schlichting waves correspond to solutions of the Orr-Sommerfeld equation (1.7). The authors construct the slow and fast modes of the homogeneous Orr-Sommerfeld operator via a modified Rayleigh-Airy iteration, with the key new ingredient being the replacement of the singular Rayleigh operator by Ray_{c+i|ε|α^{-3/2}}. They then solve the dispersion relation (5.3) and obtain Theorem 1.1: for every α in (Aν^{1/8}, Bν^{1/12}) there is a phase speed c(α) and a W^{2,∞} T-S wave, with Im(c) crossing from negative to positive at lower and upper branch points and with α Im(c) ∼ ν^{1/4} in the interior. Much of the technical machinery, especially the Rayleigh and Airy estimates in Sections 2 and 3, is quoted from the authors' earlier paper [27].","tokens_in":64527,"tokens_out":10895,"duration_ms":101118,"significance":"If the theorem is valid, it gives the first rigorous confirmation of the neutral stable curve for this problem, complementing the construction of unstable T-S waves by Grenier, Guo, and Nguyen. The result is substantive: it identifies the asymptotic locations α ∼ ν^{1/8} and α ∼ ν^{1/12}, gives the growth-rate scaling α Im(c) ∼ ν^{1/4}, and removes the analyticity assumption of [15] through a modified Langer transformation. The derivation is not fitted: c0 is a fixed regularization parameter and the sign changes come from the dispersion relation, so the theorem is explicit and falsifiable. The proof is long and depends heavily on estimates quoted from the authors' earlier paper [27], which limits the self-containedness of the manuscript.","major_comments":[{"comment":"The proof of the lower-branch crossing is incomplete. After reducing to F(−κη(0)) + O(ν^{1/16}) = u′(0)c/α, the text invokes “the property of F(z) (see Fig 3.2 in [26])” to find z0 ∈ [2, 2.5] with Im F(z0) = 0 and Im F(z) < 0 for z ∈ [2, z0). It then asserts that a perturbation z1 with |Im z1| ≤ Cν^{1/16} yields Im(F(z0+z1)+O(ν^{1/16})) < −(1/8)ν^{1/16}. This requires a quantitative transversality statement, e.g. a lower bound on |Im F′(z0)| or at least a known finite order of vanishing at z0. A textbook figure gives the qualitative sign only; if the derivative at z0 is zero or too small, the O(ν^{1/16}) error can flip the sign and no A0 with Im(c) < 0 follows. Since Theorem 1.1's first bullet depends exactly on this point, Proposition 5.1 needs a precise lemma giving the definition of F and a quantitative lower bound on |Im F| on a ν-independent interval to the left of z0.","section":"Section 5, Proposition 5.1, Case 1, Eq. (5.9)"},{"comment":"The construction of the slow and fast modes, and hence the dispersion relation, rests on a large body of estimates quoted without proof from [27]: Proposition 2.1, Lemmas 2.2 and 2.4, Lemmas 3.1–3.4, Propositions 3.5–3.7, Theorems 3.8–3.10, and Lemma A.5. Because [27] is co-authored by two of the present authors, this is not circular in the sense of fitting a target quantity, but it makes the central claim difficult to verify from the manuscript alone. The key assumption behind (1.12) — that Ray_{c+i|ε|α^{-3/2}} approximates the Orr-Sommerfeld operator to the order required by the iteration — is justified only through these quoted estimates. The paper should either state the needed estimates as explicit hypotheses or provide proofs; at minimum it should explicitly verify that the hypotheses of [27], in particular the c0 range |ε|^{2/3} ≪ c0 ≪ |ε|^{1/3}, are satisfied with the choice c0 = |ε|α^{-3/2} throughout the interval (Aν^{1/8}, Bν^{1/12}).","section":"Sections 2–3 and Theorem 4.1"}],"minor_comments":[{"comment":"There are several spelling errors: “Naiver” before (1.3), “Nyugen” in the discussion of [15], “exsit” in Theorem 1.1, and “adn” in reference [3].","section":"Throughout"},{"comment":"In the displayed partial derivatives, the third line repeats ∂crFi; it should presumably be ∂ciFr, and the fourth line should be ∂ciFi.","section":"Proposition 5.2"},{"comment":"The Tietjens function F(z) is not defined in the paper; the proof should provide its formula or a precise reference with a quantitative statement, rather than only citing a figure.","section":"Eq. (5.9)"},{"comment":"The notation H1, H2 is used both for the parameter sets in (2.5) and for the source-term functions in (4.12), which is confusing but harmless.","section":"Section 4.2"},{"comment":"The strict inequality A0 > A appears in Theorem 1.1 but only A0 is mentioned in Proposition 5.1; the statements should be made consistent.","section":"Theorem 1.1 and Proposition 5.1"}],"recommendation":"major_revision","confidential_remarks":"The novelty over [27] and [15] is real but essentially consists of completing the neutral-curve crossing within the same technical framework. The lower-branch Tietjens issue is a genuine missing quantitative argument and should be resolved before acceptance; it is not a question of physical plausibility but of rigorous proof. The editor may also wish to consider whether the heavy quotation of estimates from the authors' own paper [27] is acceptable without a more self-contained presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a serious paper with a genuine new result, and the one load-bearing soft spot is a missing quantitative lemma, not a wrong idea.\n\nWhat’s new: it constructs stable and neutral T-S waves on both branches of the neutral curve, removes the analyticity assumption from [15] via a modified Langer transformation, and gets the interior growth rate α Im(c) ~ ν^{1/4}. That is a real advance over Grenier–Guo–Nguyen and over Bian–Grenier. The modified Rayleigh–Airy iteration is laid out in detail, Sections 2–4 are coherent, and the dispersion-relation framework is the right way to organize the problem. The paper is honest about what is borrowed and what is added.\n\nThe soft spot: the lower-branch crossing in Proposition 5.1, Case 1, rests on a figure-based property of the Tietjens function. The argument needs a quantitative statement: near the zero z0 of Im F, the curve must leave the real axis with a slope large enough to dominate the O(ν^{1/16}) error. The paper neither states nor proves that. The upper branch uses explicit inequalities and is fine. This is a genuine gap in the proof as written, but it is fixable — the Tietjens function is classical and a short lemma with a derivative bound should settle it. I do not think the main theorem is wrong; I think a referee should demand the lemma before publication.\n\nOn the reliance on [27]: heavy, but not circular. Those are published, peer-reviewed lemmas. Still, since two authors co-authored [27], a referee should spot-check that the quoted estimates cover the exact parameter ranges used here, especially near ci = 0.\n\nWho this is for: anyone working on hydrodynamic stability or boundary layer transition. It deserves a serious referee. I would send it to review and ask for a quantitative proof of the Tietjens property, plus a careful statement of which [27] estimates are being quoted.","headline":"The neutral-curve construction is real progress, and the one real gap is a missing quantitative lemma on the Tietjens function on the lower branch, which looks fixable.","tokens_in":64948,"tokens_out":1774,"would_cite":true,"duration_ms":20321,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76E05","35Q30","76D10","35B35"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that Tollmien-Schlichting waves exist throughout a neighborhood of the neutral stability curve, with stable waves below, growing waves inside, and exact neutral modes at the two branch crossings.","keywords":["Tollmien-Schlichting waves","neutral stability curve","Orr-Sommerfeld equation","modified Rayleigh operator","Rayleigh-Airy iteration","triple-deck theory","boundary layer transition","linearized Navier-Stokes"],"falsifier":"A direct numerical solution of the Orr-Sommerfeld eigenvalue problem for a profile satisfying (1.8), at small viscosity, that computes Im(c(α)) through the whole range ($Aν^{{1/8}}$, $Bν^{{1/12}}$) would falsify the theorem if it did not show two sign changes near those scalings or if the interior growth rate α Im(c) did not scale like $ν^{{1/4}}$.","tokens_in":63870,"feed_emoji":"🌊","tokens_out":7478,"duration_ms":67766,"temperature":0.7,"pith_summary":"Boundary-layer flow over a flat plate is known to develop Tollmien-Schlichting (T-S) waves, and the physical theory predicts that stability is lost and regained on a two-branch neutral curve. Previous rigorous work constructed the unstable T-S waves away from the neutral curve. This paper establishes the neutral curve itself: for small viscosity, it constructs T-S wave solutions for every wave number α between a lower branch α ∼ $ν^{{1/8}}$ and an upper branch α ∼ $ν^{{1/12}}$, with the temporal growth rate Im(c) negative below the lower branch, positive in the interior, and negative above the upper branch, crossing zero exactly on the two branches. This confirms the transition picture that T-S waves are responsible for the first stage of laminar-turbulent transition.","feed_headline":"Tollmien-Schlichting waves proven to cross the neutral curve","feed_subtitle":"A modified Rayleigh-Airy construction places the stable-to-unstable transition at α ∼ ν^{1/8} and α ∼ ν^{1/12}.","key_machinery":"The load-bearing object is the modified Rayleigh operator Ray_ĉ[φ] = (u−ĉ)(∂$_Y^{2}$−$α^{2}$)φ − u''φ with ĉ = c + i|ε| $α^{{-3/2}}$, a regularized version of the singular Rayleigh operator used on the neutral curve. The imaginary shift c0 = |ε| $α^{{-3/2}}$ encodes the idea that viscous diffusion regularizes the critical layer and effectively pushes the eigenvalue left in the complex plane. Around this operator the paper sets up a modified Rayleigh-Airy iteration that builds two independent solutions of the Orr-Sommerfeld equation, the slow mode φ_s and fast mode φ_f, and then solves the dispersion relation φ_s(0)/∂_Y φ_s(0) = φ_f(0)/∂_Y φ_f(0). This dispersion relation, valid for ci ≤ 0, is what produces the neutral crossings.","core_discovery":"The central result, Theorem 1.1, is that for shear profiles satisfying the structural assumptions (1.8), for every sufficiently small viscosity ν and every α in ($Aν^{{1/8}}$, $Bν^{{1/12}}$), there is a complex phase speed c(α) and a $W^{{2,∞}}$ solution of the linearized Navier-Stokes equations in the Tollmien-Schlichting form. The imaginary part of c changes sign twice: it is negative at the lower-branch edge, zero at an intermediate value A_c $ν^{{1/8}}$, positive in the interior, zero again at B_c $ν^{{1/12}}$, and negative above the upper branch; in the interior the growth rate satisfies α Im(c) ∼ $ν^{{1/4}}$. The paper interprets the two zeros as the lower and upper branches of the neutral stability curve, and identifies the lower branch with the competition between viscous diffusion and sublayer instability and the upper branch with the stabilizing curvature u''(0)<0 of the background profile.","pith_inferences":["If the modified-Rayleigh shift |ε|α^{-3/2} is the right regularization generally, the same device should locate neutral curves for other high-Reynolds shear flows, including channel and pipe flows, where the analogous singular Rayleigh problem appears.","The paper notes that for the Blasius profile the physical upper branch is α_up ∼ ν^{1/20}; extending the transition argument to that branch would require additional estimates beyond the present theorem, and the method developed here is a plausible starting point.","An immediate testable extension is to compare the predicted dispersion relation, expressed through the Tietjens-type function F(−κη(0)), with direct numerical solution of the full Orr-Sommerfeld equation at finite but small ν; the crossing locations and the ν^{1/4} growth law are quantitative enough to check numerically."],"forward_implications":["The lower and upper branches of the neutral curve are located at α ∼ ν^{1/8} and α ∼ ν^{1/12}, and on each branch there is a neutral T-S wave with Im(c)=0.","For wave numbers strictly inside the two branches the T-S waves are unstable, with growth rate α Im(c) of order ν^{1/4}; outside the branches they are damped.","The construction admits shear profiles without analyticity; the only structural conditions are the ones in (1.8), and u''(0)<0 is not needed near the lower branch, which allows the Blasius profile there.","The transition of linear stability near the lower branch is governed by diffusion versus sublayer instability, while near the upper branch it is governed by the curvature of the background flow.","The neutral stable curve separates stable from unstable Tollmien-Schlichting modes exactly at the predicted scalings, matching the physical neutral curve of Tollmien, Schlichting and Lin."],"supporting_citations":[{"why":"Supplies the Rayleigh-Airy iteration framework and the prior construction of unstable T-S waves that this paper extends to the neutral curve.","marker":"[15]"},{"why":"Provides the rigorous estimates for the Rayleigh and Airy equations, including the modified Langer transformation, on which the iteration and dispersion relation rest.","marker":"[27]"},{"why":"Contains the asymptotic neutral-curve calculations of Tollmien, Schlichting and Lin that the theorem is designed to confirm.","marker":"[26]"},{"why":"Earlier progress on the neutral curve for analytic basic flows, which this paper goes beyond by treating non-analytic profiles and proving the full transition.","marker":"[4]"},{"why":"The Airy-function asymptotics used for the modified Langer transformation and the boundary-value expansions in the dispersion relation.","marker":"[29]"}],"fun_headline_variants":["Neutral stability curve for Tollmien-Schlichting waves confirmed","T-S waves near neutral curve: existence proven","Triple-deck method proves neutral stability curve","Boundary layer stability: neutral curve existence established","Neutral curve existence for Tollmien-Schlichting waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the claim that near the neutral curve the Orr-Sommerfeld operator is accurately approximated, in the upper and main decks, by the modified Rayleigh operator with the specific imaginary shift |ε|$α^{{-3/2}}$, an approximation justified heuristically and whose rigorous estimates are quoted from an earlier paper rather than proved here.","fun_headline_variants_meta":{"raw":{"variants":["Neutral stability curve for Tollmien-Schlichting waves confirmed","T-S waves near neutral curve: existence proven","Triple-deck method proves neutral stability curve","Boundary layer stability: neutral curve existence established","Neutral curve existence for Tollmien-Schlichting waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000858,"raw_usage":{"total_tokens":3737,"prompt_tokens":968,"completion_tokens":2769,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":584,"completion_tokens_details":{"reasoning_tokens":2691}},"tokens_in":584,"tokens_out":2769,"duration_ms":19544,"temperature":1.0,"reasoning_tokens":2691,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T12:45:53.979349+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical solution of the Orr-Sommerfeld eigenvalue problem for a profile satisfying (1.8), at small viscosity, that computes Im(c(α)) through the whole range ($Aν^{{1/8}}$, $Bν^{{1/12}}$) would falsify the theorem if it did not show two sign changes near those scalings or if the interior growth rate α Im(c) did not scale like $ν^{{1/4}}$.","supporting_citations":[{"cited_title":"Grenier, Y","cited_arxiv_id":null,"evidence_quote":"Supplies the Rayleigh-Airy iteration framework and the prior construction of unstable T-S waves that this paper extends to the neutral curve."},{"cited_title":"Masmoudi, Y","cited_arxiv_id":null,"evidence_quote":"Provides the rigorous estimates for the Rayleigh and Airy equations, including the modified Langer transformation, on which the iteration and dispersion relation rest."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the asymptotic neutral-curve calculations of Tollmien, Schlichting and Lin that the theorem is designed to confirm."},{"cited_title":"Asymptotic behaviour of solutions of linearized Navier Stokes equations in the long waves regime","cited_arxiv_id":"2312.16938","evidence_quote":"Earlier progress on the neutral curve for analytic basic flows, which this paper goes beyond by treating non-analytic profiles and proving the full transition."},{"cited_title":"Olver, D","cited_arxiv_id":null,"evidence_quote":"The Airy-function asymptotics used for the modified Langer transformation and the boundary-value expansions in the dispersion relation."}],"review_version":1}