{"id":"6a4227b7-ab47-4add-9c72-f733b7a9e1ed","arxiv_id":"2411.15440","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"A product of height-averaged damping and forcing is introduced as a universal ranking parameter for capillary rise, and a genetic algorithm using it proposes channel shapes with faster modeled wicking.","lead":"This paper models capillary rise in wavy channels as a damped linear system and introduces a single parameter, the product of damping and forcing, to rank how fast different channels wick. It then uses that parameter in a genetic algorithm to propose channel shapes that rise faster than earlier designs.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zeta = <D>·<F> is not a sufficient statistic for Eq. (4): even constant-coefficient examples with equal zeta have different rise curves, and lower-zeta cases can have larger initial slopes, so the monotonic ranking and Eq. (6) objective are unsupported.","rationale":"The reader identified the same load-bearing premise: the collapse of D(h),F(h) to height-averaged values whose product orders dynamics. My analysis strengthens this: even in the paper's own constant-coefficient toy model, the product is not an order-preserving statistic. Consequently, the asserted monotonicity of zeta is not established by Fig. 2, and the optimization objective in Eq. (6) may maximize a proxy that does not correspond to faster rise. This is an internal inconsistency of the central claim, not merely a disagreement with prior literature. I nevertheless note the paper's useful pieces: the geometric parametrization is flexible, the numerical solution of Eq. (4) is a legitimate way to explore the design space, and a scalar objective is a sensible idea for stochastic search. The flaw is in replacing the true objective by zeta without proof. The reader's REJECT verdict is appropriate; I would keep it unchanged, though CONDITIONAL would be defensible if the authors supply the missing comparison, e.g., direct optimization against rise time and a monotonicity check over a dense grid.","tokens_in":6718,"tokens_out":3635,"duration_ms":32537,"concrete_test":"Run the constant-coefficient check: solve dh/dt = -D h + F for (D,F) = (10,10), (1,100), (1000,100), (1,1000) and compare time to reach a series of target heights. If equal-zeta curves differ or a lower-zeta case reaches any target sooner, the universal ranking claim is refuted. Then rerun the genetic algorithm of Eq. (6) with the direct objective 'minimize time to reach h*' for several h*, holding the same design space and contact angle; if the resulting optimal geometries differ from the zeta-optimized one, or the zeta-optimized geometry is not uniformly faster, the proxy objective is invalid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing claim is that zeta = <D> · <F> is a universal monotonic descriptor and that maximizing it yields faster rise. The paper's motivating ODE dh/dt = -D h + F undermines this: the solution is h(t) = F/D (1 - e^{-Dt}), which depends on D and F separately, not only on their product. Two cases with identical zeta = D F can differ wildly: (D,F) = (10,10) and (1,100) both have zeta = 100, but their equilibrium heights are 1 and 100 and their relaxation times differ by 10x. Even worse for monotonicity, (D,F) = (1000,100) has zeta = 10^5 but initial slope h'(0) = 100, while (D,F) = (1,1000) has zeta = 10^3 and initial slope h'(0) = 1000: the lower-zeta system rises faster initially and reaches moderate heights sooner. Thus the premise 'larger zeta means faster rise' is false in exactly the constant case used as motivation. For the actual h-dependent D(h),F(h), replacing them by height averages is therefore not justified by the paper's argument; Fig. 2 is presented as evidence but only demonstrates a correlation selected post hoc, not a monotonic relationship. Since Eq. (6) maximizes this proxy scalar, the optimized channel and the claimed ~30% improvement are not shown to be optimal for any concrete objective such as time to reach a target height. This is internal to the paper's own model, not a disagreement with external consensus.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a one-dimensional energy-balance model for capillary rise in axisymmetric channels with slowly varying radius, written in compact form as dh/dt = -D(h)h + F(h). It introduces a height-averaged product zeta = <D><F> as a universal parameter that is claimed to order rise dynamics from wicking to anti-wicking, and it uses a genetic algorithm with objective J = zeta to optimize channel geometries. The optimized geometries are compared with the prior designs of Gorce et al. [7] and Figliuzzi and Buie [15], and the paper reports faster rise and a roughly 30% increase in zeta. The central claim is that this single parameter characterizes rise dynamics across geometric, wetting, and pressure conditions.","tokens_in":6992,"tokens_out":6472,"duration_ms":60595,"significance":"If valid, a single parameter that ranks capillary-rise performance across geometry, wettability, and applied pressure would be practically useful and would unify disparate design rules. The paper has real strengths: Eq. (4) is derived from an energy-balance formulation rather than fitted ad hoc, the model is benchmarked against external geometries from [7] and [15], and the optimized shapes are concrete falsifiable predictions. However, the central universal-parameter claim is not supported by the evidence in the manuscript. The motivating constant-coefficient system itself shows that the product D*F is not a sufficient statistic for the rise curve, and no proof or controlled numerical experiment is given that replacing D(h) and F(h) by height averages preserves the ordering of rise dynamics. Because the optimization objective in Eq. (6) is exactly this product, the reported design improvements are not tied to a well-defined rise-height or rise-time objective. The core contribution therefore needs substantial reformulation, not just local revision.","major_comments":[{"comment":"The claim that eta = D*F controls rise behavior is false for the very equation used to motivate it. For dh/dt = -D h + F with h(0) = 0, the solution is h(t) = (F/D)(1 - e^{-Dt}), so the equilibrium height is F/D and the relaxation time is 1/D; these are independent pieces of information, not functions of eta only. Equal-eta pairs such as (D,F) = (10,10) and (1,100) have the same eta = 100 but equilibrium heights 1 and 100 and relaxation times 0.1 and 1. Lower-eta cases can also rise faster initially: (D,F) = (1,1000) has eta = 10^3 and initial slope 1000, while (D,F) = (1000,100) has eta = 10^5 and initial slope 100. The text's statements that higher eta gives faster initial rise and earlier equilibration are therefore contradicted by the manuscript's own model, and this undermines the analogy used to justify zeta.","section":"After Eq. (5), constant-coefficient motivation"},{"comment":"The step from the constant-coefficient discussion to the height-averaged product zeta = <D><F> is not justified. Equation (4) has h-dependent D and F; replacing them by their averages changes the solution, and no argument is given to show that geometries with larger zeta rise faster at every time or even at a fixed target height. Fig. 2(a) is a colored ensemble of trajectories, not a controlled test: the curves vary simultaneously in D(h), F(h), and zeta, so the visual correlation does not establish the claimed monotonic relation. The statement that 'the monotonic relationship between rise height and zeta across the entire design space' is observed would need at minimum a quantitative measure, such as time to a fixed height as a function of zeta with independent variation of D and F, before it can support optimization.","section":"Paragraph defining zeta after Eq. (5)"},{"comment":"Because the cost function J = <D><F> is the very quantity whose validity as a ranking scalar is unestablished, the optimization results do not show that the resulting geometries achieve faster rise for any concrete performance measure. The text reports that the optimized curves 'surpass the heights predicted by [7,15]' and a '~30% increase', but the 30% figure refers to zeta, not to a rise-height or rise-time improvement. In the constant-coefficient case, maximizing zeta can select a large-F/large-D geometry that equilibrates quickly at a low height, so the design recommended by Eq. (6) is not necessarily the one that reaches a prescribed height first. The claim that these geometries are optimal is therefore not supported by the stated objective.","section":"Eq. (6), optimization objective"},{"comment":"Equation (2) as printed contains delta-dot-E on both sides and appears dimensionally inconsistent; as the stated basis for Eq. (4), it must be rewritten with all terms explicitly defined. The gravitational potential derivative should presumably involve dh/dt, but none appears in the first term as written. Without a corrected statement of this balance, the reduction to Eq. (4) cannot be checked.","section":"Eq. (2)"}],"minor_comments":[{"comment":"The coefficients P_A and P_B in the pressure ramp P = P_A + P_B t are introduced but never given values or units, so the hydrophobic-pressure cases in Fig. 2 are not reproducible.","section":"Pressure ramp definition, after Eq. (4)"},{"comment":"The phrase 'reduced effective viscosity' should be defined quantitatively; the non-dimensional mu-bar is an integral measure of channel resistance, not a physical viscosity, and the text should state exactly what is plotted.","section":"Fig. 2 caption"},{"comment":"The calibration of the geometric parameters R_hat, R_tilde, and lambda to the power-law profile is not described quantitatively; please report the fitted values and the fitting procedure so the benchmark comparison is reproducible.","section":"Calibration paragraph, near R(h) = (1-h/l)^(5/6)"},{"comment":"The manuscript refers to Supplementary Material for numerical validation, but the validation details are not included in the text; either include the supplementary material in the submission or summarize the validation in the paper.","section":"Reference [16]"}],"recommendation":"reject","confidential_remarks":"The central universal-parameter claim fails on the manuscript's own constant-coefficient example, and the averaging step from D(h), F(h) to zeta is asserted rather than proved. Since the optimization and the claimed advances over [7,15] depend entirely on that parameter, I do not see a local fix that preserves the paper's stated contribution; a substantial reformulation would be required. The manuscript would also need a corrected Eq. (2) and quantitative reporting of the optimization objective before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper's central parameter zeta = <D>·<F> does not do what it claims. The constant-coefficient equation used as motivation gives h(t) = F/D (1 - e^{-Dt}), which depends on D and F separately, not only on their product. Two cases with identical zeta can have equilibrium heights a factor of 100 apart, and a lower-zeta system can have a larger initial slope. So the monotonic ranking and the GA objective in Eq. (6) are unsupported. This is a load-bearing flaw, not a nit.\n\nWhat's genuinely useful: the paper sets up a compact first-order framework dh/dt = -D(h)h + F(h) for axisymmetric channels, and it covers a wide design space including hydrophilic and hydrophobic cases with external pressure. The idea of using a scalar objective for genetic optimization of channel shape is attractive, and benchmarking against the geometries of [7] and [15] is the right way to calibrate. The optimized shapes do produce faster modeled rise in the full equation, which suggests the search itself might be fine; the problem is the claimed justification via zeta.\n\nWhere it falls apart: the leap from constant coefficients to height-averaged coefficients is never justified. The paper asserts monotonicity from Fig. 2 without a proof, and the stress-test counterexample shows that even for the simplest case, higher product does not imply faster rise. The averaging step could easily reorder cases. Also, the derivation from Eqs. (1)-(2) to (5) is too compressed; there are apparent dimensional and sign issues (the OCR doesn't help, but the published text reads garbled in places). No experimental validation is given, and no code or data are released, which makes the optimization hard to reproduce.\n\nBottom line: the framework is standard Lucas-Washburn with a new averaging gimmick, and the gimmick fails on inspection. The paper is not ready for peer review as written. A serious revision would need to either prove a ranking theorem for zeta under the actual h-dependent coefficients, or replace the objective with something concretely tied to rise time (e.g., time to reach 90% of equilibrium). If the authors did that, the optimization study might be worth revisiting. As is, I'd desk-reject.","headline":"Zeta = <D>·<F> is not a sufficient statistic, so the paper's central universal-parameter claim and optimization objective are unsupported; desk-reject.","tokens_in":7593,"tokens_out":3501,"would_cite":false,"duration_ms":31440,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single parameter, the product of damping and forcing, orders capillary rise dynamics across channel geometries.","keywords":["capillary rise","wicking","anti-wicking","damped system","genetic algorithm","axisymmetric channels","lubrication approximation","microfluidics"],"falsifier":"Simulate or measure, for a pair of channels deliberately constructed to have equal $\\zeta$ but strongly different $D(h)$ and $F(h)$ profiles, the full rise curves from Eq. (4) or from experiment; if the curves differ in ordering or equilibration, the averaging collapse fails and the genetic-algorithm objective in Eq. (6) is not a faithful proxy for rise speed.","tokens_in":6408,"feed_emoji":"💧","tokens_out":6113,"duration_ms":50366,"temperature":0.7,"pith_summary":"Capillary rise in a channel is usually thought to depend separately on geometry, wettability, and applied pressure. This paper argues that for axisymmetric channels with gradually varying radius, the front height follows a linear damped-system equation, and that a single number—the height-averaged product of damping and forcing, $\\zeta = \\langle D\\rangle \\langle F\\rangle$—ranks the whole spectrum from fast wicking to strongly slowed (anti-wicking) rise. Because higher $\\zeta$ consistently predicts faster initial rise and earlier equilibration in the solutions the authors compute, they use $\\zeta$ as the objective in a genetic algorithm that searches over channel shapes. The optimized shapes are reported to rise faster and to enter the viscous regime earlier than the two previous best-known channel geometries, with $\\zeta$ about 30% higher under the same normalization.","feed_headline":"One number ranks capillary rise from wicking to anti-wicking","feed_subtitle":"The product of damping and forcing predicts rise speed and guides an optimizer to faster channel shapes.","key_machinery":"The load-bearing identity is the damped-system form of the governing equation, $\\frac{d\\bar{h}}{d\\bar{t}} = -D(\\bar{h})\\bar{h} + F(\\bar{h})$, where the damping coefficient $D$ gathers viscous resistance, gravitational effects, and geometric slope, and the forcing $F$ gathers capillary driving plus any external pressure ramp. The new object is $\\zeta = \\langle D\\rangle \\cdot \\langle F\\rangle$, the product of the height-averaged coefficients, introduced by analogy with $\\eta = D\\cdot F$ for the constant-coefficient equation. The genetic algorithm in Eq. (6) maximizes $J(x) = \\zeta$ over the shape parameters of the radius profile $R(x) = \\hat{R} + \\tilde{R}\\cos(n\\pi/2 - 2\\pi x/\\lambda)$, and the monotonic $\\zeta$–height relationship is what allows the optimizer to sidestep solving the full integro-differential system at every step.","core_discovery":"The paper claims that the dimensionless front dynamics reduce to $\\frac{d\\bar{h}}{d\\bar{t}} = -D(\\bar{h})\\bar{h} + F(\\bar{h})$, with $D$ and $F$ expressed in terms of channel radius, slope, contact angle, density contrast, viscosity, and external pressure. The central discovery is that the height-averaged product $\\zeta = \\langle D\\rangle \\langle F\\rangle$ behaves like the constant-coefficient product $\\eta = D\\cdot F$ in a linear damped system: the authors find a monotonic relation between $\\zeta$ and rise height across hydrophilic and hydrophobic channels, with and without external pressure. Using $\\zeta$ as the fitness function $J(x) = \\langle D(x)\\rangle \\cdot \\langle F(x)\\rangle$ in a genetic algorithm, the paper reports new channel geometries whose modeled rise exceeds that of the power-law optimal channel of Gorce et al. and the polynomial optimal channel of Figliuzzi and Buie. The conclusion is that a single scalar can characterize both wicking enhancement and inhibition in the same framework.","pith_inferences":["The paper's constant-coefficient argument shows that the product $\\eta = D\\cdot F$ matters for the analogue, but it does not prove that the height-averaged product $\\zeta$ preserves the ordering for the full variable-coefficient problem; that is the assumption a reader should test.","Since $\\zeta$ combines two quantities that enter the solution separately (equilibrium height is $F/D$ and relaxation rate is $D$), designs with identical $\\zeta$ could have different balances of equilibrium and speed; the optimization would then be selecting a proxy without separating these effects.","One could test the framework further by fabricating two channels with equal $\\zeta$ but different $D$ and $F$ profiles and measuring whether the rise curves collapse, which would confirm averaging is legitimate beyond the simulation cases shown.","The radius parametrization is like a truncated Fourier series; extending it with more terms may let the same optimizer find even faster shapes, but the gain would be bounded by the validity of the lubrication approximation."],"forward_implications":["If $\\zeta$ is truly universal, then channel design for either fast wicking or deliberate flow inhibition reduces to maximizing or minimizing a single scalar, rather than solving the full integro-differential equation.","The optimized geometries produced by the genetic algorithm raise the modeled liquid height beyond the previously reported power-law and polynomial optimal channels, implying the prior shapes were not global optima in this design space.","The damped-system form suggests a direct analogy with electrical RC circuits, where the product of damping and forcing plays the role of a dynamic response parameter, potentially enabling circuit-inspired intuition for microfluidic transport.","Because $\\zeta$ also tracks the onset of the viscous regime, it can be used to predict when the $h \\propto \\sqrt{t}$ behavior begins, which matters for timing-sensitive applications such as paper-based diagnostics."],"supporting_citations":[{"why":"It supplies the Lucas-Washburn equation, the classical capillary-rise model that this damped-system formulation extends.","marker":"[1]"},{"why":"It supplies the Washburn version of the capillary-rise equation and the classical $h \\propto \\sqrt{t}$ scaling that the paper generalizes.","marker":"[2]"},{"why":"It provides the power-law optimal channel geometry used to fit the radius model and is the main baseline that the optimized designs are claimed to surpass.","marker":"[7]"},{"why":"It demonstrates that capillary geometry can alter the classical power-law exponents, motivating the search for a parameter that captures geometry in rise dynamics.","marker":"[14]"},{"why":"It provides the polynomial-optimal channel geometry and the prior optimization approach that this work compares against and claims to outperform.","marker":"[15]"},{"why":"It supplies the numerical method and validation for the solution of the governing equation, supporting the computed rise curves.","marker":"[16]"},{"why":"It supplies the genetic-algorithm method used in the stochastic optimizer that maximizes the proposed parameter.","marker":"[17]"}],"fun_headline_variants":["A scalar product captures wicking and anti-wicking in capillary rise","Capillary dynamics: one number from damping and forcing","Wicking to anti-wicking: a universal single-scalar framework","Optimize channel geometry with the damping-forcing product"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that replacing the height-dependent damping and forcing with their averaged values keeps the ranking of rise speeds intact, so that the single product $\\zeta = \\langle D\\rangle \\langle F\\rangle$ is enough to compare channels; if two channels with the same $\\zeta$ but very different $D$ and $F$ profiles rise at different rates, the framework loses its predictive power.","fun_headline_variants_meta":{"raw":{"variants":["A scalar product captures wicking and anti-wicking in capillary rise","Capillary dynamics: one number from damping and forcing","Wicking to anti-wicking: a universal single-scalar framework","Optimize channel geometry with the damping-forcing product"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001503,"raw_usage":{"total_tokens":6007,"prompt_tokens":900,"completion_tokens":5107,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":516,"completion_tokens_details":{"reasoning_tokens":5037}},"tokens_in":516,"tokens_out":5107,"duration_ms":30503,"temperature":1.0,"reasoning_tokens":5037,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:18:27.992502+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate or measure, for a pair of channels deliberately constructed to have equal $\\zeta$ but strongly different $D(h)$ and $F(h)$ profiles, the full rise curves from Eq. (4) or from experiment; if the curves differ in ordering or equilibration, the averaging collapse fails and the genetic-algorithm objective in Eq. (6) is not a faithful proxy for rise speed.","supporting_citations":[{"cited_title":"Lucas, Kolloid-Zeitschrift 23, 15 (1918)","cited_arxiv_id":null,"evidence_quote":"It supplies the Lucas-Washburn equation, the classical capillary-rise model that this damped-system formulation extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the Washburn version of the capillary-rise equation and the classical $h \\propto \\sqrt{t}$ scaling that the paper generalizes."},{"cited_title":"Gorce, I","cited_arxiv_id":null,"evidence_quote":"It provides the power-law optimal channel geometry used to fit the radius model and is the main baseline that the optimized designs are claimed to surpass."},{"cited_title":"Reyssat, L","cited_arxiv_id":null,"evidence_quote":"It demonstrates that capillary geometry can alter the classical power-law exponents, motivating the search for a parameter that captures geometry in rise dynamics."},{"cited_title":"Figliuzzi and C","cited_arxiv_id":null,"evidence_quote":"It provides the polynomial-optimal channel geometry and the prior optimization approach that this work compares against and claims to outperform."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the numerical method and validation for the solution of the governing equation, supporting the computed rise curves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the genetic-algorithm method used in the stochastic optimizer that maximizes the proposed parameter."}],"review_version":1}