{"id":"a51aba75-db15-483d-9401-8c446ad99f35","arxiv_id":"1908.03705","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The evaporation lifetime of a droplet in a two-dimensional array scales with the ratio of droplet top-view area to its Voronoi cell area, independent of array symmetry and substrate.","lead":"This paper finds that a droplet's evaporation lifetime in any two-dimensional array depends mainly on the empty area around it, measured by a Voronoi cell, and not on the exact arrangement of neighboring droplets. It also shows that closely packed asymmetric arrays create directional internal flows that can be used to pattern deposited particles.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The universal lifetime scaling (Eqs. 7-8) depends on two empirical constants, A and β; if these are fitted to the same data, the Fig. 3 collapse does not test the model.","rationale":"The reader's formal 'weakest_assumption' points to the uniform-concentration/vertical-diffusion idealization in Eq. (2). That is a reasonable modeling concern, but the paper's own stated universality is equally vulnerable to the calibration status of A and β. The reader's rationale does mention that A and β are not derived, so there is partial overlap, but the stress-test should focus on the single most load-bearing issue: the predictive claim 'uniquely described by the Voronoi tessellation area' cannot be evaluated unless the constants in Eqs. (7)-(8) are independent of the data they explain. The proposed cross-validation test directly settles this: if the constants are universal, held-out predictions should match; if they were fitted, they will not. I also note the paper's genuine strengths: a clean geometric parameter (Voronoi area), a simple two-step diffusion model, and experimental coverage across substrates, asymmetries, and particle loading. These make the paper worth conditional publication, but the current draft does not yet establish that the model predicts rather than rationalizes the collapse.","tokens_in":8461,"tokens_out":5884,"duration_ms":66185,"concrete_test":"Perform a leave-one-substrate-out cross-validation. Fit A and β (or the combined prefactor) to the experimental t_c/t_uc data from four of the five substrates; use those constants to predict the held-out substrate's lifetimes from its A_df/A_v values using Eqs. (7)-(8). Report residuals against the experimental scatter. If the held-out predictions fall outside the within-substrate repeatability, A and β are not universal and the collapse is a fit artifact; if they agree, the predictive claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations (7)-(8) predict t_c/t_uc = 1 + C(θ) A_df/A_v with C(θ) = 2Aβf(θ) (hydrophilic) or 2Aβf(θ)/cosec²(θ) (hydrophobic). The two prefactors A and β are not derived. A is introduced as a constant '~O(1)' for the average unconfined decay rate, and β enters as L_a = βR_ci with 'β~4' stated from Supplementary Figure S5. The paper does not demonstrate that A and β are fixed a priori from independent measurements; if they are calibrated to the same experimental collapse, the apparent universality in Fig. 3 is a two-parameter fit. The central claim—that lifetime is uniquely determined by the Voronoi cell area—is only a prediction if A and β are universal, substrate-independent constants. This is especially load-bearing because the theoretical bounds in Fig. 3(a) are themselves generated by these constants; a post-hoc choice of A and β could reproduce the data for any monotone trend. The reader's concern about lateral vapor exchange (Eq. 2) is related but secondary: even under the vertical-diffusion model, the predictive content is controlled by the calibration status of A and β.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a generalized theoretical description of the evaporation lifetime of droplets arranged in arbitrary two-dimensional arrays. The central idea is to quantify vapor-mediated confinement by the Voronoi tessellation area A_v of each droplet. Evaporation is modeled as two serial diffusion steps: first from the droplet surface to a vapor-rich confinement region, and then vertically through the cell area A_v to the ambient. The resulting scaling, Eqs. (7)-(8), predicts t_c/t_uc as a linear function of A_df/A_v with a slope that depends only on the initial contact angle through known functions and the constants A and beta. The authors validate the prediction with experiments on glass, PDMS, GDL, and several other substrates, for symmetric and asymmetric arrays, and for nanoparticle-laden droplets. They further report that asymmetric confinement creates directional internal flows that can be exploited for patterned particle deposition.","tokens_in":8717,"tokens_out":4636,"duration_ms":48123,"significance":"If the proposed scaling holds, it provides a remarkably simple and potentially universal predictive tool for droplet lifetimes in arbitrary 2D arrays, covering different substrates, array asymmetries, and particle loadings with a single geometric parameter, the Voronoi cell area. The companion observation that local evaporation flux asymmetries can be tuned independently of global lifetime is also of practical interest for particle-patterning applications. The paper combines a clear physical picture, a compact analytical result, and a broad experimental data set, including several substrate chemistries. The strength of the claim, however, rests on the status of the two constants A and beta, since they enter the prefactor of the central scaling and currently are not derived from first principles or independently calibrated in the main text.","major_comments":[{"comment":"The prefactor C(theta) in Eqs. (7)-(8) contains two constants, A and beta, whose values are not established a priori. A is introduced as '~O(1)' immediately before Eq. (1), and beta~4 is stated to follow from Supplementary Figure S5. If A and beta are inferred from the same t_c/t_uc data plotted in Fig. 3, the reported collapse is a two-parameter fit rather than an independent test of the Voronoi-area universality. The authors should either derive A and beta from independent measurements (e.g., unconfined single-droplet decay curves and direct vapor-concentration profiling), or provide a sensitivity analysis showing that a single a priori (A,beta) pair simultaneously collapses all substrates, array geometries, and particle loadings.","section":"Theory, Eqs. (7)-(8)"},{"comment":"The vertical-diffusion model in Eq. (2) neglects lateral vapor exchange between neighboring Voronoi cells and in-plane concentration gradients within a cell. The central claim that the lifetime depends only on A_df/A_v, and not on the asymmetry parameter L_max/L_min, would break if lateral exchange contributes appreciably. The manuscript should quantify this contribution, for example by comparing the model with full three-dimensional diffusion simulations for the most asymmetric and most sparse cases tested, or by an explicit estimate of the ratio of lateral to vertical transport (e.g., a geometric factor involving A_v^{1/2}/L_a) below which lateral exchange is negligible.","section":"Eq. (2) and Fig. 3"},{"comment":"The step from the ratio of average volumetric decay rates to the ratio of lifetimes implicitly assumes that the initial droplet volume is the same for confined and unconfined droplets and that both evaporate to zero volume. More importantly, the model uses only the initial contact radius and angle, R_ci and theta_ci, in the rate expressions (1) and (7)-(8). For droplets evaporating in constant-angle, constant-radius, or stick-slip modes, R and theta vary with time, and f(theta) is a function of the instantaneous angle. The authors should state which evaporation mode their experiments follow and justify that A and beta, together with the initial-angle substitution, account for this time dependence.","section":"Eqs. (5)-(6)"}],"minor_comments":[{"comment":"The caption reads 'Area of the Voronoi cell, A_v quantifies the of vapor accumulation'; it should read 'quantifies the extent of vapor accumulation'.","section":"Figure 1 caption"},{"comment":"The notation for the accumulation length is inconsistent: the text introduces \\bar{L_a}, while Eq. (3) and subsequent equations use L_a without the overbar. Please define the symbol once and use it consistently.","section":"Throughout"},{"comment":"The statement that 'raw data shows variation of two times' is ambiguous: it could mean two repeated experiments, a factor-of-two spread, or two orders of magnitude. Please clarify and report the number of repeats and the corresponding error bars in Fig. 3.","section":"Fig. 3 discussion"},{"comment":"The reference list is inconsistent: 'Toledano, P., Mettout, B., Aroyo, M., & Mato, J. P.' appears in the text as both 'Toledano' and 'Tolédano'; please standardize the spelling.","section":"References"},{"comment":"The intermediate denominator A_v cosec^2(theta_ci) is algebraically correct but easy to misread; inserting A_df explicitly in the final line would help the reader follow the simplification.","section":"Eq. (8)"},{"comment":"The text says a '5 X 5 two-dimensional droplet array' with the center droplet 'surrounded by six adjacent droplets'; this implies a hexagonal packing rather than a square grid. Please clarify the actual arrangement used.","section":"Experimental methods"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely problem and the experimental data set is valuable. The main risk is that the two constants A and beta are not independently fixed, which would reduce the claimed universality to a fitting exercise. If the authors can supply an independent calibration or a convincing sensitivity analysis, the paper could become publishable. The journal should also consider encouraging the authors to deposit raw experimental data and analysis scripts to support the collapse claimed in Fig. 3."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here’s what I’d want you to know before reading: the paper has a genuinely useful idea, and it’s backed by a fair amount of experiment, but the headline universal prediction is not as clean as it looks because two constants in the final formula are not fixed before the data.\n\nWhat is new: nobody has used Voronoi tessellation to collapse evaporation lifetimes for arbitrary 2D droplet arrays. The authors show that for the center droplet in a 5x5 array, the lifetime ratio depends on the ratio of droplet top-view area to Voronoi cell area, and is insensitive to positional asymmetry. They demonstrate this over multiple substrates, contact angles from 35° to 135°, and with nanoparticle-laden droplets. The internal-flow directionality and deposit patterning for asymmetric arrays is a nice secondary result, even though it is not the main claim.\n\nWhat they do well: the derivation from Fick’s law is short and readable, and it reduces cleanly to earlier linear-array and single-droplet results. The comparison with Carrier et al. is fair, and the physics story—vapor fills the Voronoi cell, then escapes vertically—is transparent enough to test.\n\nSoft spots: the final lifetime formula, Eqs. (7)-(8), depends on A, stated as O(1), and β, stated as ~4 from Supplementary Fig. S5. It is not shown that either is fixed a priori from an independent measurement. If A and β are chosen after seeing the data, the agreement in Fig. 3 is a two-parameter fit and the collapse does not by itself validate the universal scaling. This is the load-bearing issue. The reader’s concern about lateral vapor exchange between cells is secondary; even if the vertical-diffusion model is right, the predictive content is controlled by the calibration status of A and β. There is also a minor notational issue in Eq. (3): A appears in the first-step average rate but not in the diffusive flux, so its survival into the final constant is more a modeling convention than a derived quantity.\n\nWho it is for: droplet evaporation and soft-matter/fluid-dynamics researchers, and people designing drying protocols for inkjet, spray coating, or colloid patterning. It is a solid within-subfield contribution, not a field-reshaping result.\n\nRecommendation: send it to peer review. A good referee should ask the authors to state exactly how A and β are obtained, and ideally to perform an out-of-sample test in which the constants are fixed on one subset and the prediction is checked on another. The Voronoi scaling deserves to be in the literature, but the universality claim should be advertised only for the regime in which the constants are validated.","headline":"A genuinely useful Voronoi-based lifetime scaling for arbitrary 2D droplet arrays, but the universal prediction is only as strong as the provenance of the two constants A and β.","tokens_in":9239,"tokens_out":4204,"would_cite":true,"duration_ms":49098,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A droplet's evaporation lifetime in a two-dimensional array depends only on the area of its Voronoi cell, not on the positions of its neighbors.","keywords":["cooperative evaporation","sessile droplets","two-dimensional droplet arrays","Voronoi tessellation","evaporation lifetime","vapor confinement","particle deposition","internal flow symmetry breaking"],"falsifier":"Print two droplet arrays whose center droplets have identical $A_{df}/A_v$ but very different cell shapes—one compact hexagonal cell and one long, thin cell—and measure the center-droplet lifetimes; if they differ beyond experimental scatter, the Voronoi-area-only description fails. A complementary check is a 3D diffusion simulation of both geometries that includes lateral vapor flux between cells.","tokens_in":8248,"feed_emoji":"💧","tokens_out":7951,"duration_ms":73747,"temperature":0.7,"pith_summary":"This paper claims that in an arbitrarily arranged two-dimensional array of sessile droplets, any droplet's evaporation lifetime is controlled by a single geometric quantity: the area of its Voronoi cell, the patch of the plane closer to that droplet than to any other. The authors derive a linear scaling, $t_c/t_{uc} = 1 + C(\\theta)\\,A_{df}/A_v$, in which the confined-to-unconfined lifetime ratio grows with the ratio of the droplet's top-view area $A_{df}$ to the Voronoi cell area $A_v$, and the coefficient $C(\\theta)$ depends only on the initial contact angle. They report experimental collapse onto this line for water and nanoparticle-laden droplets on substrates from hydrophilic glass to superhydrophobic GDL, for both symmetric and asymmetric arrays. If correct, this replaces the many-body vapor-confinement problem with a simple geometric partition, letting engineers predict drying times of printed or sprayed droplet patterns from an image. It also reveals that array asymmetry, while leaving the global lifetime unchanged, redirects internal flows and can steer where dissolved particles deposit.","feed_headline":"Voronoi area alone fixes droplet lifetime in 2D arrays","feed_subtitle":"One geometric ratio—droplet area over Voronoi cell—predicts evaporation time across substrates and asymmetry.","key_machinery":"The central object is the Voronoi tessellation of the droplet array, which partitions the plane into cells $A_v$; each cell is the spatial footprint of vapor confinement belonging to one droplet. The load-bearing quantity is the confinement ratio $A_{df}/A_v$, which enters linearly in the lifetime scaling and absorbs all geometric information about the array. The mechanism is a two-resistance diffusion model: Fickian escape from the droplet surface to the uniform confinement concentration $c'_\\infty$ is coupled to one-dimensional vertical diffusion through the cell area to the ambient, and the product of the two resistances produces the coefficient $2A\\beta f(\\theta_{ci})$ (with an additional $\\mathrm{cosec}^2\\theta_{ci}$ factor for hydrophobic substrates). The accumulation length $\\bar{L}_a$, the distance over which the vapor relaxes to the ambient, is scaled as $\\beta R_{ci}$ with $\\beta \\sim 4$, fixing the theory with no free fitting parameters. Transient vapor-front merging is neglected because its timescale, $10^{-2}$–$10^{-1}$ s, is far shorter than the droplet lifetime, $\\sim 10^3$ s.","core_discovery":"The paper's central claim is that cooperative evaporation in a 2D droplet array is universal in a specific sense: the lifetime of a droplet relative to an isolated one is insensitive to the relative positions of surrounding droplets and to the asymmetry of the array, and is uniquely determined by the Voronoi tessellation area $A_v$. The supporting model treats evaporation in two steps: vapor first diffuses from the droplet surface to a quasi-uniform 'confinement' region with concentration $c'_\\infty$, then escapes vertically through the cell cross-section $A_v$ to the ambient. Equating these fluxes yields the predictive formulas of Eqs. (7)-(8), $t_c/t_{uc} = 1 + [2A\\beta f(\\theta_{ci})]\\,A_{df}/A_v$ for hydrophilic substrates and the corresponding cosecant-weighted form for hydrophobic ones, with $\\beta \\approx 4$ and $f(\\theta)$ the standard single-droplet function. The paper reports that measurements across substrates, particle loadings, and symmetric versus asymmetric layouts collapse onto this trend, and that the earlier flat-superdrop model overshoots to a factor of two in sparse arrays while the new model agrees with experiments. A companion result is that high confinement combined with geometric asymmetry breaks the radial symmetry of local evaporation flux, generating a unidirectional internal flow and preferential deposition of particles at the less-confined side of the droplet.","pith_inferences":["The same Voronoi reduction suggests a mean-field structure, so the scaling might extend to polydisperse or three-dimensionally arranged droplets through weighted (power) Voronoi diagrams; the paper does not test this.","Because global lifetime and internal flow direction are decoupled, one could deliberately distort an array to pattern deposits while holding all droplet lifetimes fixed, which is a design strategy the paper hints at but does not fully develop.","A boundary of the model is likely at very low confinement (large, elongated cells), where lateral vapor exchange between cells should become comparable to vertical escape; locating where the linear law breaks is a natural next test.","A full 3D diffusion simulation with the same droplet positions could independently verify whether the lifetime collapse on $A_{df}/A_v$ is exact or only approximate; if substantial deviations appear, a correction term based on cell shape would be needed."],"forward_implications":["For any 2D droplet array, a top-view image plus the initial contact angle is enough to predict each droplet's evaporation lifetime: compute the Voronoi cell area and read $t_c/t_{uc}$ from the linear law.","The universal scaling extends to disordered and asymmetric arrays and to nanoparticle-laden droplets (0.5–5 wt% tested), so pattern-level evaporation control does not require resolving pairwise interactions.","Sparse arrays are handled correctly, where the earlier superdrop approximation breaks down and overshoots the lifetime ratio to 2.","Asymmetry can be used as a design input: at high confinement, an asymmetric array creates a unidirectional internal flow and moves particles preferentially to the less-confined side, without changing the global drying time.","Substrate wettability enters only through known functions of $\\theta_{ci}$, so the same formula is expected to hold for any droplet-substrate combination within the model's validity."],"supporting_citations":[{"why":"Supplies the Voronoi tessellation machinery used to partition the array into confinement cells.","marker":"Aurenhammer (1991)"},{"why":"Provides the flat-superdrop collective-evaporation model that the new scaling corrects and outperforms for sparse arrays.","marker":"Carrier et al. (2016)"},{"why":"Prior universal evaporation model for ordered droplet arrays that the present formulation generalizes to arbitrary arrangements.","marker":"Hatte et al. (2019)"},{"why":"Introduces the confined vapor-concentration picture and accumulation length that underlie the two-step diffusion equations.","marker":"Bansal et al. (2017b)"},{"why":"Gives the single-droplet evaporation-rate function $f(\\theta)$ used as the unconfined baseline.","marker":"Picknett & Bexon (1977)"},{"why":"Demonstrates experimentally that interacting droplets lengthen the central droplet's lifetime, motivating the confinement scaling.","marker":"Laghezza et al. (2016)"}],"fun_headline_variants":["Voronoi area alone sets droplet lifetime in 2D arrays","One geometric ratio predicts evaporation in droplet arrays","Droplet lifetimes in arrays collapse to a Voronoi rule","Cooperative evaporation universal: Voronoi area is key"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes that after a short transient, the vapor in each Voronoi cell is well mixed and escapes only straight upward through the cell's area, so neighboring cells do not trade vapor sideways.","fun_headline_variants_meta":{"raw":{"variants":["Voronoi area alone sets droplet lifetime in 2D arrays","One geometric ratio predicts evaporation in droplet arrays","Droplet lifetimes in arrays collapse to a Voronoi rule","Cooperative evaporation universal: Voronoi area is key"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00035,"raw_usage":{"total_tokens":1919,"prompt_tokens":959,"completion_tokens":960,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":575,"completion_tokens_details":{"reasoning_tokens":892}},"tokens_in":575,"tokens_out":960,"duration_ms":9766,"temperature":1.0,"reasoning_tokens":892,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:05:18.072240+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Print two droplet arrays whose center droplets have identical $A_{df}/A_v$ but very different cell shapes—one compact hexagonal cell and one long, thin cell—and measure the center-droplet lifetimes; if they differ beyond experimental scatter, the Voronoi-area-only description fails. A complementary check is a 3D diffusion simulation of both geometries that includes lateral vapor flux between cells.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Voronoi tessellation machinery used to partition the array into confinement cells."}],"review_version":1}