{"id":"6dbd8e82-4224-4d7b-a8f9-31ea072a955a","arxiv_id":"1908.07971","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For 2D drops pinned on sharp edges, stability under constant-pressure and constant-volume perturbations is set by turning points in the pressure-area plane, giving size-dependent critical pinning angles.","lead":"This paper maps how the stability of two-dimensional liquid drops pinned on sharp edges depends on drop size, measured by the Bond number. It shows that closed drops become more stable as they shrink while open drops become less stable, which matters for designing superhydrophobic surfaces.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability boundaries rest on unverified absence of pre-fold eigenvalue crossings; direct second-variation check needed before accepting the P0-a turning-point map.","rationale":"The reader correctly identified the reliance on Maddocks' turning-point theorem without a direct second-variation check as the weakest assumption. I agree that this is the load-bearing point: the paper's quantitative stability boundaries, and especially the abstract's statement about size-dependent stability, are only as reliable as the theorem's applicability to the specific branch. My stress-test sharpens the reader's concern in two ways. First, the theorem requires that no other eigenvalue of the full Hessian vanish before the fold; the paper provides no spectral check. Second, the half-domain formulation only captures symmetric perturbations, while the phrase 'planar pinned perturbations' does not by itself exclude antisymmetric in-plane modes; these are not discussed. The paper does explicitly acknowledge the exclusion of longitudinal Rayleigh-Plateau modes, which is an honest limitation, but the planar antisymmetric case is not acknowledged. Because the missing check is concrete and numerical, and because the central claims are otherwise coherent and consistent with known Bo=0 results, I would not reject the paper. I would recommend acceptance conditional on the eigenvalue verification described above. If that verification passes, the turning-point argument is sound and the verdict should be ACCEPT; if it fails, the stability map would need to be recomputed.","tokens_in":8224,"tokens_out":15404,"duration_ms":173826,"concrete_test":"For a representative Bond number, e.g. Bo=1, compute the linear stability of the equilibrium branch over the full range of pinning angles by discretizing the second variation of the energy functional (1). Include both symmetric perturbations satisfying ζ_dot(sa)=0 and antisymmetric perturbations satisfying ζ(sa)=0 at the symmetry axis, with pinned endpoints at the edges, and track the smallest eigenvalue of each symmetry class as a function of θ0. Verify that the first zero crossing coincides with the green constant-pressure fold (θ0≈240° for Bo=1) and the yellow constant-volume fold (θ0≈266° for Bo=1). If an eigenvalue crosses before these folds, or if an antisymmetric mode gives a lower threshold, the stability map and the associated conclusions in the abstract need revision; if no earlier crossing exists, the turning-point argument is confirmed.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that stability changes exactly at extrema of P0(a) and a(P0), via Maddocks' turning-point theorem as applied in Section II.B. The theorem is invoked without verifying its hypotheses for this infinite-dimensional, free-endpoint problem: no direct computation of the second variation is reported, and only planar, pinned perturbations are considered. If any symmetric or antisymmetric eigenmode crosses zero before the first fold in the P0-a branch, the critical angles in Fig. 3 and the abstract's size-stability statement would be wrong. The paper explicitly excludes longitudinal Rayleigh-Plateau modes, and the half-domain formulation with the symmetry condition ζ_dot(sa)=0 means antisymmetric planar modes are not analyzed; such modes could, in principle, become unstable earlier. Maddocks' theorem identifies stability changes at folds only when no other degeneracy occurs, so the load-bearing condition is that the lowest eigenvalue of the full Hessian vanishes at the folds and nowhere earlier. This condition is plausible for the branches shown, but it is not checked, and the branch ends at singular limits (θ0→0 and θ0→360) where the zero-area configuration may carry additional zero modes.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the equilibrium and stability of two-dimensional liquid drops pinned on sharp edges under gravity. The authors introduce a variational formulation with a pressure eigenvalue P0, solve the Euler-Lagrange equations for a range of Bond numbers, and use turning-point theory (Maddocks' theorem) to claim that stability changes occur at extrema of P0(a) for open (constant-pressure) drops and at extrema of a(P0) for fixed-volume drops. They construct a stability map in the Bond-number/pinning-angle plane, apply the results to a multiscale model of superhydrophobic substrates, and draw qualitative conclusions about how drop size affects stability.","tokens_in":8397,"tokens_out":14827,"duration_ms":134432,"significance":"If the stability boundaries are correct, the paper offers an elegant, parameter-free explanation of how drop size (Bond number) controls the stability of pinned interfaces, which is directly relevant to superhydrophobicity and droplet microfluidics. The variational turning-point approach is a standard and powerful tool, and the multiscale discussion makes useful contact with experimentally observed wetting transitions. However, the central stability map rests on assumptions about the second variation that are not verified, and some of the abstract's size-stability claims are not derived explicitly. These issues temper the otherwise significant contribution.","major_comments":[{"comment":"The stability analysis is performed on the half-domain with the symmetry boundary condition ζ_dot(sa)=0, which restricts the admissible perturbations to those symmetric about the vertical axis. Antisymmetric planar perturbations (which would have a nonzero vertical component at the symmetry plane) are not captured by this formulation. The turning-point method identifies stability changes only at folds of the P0-a branch; a symmetry-breaking bifurcation or any other eigenmode crossing zero at a point that is not a fold would not be detected. Since the stability boundaries in Fig. 3 are the central result of the paper, the authors must either verify directly (e.g., by discretizing the second variation and checking positive definiteness) that no eigenmode of the full planar problem becomes unstable before the folds, or explicitly restrict and justify the analysis to symmetric perturbations.","section":"Section II.B, Eq. (2e)"},{"comment":"The statements \"Drops with a fixed volume become more stable as they shrink in size\" and \"open drops... are less stable as their associated Bond number decreases\" are not derived in the manuscript. The stability map in Fig. 3 gives critical angles as functions of Bo, but it does not track drops of fixed physical volume as size changes; the dimensionless area a = A/(2Wd²) varies when Wd changes at fixed A. As presented, the claims appear inconsistent with the monotonic rise of the stability boundaries with Bo seen in Fig. 2 (e.g., the yellow constant-volume limit increases from θ0=90° at Bo=0 to θ0=180° at large Bo). The paper should provide the missing iso-volume analysis or qualify these claims.","section":"Abstract and Section III"}],"minor_comments":[{"comment":"The statement that P0 is \"a one-to-one function of the pinning angle\" is inaccurate when P0(θ0) has extrema, as in Fig. 2; the authors likely mean \"single-valued.\" Please rephrase.","section":"Section I.A"},{"comment":"The sentence \"the stability can be studied in the preferred plane P0 − (−E_P0)\" is confusing; it would be clearer to say \"the P0-a plane, where E_P0 = -a.\"","section":"Section II.B"},{"comment":"The paper does not describe the numerical method used to solve the boundary value problem (2) and generate the equilibrium branches; adding a brief description would help reproducibility.","section":"Figures 2 and 4"},{"comment":"There are minor typographical and reference issues, including a duplicated journal name in reference [14] and missing apostrophes in \"drop's\" in a few places.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The main concern is the unverified absence of pre-fold eigenmode crossings and the restriction to symmetric perturbations; this directly affects the reliability of the stability map, which is the paper's central claim. The abstract's size-stability statements also need to be reconciled with the presented stability boundaries. These issues are fixable through additional numerical checks and more careful framing, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThis is a clean, useful paper on the equilibrium and stability of two-dimensional drops pinned on sharp edges. The main new result is the stability map in the Bond-number–pinning-angle plane, with constant-pressure and constant-volume stability limits separated. At Bo=0 they recover Speth and Lauga's result, which is a good check. The method is the standard turning-point argument: stability changes at folds of the pressure-area relation. They invoke Maddocks' theorem rather than proving the second variation directly, and that is the biggest soft spot. The theorem is not a black box; you need to know that no other eigenvalue crosses before the fold. The paper does not verify that. It also explicitly excludes longitudinal Rayleigh-Plateau modes and, because they use the half-domain with a symmetry condition, antisymmetric planar modes are not analyzed. The authors are transparent about the first exclusion but silent on the second. That said, I think the stability map is probably right for the class of perturbations considered. The turning-point method is well established in this community, and the Bo=0 agreement is reassuring. It would be better if they had reported a direct check of the second variation for a few representative branches, or provided code/data. But for a theoretical paper, this is a reasonable level of rigor. The multiscale discussion in Section III is more qualitative, but it ties the results to superhydrophobicity in a way that is suggestive rather than rigorous. The citation pattern is fine. Overall, I'd send this to a referee who knows capillary stability. It deserves a serious look, not a desk reject. If the referee presses on the Maddocks hypothesis, the authors may need to add a few numerical eigenvalue calculations, but that's a revision, not a rejection.","headline":"A solid stability map for 2D pinned drops, with the main caveat that the turning-point theorem is invoked rather than verified.","tokens_in":8930,"tokens_out":4894,"would_cite":false,"duration_ms":46466,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.55.D-"],"model":"deepseek-v4-flash","headline":"The stability of two-dimensional drops pinned on sharp edges is controlled by the Bond number, with critical pinning angles given by extrema of the pressure and area curves.","keywords":["pinned drops","Bond number","stability","turning-point methods","superhydrophobicity","Cassie-Baxter states","pressure eigenvalue","two-dimensional drops"],"falsifier":"For one Bond number, say $Bo=1$, compute the full spectrum of the second variation of the energy along the equilibrium branch and locate the first angle at which an eigenvalue crosses zero; if that angle differs from the extremum of $P_0(a)$ or $a(P_0)$, the turning-point criterion misses an instability. Alternatively, measure depinning angles of pressure-controlled drops pinned on a sharp-edged aperture of known width and compare them with the predicted critical angles.","tokens_in":8027,"feed_emoji":"💧","tokens_out":7284,"duration_ms":68280,"temperature":0.7,"pith_summary":"This paper establishes that the stability of two-dimensional drops pinned on sharp edges is controlled by the Bond number, the ratio of drop size to capillary length. Equilibrium shapes exist for a broad range of pinning angles, but stability restricts the realizable ones: closed drops with fixed volume become more stable as they shrink, while open drops connected to a reservoir become less stable as the Bond number decreases. The critical angles at which stability is lost are the extrema of the pressure eigenvalue $P_0(a)$ for open drops and of the area $a(P_0)$ for closed drops. These results give a quantitative stability map for the non-wetting Cassie-Baxter states that underpin superhydrophobic surfaces and for evaporating drops on micropatterned substrates.","feed_headline":"Smaller pinned drops are more stable, unless open to a reservoir","feed_subtitle":"Critical pinning angles follow extrema of pressure and area curves, giving a stability map for superhydrophobic surfaces.","key_machinery":"The central object is the pressure eigenvalue $P_0$, the pressure at zero height scaled by $\\sigma/W_d$, together with the scaled cross-sectional area $a$; both are functions of the pinning angle $\\theta_0$ and the Bond number $Bo$. Stability is read off from the curves $P_0(a)$ and $a(P_0)$: by Maddocks' turning-point theorem, unconstrained extremals change stability at folds of the solution branches, that is, at local extrema of $P_0(a)$, and constant-volume extremals change stability at local extrema of $a(P_0)$. The theorem replaces a direct check of the second variation of the energy functional with a geometric condition on computed equilibrium branches.","core_discovery":"Working in a planar arc-length formulation with the pressure eigenvalue $P_0$ as an unknown, the paper computes equilibrium branches for drops pinned on two sharp edges and locates their stability changes using turning-point arguments. It finds that unconstrained (constant-pressure) drops lose stability at extrema of $P_0$ as a function of area $a$, whereas constant-volume drops lose stability at extrema of $a$ as a function of $P_0$. The resulting stability diagram in the Bond-number--pinning-angle plane has two boundaries: below the constant-pressure limit the drop is stable to both perturbation classes; between the constant-pressure and constant-volume limits it is stable only to volume-preserving perturbations; above both it is unstable. At zero Bond number the two limits meet at $\\theta_0 = 90^\\circ$, matching earlier theory and experiment, and at large Bond number they approach the common limit $\\theta_0 = 180^\\circ$, a semi-infinite liquid layer. The same fold criterion applies to inverted (hanging) drops with negative pinning angles.","pith_inferences":["A direct extension of the fold criterion is that axisymmetric three-dimensional pinned drops should have stability boundaries set by extrema of $P_0(V)$ and $V(P_0)$, so the qualitative size trends found here would likely carry over to spherical-cap-like drops.","A testable experimental check would be to measure, at fixed Bond number, the pinning angle at which a pressure-controlled pinned drop depins and compare it with the extremum of $P_0(a)$; a mismatch would signal that a mode excluded by the planar analysis is active.","Because longitudinal Rayleigh-Plateau modes are excluded, the stability diagram should be read as the boundary for planar pinned perturbations; for drops that are long in the third dimension, three-dimensional bulging modes would set a lower practical stability limit.","The multiscale pressure-matching argument suggests a design rule: making the microfeature spacing much smaller than the drop radius keeps the underside interface flat and stable, so the robustness of superhydrophobic surfaces should scale with the ratio of the two Bond numbers."],"forward_implications":["On a superhydrophobic substrate, the Gibbs pinning condition is necessary but not sufficient: a drop can be in mechanical equilibrium yet unstable, so the attainable non-wetting states are those inside the stability region of the $Bo$--$\\theta_0$ map.","A shrinking closed drop becomes more stable, so an evaporating drop on a micropatterned surface should remain non-wetting until its size approaches the microstructure scale, at which point the interface underneath develops curvature and the drop transitions to a wetted state.","An open drop fed by a reservoir behaves oppositely: lowering the Bond number narrows the range of stable pinning angles, so pressure-controlled drops destabilize more easily at small scales.","When the drop and microstructure scales are widely separated, the micro-interface under the drop is nearly flat because its pressure eigenvalue must match that of the macro-interface; this links the two scales and explains why hierarchical micropatterning widens the non-wetting regime.","For $Bo \\to 0$ both stability limits converge to $\\theta_0 = 90^\\circ$, reproducing prior theoretical and experimental results, and for large $Bo$ they converge to $\\theta_0 = 180^\\circ$, where constant-pressure and constant-volume perturbations coincide."],"supporting_citations":[{"why":"Introduces the turning-point method that the paper uses to locate stability changes.","marker":"[4]"},{"why":"Supplies the core theorem that stability changes occur at folds of the solution branches in the $P_0$--$a$ plane.","marker":"[5]"},{"why":"Provides the modern discussion of unconstrained versus constrained stability that frames the two perturbation classes.","marker":"[6]"},{"why":"Gives the prior theoretical result that the zero-Bond-number critical angle is $90^\\circ$, used as a benchmark.","marker":"[7]"},{"why":"Provides the experimental observation of the same zero-Bond-number limit.","marker":"[8]"},{"why":"Motivates the constant-pressure stability scenario through experiments with liquid supplied through tubing.","marker":"[9]"},{"why":"Reports the nearly flat micro-interface under drops that the multiscale model reproduces.","marker":"[14]"},{"why":"Documents the evaporating-drop transition to wetting that the paper explains by size-dependent stability.","marker":"[18]"},{"why":"Establishes the robustness benefit of hierarchical micropatterning that the pressure-matching mechanism accounts for.","marker":"[23]"},{"why":"Analyzes the longitudinal Rayleigh-Plateau instability that the planar analysis deliberately excludes.","marker":"[19]"}],"fun_headline_variants":["Pinned drop stability hinges on pressure and area folds","Smaller sealed drops stable, open drops not at low Bond number","Stability map for pinned drops from extrema of pressure and area","Fold points set stability limits for two-dimensional pinned drops","Open and closed drops follow opposite stability rules"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that a pinned drop loses stability precisely when the pressure-versus-area relation turns over, and no other kind of disturbance, such as a bulge along the third dimension, goes unstable before that point.","fun_headline_variants_meta":{"raw":{"variants":["Pinned drop stability hinges on pressure and area folds","Smaller sealed drops stable, open drops not at low Bond number","Stability map for pinned drops from extrema of pressure and area","Fold points set stability limits for two-dimensional pinned drops","Open and closed drops follow opposite stability rules"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000593,"raw_usage":{"total_tokens":2777,"prompt_tokens":945,"completion_tokens":1832,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":1762}},"tokens_in":561,"tokens_out":1832,"duration_ms":12793,"temperature":1.0,"reasoning_tokens":1762,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:51:49.137023+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For one Bond number, say $Bo=1$, compute the full spectrum of the second variation of the energy along the equilibrium branch and locate the first angle at which an eigenvalue crosses zero; if that angle differs from the extremum of $P_0(a)$ or $a(P_0)$, the turning-point criterion misses an instability. Alternatively, measure depinning angles of pressure-controlled drops pinned on a sharp-edged aperture of known width and compare them with the predicted critical angles.","supporting_citations":[{"cited_title":"Bhushan and Y","cited_arxiv_id":null,"evidence_quote":"Introduces the turning-point method that the paper uses to locate stability changes."},{"cited_title":"Lenz and R","cited_arxiv_id":null,"evidence_quote":"Supplies the core theorem that stability changes occur at folds of the solution branches in the $P_0$--$a$ plane."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the prior theoretical result that the zero-Bond-number critical angle is $90^\\circ$, used as a benchmark."},{"cited_title":"Bostwick and P","cited_arxiv_id":null,"evidence_quote":"Provides the experimental observation of the same zero-Bond-number limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Motivates the constant-pressure stability scenario through experiments with liquid supplied through tubing."},{"cited_title":"Hensel, A","cited_arxiv_id":null,"evidence_quote":"Reports the nearly flat micro-interface under drops that the multiscale model reproduces."},{"cited_title":"Josserand and S","cited_arxiv_id":null,"evidence_quote":"Documents the evaporating-drop transition to wetting that the paper explains by size-dependent stability."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the robustness benefit of hierarchical micropatterning that the pressure-matching mechanism accounts for."},{"cited_title":"Zhang, Q","cited_arxiv_id":null,"evidence_quote":"Analyzes the longitudinal Rayleigh-Plateau instability that the planar analysis deliberately excludes."}],"review_version":1}