{"id":"7378c46b-c5c8-498c-b63f-0b9b1e890816","arxiv_id":"2605.31380","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"DNS in quasi-1D Rayleigh-Bénard convection identifies subcritical transition to turbulence with three distinct hysteresis loops, showing linear stability of the steady state but finite-amplitude triggering of turbulence.","lead":"Direct numerical simulations of Rayleigh-Bénard convection in a narrow quasi-one-dimensional box reveal three coexisting states—steady convection, oscillatory chaos, and intermittent turbulence—with abrupt jumps and multiple hysteresis loops in heat and momentum transport. A smart generalist might read it because the work challenges the standard picture of how buoyancy-driven flows become turbulent and hints at shared mechanisms with shear-driven turbulence.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"Quasi-1D confinement may induce apparent subcriticality via mode suppression not present in wider domains","rationale":"The reader's weakest assumption is precisely the load-bearing point; the full text would need to contain explicit wider-domain controls or linear-stability analysis to remove the concern. Absent those, the verdict moves from UNVERDICTED to CONDITIONAL pending verification that the subcritical route is not confinement-induced.","tokens_in":1735,"tokens_out":319,"duration_ms":18867,"concrete_test":"Repeat the Ra sweep at the reported subcritical window using identical numerics but with horizontal aspect ratio doubled; verify whether the steady state remains linearly stable under small random perturbations and whether the three hysteresis loops survive. A change to supercritical transition or appearance of linear instability falsifies the claim that the observed subcriticality is geometry-independent.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim requires that the reported linear stability of the steady convective state (to infinitesimal perturbations) together with its finite-amplitude transition to intermittent turbulence reflects intrinsic buoyancy-driven dynamics rather than an artifact of the narrow geometry. In standard Rayleigh-Bénard convection the primary bifurcation is supercritical; the quasi-1D setup can eliminate transverse modes whose linear growth would otherwise destroy the steady state, potentially creating artificial hysteresis loops and subcritical behavior. The abstract explicitly ties the result to both the quasi-1D domain and the static/quasi-static DNS, so any untested dependence on aspect ratio directly undermines the assertion of a general subcritical route.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript presents static and quasi-static direct numerical simulations of Rayleigh-Bénard convection in a quasi-one-dimensional domain. It reports the coexistence of steady convection, oscillatory chaos, and intermittent turbulence within a narrow Ra range, with abrupt jumps in Nu and Re, three distinct hysteresis loops (normal, reverse, anomalous), and a linearly stable steady state that transitions to turbulence under finite-amplitude perturbations, establishing a subcritical route contrary to the prevailing supercritical view.","tokens_in":1876,"tokens_out":442,"duration_ms":19104,"significance":"If the subcritical transition and hysteresis are shown to be intrinsic rather than confinement-induced, the result would be significant for fluid dynamics, providing the first clear demonstration of subcriticality in buoyancy-driven turbulence and supporting a unified framework with shear-driven flows. The use of both static and quasi-static DNS to track the states is a methodological strength.","major_comments":[{"comment":"Abstract: the claim of a general 'subcritical route' for buoyancy-driven turbulence is load-bearing on the assertion that the observed linear stability and finite-amplitude transition are not artifacts of quasi-1D mode suppression. The abstract explicitly links the result to the narrow domain and the DNS protocols, yet no tests varying the transverse aspect ratio are referenced to rule out elimination of linearly unstable transverse modes that would destroy the steady state in wider geometries.","section":"Abstract"},{"comment":"The central evidence for subcriticality (linear stability to infinitesimal perturbations combined with transition under finite-amplitude disturbances) requires explicit verification that the steady convective base state remains linearly stable only because of the confinement; without such checks or comparisons to wider domains, the distinction from standard supercritical primary bifurcation in Rayleigh-Bénard convection cannot be established as intrinsic.","section":"Abstract"}],"minor_comments":[{"comment":"The abstract mentions 'static and quasi-static DNS' but provides no details on grid resolution, time-stepping, or perturbation protocols; these should be added to allow reproducibility assessment.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading and constructive comments. We respond point-by-point to the major comments, emphasizing that our study is confined to the quasi-one-dimensional setup as stated in the title, abstract, and manuscript.","responses":[{"response":"Our work is explicitly limited to a quasi-one-dimensional domain, as indicated throughout the manuscript including the abstract, which links the findings to the narrow domain and the static/quasi-static DNS protocols. We demonstrate that, within this confined geometry, the steady convective state is linearly stable to infinitesimal perturbations yet transitions to intermittent turbulence under finite-amplitude disturbances, establishing subcriticality in this setting. We make no claim that the behavior is independent of confinement or occurs in wider domains; the result shows that a subcritical route is possible under quasi-1D confinement. Tests varying the transverse aspect ratio are not included, as they lie outside the scope of the present study.","revision_made":"no","referee_comment":"[Abstract] Abstract: the claim of a general 'subcritical route' for buoyancy-driven turbulence is load-bearing on the assertion that the observed linear stability and finite-amplitude transition are not artifacts of quasi-1D mode suppression. The abstract explicitly links the result to the narrow domain and the DNS protocols, yet no tests varying the transverse aspect ratio are referenced to rule out elimination of linearly unstable transverse modes that would destroy the steady state in wider geometries."},{"response":"The evidence for subcriticality consists of the verified linear stability of the steady state to infinitesimal perturbations together with its transition under finite-amplitude disturbances, as obtained from the DNS in the quasi-1D domain. We do not claim this linear stability holds only because of confinement in a manner that makes subcriticality intrinsic to all buoyancy-driven flows, nor do we assert a distinction from the standard supercritical bifurcation outside the confined geometry. The abstract and text present the subcritical route as occurring under the specified quasi-1D confinement, contrary to the prevailing view for standard cases. Direct comparisons to wider domains are not performed in this work.","revision_made":"no","referee_comment":"[Abstract] The central evidence for subcriticality (linear stability to infinitesimal perturbations combined with transition under finite-amplitude disturbances) requires explicit verification that the steady convective base state remains linearly stable only because of the confinement; without such checks or comparisons to wider domains, the distinction from standard supercritical primary bifurcation in Rayleigh-Bénard convection cannot be established as intrinsic."}],"tokens_in":1371,"tokens_out":532,"duration_ms":23657,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point is that this work finds a subcritical route to turbulence in buoyancy-driven flow, with three coexisting states and multiple hysteresis loops in a quasi-1D domain.\n\nStatic and quasi-static DNS show steady convection, oscillatory chaos, and intermittent turbulence in a narrow Ra range. Transitions produce jumps in Nu and Re, the steady state is linearly stable to small disturbances but responds to finite-amplitude ones, and three hysteresis loops appear in the Nu-Ra plane.\n\nThe paper does a clear job documenting these states and the global transport changes. It supplies a concrete numerical example that differs from the usual supercritical picture.\n\nThe soft spot is the quasi-1D confinement itself. Standard Rayleigh-Benard convection has a supercritical primary bifurcation, and narrowing the domain can remove transverse modes that would otherwise destabilize the base flow linearly. The abstract ties the subcritical behavior to this geometry, yet without aspect-ratio checks or explicit mode analysis the result risks being an artifact rather than a general feature of buoyancy-driven flows.\n\nThis is for people working on convection transitions and hysteresis. It deserves peer review so referees can examine the resolution, perturbation methods, and whether the subcriticality survives in wider domains.","headline":"The paper reports subcritical turbulence onset with hysteresis in quasi-1D Rayleigh-Benard convection, but the narrow geometry is likely creating the effect through mode suppression.","tokens_in":2363,"tokens_out":324,"would_cite":false,"duration_ms":16995,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Buoyancy-driven convection transitions to turbulence subcritically through finite disturbances and multiple hysteresis loops.","keywords":["Rayleigh-Bénard convection","subcritical transition","buoyancy-driven turbulence","hysteresis loops","direct numerical simulation","finite-amplitude instability","quasi-one-dimensional confinement"],"falsifier":"An experiment or simulation in the same quasi-one-dimensional domain that finds only supercritical onset, no hysteresis, and no stable coexistence of the three states would falsify the subcritical claim.","tokens_in":2631,"feed_emoji":"","tokens_out":686,"duration_ms":14479,"temperature":0.7,"pith_summary":"The paper establishes that in a quasi-one-dimensional Rayleigh-Bénard setup, the flow can pass from steady convection to intermittent turbulence without the gradual growth of small disturbances that defines a supercritical transition. Instead, the steady state remains linearly stable yet yields abruptly to turbulence once a finite-amplitude kick is applied, producing three coexisting states and three distinct hysteresis loops in the global transport quantities. A sympathetic reader would care because this route mirrors the well-known subcritical paths in shear-driven flows and therefore supplies the first concrete bridge toward a single description of how ordered convection gives way to turbulence in buoyancy-driven systems.","feed_headline":"Buoyancy-driven flows reach turbulence subcritically","feed_subtitle":"Simulations reveal steady convection jumps to chaos and turbulence via finite disturbances, producing three hysteresis loops in transport qu","key_machinery":"The quasi-one-dimensional confinement that isolates a narrow band of Rayleigh numbers where steady convection, oscillatory chaos, and intermittent turbulence coexist and exhibit normal, reverse, and anomalous hysteresis loops.","core_discovery":"Static and quasi-static direct numerical simulations identify a narrow Rayleigh-number interval containing three coexisting states—steady convection, oscillatory chaos, and intermittent turbulence—linked by abrupt jumps and pronounced hysteresis in both Nusselt and Reynolds numbers. The steady convection state resists infinitesimal perturbations yet becomes unstable to finite-amplitude disturbances, furnishing the defining signature of a subcritical transition. This observation directly contradicts the longstanding supposition that buoyancy-driven turbulence always onsets supercritically.","pith_inferences":["The same subcritical mechanism may appear in other geometries once the confinement aspect ratio is reduced sufficiently to suppress secondary instabilities.","Engineering models of natural convection could incorporate finite-amplitude thresholds rather than relying solely on linear stability criteria.","Varying the Prandtl number across the identified window would test whether the three hysteresis loops persist or merge."],"forward_implications":["The Nusselt and Reynolds numbers exhibit discontinuous jumps accompanied by three separate hysteresis loops when the Rayleigh number is varied.","Steady convection persists under infinitesimal noise but collapses to intermittent turbulence once a sufficient finite disturbance is introduced.","The subcritical route supplies a common instability mechanism that can be compared directly with those already established for shear-driven flows.","The narrow parameter window of multistability can be used to study switching between ordered and turbulent states under controlled forcing."],"fun_headline_variants":["Subcritical turbulence in buoyancy-driven convection flows","Multiple hysteresis loops in subcritical Rayleigh-Benard flow","Steady convection transitions via finite amplitude disturbances","Three coexisting states in subcritical buoyancy transition"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The linear stability of the steady state and its response to finite-amplitude disturbances are physical properties of the buoyancy-driven flow rather than artifacts introduced by the confinement geometry or the chosen numerical methods.","fun_headline_variants_meta":{"raw":{"variants":["Subcritical turbulence in buoyancy-driven convection flows","Multiple hysteresis loops in subcritical Rayleigh-Benard flow","Steady convection transitions via finite amplitude disturbances","Three coexisting states in subcritical buoyancy transition"]},"model":"grok-4.3","cost_usd":0.009276,"raw_usage":{"total_tokens":4150,"prompt_tokens":664,"num_sources_used":0,"completion_tokens":51,"cost_in_usd_ticks":92762000,"prompt_tokens_details":{"text_tokens":664,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3435,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":664,"tokens_out":51,"duration_ms":24156,"temperature":1.0,"reasoning_tokens":3435,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-28T20:47:57.868720+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"An experiment or simulation in the same quasi-one-dimensional domain that finds only supercritical onset, no hysteresis, and no stable coexistence of the three states would falsify the subcritical claim.","supporting_citations":[],"review_version":1}