{"id":"59bb5309-bec9-4aa0-8420-eccaced0c53a","arxiv_id":"2501.07664","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a 2D hyperelastic solid, two pressurized cylindrical cavities attract at low pressure but repel beyond a critical separation at high pressure, a behavior linear elasticity cannot produce.","lead":"This computational study shows that two pressurized cylindrical cavities in a soft hyperelastic material attract each other at low pressures but can repel at higher pressures beyond a critical separation. The transition, which does not exist in linear elasticity, could help engineers control cavity positions in soft materials using pressure and material stiffness.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Repulsive regime may be a finite-domain artifact: the traction-free outer boundary at 50R is never tested by a nonlinear domain-size convergence study.","rationale":"The reader correctly identified both the undocumented curve fit and the finite outer boundary as concerns. I focus on the finite boundary as the single most load-bearing issue because it threatens the qualitative claim itself, not just the quantitative extraction of the transition. The fit affects the location of the zero crossing and the reported critical separations, but if the raw energy data in Appendix E already show a maximum, the existence of a sign change in the driving force is not created by the fit. By contrast, the truncated traction-free boundary at 50R is a modeling assumption that is not validated for the high-pressure, large-deformation regime: the deformed cavities are only about 10 radii from the boundary at the largest separations, and the boundary-induced contribution to the interaction force can be of the same order as the true two-cavity interaction at large η/R. The manuscript's mesh convergence study varies only element size at a fixed outer radius, and the linear-elastic validation cannot expose this effect because it is negligible in the small-deformation regime. The proposed domain-size convergence test would directly settle whether the repulsive regime is intrinsic or an artifact of the finite computational domain. I therefore keep the reader's CONDITIONAL verdict unchanged: the paper is promising but incomplete without this check. I credit the paper for the single-cavity analytical validation, the linear-elastic series validation, and the careful mesh refinement, all of which are genuine supporting evidence, but none of them addresses the missing infinite-domain limit at high pressure.","tokens_in":11645,"tokens_out":17050,"duration_ms":184143,"concrete_test":"For the neo-Hookean model at P/μ=1.5 and 1.75, repeat the two-cavity FE calculation with outer boundary radii R_o/R = 50, 100, 200, and 500, keeping the mesh size near the cavities the same (e.g., element size min/R=0.02 as in Mesh 3). From the raw, unfitted potential energies, compute dΠ/dη by finite differences over η/R points bracketing the suspected maximum, or with the same fitting procedure but a documented functional form and residual analysis. Record the maximum positive driving force and the critical separation as functions of R_o/R, and check whether Π(η/R=20) approaches the infinite single-cavity value as R_o grows. If the positive-driving-force region persists and converges with increasing R_o, the repulsion is real; if it shrinks or vanishes, the finite traction-free boundary at 50R is the cause.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.4 fixes the outer boundary at radius 5 with R=0.1, i.e., R_o/R=50, and the infinite-domain claim is never checked against a larger domain in the nonlinear regime. This matters because the central observable, F=-dΠ/dη, changes sign only at high pressure, where the deformation is large: the infinite-medium single-cavity solution (Eq. B.43) gives λ_a≈3.6 at P/μ=1.75, so at η/R=20 each cavity surface is only about 10 deformed radii from the traction-free outer boundary. As η/R increases from 2.2 to 20, the cavity centers move from ±1.1R to ±10R, i.e., each cavity approaches the free boundary by roughly 18% of its initial standoff. A free boundary is softer than the infinite medium and tends to attract a pressurized cavity, so this 'image' effect adds a contribution to Π that decreases with η and hence a positive contribution to F. That is precisely the direction of the claimed repulsion. The mesh refinement study in Appendix C fixes R_o and therefore cannot detect this; the linear-elastic validation is insensitive to it because the effect is O((R/R_o)^2), while at high pressure the expanded cavities make the effective standoff only about 10 deformed radii. The non-monotonic energy seen in Appendix E could therefore be a finite-domain artifact rather than an intrinsic two-cavity interaction in an infinite medium.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two equal pressurized cylindrical cavities in a 2D infinite hyperelastic medium under plane strain, using ABAQUS finite elements for neo-Hookean, Mooney-Rivlin, and Arruda-Boyce materials. The interaction is quantified by the driving force F = -dΠ/dη, where Π is the total potential energy and η is the center-to-center separation. The authors validate their computational pipeline against a corrected bipolar-coordinate linear-elastic series solution and against an analytical single-cavity neo-Hookean solution, and they report mesh-convergence studies and raw energy data. Their central claim is that at low pressure-shear modulus ratios the interaction is always attractive, while at higher ratios the energy becomes non-monotonic: below a critical separation the force is attractive and above it repulsive, a regime absent in linear elasticity. The paper also reports that Mooney-Rivlin results are nearly independent of the stiffening parameter α, while Arruda-Boyce results depend on λ_m, with smaller λ_m giving a smaller critical separation.","tokens_in":11862,"tokens_out":4205,"duration_ms":42050,"significance":"If the repulsive regime is genuine, the result is counterintuitive and relevant for soft-material design, phase-separation problems, and geomechanics. The paper has clear strengths: it validates the numerical method against an external linear-elastic series solution and a single-cavity analytical benchmark, reports raw potential-energy data in Appendix E, and scans material parameters rather than fitting them to obtain the claimed behavior. The central limitation is that the main observable is obtained by differentiating an unspecified smooth-curve fit to discrete energy values, and the infinite-domain claim is never checked by a nonlinear outer-boundary convergence study. Because both issues bear directly on whether the repulsive regime is a property of the infinite hyperelastic medium or an artifact of the numerical procedure, the result is defensible only after those points are addressed.","major_comments":[{"comment":"The driving force F = -dΠ/dη is computed from a smooth curve fitted to discrete finite-element energy values, but the paper does not state the fitting function, the smoothing procedure, or the sensitivity of the derivative to the fit. Since the existence of the repulsive regime and the value of (η/R)_critical are defined by a sign change of this numerical derivative, the reader cannot tell whether the transition is a feature of the underlying FE solution or an artifact of the chosen interpolant. Please report the fitting form, fit residuals for each P/μ, and demonstrate that the sign change and critical separation are insensitive to the fitting procedure.","section":"Section 2.4, Figures 4–5"},{"comment":"The infinite-medium claim is not checked for outer-boundary effects in the nonlinear regime. All simulations fix R_o/R = 50, and the mesh-convergence study in Appendix C varies mesh density at this fixed outer radius. At P/μ = 1.75, Eq. (B.43) gives λ_a ≈ 3.6, so at η/R = 20 the cavity wall is only about 10 deformed radii from the traction-free outer boundary; moreover, as η/R increases from 2.2 to 20, the cavities move from ±1.1R to ±10R, approaching that boundary by roughly 18% of its initial standoff. A traction-free boundary is softer than the infinite medium and attracts a pressurized cavity, producing a contribution to F that is positive—the same direction as the claimed repulsion. A domain-size study with R_o/R = 100 and 200 at high P/μ and large η/R is required to rule out a finite-domain artifact; the linear-elastic validation cannot detect this effect because it is O((R/R_o)^2) at low pressure and small deformation.","section":"Section 2.4, Appendix C"}],"minor_comments":[{"comment":"The sentence 'Effect of strain stiffening on these interactions are also analyzed' has a subject-verb agreement error.","section":"Abstract and Section 4"},{"comment":"The caption does not identify which curve corresponds to which P/μ value, making it difficult to connect the plotted potential-energy and driving-force curves to the pressure levels discussed in the text.","section":"Figure 4"},{"comment":"The text says the corrected solution and its MATLAB implementation are included, but no MATLAB code appears in the appendix; only the analytical formulas are given.","section":"Appendix A"},{"comment":"The notation mixing summation and integration symbols (χ0X, πX) is confusing; please clarify that this is a numerical quadrature rule.","section":"Equation (26)"},{"comment":"The axes are unlabeled; the plotted 'Area of deformed cavity' should state the normalization and whether it is the deformed area or the area change relative to the undeformed cavity.","section":"Figure 12"},{"comment":"The statement that 'mesh independence was found to hold across pressures and models' is not supported by the displayed results, which show only a neo-Hookean model at P/μ = 1.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The two major comments are load-bearing: the fitting procedure is the only route to the driving force, and the finite-domain concern directly targets the repulsive regime. If the authors can show that the sign change persists with a larger outer radius and is robust to the fitting choice, the result would be a valuable contribution; otherwise the central claim may be an artifact. The paper would also benefit from making the raw data and fitting code available."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the Saeedi–Kothari paper on pressurized cavity pairs in hyperelastic media. The headline finding—a pressure-dependent transition from attraction to repulsion—is genuinely new as far as the cited literature goes, and the paper is cleanly executed: the linear-elastic bipolar validation, the single-cavity neo-Hookean check, and the mesh convergence study all look solid. The authors are upfront that the linear problem is always attractive and that the repulsive regime only appears in nonlinear models at P/μ ≳ 1.\n\nThe soft spot is the finite outer boundary. They fix R_o/R = 50 and never test a larger domain in the nonlinear regime. That matters more than it looks: at P/μ = 1.75 the single-cavity solution gives a deformed radius of ~3.6R, so the traction-free outer boundary is only about 13 deformed radii from the cavity surface at the largest separation they simulate. A free boundary is more compliant than the infinite medium and will attract a pressurized cavity, adding a positive contribution to the driving force—the same sign as the claimed repulsion. The effect grows with η because the cavity approaches the boundary, and it grows with pressure because the cavity expands. The linear validation can't catch it because cavities stay small; the mesh convergence can't catch it because R_o is fixed. This is a real hole in the evidence.\n\nA second, smaller issue: the driving force is the derivative of a smooth-curve fit to discrete energy data, but the fitting function and sensitivity are not reported. The raw energy data in Appendix E shows the same non-monotonic shape, so this is probably not the main risk, but it lowers reproducibility.\n\nOverall: the paper is worth a serious referee. The question is important for soft-matter design, and the computational setup is competent. I would ask for a domain-size convergence study at P/μ = 1.5 and 1.75 (e.g., R_o/R = 100, 200) to show that the repulsive regime survives, and for a transparent description of the fitting. If the repulsion vanishes in the infinite-domain limit, the central claim is an artifact; if it persists, this is a useful result.\n\nRecommendation: send to peer review with a request for major revision on exactly this point.","headline":"The paper's new repulsive regime is worth taking seriously, but the finite outer boundary at 50R makes the finding unproven until a nonlinear domain-size convergence check is done.","tokens_in":12416,"tokens_out":5112,"would_cite":false,"duration_ms":48570,"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":"Two pressurized cavities in a hyperelastic medium switch from attraction to repulsion once the pressure-to-shear-modulus ratio exceeds a critical value.","keywords":["hyperelasticity","pressurized cavities","cavity interaction energy","driving force","neo-Hookean material","Mooney-Rivlin material","Arruda-Boyce material","strain stiffening"],"falsifier":"Recompute the neo-Hookean two-cavity energy at $P/\\mu=1.5$ with the separation sampled every $\\Delta(\\eta/R)=0.01$ around the reported $(\\eta/R)_{\\rm critical}$, and take the driving force directly from finite differences of the raw, unsmoothed finite-element energies; if the derivative never changes sign, the repulsion is an artifact of the smoothing.","tokens_in":11384,"feed_emoji":"🧲","tokens_out":7919,"duration_ms":72617,"temperature":0.7,"pith_summary":"The paper establishes that two equal pressurized cylindrical cavities in a 2D hyperelastic medium can repel each other, something linear elasticity says never happens. Below a pressure-to-shear-modulus ratio of about one, the interaction is always attractive; above that ratio, the energy landscape becomes non-monotonic, so cavities closer than a critical separation attract while wider-separated cavities repel. The same transition appears in neo-Hookean, Mooney-Rivlin, and Arruda-Boyce materials, with strain stiffening in the Arruda-Boyce model moving the critical separation. If correct, this gives a quantitative lever for arranging or stabilizing cavities in soft materials and underground structures by choosing pressure and material parameters rather than only geometry.","feed_headline":"Pressurized cavity pairs flip from attraction to repulsion","feed_subtitle":"Above a critical pressure, soft elastic materials push widely spaced cavities apart—a tunable knob for bubble and tunnel control.","key_machinery":"The central object is the configurational driving force $F = -d\\Pi_{\\rm eqm}/d\\eta$, the negative derivative of the equilibrium total potential energy with respect to the separation $\\eta$ between cavity centers. The authors compute $\\Pi_{\\rm eqm}$ from finite-element solutions as the strain energy minus the pressure-volume work of the two cavities, fit a smooth curve through the discrete energy values, and read the nature of the interaction from the sign of the derivative. The nonlinear constitutive models are what let the energy landscape acquire a maximum at high pressure, while the linear-elastic landscape stays monotone; strain stiffening, especially in the Arruda-Boyce model, shifts the location of that maximum.","core_discovery":"For two equal pressurized circular cavities in an infinite, incompressible hyperelastic medium under plane strain, the total potential energy as a function of center-to-center separation is monotonically decreasing at low pressures ($P/\\mu \\lesssim 1$), so the configurational driving force $F = -d\\Pi_{\\rm eqm}/d\\eta$ is always negative and the cavities attract. At higher pressures ($P/\\mu \\gtrsim 1$), the energy curve develops a maximum at a critical separation $(\\eta/R)_{\\rm critical}$; the driving force is negative for closer spacings and positive for wider ones, meaning sufficiently separated cavities repel. The critical separation shifts with pressure and with strain-stiffening parameters, and the transition appears in neo-Hookean, Mooney-Rivlin, and Arruda-Boyce models, with Arruda-Boyce showing the strongest dependence on its limiting-stretch parameter.","pith_inferences":["If the transition survives direct differentiation of the raw energies, the energy maximum implies an unstable separatrix; in multi-cavity arrays this could organize cavities into repulsive lattices that resist coalescence without walls.","A natural extension beyond the paper is to spherical cavities in 3D, where the same energy-maximum mechanism may regulate cavitation damage and phase-separated droplet spacing in soft solids.","The finite outer radius used to approximate infinity should be tested against larger domains at high $P/\\mu$; if the energy peak shifts or disappears with domain size, the infinite-domain interpretation needs revisiting.","The exact value of the reported critical separation is the softest quantitative output because it comes from a smoothed derivative; direct finite differences of the raw finite-element energies would settle how much of the phase diagram is robust."],"forward_implications":["Below $P/\\mu \\simeq 1$ the two-cavity system behaves like the linear-elastic case: attraction at every separation, so no equilibrium spacing exists.","Above the threshold there is an unstable equilibrium at $(\\eta/R)_{\\rm critical}$; cavities perturbed inward attract and coalesce, while cavities perturbed outward separate further.","The critical separation depends on pressure: higher pressure moves the transition, and in Arruda-Boyce materials a smaller limiting stretch $\\lambda_m$ lowers the critical separation.","Because all three constitutive models show the transition, repulsion at large separation is a generic consequence of nonlinear elasticity rather than a quirk of one material law.","The energy-based driving-force formulation gives a direct design target: match the critical separation to a desired tunnel or bubble spacing by choosing the pressure ratio."],"supporting_citations":[{"why":"Introduces the bipolar-coordinate solution technique used to validate the linear-elastic limit.","marker":"[8]"},{"why":"Supplies the series solution for stresses around two equal circular holes that the paper corrects and extends.","marker":"[11]"},{"why":"Gives the earlier interaction-energy calculation whose repulsion conclusion this paper shows to be wrong, motivating the energy-based correction.","marker":"[14]"},{"why":"Provides the neo-Hookean strain-energy form used in the nonlinear simulations.","marker":"[15]"},{"why":"Provides the Mooney-Rivlin model used as a second nonlinearity benchmark.","marker":"[16]"},{"why":"Supplies the Arruda-Boyce model and its consistency condition, the strain-stiffening law tested here.","marker":"[17]"},{"why":"Parameterizes strain stiffening in the Mooney-Rivlin model through the parameter $\\alpha$.","marker":"[24]"},{"why":"Supplies the finite-element solver used to generate all the nonlinear potential-energy data.","marker":"[25]"}],"fun_headline_variants":["Cavity pairs: pressure flips attraction to repulsion","Soft solids: pressure switches cavity interaction sign","A critical pressure turns cavity attraction repulsive","Pressure tunes cavity pairs from pull to push","Hyperelastic cavities: repulsion beyond critical pressure"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The sign change in the driving force is read from the derivative of a smooth curve fitted through discrete finite-element energy values, so the entire repulsive regime depends on that fitted curve faithfully representing the true energy landscape.","fun_headline_variants_meta":{"raw":{"variants":["Cavity pairs: pressure flips attraction to repulsion","Soft solids: pressure switches cavity interaction sign","A critical pressure turns cavity attraction repulsive","Pressure tunes cavity pairs from pull to push","Hyperelastic cavities: repulsion beyond critical pressure"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000153,"raw_usage":{"total_tokens":1161,"prompt_tokens":852,"completion_tokens":309,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":468,"completion_tokens_details":{"reasoning_tokens":238}},"tokens_in":468,"tokens_out":309,"duration_ms":3703,"temperature":1.0,"reasoning_tokens":238,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:39:21.995734+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the neo-Hookean two-cavity energy at $P/\\mu=1.5$ with the separation sampled every $\\Delta(\\eta/R)=0.01$ around the reported $(\\eta/R)_{\\rm critical}$, and take the driving force directly from finite differences of the raw, unsmoothed finite-element energies; if the derivative never changes sign, the repulsion is an artifact of the smoothing.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the bipolar-coordinate solution technique used to validate the linear-elastic limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the series solution for stresses around two equal circular holes that the paper corrects and extends."},{"cited_title":"Davanas, Analysis of elastic interactions between holes, Journal of Materials Science 27 (1992) 1589–1598","cited_arxiv_id":null,"evidence_quote":"Gives the earlier interaction-energy calculation whose repulsion conclusion this paper shows to be wrong, motivating the energy-based correction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the neo-Hookean strain-energy form used in the nonlinear simulations."},{"cited_title":"Mooney, A theory of large elastic deformation, Journal of applied physics 11 (9) (1940) 582–592","cited_arxiv_id":null,"evidence_quote":"Provides the Mooney-Rivlin model used as a second nonlinearity benchmark."},{"cited_title":"Arruda, M","cited_arxiv_id":null,"evidence_quote":"Supplies the Arruda-Boyce model and its consistency condition, the strain-stiffening law tested here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Parameterizes strain stiffening in the Mooney-Rivlin model through the parameter $\\alpha$."},{"cited_title":"Syst` emes, Abaqus 2022, simulia, Dassault Syst` emes, Providence, RI, USA (2022)","cited_arxiv_id":null,"evidence_quote":"Supplies the finite-element solver used to generate all the nonlinear potential-energy data."}],"review_version":1}