{"id":"96c64031-8bbc-4cea-ab7f-511a1fe723e0","arxiv_id":"2607.28938","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Zero-mean thermal forcing over a patterned wall drifts a liquid-vapor interface with speed proportional to wall amplitude, temperature swing, and phase, with spinodal wetting giving bounded gain.","lead":"This paper shows how a liquid-vapor interface drifting over a patterned, unevenly heated wall can move in one direction even when the heating averages to zero: the motion is set by the phase overlap between the wall pattern and the thermal drive. It then ranks which wetting conditions amplify this effect and argues that every observed drift must be balanced by an equal and opposite momentum reservoir.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Drift law is asserted, not derived: the response field M and prefactor M_h are never connected to the diffuse-interface free energy (Eq. 9), so the central quantitative claim is a phenomenological ansatz.","rationale":"The reader's verdict (CONDITIONAL) is appropriate, but the weakest assumption they identified (Section IV source calculations) is not the most load-bearing gap. The drift law itself, which is the paper's strongest quantitative claim, lacks a derivation: M is never defined in the diffuse-interface model, and M_h is not computed. This is a more fundamental problem because the bounded-gain analysis presupposes the drift law. The exact nulls and phase-selection rules are symmetry consequences and remain credible, so the paper is not worthless; it is an incomplete derivation. The requested concrete test would settle whether the drift law is a genuine consequence of the model or an imposed ansatz. If the test succeeds, CONDITIONAL remains; if it fails, the verdict should move to REJECT. Since the test is not yet done, the reader's CONDITIONAL verdict is unchanged.","tokens_in":9619,"tokens_out":5745,"duration_ms":57842,"concrete_test":"Derive the linearized response of the diffuse-interface model (9) with h(y) = h0 + ε_h sin(qy+φ) to first order in ε_h, solve the coupled phase-field and Darcy–Brinkman equations, and show that the mode-q velocity component is exactly (1/2) ε_h M_h γ_T q Θ sinφ with an explicit algebraic expression for M_h in terms of A, κ, c, h0, and the Brinkman friction. If no such expression exists, the central drift law is not a prediction of the stated model. Alternatively, run the reported numerical code with M defined directly from the ϕ field (not as an ad hoc input) and report the measured M_h from the drift amplitude; if the measured value differs by more than a few percent from the fitted value used to produce Fig. 1, the numerical support is circular.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the covariance between a response field M(y) and a zero-mean drive f(y) actually produces a physical drift velocity U in the diffuse-interface model. Section III defines the covariance abstractly, introduces the free energy (9), and then states Eq. (14) as the leading drift law without derivation. M(y) is never specified in terms of the order parameter ϕ or the wall field h(y); M_h = ∂M/∂h|_{h0} is not computed; and the step from covariance to U is not shown. The numerical checks in Fig. 1 and Table I verify the scaling and phase law, but if the simulation defines M implicitly or fits the prefactor, the test is circular. The nulls and sign reversals follow from the product form and are robust, but the quantitative proportionality U ∝ ε_h Θ sinφ is load-bearing and unsupported. Section IV's bounded-gain analysis inherits this issue because the gain factor multiplies an unproven prefactor; the later 'source calculations' (Eqs. 20–21) are also unattributed and not derived from (9).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a conservative continuum accounting scheme in which zero-mean tangential forcing can generate directed interfacial drift through the covariance between a local response field M(y) and a zero-mean drive f(y). In a diffuse-interface realization with sinusoidal wall heterogeneity and Marangoni thermal drive, it claims the drift law U ∝ ε_h Θ sin φ with exact nulls, sign reversals, and a momentum ledger that identifies the compensating reservoir. It then analyzes wetting susceptibility as a bounded gain factor: first-order spinodal conditions give finite amplification, critical wetting saturates, bulk criticality suppresses, and thermal noise caps the gain. Reduced numerical models are used to check the phase law, spinodal optimum, momentum-budget closure, chirality, and grid convergence.","tokens_in":9929,"tokens_out":4500,"duration_ms":47206,"significance":"If the central drift law were derived from the diffuse-interface free energy, the paper would offer a useful symmetry-constrained selection rule and a disciplined accounting framework for rectified interfacial transport. The covariance identity, exact nulls, sign reversal, and chiral selection rule are clean and testable, and the archived code and reproducibility measures are strengths. However, the quantitative law and the bounded-gain scalings are currently asserted rather than derived, and the numerical checks are largely self-consistency tests. The potential significance is real but contingent on closing these gaps.","major_comments":[{"comment":"The drift law is not derived. Eq. (13) postulates M[h]=M0+ε_h M_h sin(qy+φ) without specifying M in terms of the order parameter, the wall field, or the free energy (9). Eq. (14) states 'perturbation theory then gives' U=..., but neither the linear response of M to h nor the conversion of covariance into a velocity is shown. Since M_h is never connected to Eq. (9), the central quantitative claim is an ansatz, not a prediction of the diffuse-interface model. Please derive M_h from the model or explicitly label Eq. (14) as a phenomenological constitutive assumption.","section":"Sec. III, Eqs. (13)-(14)"},{"comment":"The 'source calculations' leading to K∼t^{3ν-1} and χmax∝(k_B T_abs)^{-1/3} are not presented. The text lists ingredients for Eq. (20) but does not show the derivation; Eq. (21) is justified only by the assumed saddle-node normal form (18), which is itself asserted. These scaling laws are load-bearing for the bounded-gain conclusion, so they must be derived from Eq. (9) or clearly identified as external results with explicit references.","section":"Sec. IV, Eqs. (20)-(21)"},{"comment":"The numerical verification is circular as reported. R²=1.0 and a maximum prediction error of 4.24×10^{-21} indicate that the Brinkman solver is reproducing the same linear relation that was inserted into the model, not testing an independent prediction. Please distinguish consistency checks from validation, report the residuals on an independent calculation or against a separately derived solution, and clarify whether the solver was solving equations identical to the assumed drift law.","section":"Appendix A, Table I"},{"comment":"M_h and D are never computed or bounded; they remain free parameters. The scaling laws and gain discussion involve products of these unknowns. Even if the nulls and phase reversals are robust, the predicted amplitude is not testable until M_h and D are estimated from the model or constrained by independent data. Please provide expressions or bounds for these quantities.","section":"Sec. III-IV, free parameters"}],"minor_comments":[{"comment":"The ledger includes EM, radiation, and '···' terms that are never used in the paper. Clarify that these are placeholders and how they would be closed, or drop them to avoid implying a completeness that is not demonstrated.","section":"Sec. II, Eq. (4)"},{"comment":"The axis labels use symbols (V_drift, ε_h, φ) without defining dimensionless units. Please state the normalization and parameter values in the captions.","section":"Fig. 1 and Fig. 2"},{"comment":"The phrase 'source calculations' is vague; either present those calculations in a supplementary appendix or cite a specific reference for each result.","section":"Sec. IV"},{"comment":"The paper would benefit from a glossary defining M, M_h, ε_h, Θ, φ, and K_X early in Section III, since the readability of the later scaling arguments depends on these symbols.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper has a useful conceptual core, but the main technical claim (Eq. 14) is not derived, and the numerical tables overstate verification. I would be open to a revision that supplies the missing derivations and re-scopes the numerical results as consistency checks rather than independent tests. The topic is within the journal's scope, and the selection-rule framework is potentially valuable, but the current form is not sufficient for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper is a clean synthesis, and the honesty about scope is real. The announced drift law U∝ε_h Θ sinφ is the sort of thing you'd expect from a symmetry argument, and the numerical checks confirm it, but the derivation connecting the covariance to a physical velocity is skipped; M_h is introduced as a response coefficient and never computed from the free energy of Eq. (9). That is the one load-bearing gap.\n\nWhat is actually new: the wetting-susceptibility gain ranking (spinodal amplification, critical-wetting saturation, bulk-criticality suppression) applied to Marangoni rectification. The explicit momentum ledger, treating prior work as constraints, is a useful discipline. The numerical appendix is reproducible, code archived, and the budget closure to 1e-16 is solid. The paper earns credit for not overclaiming: it repeatedly calls itself an accounting scheme and lays out falsification protocols.\n\nSoft spots: the prefactor. The paper says 'perturbation theory then gives' Eq. (14) but the derivation is not shown. To be fair, the phase law and nulls follow from the sinusoidal form and symmetry; the prefactor is what's unproven. The stress-test also flags the Section IV scaling laws as allegedly unattributed, but on reading, the exponents are sketched in the text (e.g., they derive K~t^{3ν−1} from scaling of γ, qξ, and I). The one-third thermal cap is given with a short explanation, though the constants like n are not specified. So that criticism is overblown; the real gap is M_h (and D, the drive that reaches the interface) as free parameters.\n\nDoes the central argument hold up? As an accounting framework, yes. As a quantitative prediction, no, unless M_h is computed from the model or measured. The paper's own scope statement tempers this, so it's a moderate limitation, not a fatal one.\n\nWho's it for: people working on thermocapillary pumping, patterned wetting, or soft-matter ratchets. Good reading group material. It deserves a serious referee: it's coherent, testable, and provocative in a useful way. The referee should ask for either a derivation of M_h from Eq. (9) or an explicit statement that M_h is a phenomenological coefficient fitted to the reduced numerics, plus a bit more detail on the thermal-cap normalization.\n\nRecommendation: send to review; conditional accept after a revision that either derives or demotes the prefactor.","headline":"A transparent accounting framework for interfacial rectification, with robust selection rules but a quantitative drift law whose amplitude coefficient is never derived.","tokens_in":10401,"tokens_out":4175,"would_cite":false,"duration_ms":41448,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["68.03.-g","47.55.-t"],"model":"deepseek-v4-flash","headline":"The paper argues that directed transport of liquid–vapor interfaces under zero-mean forcing is controlled by the covariance between a local response field and the drive, yielding the drift law U∝ε_hΘ sinφ with exact symmetry nulls, and that","keywords":["momentum rectification","covariance","Marangoni drift","wetting transitions","liquid-vapor interface","diffuse-interface model","zero-mean forcing","bounded gain"],"falsifier":"Measure the drift speed of a liquid film over a substrate with sinusoidal wetting heterogeneity h(y)=h0+ε_h sin(qy+φ) and sinusoidal thermal profile T(y)=Θ sin(qy). The law U∝ sin φ requires exact nulls at φ=0 and π and sign reversal across φ=0; any off-set drift at null phase, or a drift that does not vanish as Θ→0, would falsify the covariance mechanism.","tokens_in":9528,"feed_emoji":"💧","tokens_out":4008,"duration_ms":36176,"temperature":0.7,"pith_summary":"This paper seeks to establish that zero-mean tangential forcing can still drive a liquid–vapor interface through a covariance effect: the overlap between a structured wall response and a structured thermal drive, not the mean drive, determines the drift. In a minimal diffuse-interface model, the leading drift is U∝ε_hΘ sinφ — proportional to wall heterogeneity, thermal amplitude, and the sine of their phase shift — with exact nulls when either factor vanishes or the phase inverts. Wetting susceptibility acts as a bounded gain: first-order spinodal conditions amplify, critical wetting saturates because the interface unbinds from the drive, and bulk criticality suppresses the channel as the interface dissolves. The paper frames itself as a symmetry-constrained momentum accounting scheme, identifying compensating momentum reservoirs for every rectified channel.","feed_headline":"Covariance, not mean force, drives liquid–vapor drift","feed_subtitle":"A zero-mean thermal pattern above a patterned wall yields U∝εΘ sinφ; wetting physics bounds the gain.","key_machinery":"The central object is the covariance ⟨Mf⟩−⟨M⟩⟨f⟩ between a local response field M(y) (mobility, wetting strength, slip, compliance) and a zero-mean tangential drive f(y). In a diffuse-interface free energy F=∫[A/4(φ^2−1)^2+κ/2|∇φ|^2] + Robin wall term, the covariance yields a Marangoni traction f_M=γ_T(I−nn)∇T δ_ξ; the drift law follows from perturbation theory with M[h(y)]=M0+ε_h M_h sin(qy+φ). The susceptibility χ_q=D/(γ q^2+K_X) links covariance to gain, with K_X from the saddle-node normal form V(ℓ)=ℓ^3/3−Δℓ.","core_discovery":"The central discovery is that covariance-driven rectification gives a phase-selective drift law: for a Marangoni-driven liquid–vapor interface above a weakly structured wall with wetting-strength modulation h(y)=h0+ε_h sin(qy+φ) and zero-mean temperature variation T(y)=Θ sin(qy), the drift is U=(1/2)ε_h M_h γ_T q Θ sin φ, with exact nulls at ε_h=0, Θ=0, and sign reversal under φ→−φ or Θ→−Θ. Wetting susceptibility near a first-order spinodal diverges as χ∝Δ^{-1/2}, providing a bounded amplification regime, while critical wetting saturates because the unbinding interface loses its short-range drive, and bulk criticality suppresses the channel as K∼t^{3ν−1}.","pith_inferences":["If the covariance law is generic, similar rectification could arise in other driven soft-interface contexts (e.g., chemical patterning, electric fields) where a local response field correlates with a zero-mean drive; the paper only sketches these reservoirs.","The saddle-node normal form assumption may understate the role of non-mean-field fluctuation corrections in real water; testing the predicted U∝ε_hΘ sinφ at finite temperature and near contact-line depinning could reveal whether the bounded-gain result survives beyond the reduced model.","The paper's accounting scheme might be applied to biomolecular or colloidal transport where momentum reservoirs are often ignored; the budget-closure discipline is a transferable methodology.","A direct experimental falsifier: measure drift vs phase φ on a patterned substrate with controlled thermal wave; the law requires pure sine with nulls — any reproducible deviation (e.g., quadratic component) would falsify the mechanism."],"forward_implications":["If this is correct, a zero-mean thermal pattern above a heterogeneous wall can produce a predictable directed interface drift whose sign is controlled by the phase shift, enabling phase-selective microfluidic pumping.","The nulls (ε_h=0, Θ=0, phase inversion) give explicit falsification handles: no drift without heterogeneity or drive, and sign reversal under inversion.","Wetting transitions organize a gain ladder: first-order spinodal is the exploitable regime, critical wetting saturates, bulk criticality suppresses; so design should target spinodal-like conditions while maintaining wall coupling.","Momentum budget closure mandates that any measured drift be matched by a compensating reservoir (vapor, waves, wall stress, etc.), which diagnoses hidden channels.","Transverse chirality rule ⟨f_z⟩∝c(b^2−a^2)cos(φ1+φ2) offers a 3D test with mirror reversal."],"fun_headline_variants":["Covariance rectifies zero-mean forces at liquid-vapor walls","Drift law U∝εΘ sinφ emerges from covariance","Wetting susceptibility controls interfacial drift gain","Zero-mean drive, directed flow via covariance"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The bounded-gain and criticality suppression conclusions rest on the asserted saddle-node normal form for the binding potential near a first-order spinodal and on the source calculations for the thermal cap and K∼t^{3ν−1}; if those are not valid, the useful spinodal amplification is unsupported.","fun_headline_variants_meta":{"raw":{"variants":["Covariance rectifies zero-mean forces at liquid-vapor walls","Drift law U∝εΘ sinφ emerges from covariance","Wetting susceptibility controls interfacial drift gain","Zero-mean drive, directed flow via covariance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000621,"raw_usage":{"total_tokens":2777,"prompt_tokens":866,"completion_tokens":1911,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":1854}},"tokens_in":610,"tokens_out":1911,"duration_ms":14548,"temperature":1.0,"reasoning_tokens":1854,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T16:51:33.514613+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the drift speed of a liquid film over a substrate with sinusoidal wetting heterogeneity h(y)=h0+ε_h sin(qy+φ) and sinusoidal thermal profile T(y)=Θ sin(qy). The law U∝ sin φ requires exact nulls at φ=0 and π and sign reversal across φ=0; any off-set drift at null phase, or a drift that does not vanish as Θ→0, would falsify the covariance mechanism.","supporting_citations":[],"review_version":1}