{"id":"69353881-1330-426a-a667-7c3bc87e2b09","arxiv_id":"2507.00774","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A pair of water drops on a superhydrophobic substrate evaporates slower than a lone drop due to vapor shielding, and a diffusion-plus-cooling-plus-convection model approximates the heated case.","lead":"Drops placed side by side on a heated superhydrophobic surface evaporate more slowly than an isolated drop, because vapor from one drop shields the other, extending lifetime 1.6 times at room temperature and 1.2 times at 50°C. The paper checks whether standard diffusion models, upgraded with cooling and convection corrections, can predict these dynamics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'comprehensive model accurately captures' claim rests on a single half-volume match produced by two correlations explicitly outside their validity ranges; a full-curve error check is needed before the model claim can be accepted.","rationale":"The experimental observation of vapor shielding is solid: the lifetime ratio changes from 1.6 at room temperature to 1.2 at 50 C in the expected direction, and the repeated measurements give consistent morphology data. The stress test therefore targets the modeling claim. Eqs. (14)-(15) multiply a diffusion-limited two-drop rate by two empirical correlations, and the authors themselves state that both are outside the experimental regime: K(E,theta) is for an isolated drop under saturated conditions, and the Kelly-Zion correlation is for a single low-contact-angle drop. Since no parameter is fitted, the single 14% half-time agreement at 50 C is weak evidence for 'accurately captures', especially because the cooling factor reduces the rate while the convection factor increases it, so errors can cancel. The proposed full-curve error check uses data already present in Fig. 10 and would settle whether the agreement is a genuine prediction of the temporal dynamics or an artifact of evaluating only one time point. This does not change the reader's CONDITIONAL verdict: the experimental contribution remains valuable, but the modeling claim needs to be stated more cautiously until a full-curve or multi-temperature validation is provided.","tokens_in":18300,"tokens_out":8871,"duration_ms":105360,"concrete_test":"Recompute the D_f+Ec+Cv prediction (Eqs. 14-15) against the experimental V/V0(t) curves in Fig. 10 at t/t_f = 0.25, 0.5, and 0.75, and report the maximum and mean absolute error over the full measured range (0 < t/t_f < 0.8) for both the isolated and two-drop cases at Ts = 50 C. If the full-curve error exceeds about 20% or is strongly non-uniform in time, the claim that the model 'accurately captures' the dynamics is not supported and should be tempered.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The experimental core—paired drops on a superhydrophobic substrate live 1.6x (27 C) and 1.2x (50 C) longer than isolated drops—is well supported by repeated shadowgraphy measurements. The load-bearing weakness is in the theoretical half of the central claim, namely that Eqs. (14)-(15) 'accurately capture' the elevated-temperature dynamics. Section III A, in the paragraph after Eq. (15), states explicitly that the evaporative-cooling factor K(E,theta) from Shen et al. [69] applies to an isolated drop under saturated conditions, and that the convection enhancement 0.31 Gr^0.216 from Kelly-Zion et al. [66] corresponds to a single-drop system with low contact angles. Both factors are therefore outside the two-drop, superhydrophobic, RH=16% regime of the experiments. The model has no adjustable parameters, and its only quantitative support reported for the two-drop case is the 14% error in the time to reach V/V0=0.5 at 50 C (Fig. 10d). No full-curve error metric is given. Because K=0.54 suppresses evaporation while the Gr term enhances it, the half-time match can arise from compensation between two out-of-range corrections rather than from a faithful representation of the physics. The paper itself flags this limitation, yet the abstract and conclusion still claim the model 'accurately captures' the dynamics; that claim is not established by the presented evidence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports an experimental and theoretical study of sessile water droplets (V0 ≈ 4 µl) evaporating on a micro-nano textured superhydrophobic aluminum substrate at Ts = 27 °C and 50 °C, RH = 16%, comparing an isolated droplet with a two-drop configuration. The main experimental findings are that droplets in the two-drop configuration evaporate more slowly due to vapor shielding, with lifetimes 1.6 times (27 °C) and 1.2 times (50 °C) those of an isolated droplet, and that the evaporation mode differs with temperature and configuration. The theoretical part combines a diffusion-based model (Popov) with an evaporative-cooling correction K(E,θ) from Shen et al. and a natural-convection enhancement Er = 0.31Gr^0.216 from Kelly-Zion et al., claiming this combined model 'accurately captures' the elevated-temperature dynamics while diffusion alone suffices at room temperature.","tokens_in":18566,"tokens_out":4798,"duration_ms":51714,"significance":"If substantiated, the experimental results provide a useful and novel dataset on multi-droplet evaporation on superhydrophobic surfaces at elevated temperatures, a regime that has received little attention. The repeated shadowgraphy measurements with uncertainty bars, the clear reporting of lifetimes and evaporation modes, and the absence of fitted parameters in the model are strengths. The experimental core—paired droplets live 1.6× and 1.2× longer than isolated droplets—is well supported. However, the theoretical modeling claim as stated is not established: the 'accurate' combined model rests on two empirical correlations explicitly outside their validity range and is supported only by a single half-volume-time match at 50 °C, raising the risk of error compensation.","major_comments":[{"comment":"The claim that the combined D_f+Ec+Cv model 'accurately captures' the elevated-temperature dynamics is not supported by the evidence presented. The only quantitative support in the two-drop case is a 14% overprediction of the half-volume time at Ts = 50 °C (Fig. 10d), and no full-curve error metric (e.g., L2 relative error in V/V0(t)) is reported. The text itself states that K(E,θ) is derived for an isolated drop under saturated conditions and that the convection correlation 0.31Gr^0.216 corresponds to a single-drop system with low contact angles; both are outside the present superhydrophobic, RH=16%, paired-drop regime. Since K=0.54 suppresses evaporation while the Gr term enhances it, the observed agreement could arise from compensation of two out-of-range corrections. Please either soften the claim to 'improves agreement' or strengthen it with a sensitivity analysis for K and Er and a full-curve error quantification.","section":"§III A, Eq. (15) and Fig. 10(d); Abstract and Conclusion"},{"comment":"The geometric ratios used for the two-drop model are internally inconsistent with the reported experimental dimensions. The text reports Lc/Rc ≈ 5.46 and Le/Rc = 0.52, but the measured initial droplet diameter is d0 ≈ 2.1 mm (Rc ≈ 1.05 mm) with Lc ≈ 2.46 mm and Le ≈ 0.22 mm, which gives Lc/Rc ≈ 2.3 and Le/Rc ≈ 0.2. Since the dimensionless concentration field φ in Eq. (6) depends on Rc, h, and Lc, the reported 18% underprediction at room temperature and 14% overprediction at 50 °C could be influenced by using a different value of Rc in the calculations than in the experiments. Please state the exact values of Rc, h, Lc, and Le used in Eqs. (9) and (15) and verify their consistency with the measured data.","section":"§III A, after Eq. (3)"},{"comment":"For the isolated drop at Ts = 50 °C, the paper states that the combined model 'agrees well' with experiments but reports no quantitative error. The preceding D_f+Ec model overpredicts the half-volume time by 16%, and the addition of the convection enhancement changes this to an unreported value. Please report the half-time error and, ideally, a full-curve error metric for the isolated-drop case at 50 °C; without this, the claim that Eqs. (14)–(15) accurately capture the isolated-drop dynamics is not quantitatively established.","section":"§III A and Fig. 10(c)"}],"minor_comments":[{"comment":"There are several typographical errors: 'theoreticaly' (start of §III A), 'dimater' (Fig. 2 caption), 'at at Ts' (Fig. S3 caption), 'op surface' (Fig. 2 caption), and 'boemite' (Fig. 3 caption, should be 'boehmite').","section":"Throughout"},{"comment":"The phrase 'both isolated and single-drop systems' appears to be a slip; it should presumably be 'isolated and two-drop systems'.","section":"§III A, after Eq. (15)"},{"comment":"The sentence 'Changing this θ value in the range 155° ± 5°' is inconsistent with the just-stated calculation value θ = 150°; the intended range is likely 150° ± 5°.","section":"§III A, after Eq. (3)"},{"comment":"The abstract and conclusion use the phrase 'accurately captures' for the combined model, while the conclusion's own wording later softens to 'improves agreement'; aligning these statements would more accurately reflect the evidence.","section":"Abstract and Conclusion"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid experimental paper whose headline numbers—1.6x lifetime at 27°C and 1.2x at 50°C—are well supported by repeated shadowgraphy runs. The modeling half overreaches: the \"comprehensively captures\" language in the abstract and conclusion is not backed by the reported metrics.\n\nWhat's new: first experimental study of two sessile drops on a micro-nano textured superhydrophobic substrate at room and elevated temperature, with careful side-by-side imaging. The lifetime ratios and the vapor-shielding interpretation are plausible. The mode analysis (CCA, mixed, stick-slip) is consistent with prior work, and the supplementary repeats show the scatter. The comparison uses established models with no fitted constants, so the experimental results stand independently of the theoretical section.\n\nSoft spots: the model leans on two correlations the authors themselves admit are outside their validity range—the cooling factor K(E,θ) from Shen et al. (isolated drop, saturated conditions) and the convection enhancement 0.31Gr^0.216 from Kelly-Zion et al. (single, low-contact-angle drops). For the two-drop 50°C case, the combined model overpredicts the half-evaporation time by 14%. That is a single time point, not a full-curve metric, and with K suppressing while the Gr term enhances, the agreement could reflect compensating errors. The paper flags these limitations in prose but then repeats \"accurately captures\" in the abstract and conclusion; that claim should be softened. Also, the model assumes equal-size symmetric evaporation with J = J_iso/(1+φ), while the paper elsewhere emphasizes asymmetric evaporation; that inconsistency is not addressed. Minor: Sec. III A says \"isolated and single-drop systems\" where it presumably means the two-drop configuration.\n\nNone of this undercuts the experimental core. The lifetime ratios and mode transitions are useful data for spray cooling and coating applications, and the paper honestly reports its own limitations in the body.\n\nFor peer review: yes, send it out. A referee should ask for a full-curve error metric, a more careful statement of what the model does and does not capture, and a discussion of the symmetric-evaporation assumption. With those revisions the paper would be solid. I'd cite it for the lifetime ratios if I worked on multi-drop evaporation.","headline":"Useful experimental dataset on two-drop evaporation on superhydrophobic surfaces at elevated temperature; the modeling claim is stronger than the evidence supports.","tokens_in":19110,"tokens_out":2212,"would_cite":true,"duration_ms":26125,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two closely spaced droplets on a superhydrophobic surface evaporate more slowly than an isolated droplet—1.6 times slower at 27 °C and 1.2 times slower at 50 °C—because vapor shielding raises the local humidity between them.","keywords":["Sessile droplet evaporation","Superhydrophobic substrate","Vapor shielding","Two-drop configuration","Evaporative cooling","Natural convection","Contact angle dynamics","Shadowgraphy imaging"],"falsifier":"Measure the interfacial temperature of the two drops during evaporation with infrared thermography at $T_s = 50^\\circ$C: the model's cooling factor predicts a specific suppression of the interface temperature relative to the substrate, and if a high-contact-angle superhydrophobic pair does not show that suppression, the 14% agreement at half-volume would have to be considered coincidental.","tokens_in":18020,"feed_emoji":"💧","tokens_out":12345,"duration_ms":116849,"temperature":0.7,"pith_summary":"This paper investigates how two water droplets sitting side by side on a water-repellent micro-nano textured superhydrophobic aluminum surface evaporate, compared with an isolated droplet, at 27 °C and 50 °C. It claims that the two-drop configuration evaporates more slowly because vapor from each droplet accumulates between them—a vapor-shielding effect—producing asymmetric evaporation and longer lifetimes: about 1.6 times longer at room temperature and 1.2 times longer at 50 °C. The paper further claims that a purely diffusion-based model captures evaporation at room temperature, while at 50 °C the model must also include evaporative cooling and buoyancy-driven natural convection to track the measured volume. The practical stakes: droplet arrays on superhydrophobic surfaces evaporate at different rates than isolated droplets, and the difference depends strongly on substrate temperature.","feed_headline":"Twin droplets evaporate up to 1.6x slower than a lone drop","feed_subtitle":"At 50°C the slowdown shrinks to 1.2x — and diffusion alone stops predicting the evaporation rate.","key_machinery":"The load-bearing object is the vapor-shielding factor $1/(1+\\varphi)$, where $\\varphi$ is the dimensionless vapor-concentration field of the isolated spherical cap (given analytically by Eqs. (6)--(8)); the pair model divides the isolated evaporation rate by this factor, assuming equal drops and symmetric evaporation. At elevated temperature, two empirical corrections carry the argument: the evaporative-cooling factor $K(E,\\theta)$ (with $E = 0.19$ and $\\theta = 150^\\circ$, giving $K = 0.54$) and the natural-convection enhancement $E_r = 0.31 Gr^{0.216}$ based on the Grashof number. These are combined with the spherical-cap diffusion solution $f(\\theta)$ into Eqs. (14)--(15), which are integrated to predict $V(t)$.","core_discovery":"On a superhydrophobic substrate, two droplets placed side by side with edge-to-edge gap $L_e \\le 0.37$ mm live significantly longer than a single drop. At $T_s = 27^\\circ$C the pair's average lifetime is 2050 s versus 1269 s for the isolated drop (1.6$\\times$); at $T_s = 50^\\circ$C it is 593 s versus 494 s (1.2$\\times$). The paper attributes this to vapor shielding: the inner sides of the pair see a higher local vapor concentration, so the evaporation flux is asymmetric and the total rate is reduced. It supports this with a sequence of theoretical models: Eq. (3) (diffusion only) suffices at room temperature; at 50 $^\\circ$C the isolated drop and the pair require the diffusion equation multiplied by the evaporative-cooling factor $K(E,\\theta) = 0.54$ and by $1 + 0.31 Gr^{0.216}$ for natural convection, giving Eq. (14) for a single drop and Eq. (15) for the pair. The combined model overestimates the half-volume time at 50 $^\\circ$C for the two-drop system by only 14%, whereas diffusion alone underpredicts it by 31%.","pith_inferences":["The paper does not vary the edge-to-edge gap $L_e$ systematically; a testable extension is to measure pair lifetime against $L_e/R_c$ and compare with the $1/(1+\\varphi)$ prediction, which should show shielding decaying smoothly as the drops separate.","The model assumes equal-sized drops with symmetric evaporation; if one drop is smaller, the shielding field is asymmetric and the smaller drop should be shielded more strongly, a prediction that could be checked by dispensing unequal volumes.","If the trend extrapolates, dense droplet arrays on superhydrophobic surfaces at near-room temperature will show much longer collective lifetimes, while at elevated temperature natural convection short-circuits the shielding—a consideration for cooling and anti-icing applications.","Because both borrowed correlations were derived for single drops, the 14% error at 50 °C for the pair is not strong evidence by itself; a direct test of the two corrections on a superhydrophobic pair would separate mechanism from curve-fitting."],"forward_implications":["At 27 °C a paired droplet takes about 1.6 times as long to evaporate as an isolated one; at 50 °C the ratio drops to about 1.2, so heating a superhydrophobic substrate weakens vapor shielding.","At room temperature the diffusion-only model (Eq. 3) is sufficient; at 50 °C, diffusion alone underpredicts the half-volume time by 31% for the pair, while the full model with evaporative cooling and convection overpredicts it by only 14%.","In the pair, the evaporation rate converges to the isolated-drop rate after $t/t_{f,iso} \\approx 0.7$ at 27 °C and $\\approx 0.4$ at 50 °C, meaning the late-stage pair behaves like two independent drops.","Both isolated and paired droplets shift from mostly constant-contact-angle evaporation at room temperature to mixed-mode evaporation at 50 °C, with stick-slip events in both cases."],"supporting_citations":[{"why":"It supplies the multi-droplet shielding relation $J = J_{\\rm iso}/(1+\\varphi)$ and the analytical concentration field used for the pair.","marker":"[17]"},{"why":"It validates the multi-droplet model for center-to-center spacings with $L_c/R_c > 3$, the regime of these experiments.","marker":"[16]"},{"why":"It provides the evaporative-cooling correction $K(E,\\theta)$ that scales the diffusion rate at 50 °C.","marker":"[69]"},{"why":"It provides the empirical natural-convection enhancement $E_r = 0.31 Gr^{0.216}$ used at elevated temperature.","marker":"[66]"},{"why":"It supplies the heated hydrophobic/superhydrophobic diffusion-evaporation framework that the paper extends.","marker":"[56]"},{"why":"It gives the spherical-cap diffusion solution $f(\\theta)$ used in the single- and two-drop rate equations.","marker":"[54]"},{"why":"It supports the constant-contact-angle evaporation mode assumed for superhydrophobic droplets.","marker":"[55]"},{"why":"It provides prior multi-drop evaporation experiments on heated hydrophilic and hydrophobic substrates that the paper compares against.","marker":"[20]"},{"why":"It documents room-temperature collective evaporation on superhydrophobic surfaces, the experimental baseline that the paper extends to elevated temperature.","marker":"[22]"}],"fun_headline_variants":["Droplet pairs survive 1.6x longer via vapor shielding","Vapor shielding: twin drops live longer on superhydrophobic surfaces","At room temp, paired droplets last 1.6x longer than singles","Heated surfaces: droplet pairs still evaporate slower","Twin droplets: vapor shielding extends lifetime by 1.6x"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The elevated-temperature model imports an evaporative-cooling correction and a buoyant-convection enhancement that were measured for isolated, low-contact-angle drops under saturated conditions, and applies them unchanged to each member of a closely spaced pair on a superhydrophobic surface, while also assuming the two drops are equal-sized and evaporate symmetrically.","fun_headline_variants_meta":{"raw":{"variants":["Droplet pairs survive 1.6x longer via vapor shielding","Vapor shielding: twin drops live longer on superhydrophobic surfaces","At room temp, paired droplets last 1.6x longer than singles","Heated surfaces: droplet pairs still evaporate slower","Twin droplets: vapor shielding extends lifetime by 1.6x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000336,"raw_usage":{"total_tokens":1920,"prompt_tokens":1064,"completion_tokens":856,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":680,"completion_tokens_details":{"reasoning_tokens":763}},"tokens_in":680,"tokens_out":856,"duration_ms":9919,"temperature":1.0,"reasoning_tokens":763,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:07:31.784481+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the interfacial temperature of the two drops during evaporation with infrared thermography at $T_s = 50^\\circ$C: the model's cooling factor predicts a specific suppression of the interface temperature relative to the substrate, and if a high-contact-angle superhydrophobic pair does not show that suppression, the 14% agreement at half-volume would have to be considered coincidental.","supporting_citations":[{"cited_title":"D.; Stone, H","cited_arxiv_id":null,"evidence_quote":"It supplies the multi-droplet shielding relation $J = J_{\\rm iso}/(1+\\varphi)$ and the analytical concentration field used for the pair."},{"cited_title":"J.; Ouali, F","cited_arxiv_id":null,"evidence_quote":"It validates the multi-droplet model for center-to-center spacings with $L_c/R_c > 3$, the regime of these experiments."},{"cited_title":"Numerical and theoretical analysis of fast evaporating sessile droplets with coupled fields","cited_arxiv_id":null,"evidence_quote":"It provides the evaporative-cooling correction $K(E,\\theta)$ that scales the diffusion rate at 50 °C."},{"cited_title":"L.; Pursell, C","cited_arxiv_id":null,"evidence_quote":"It provides the empirical natural-convection enhancement $E_r = 0.31 Gr^{0.216}$ used at elevated temperature."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the heated hydrophobic/superhydrophobic diffusion-evaporation framework that the paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the spherical-cap diffusion solution $f(\\theta)$ used in the single- and two-drop rate equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supports the constant-contact-angle evaporation mode assumed for superhydrophobic droplets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides prior multi-drop evaporation experiments on heated hydrophilic and hydrophobic substrates that the paper compares against."},{"cited_title":"Collective behavior of evaporating droplets on superhydrophobic surfaces","cited_arxiv_id":null,"evidence_quote":"It documents room-temperature collective evaporation on superhydrophobic surfaces, the experimental baseline that the paper extends to elevated temperature."}],"review_version":1}