{"id":"44dc16ee-d06f-4f37-96c4-150ec704b0ef","arxiv_id":"1908.03768","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A practical polydisperse Morse-Witten foam simulation matches direct calculations and yields scaling laws for excess energy and coordination near jamming.","lead":"The authors turn the Morse-Witten theory of bubble contacts into a working simulation of wet 2D foams with mixed bubble sizes. The method reaches high liquid fractions where older foam software struggles, and it reproduces known results while suggesting new scaling behavior near jamming.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (10)'s linear-plus-quadratic superposition of single-contact Morse-Witten solutions is validated only for two bubbles; no test exercises a bubble with three or more contacts, so f^2 errors in disordered multi-contact foams remain unquantified.","rationale":"The reader's weakest assumption identifies exactly the load-bearing concern: Eq. (10) generalizes a two-bubble quadratic correction to multi-contact situations without a direct check. I agree with that assessment. The paper does provide genuine independent support where it exists: the two-bubble test in Figure 3 and the match with Princen's exact monodisperse result in Figure 6(a) are meaningful, and the convergence above phi = 0.12 is a useful practical contribution. However, neither test exercises the specific superposition that makes the polydisperse disordered case work. The Princen case is symmetric enough that many cross-contact terms may cancel or be uniform, and the Plat comparison on a single ten-bubble system is an aggregate test in which the contact-network shift of about Delta phi = 0.01 could hide systematic errors. The proposed concrete test isolates the assumption by comparing Eq. (10) against a full nonlinear calculation for a three-contact bubble, including polydisperse radii and unequal forces. If that test passes at the force levels reached near phi = 0.12, the central claim is adequately supported. If it fails, the conditional verdict should be revised toward rejection or the claim narrowed to a more restricted liquid-fraction range. Since the reader already assigned CONDITIONAL with moderate confidence, this stress test does not change the verdict; it reinforces the specific condition under which the paper should be accepted.","tokens_in":9878,"tokens_out":9003,"duration_ms":87283,"concrete_test":"Solve the full nonlinear Young-Laplace equilibrium shape of a single 2D bubble with three contacts arranged at 120 degrees with equal dimensionless force f, using a direct numerical solver such as Surface Evolver or a boundary-integral method, with contacts represented as finite flat segments. Extract the contact-point displacement x_ij along each contact direction and compare it with Eq. (10) for f = 0.05, 0.17, 0.3, and 0.5, where f = 0.17 is a rough estimate of the force at phi = 0.12 obtained from Eq. (19). Repeat with unequal forces and with one neighbour of radius R_j = 0.7 R_i to mimic polydispersity. If the error exceeds the approximately 2% benchmark of Figure 3 at the same f, Eq. (10)'s multi-contact superposition is not validated and the headline claim would need to be narrowed; if the error remains below that benchmark, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the implementation produces polydisperse 2D foam structures consistent with direct calculations. The most load-bearing step is Eq. (10), which constructs the deformation x_ij at contact j of bubble i as a linear superposition of single-contact solutions plus a pairwise quadratic correction. The quadratic correction is derived in Section 2.2 for an isolated two-bubble contact, and Figure 3 validates that pair formula to under 2% error up to f = 0.5. But in a disordered foam each bubble has three or more contacts, and Eq. (10) assumes that the only O(f^2) contribution to x_ij comes from F_ij^2, with no cross-contact quadratic terms and no modification of the single-contact g(theta) pattern by other forces. The authors explicitly acknowledge in Section 2.2 that 'the linearised theory contains errors of order f^2 from the outset which we do not claim to eliminate.' The existing validation does not isolate this assumption: the Princen comparison in Figure 6(a) is monodisperse and highly symmetric, and the Plat comparison is a single 10-bubble system in which many effects vary simultaneously, with contact changes compared only after a phi shift. Thus the accuracy of Eq. (10) for a general multi-contact polydisperse bubble is the unvalidated link between the pair-level theory and the headline foam-scale claim. If cross-contact f^2 terms are comparable to the retained F_ij^2 term, the simulated structures, energies, and scaling exponents could shift, especially at phi = 0.12 where forces are largest.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an implementation of Morse-Witten theory for polydisperse two-dimensional foams. It derives a two-bubble force-deformation relation (Eq. 7) with an O(f^2) correction, extends it to multi-contact bubbles via Eq. (10), and introduces an iterative fixed-point scheme that updates forces, deformations, and centre positions while monitoring the contact network. Validation includes a two-bubble separation test (relative error <2% for f <= 0.5), a monodisperse hexagonal case that reproduces Princen's exact excess energy, and a 10-bubble comparison with the Plat simulator in which the sequence of contact changes is nearly reproduced after a liquid-fraction shift. For disordered 100-bubble systems the paper reports epsilon ~ Delta-phi^2.2 and Delta-Z ~ Delta-phi^0.52, and it outlines a 3D extension.","tokens_in":10213,"tokens_out":9424,"duration_ms":94961,"significance":"If the multi-contact approximation is accepted, the paper fills a genuine gap: existing 2D foam simulators such as Plat struggle near the wet limit, while the proposed method is reported to converge for liquid fractions down to about 0.12. The two-bubble test and the exact Princen comparison are useful quantitative anchors, and the algorithm is described in enough detail to be reimplemented. The reported scaling laws are relevant to the debate on foam jamming, and the explicit acknowledgement of the theory's O(f^2) limitations is a strength. However, the significance is limited by the thin multi-contact validation and the absence of uncertainty estimates on the scaling exponents, which is why the manuscript needs further work.","major_comments":[{"comment":"The load-bearing step from the validated two-bubble law to a foam is the superposition in Eq. (10), whose only O(f^2) term is proportional to F_ij^2 and contains no cross-contact quadratic terms. The numerical check in Fig. 3 validates Eq. (9) for an isolated pair over 0 < f < 0.5, but no controlled test exercises a bubble with three or more contacts at comparable forces. The 10-bubble Plat comparison in Section 5 is a single disordered configuration and is reported only as a sequence of contact changes with a systematic shift Delta-phi ~ 0.01, so it cannot quantify the omitted cross-contact O(f^2) errors. The authors' own statement in Section 2.2 that 'the linearised theory contains errors of order f^2 from the outset' is a relevant limitation, but the magnitude of the multi-contact error remains unquantified, and this is the key link between the pair-level theory and the paper's central claim. A direct comparison for a small cluster, such as a three-bubble triangle, against Surface Evolver or Plat would provide a concrete bound on this error.","section":"Section 2.2, Eq. (10), Fig. 3"},{"comment":"The reported scaling laws epsilon ~ Delta-phi^2.2 and Delta-Z ~ Delta-phi^0.52 are presented without uncertainties or a statement of the fitting procedure. The phi_c values used to define Delta-phi are obtained by fitting the lowest eight points above an arbitrary threshold of 10^-4, with no sensitivity analysis; different choices of the fitting window and threshold could change the exponents, which are the main quantitative results of the disordered-foam simulations. The authors should report error bars, for example from a bootstrap over the 1000 systems, and show the stability of the exponents with respect to the fitting window.","section":"Section 5, Figs. 6(b) and 7"},{"comment":"The text says the simulations were run for liquid fractions from 0.18 to 0.12 and that the expected critical liquid fraction is phi_c ~ 0.16, but then reports the fitted average as 0.843 +/- 0.003. This is internally inconsistent: 0.843 is a packing fraction, not a liquid fraction, and if it were used in Delta-phi = phi_c - phi the plotted range in Fig. 6(b) would not be accessible. The manuscript must clarify which convention is used and correct this value, because the scaling analysis depends on the definition of Delta-phi.","section":"Section 5, phi_c determination"},{"comment":"The validation against Plat is limited to a single 10-bubble system and to the sequence of contact changes, compared only after a shift in liquid fraction of roughly Delta-phi ~ 0.01. There is no quantitative comparison of bubble positions, coordination numbers, or energies. A few small-system comparisons with a defined error metric would substantially strengthen the claim that the Morse-Witten implementation reproduces direct calculations for polydisperse foams.","section":"Section 5, Fig. 5"}],"minor_comments":[{"comment":"The caption refers to 'Eqn (??)'; this should be Equation (9).","section":"Fig. 3 caption"},{"comment":"The phrase 'termsoforder f 2 orhigherneedtobeconsidered' is missing spaces and should be typeset as 'terms of order f^2 or higher need to be considered'.","section":"Section 2.2, after Eq. (9)"},{"comment":"The update F^(n+1) = a F^(n+1) + (1-a) F^(n) uses the same superscript for the new force and the mixed force; introducing an intermediate symbol for the undamped update would remove the ambiguity.","section":"Eq. (14)"},{"comment":"The symbol N is used both for the number of bubbles and for the number of contacts in the sum over k; a different index or letter would avoid confusion.","section":"Eq. (10) and Section 3.1"}],"recommendation":"major_revision","confidential_remarks":"The core two-bubble derivation and the exact Princen comparison appear sound, and the 0.843 value is most plausibly a typo or a packing-fraction convention, but it must be corrected before publication. The main risks are the unquantified multi-contact O(f^2) error in Eq. (10) and the missing uncertainty analysis for the scaling exponents. I see no attribution or scope problem; the paper is a reasonable methods contribution to cond-mat.soft."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about arXiv:1908.03768. First, it is the first working implementation of Morse-Witten theory for polydisperse 2D foams, with a contact-network update scheme that handles topological changes. That fills a genuine gap in wet-limit foam simulation. Second, the central force-deformation relation for multibubble contacts, Eq. (10), is only tested for a single pair of bubbles. The multi-contact superposition is an assumption, and the authors say so. That is the soft spot to focus on.\n\nWhat the paper does well: it gives a clean derivation of the two-bubble contact law with a quadratic correction, checks it numerically against the linearized profiles to under 2% error up to f=0.5, and fixes a small error in the authors' earlier paper [5] (the f^2 coefficient). The monodisperse energy matches Princen's exact result in the wet limit, and the 10-bubble comparison with Plat reproduces the same contact-change sequence, just shifted by about Δφ=0.01. The algorithm also converges at liquid fractions above 0.12, where most Plat runs fail. For a simulation paper, that is real evidence.\n\nThe soft spots are mostly around the multi-contact step. Eq. (10) superposes single-contact linearized solutions and adds a pairwise quadratic correction, but no test exercises a bubble with three or more contacts. In a disordered foam that is the normal case, so the f^2 error could be larger than the two-bubble test suggests. The authors acknowledge the linearized theory has f^2 errors 'from the outset,' but they do not quantify them for multi-contact configurations. The Princen comparison is monodisperse and symmetric; the Plat comparison is a single 10-bubble system. So the headline scaling results (ε∝Δφ^2.2, ΔZ∝Δφ^0.52) rest on the unvalidated multi-contact assumption. Also, those exponents have no error bars, and φc is obtained from a per-system fit without sensitivity analysis. Missing code limits reproducibility. These are fixable for a revised version.\n\nNone of this kills the paper. The central claim is stated carefully: the structures are consistent with direct calculations within the limitations of the theory. That claim holds, with the caveat that 'limitations' include the unquantified multi-contact error. The paper deserves a serious referee. I would send it to review, and ask for one additional test that isolates a three-contact or four-contact configuration, plus error bars on the exponents and, ideally, code release.\n\nWho is this for: people working on wet foam or emulsion simulation, and anyone interested in jamming of soft disks. I would not cite it in my own work unless I move into foam simulation, but I would put it in front of the reading group.","headline":"First real polydisperse Morse-Witten foam simulator, with solid pair-level checks but a multi-contact superposition that is tested only in two-bubble geometry; worth refereeing after some extra validation.","tokens_in":10760,"tokens_out":2405,"would_cite":false,"duration_ms":24261,"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":"This paper extends Morse-Witten bubble-contact theory to polydisperse foams and turns it into a working 2D simulation for the wet limit.","keywords":["Morse-Witten theory","wet foams","polydisperse bubbles","two-dimensional foams","jamming transition","force networks","liquid fraction","foam simulation"],"falsifier":"A direct numerical solution of the Young-Laplace equation for two unequal bubbles pressed together across the reported force range $0<f<0.5$ would test Equation (9): if the center-to-center separation deviates by more than the claimed $\\sim2\\%$, or if the deviation grows with the number of contacts in a larger foam, the superposition assumption is not safe.","tokens_in":9693,"feed_emoji":"🫧","tokens_out":8465,"duration_ms":79859,"temperature":0.7,"pith_summary":"Morse and Witten's theory gives analytic formulas for how a bubble deforms under a contact force, but until now it has been applied only to monodisperse or nearly monodisperse foams. This paper shows how to include bubbles of arbitrary different sizes, adding a quadratic-in-force correction to the deformation law for a contact between two unequal bubbles. It then packages this into an iterative simulation that follows the changing contact network of a disordered foam. The resulting equilibrated 2D structures agree with direct calculations for liquid fractions above about 0.12, a regime in which at least 80% of standard geometric foam simulations fail to converge. The same reasoning is carried through to three dimensions, where the path to a future implementation is laid out.","feed_headline":"Morse-Witten theory reaches polydisperse wet 2D foams","feed_subtitle":"New simulation matches direct calculations and gives jamming scalings in a regime where other schemes fail to converge.","key_machinery":"The load-bearing object is the generalized force-deformation relation (Equation 10), which superposes linearized single-contact Morse-Witten deformation profiles from all contacts of a bubble and adds a quadratic correction that encodes polydispersity. The single-contact profile comes from solving the linearized Young-Laplace equation for a 2D bubble pressed by a point force while keeping the centroid fixed, giving the deformation function $g(\\theta)=(\\pi-\\theta)\\sin\\theta - \\cos\\theta/2 - 1$. This relation turns a foam into a central-force network on centroid coordinates, and the paper surrounds it with an iteration that updates deformations, forces, and positions while monitoring overlaps and removing spurious negative-force contacts.","core_discovery":"The paper claims that the Morse-Witten description of bubble contacts can be made genuinely polydisperse: when two bubbles of radii $R_1$ and $R_2$ meet, the inward displacement of each contact point is not the linear Morse-Witten result but includes a term of order $f^2$ involving $2 + R_1/R_2 + R_2/R_1$, and this corrected law (Equation 10) can be applied to every contact of every bubble. The resulting system, in which bubbles are represented by centroid positions, contact forces, and contact deformations, satisfies force-deformation consistency, deformation-displacement consistency, action-reaction, and force equilibrium through an iterative scheme that adds and removes contacts as the foam evolves. The authors find that their equilibrated structures reproduce contact changes seen in direct calculations, match the exact result for ordered monodisperse wet foams, and give a critical liquid fraction of $\\varphi_c \\simeq 0.843$, consistent with published values. Near jamming they report an excess energy scaling $\\varepsilon \\propto \\Delta\\varphi^{2.2}$ and an excess coordination number scaling $\\Delta Z \\propto \\Delta\\varphi^{0.52}$.","pith_inferences":["One testable extension is to seed a two-bubble calculation with unequal radii and measure whether the predicted force-separation law, Equation (7), holds for polydispersities up to $\\Delta R/R_0 = 0.45$; the paper reports deviations under 2% only in its two-bubble test.","The discrepancy between the $\\Delta Z \\propto \\Delta\\varphi^{0.52}$ scaling found here and the linear scaling reported from the geometric simulator suggests the linear law may be an artefact of contact detection; comparing both methods on identical initial packings would resolve it.","The force networks generated by this method could be used to study the tail of the contact-force distribution in polydisperse foams, a question this paper only begins with its preliminary distribution."],"forward_implications":["Foam simulations can now be run reliably in the wet limit, at liquid fractions above about 0.12, where the standard geometric simulation approach most often fails to converge.","For ordered monodisperse foams the model reproduces the exact excess-energy result, so the scheme can serve as a validated baseline for wet foam mechanics.","Near jamming the simulation gives specific scalings, excess energy $\\varepsilon \\propto \\Delta\\varphi^{2.2}$ and excess coordination number $\\Delta Z \\propto \\Delta\\varphi^{0.52}$, plus a critical liquid fraction consistent with published values.","The same machinery carries over to three dimensions; the paper writes down the 3D force-deformation relation, including its logarithmic term, as the foundation of a future implementation."],"supporting_citations":[{"why":"supplies the original single-contact elastic theory that this paper extends to polydisperse foams.","marker":"[1]"},{"why":"provides the 2D single-bubble deformation function $g(\\theta)$ and the linear contact law used as the starting point.","marker":"[5]"},{"why":"shows an earlier application of Morse-Witten theory to crystalline emulsions, the restricted setting this paper generalises.","marker":"[3]"},{"why":"introduces the iterative treatment of many-body contact forces and provides the monodisperse comparison.","marker":"[4]"},{"why":"supplies the reference 2D foam simulator that generates test packings and provides the comparison structures.","marker":"[10]"},{"why":"gives the exact excess-energy result for ordered monodisperse wet foams used to validate the model.","marker":"[17, 18]"},{"why":"reports the linear $\\Delta Z$ scaling that this paper's $\\Delta\\varphi^{0.52}$ result challenges.","marker":"[16]"},{"why":"provides the soft-disk jamming scaling $\\Delta Z \\propto \\Delta\\varphi^{0.5}$ with which the present coordination-number result agrees.","marker":"[22]"}],"fun_headline_variants":["Morse-Witten theory goes polydisperse for 2D foams","Polydisperse wet foams now match Morse-Witten theory","2D foam simulation tackles polydisperse Morse-Witten","Corrected bubble contacts: polydisperse wet 2D foam model","Polydisperse correction extends Morse-Witten to wet 2D foams"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The construction rests on adding up the small deformations each contact would cause if isolated, trusting that the errors left by this linearization remain small even when many unequal bubbles press on each other near jamming.","fun_headline_variants_meta":{"raw":{"variants":["Morse-Witten theory goes polydisperse for 2D foams","Polydisperse wet foams now match Morse-Witten theory","2D foam simulation tackles polydisperse Morse-Witten","Corrected bubble contacts: polydisperse wet 2D foam model","Polydisperse correction extends Morse-Witten to wet 2D foams"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000836,"raw_usage":{"total_tokens":3616,"prompt_tokens":886,"completion_tokens":2730,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":502,"completion_tokens_details":{"reasoning_tokens":2632}},"tokens_in":502,"tokens_out":2730,"duration_ms":20251,"temperature":1.0,"reasoning_tokens":2632,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:02:38.382682+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical solution of the Young-Laplace equation for two unequal bubbles pressed together across the reported force range $0<f<0.5$ would test Equation (9): if the center-to-center separation deviates by more than the claimed $\\sim2\\%$, or if the deviation grows with the number of contacts in a larger foam, the superposition assumption is not safe.","supporting_citations":[{"cited_title":"Morse and T","cited_arxiv_id":null,"evidence_quote":"supplies the original single-contact elastic theory that this paper extends to polydisperse foams."},{"cited_title":"Weaire, R","cited_arxiv_id":null,"evidence_quote":"provides the 2D single-bubble deformation function $g(\\theta)$ and the linear contact law used as the starting point."},{"cited_title":"Buzza and M","cited_arxiv_id":null,"evidence_quote":"shows an earlier application of Morse-Witten theory to crystalline emulsions, the restricted setting this paper generalises."},{"cited_title":"Höhler and S","cited_arxiv_id":null,"evidence_quote":"introduces the iterative treatment of many-body contact forces and provides the monodisperse comparison."},{"cited_title":"Bolton,Software PLAT: A computer code for simulating two-dimensional liquid foams, https://github.com/fbolton/plat (1996)","cited_arxiv_id":null,"evidence_quote":"supplies the reference 2D foam simulator that generates test packings and provides the comparison structures."},{"cited_title":"Winkelmann, F","cited_arxiv_id":null,"evidence_quote":"reports the linear $\\Delta Z$ scaling that this paper's $\\Delta\\varphi^{0.52}$ result challenges."},{"cited_title":"O’Hern, L.E","cited_arxiv_id":null,"evidence_quote":"provides the soft-disk jamming scaling $\\Delta Z \\propto \\Delta\\varphi^{0.5}$ with which the present coordination-number result agrees."}],"review_version":1}