{"id":"3d644237-77fd-44c3-934c-544dc46cad68","arxiv_id":"2501.16540","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A viscoelastic polymer coating on a bursting bubble makes the jet thinner and faster when the coating is thicker, but slower and broader when the polymer is more concentrated, shifting between droplet and no-droplet regimes.","lead":"This experiment shows that a thin viscoelastic polymer layer on a bursting bubble changes the jet that sprays upward: a thicker layer makes the jet faster and thinner, while a more concentrated polymer makes it slower and thicker. The results help explain how organic slime on ocean bubbles controls the number and size of droplets that become aerosols.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The De~10^-2 no-jet-drop boundary and the viscoelastic-stress interpretation depend on a literature lambda_r correlation not validated for the PEO (MW 6x10^5) used here; direct extensional rheometry on these exact solutions would settle the regime map.","rationale":"I read the manuscript independently with the central claims in mind. The directly measured trends - v_j increasing and r_j decreasing with psi_0 at fixed c; v_j decreasing and r_j increasing with c at fixed psi_0; N_d following the same pattern; drop ejection ceasing at c = 0.2 wt% - are well supported by the high-speed images, Figs. 5-9, and the collapse-time invariance in Fig. 4. These are the empirical backbone of the paper and I see no internal inconsistency in them. The load-bearing point is the conversion of concentration into Deborah number. Because De = lambda_r/t_c and lambda_r is taken from a literature correlation rather than measured, the stated boundary De~10^-2 in Fig. 10 and the sentence 'compound bubble bursting transitions into a no-jet-drop regime for De >= 10^-2' are only as secure as that correlation. Relaxation times of PEO solutions are known to depend strongly on molecular weight, concentration, and solvent quality; applying a fit obtained for a different PEO (or different MW) introduces an unquantified systematic error that propagates directly into the regime map and into the elastic-stress interpretation. The reader's conditional verdict is appropriate. I also note a smaller caveat: the manuscript says Oht is 'relatively constant' while Table 3 shows a 64% increase from 3.72e-3 to 6.12e-3; this does not overturn the viscoelastic interpretation but means residual Newtonian viscous effects are not fully excluded. Both issues are addressable without changing the measured trends, so the verdict should stay CONDITIONAL (UNCHANGED from the reader).","tokens_in":13708,"tokens_out":7341,"duration_ms":73872,"concrete_test":"Measure the apparent extensional relaxation time lambda_r of the exact PEO solutions used (MW 6x10^5 g/mol, c = 0.01, 0.02, 0.03, 0.06, 0.1, 0.2 wt%) by capillary breakup extensional rheometry (CaBER/ROJER) or an equivalent technique at the strain rates relevant to the jet; recompute De and redraw Fig. 10. If the boundary remains at De~10^-2, the regime map and the elastic-stress interpretation stand; if the De values shift by more than a factor of ~2, the boundary should be re-expressed in terms of c or another measured fluid property.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim includes a regime map in De-psi0 space with the jet-drop/no-jet-drop boundary at De~10^-2 (Fig. 10), and the conclusion that increasing polymer concentration suppresses jetting 'due to increasing viscoelastic stresses' (Section 3, Conclusions). These statements are built on the Deborah number De = lambda_r/t_c, where lambda_r is not measured in this work. The paper adopts the semi-empirical fit lambda_r = 2.707e-7 c^1.733 ms (c in ppm) from refs. [26,49] (Section 2, Materials), but the relaxation time of dilute PEO solutions is strongly molecular-weight dependent, and the manuscript does not establish that the correlation applies to the Sigma-Aldrich PEO of MW 6x10^5 g/mol used here. If the true lambda_r at c = 0.2 wt% is, say, an order of magnitude smaller, the 'no-jet-drop' cases would sit below De~10^-2 and the claimed boundary would not be a material transition. The monotonic trends of v_j, r_j, v_d, r_d, and N_d with c and psi_0 are direct measurements and would survive a revised De, so the central empirical trends are not at risk; what is at risk is the dimensionless boundary and the elastic-stress mechanism. A secondary caveat: the text states Oht is 'relatively constant' while Table 3 lists Oht from 3.72e-3 to 6.12e-3 (a 64% increase), so residual Newtonian viscous effects are not fully excluded, but this is subordinate to the unmeasured lambda_r issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports an experimental study of Worthington jet formation from bursting compound bubbles, in which a bubble of radius R0 ≈ 1.48 mm is coated with a layer of aqueous PEO solution and rises through hexadecane. The authors use synchronized high-speed imaging to measure the cavity collapse time, capillary-wave wavelength, cavity cone angle, jet velocity and radius, top-drop velocity and radius, and the number of jet drops as functions of the compound-layer volume fraction ψ0 (0–60%) and PEO concentration c (0.01–0.2 wt%). They report that at fixed c, increasing ψ0 produces faster, thinner jets and smaller, more numerous drops; at fixed ψ0, increasing c produces slower, broader jets and eventually suppresses drop ejection at c = 0.2 wt%. The concentration effect is interpreted through the Deborah number De = λr/tc and the elastocapillary number Ec, and a De–ψ0 regime map is presented with a jet-drop/no-jet-drop boundary near De ≈ 10^-2. The apparent extensional relaxation time λr is taken from the literature correlation of Rodríguez-Díaz et al. rather than measured for the specific PEO used here.","tokens_in":14056,"tokens_out":4286,"duration_ms":44897,"significance":"If the results hold, this is the first systematic experimental characterization of bubble-bursting jets for compound bubbles with a viscoelastic layer, a configuration directly relevant to marine aerosol generation by EPS-coated bubbles. The qualitative trends in jet velocity, jet radius, and drop number versus ψ0 and c are directly supported by the imaging data and appear internally consistent, and the paper reports the uncertainty in the cavity collapse time. The main risk is that the dimensionless interpretation and the De ≈ 10^-2 boundary rely on an apparent extensional relaxation time taken from a literature correlation rather than measured for the present PEO solution; this does not threaten the direct empirical trends but does affect the claimed mechanism and the regime map. The statistical reporting is also incomplete: the compound-bubble data in Figs. 7–9 lack error bars and replicate counts.","major_comments":[{"comment":"The De–ψ0 regime map and the conclusion that drop ejection is suppressed by viscoelastic stresses depend on De = λr/tc, where λr is not measured in this work but is taken from the semi-empirical fit λr = 2.707 × 10^-7 c^1.733 ms of Rodríguez-Díaz et al. (refs. 26 and 49). The manuscript does not establish that this correlation transfers to the Sigma-Aldrich PEO of MW 6 × 10^5 g/mol used here, and relaxation times of dilute PEO solutions are strongly molecular-weight dependent. If the true λr at c = 0.2 wt% were an order of magnitude smaller, the no-jet-drop cases would fall below De ≈ 10^-2 and the claimed boundary would not be a material transition. The direct measurements of vj, rj, vd, rd, and Nd are not at risk, but the dimensionless boundary and the elastic-stress mechanism are. Please add direct extensional rheometry on these exact solutions, or at minimum a quantitative sensitivity analysis (e.g., propagating an order-of-magnitude uncertainty in λr through Fig. 10 and stating whether the boundary remains identifiable).","section":"Section 2 (Materials) and Fig. 10"},{"comment":"The manuscript does not report replicate counts or error bars for the compound-bubble data points in Figs. 7, 8, and 9; error bars appear only for the bare-bubble baseline in Fig. 7. As a result, the reader cannot assess whether the reported trends (for example, the increase of vj with ψ0 or the decrease of Nd with c) are statistically significant, nor whether the differences among PEO concentrations exceed experimental scatter. Please state the number of repeated runs per condition, report standard deviations or confidence intervals for all data points, and where possible give a statistical test or at least a clear statement of the measurement uncertainty for each plotted quantity.","section":"Figs. 7–9 and Section 3"},{"comment":"The text states that Oht remains 'relatively constant' as the polymer concentration increases, but Table 3 lists Oht values from 3.72 × 10^-3 (0.01 wt%) to 6.12 × 10^-3 (0.2 wt%), a 64% increase. Because the total Ohnesorge number changes substantially across the concentration range, the attribution of the slower and broader jets at higher c solely to increasing viscoelastic stresses is not fully separated from a possible Newtonian viscous effect. This point is subordinate to the unmeasured λr issue, but it should be addressed, for example by comparing with a Newtonian fluid of matched Oht or by quantifying the expected effect of the Oht variation on jet velocity and radius.","section":"Section 3 (Jetting dynamics) and Table 3"}],"minor_comments":[{"comment":"There are duplicated verbatim paragraphs: the introductory paragraph beginning 'Here, we experimentally investigate the bursting dynamics...' appears twice, and the concluding paragraph beginning 'We believe that our findings advance the understanding...' also appears twice. Please remove the duplicates.","section":"Section 1 and Section 4"},{"comment":"The sentence 'While PEO molecules are excepted to adsorb onto the interface' contains a typo; it should read 'expected to adsorb'.","section":"Section 3 (Cavity collapse)"},{"comment":"Reference 26 and reference 49 contain LaTeX escape sequences such as 'Rodr \\'iguez-D \\'iaz', 'Ga\\~n\\'an-Calvo', and 'Cabezas' that are not rendered; these should be cleaned up to proper UTF-8 names.","section":"References 26 and 49"},{"comment":"The right panel of the regime map projects the current experiments onto the Bo–Oht space of Walls et al. (ref. 57); please clarify in the caption or text how Oht is computed for the compound bubbles and why the no-jet-drop cases nevertheless fall within the Newtonian jet-drop region with respect to the Bo–Oht criterion.","section":"Fig. 10"},{"comment":"The definition of ψ0 = 3V0/(4πR0^3) states that V0 is obtained by image analysis right before jet formation, but it is not clear whether this volume is the initial coating volume or the volume at the moment of measurement; please clarify when and how V0 is evaluated, and whether the reported ψ0 values refer to the initial or instantaneous compound-layer fraction.","section":"Section 2 (Experimental setup)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reports a substantial experimental dataset on a relevant and understudied configuration, and the central empirical trends are likely to be robust. The main issue is that the dimensionless interpretation and the regime map hinge on a literature relaxation-time correlation that is not validated for the specific PEO solution; this is fixable within the scope of the manuscript, either by direct rheometry or by softening the De-based claims. I do not see grounds for rejection, but the statistical reporting also needs improvement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on arXiv:2501.16540. The paper does something new: it extends the group's oil-coated bubble work to a viscoelastic PEO layer on a bubble in a Newtonian bulk. The main measured trends are clear and internally consistent: for fixed polymer concentration, larger compound layer fraction gives faster, thinner jets and more drops; for fixed fraction, higher concentration gives slower, thicker jets and fewer drops. The imaging looks careful, and the collapse-time insensitivity is a useful data point. The beads-on-a-string observations and the growing gap between jet and drop velocity with increasing De are plausible and support the viscoelastic-stress narrative.\n\nThe soft spots are exactly where the stress-test note lands. First, most data plots (Figs. 7-9) show error bars only for the bare-bubble baseline; there are no replicate counts for the compound cases, so we can't judge scatter. That's fixable in revision. Second, and more important, the De-psi0 regime map and the 'no-jet-drop for De >~1e-2' boundary depend entirely on an apparent extensional relaxation time taken from a literature power-law fit for PEO, not measured on their solutions (MW 6x10^5). The stress-test note is right: if the real lambda_r is an order of magnitude different at 0.2 wt%, the boundary moves and the mechanism claim weakens. The qualitative trends with c survive, but the dimensionless boundary and the elastic-stress interpretation need direct extensional rheometry on these exact solutions, or at least a sensitivity analysis and a clear caveat.\n\nA minor complaint: the text says Oht is 'relatively constant' while Table 3 shows it goes from 3.72e-3 to 6.12e-3, a 64% increase. That's not huge but it's not nothing; they should soften that claim. Also the manuscript has some duplicated paragraphs in the intro and conclusion, and a few garbled figure references. Copy editing will catch those.\n\nOverall, the central empirical claims look solid and are a real step beyond prior work. The dimensionless framework is the fragile part. This deserves a serious referee who knows extensional rheometry and bubble bursting; with the relaxation-time issue addressed, it would be a good contribution. I'd send it to review.","headline":"New experimental data on compound bubble bursting with solid main trends, but the regime map rests on an unmeasured relaxation time that needs direct rheometry or a strong caveat.","tokens_in":14597,"tokens_out":2334,"would_cite":false,"duration_ms":21895,"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":"Bursting bubbles coated with a polymer layer produce jets whose speed and drop ejection are governed by two competing controls: coating thickness and polymer strength.","keywords":["bubble bursting","Worthington jet","viscoelastic coating","jet drops","capillary wave damping","Deborah number","polyethylene oxide","aerosolization"],"falsifier":"Directly measure the extensional relaxation time of the 0.01–0.2 wt% PEO solutions (for example by capillary breakup rheometry), recompute the Deborah number for each case, and check whether all no-jet-drop cases fall above $De \\approx 10^{-2}$ and all jet-drop cases below it; testing a second polymer with the same relaxation time but different chemistry would show whether the boundary is universal.","tokens_in":13504,"feed_emoji":"🪫","tokens_out":11419,"duration_ms":97264,"temperature":0.7,"pith_summary":"This paper experimentally establishes that a thin viscoelastic polymer layer coating a bursting bubble controls the Worthington jet and the drops it throws. At a fixed polymer concentration, increasing the coating volume fraction $\\psi_0$ makes the jet faster and narrower and produces smaller, more numerous jet drops. At a fixed $\\psi_0$, increasing polymer concentration makes the jet slower and wider and reduces the number of drops, until drop ejection stops entirely near $De \\approx 10^{-2}$. The authors attribute the first trend to stronger capillary-wave damping and a smaller cavity cone angle, and the second to growing viscoelastic stresses. The result matters because ocean bubbles are often coated with extracellular polymeric substances, so it links bubble bursting to how biological material enters the atmosphere.","feed_headline":"Thicker polymer coats make bursting bubbles jet faster and thinner","feed_subtitle":"At fixed polymer strength, thicker coats shrink drops; stronger polymer solutions slow jets and suppress drops.","key_machinery":"The load-bearing objects are the compound bubble—an air core coated by a polyethylene oxide solution (a weakly viscoelastic, nearly constant-viscosity Boger fluid) sitting in a Newtonian bulk—and two dimensionless controls: the coating volume fraction $\\psi_0$ and the Deborah number $De = \\lambda_r/t_c$, where $\\lambda_r$ is the apparent extensional relaxation time and $t_c = \\sqrt{\\rho_b R_0^3/\\gamma_e}$ is the inertio-capillary time. During cavity collapse, capillary waves run along the air-coating and coating-bulk interfaces; the wave separation shortens the characteristic wavelength $L$, increasing damping, while the converging waves form a cavity cone of semi-angle $\\beta$ that sets the jet's speed and radius. The viscoelastic stresses, scaled by $De$ and the elastocapillary number $Ec$, resist the extensional stretching of the jet, slowing and widening it. The combination of these two mechanisms explains the measured jet velocity and radius trends and sets the drop/no-drop boundary at $De \\approx 10^{-2}$.","core_discovery":"The central discovery is that a bubble coated with a weakly viscoelastic polymer layer has two separable controls on its bursting jet. Increasing the coating volume fraction $\\psi_0$ at fixed polymer concentration enhances jetting: the dimensionless jet velocity $v_j/v_{ce}$ rises and the dimensionless jet radius $r_j/R_0$ falls, leveling off near $\\psi_0 \\gtrsim 30\\%$, because the thicker layer separates capillary waves onto the air-coating and coating-bulk interfaces, shortens the capillary wavelength, damps short-wavelength disturbances, and forms a smaller cavity cone angle $2\\beta$ during collapse. Increasing polymer concentration at fixed $\\psi_0$ suppresses jetting: viscoelastic stresses stretch and thicken the jet, making it slower and wider, reducing drop number, and at $c = 0.2$ wt% (where $De \\approx 2.7 \\times 10^{-2}$) suppressing jet drops entirely. The paper reports, for the first time, a regime map in the $(De,\\psi_0)$ plane with a jet-drop/no-jet-drop boundary near $De \\approx 10^{-2}$, with $\\psi_0$ having negligible influence on that transition.","pith_inferences":["A natural extension is that the thickness-driven enhancement is a geometric and capillary effect that should occur for any immiscible low-surface-tension coating, while the suppression branch depends specifically on the polymer relaxation time.","The $De \\approx 10^{-2}$ boundary is only as robust as the published power-law relaxation-time correlation used to compute $De$; direct measurement of $\\lambda_r$ for these solutions could move the boundary.","In marine settings, bubbles coated by extracellular polymeric substances could fall on either side of the boundary depending on local polymer concentration, implying that biological aerosol fluxes may respond non-monotonically to polymer loading.","A testable extension is to measure the full drop-size distribution including satellite drops: the bead-on-a-string dynamics imply the clearest signature of increasing $De$ should appear in satellite counts and filament lifetimes rather than only in the top-drop radius."],"forward_implications":["Thicker polymeric coatings on a bubble of fixed composition produce faster, thinner jets and smaller, more numerous jet drops, so coating thickness alone can amplify the aerosol emission of the coating material.","Stronger polymer solutions produce slower, thicker jets with fewer drops, and above $De \\approx 10^{-2}$ no jet drops are ejected at all.","The cavity collapse time is unchanged across all tested concentrations and coating fractions, so the coating alters jet formation without changing the collapse timing.","Because the drop/no-drop boundary is set by $De$ rather than $\\psi_0$, the polymer relaxation time, not the coating thickness, determines whether coated bubbles release drops into the air."],"supporting_citations":[{"why":"Provides the weakly-viscoelastic bubble-bursting baseline and the power-law correlation $\\lambda_r = 2.707 \\times 10^{-7} c^{1.733}$ ms used to compute the Deborah number for every concentration.","marker":"[26]"},{"why":"Establishes the compound-bubble method and the capillary-wave-separation effect: a coating shortens the capillary wavelength and increases damping, which the paper extends from oil to polymer coatings.","marker":"[45]"},{"why":"Supplies the theoretical relation between cavity cone angle and jet speed and radius used to interpret why smaller cone angles at higher $\\psi_0$ produce faster, thinner jets.","marker":"[55]"},{"why":"Shows experimentally and numerically that polymer viscosity and relaxation time control Worthington jet velocity and where drops are emitted, the baseline the compound-layer results are compared with.","marker":"[27]"},{"why":"Provides the Oldroyd-B simulations and the $De$/$Ec$ definitions used to frame viscoelastic stress effects and elastocapillary filament thinning.","marker":"[28]"},{"why":"Supplies the Oh-Bo regime map for bare bubble jet-drop production, used to show that the no-jet-drop cases fall inside the Newtonian drop region and therefore the coating changes the boundary.","marker":"[57]"}],"fun_headline_variants":["Polymer thickness speeds bubble jets, strength slows them","Bubble burst jets: thicker coats shrink drops, stronger solutions swell them","Coating thickness and viscoelasticity independently steer bubble jets","Thicker polymer layers jet faster, stronger ones jet slower","Polymer coat tuning: fast thin jets or slow thick ones"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the relaxation time of the polymer solution, which sets the Deborah number $De$ used to draw the drop/no-drop boundary, is accurately given by a published power-law fit rather than measured for these exact solutions; if that fit does not transfer, the $De$ values and the $De \\approx 10^{-2}$ boundary would shift.","fun_headline_variants_meta":{"raw":{"variants":["Polymer thickness speeds bubble jets, strength slows them","Bubble burst jets: thicker coats shrink drops, stronger solutions swell them","Coating thickness and viscoelasticity independently steer bubble jets","Thicker polymer layers jet faster, stronger ones jet slower","Polymer coat tuning: fast thin jets or slow thick ones"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000291,"raw_usage":{"total_tokens":1766,"prompt_tokens":1078,"completion_tokens":688,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":694,"completion_tokens_details":{"reasoning_tokens":603}},"tokens_in":694,"tokens_out":688,"duration_ms":8087,"temperature":1.0,"reasoning_tokens":603,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T12:29:30.008228+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly measure the extensional relaxation time of the 0.01–0.2 wt% PEO solutions (for example by capillary breakup rheometry), recompute the Deborah number for each case, and check whether all no-jet-drop cases fall above $De \\approx 10^{-2}$ and all jet-drop cases below it; testing a second polymer with the same relaxation time but different chemistry would show whether the boundary is universal.","supporting_citations":[{"cited_title":"Bubble bursting in a weakly viscoelastic liquid","cited_arxiv_id":null,"evidence_quote":"Provides the weakly-viscoelastic bubble-bursting baseline and the power-law correlation $\\lambda_r = 2.707 \\times 10^{-7} c^{1.733}$ ms used to compute the Deborah number for every concentration."},{"cited_title":"T.; Feng, J","cited_arxiv_id":null,"evidence_quote":"Establishes the compound-bubble method and the capillary-wave-separation effect: a coating shortens the capillary wavelength and increases damping, which the paper extends from oil to polymer coatings."},{"cited_title":"M.; Blanco-Rodr \\' guez, F","cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical relation between cavity cone angle and jet speed and radius used to interpret why smaller cone angles at higher $\\psi_0$ produce faster, thinner jets."},{"cited_title":"Effect of the polymer viscosity and relaxation time on the Worthington jet produced by bubble bursting in weakly viscoelastic liquids","cited_arxiv_id":null,"evidence_quote":"Shows experimentally and numerically that polymer viscosity and relaxation time control Worthington jet velocity and where drops are emitted, the baseline the compound-layer results are compared with."},{"cited_title":"Viscoelastic Worthington jets & droplets produced by bursting bubbles","cited_arxiv_id":"2408.05089","evidence_quote":"Provides the Oldroyd-B simulations and the $De$/$Ec$ definitions used to frame viscoelastic stress effects and elastocapillary filament thinning."},{"cited_title":"bead-on-a-string","cited_arxiv_id":null,"evidence_quote":"Supplies the Oh-Bo regime map for bare bubble jet-drop production, used to show that the no-jet-drop cases fall inside the Newtonian drop region and therefore the coating changes the boundary."}],"review_version":1}