{"id":"d991747b-4016-4fdf-a542-5aff7161c77b","arxiv_id":"2506.20745","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Mushroom fibrils with wide, thin caps keep pull-off strength near the ideal value over a wide range of the parameter χ, but central adhesion defects can cut that strength by 65% or more.","lead":"This paper uses finite element simulations with a Dugdale cohesive zone law to map the pull-off strength of mushroom-shaped adhesive fibrils as a function of geometry and an adhesion parameter. It shows that wide, thin caps keep strength high by initiating detachment at the fibril center, but that even tiny central defects can sharply reduce that strength.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dimensional reduction behind Eq. (10) is verified only for c=0 envelopes; the c>0 states that fix the defect-case strengths in Figs. 7–9 are never tested at a second σ0, leaving the claimed universality of the defect design curves unproven.","rationale":"The reader's weakest_assumption is the right target. I agree with it, and I want to sharpen it: the gap is not merely 'no check for defect geometries'; it is that the σ0-checks were performed on the c=0 envelope, which is the right set for defect-free pull-off points but the wrong set for the defect problem. The defect-case pull-off strengths in Figs. 7–9 are taken from c>0 equilibrium states, none of which has been recomputed at a second σ0. The linear-superposition argument for the Dugdale states makes the reduction plausible, and the high-χ agreement with the Kassir–Bregman-based Eq. (8) is genuine supporting evidence (modulo the Eq. (8)/appendix consistency issue the reader also raises). But neither closes the verification gap at intermediate χ, where Figs. 7–9 show the strongest defect sensitivity. Since the FE code and meshes already exist, the test is cheap, which is why conditional acceptance is the right verdict: the paper should not be ACCEPT until this is done, and it should not be REJECT because the concern is narrow and addressable. I weighed two other candidate concerns. (1) Mesh resolution at the χ≈8000 edge-detachment boundary for β=1.25, ξ=0.01: the appendix's validity bound is derived from a central penny-shaped-crack strip-yield model (Eqs. A1–A3), not from an edge Dugdale zone, and the localized refinement study is not specified as to which case it covered. This is a real secondary risk to the specific number 0.78σ0 at χ≈8000, but the qualitative claim (mushroom greatly outperforms the punch at high χ) is robust to it. (2) The printed Eq. (8) inconsistency: significant as a paper-quality issue but not load-bearing for the physics if the plotted dashed lines were generated from the correct expression. The dimensional-reduction gap is more fundamental because it underpins the universality of the entire design-method claim, including the extension to defects that the statistical array application depends on. Credit is due where the paper is strong: the snapshot construction is a sound way to obtain unstable post-peak branches; the c=0 σ0-independence is directly verified for two geometries that bracket both detachment modes; the comparison with Del Campo et al. is honest (6× predicted vs. 20× reported, with a plausible explanation); and the refusal to force fits onto Fig. 3 avoids overparameterization. None of these closes the c>0 gap, so the reader's CONDITIONAL verdict remains the appropriate one.","tokens_in":23596,"tokens_out":39860,"duration_ms":417505,"concrete_test":"Recompute the central-defect pull-off curve for β=1.25, ξ=0.01 and 2c/Df=0.005 (Fig. 7) at σ0=0.05 MPa and σ0=0.2 MPa, scaling δC linearly with σ0 so that χ takes the same values as in Fig. 7 (e.g., 10, 100, 1000, 8000), using the same meshes and the same c-family snapshot construction; plot σC/σ0 vs χ for the three σ0 values. If the three curves collapse to within about 2%, Eq. (10)'s dimensional reduction is confirmed in the defect regime; if they deviate by more than a few percent, the defect design maps and the LEFM comparison are conditional on the individual values of σ0 and δC, and the claimed universality of Eqs. (2) and (10) must be restricted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the universal functional form σC/σ0 = f(χ, β, ξ, R/Df, 2c/Df) in Eqs. (2) and (10), with χ = σ0Df(1−ν²)/(2πEδC). This is what turns the FE runs into design maps and motivates the array-statistics application via ref. [31]. It rests entirely on the Tang et al. [22] dimensional reduction: σ0 and δC must enter only through χ. The only direct checks are in the 'Results for fibril detachment' section: for β=1.05, ξ=0.025 and for β=1.25, ξ=0.01, the c=0 runs (cohesive traction over the whole interface) were repeated at σ0 = 0.05, 0.1, 0.2 MPa and gave identical normalized pull-off curves. Those runs determine the defect-free pull-off points (colored triangles in Fig. 2), so the defect-free reduction is reasonably supported. They do not test the c>0 snapshot family — the states with cohesive traction applied only on r ≥ c — yet all central-defect strengths in Figs. 7–9 are read from such states, because the initial defect of radius c0 makes the pull-off state a c>0 state rather than a point on the c=0 envelope. The LEFM agreement at high χ (Eq. (8) vs. the solid curves in Fig. 7) is indirect evidence only in the small-scale-yielding tail; it says nothing about the intermediate-χ range where the defect-sensitivity transitions occur. If the reduction fails in the defect regime — e.g., through finite-strain or contact nonlinearity that is more pronounced once a traction-free crack is present — then Eq. (10), the defect design curves, and the LEFM comparison lose their claimed generality, and the central message narrows to 'runs at one material point are internally consistent' rather than universal design rules.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses finite-element simulations with a Dugdale cohesive zone model to compute pull-off curves and strengths for mushroom-shaped fibrils adhered to a rigid substrate. The central claim is that, for defect-free fibrils, the normalized pull-off strength σC/σ0 depends only on the dimensionless parameter χ and the geometric ratios β = Df/D, ξ = h/Df, and R/Df (Eqs. 2–3), and that for fibrils with a central adhesion defect it additionally depends on 2c/Df (Eq. 10). The authors present full load-extension curves for representative geometries, a systematic map of pull-off strength versus χ (Fig. 3), traction and separation profiles at pull-off (Figs. 4–6), and defect-sensitivity curves with a parameter-free LEFM prediction (Eq. 8) that matches the numerics at high χ (Figs. 7–9). The results indicate that wide, thin caps promote central detachment and maintain high strength up to very large χ, but that central defects cause severe degradation at high χ. The paper concludes with design guidance, including an application to experimental PDMS fibrils and to statistical array models.","tokens_in":23964,"tokens_out":3849,"duration_ms":42247,"significance":"If the central claims hold, the paper provides practically useful design maps for bio-inspired fibrillar adhesives and a direct input for statistical strength models of fibril arrays, since Eq. (8) gives a closed-form, parameter-free LEFM prediction for the defect-dominated regime. The work's strengths include a detailed description of the finite-element protocol, explicit statements of the modeling assumptions, and direct checks of the pull-off strength's independence from σ0 and δC separately (performed for two defect-free geometries). The comparison with the LEFM prediction is a concrete, falsifiable test that the numerics pass in the small-scale-yielding tail. However, the claimed universality of the defect-case functional form (Eq. 10) is not directly verified, and the procedure for assembling unstable trajectories from separate equilibrium solutions is an assumption that should be validated. These issues are local and fixable, which is why I recommend major revision rather than rejection.","major_comments":[{"comment":"The functional form (10) extends the Tang et al. dimensional reduction to fibrils with central defects, but the invariance with respect to σ0 and δC beyond the combination χ is tested only for the defect-free (c = 0) envelopes of two geometries in the section 'Results for fibril detachment'. All c > 0 states that determine the defect strength curves in Figs. 7–9 are computed at a single value σ0 = 0.1 MPa. If the reduction fails once a traction-free crack is present, Eq. (10) and the defect design curves lose their claimed generality. I request a repeat of at least one defect case (for example β = 1.25, ξ = 0.01, 2c/Df = 0.005) with σ0 = 0.05 MPa and σ0 = 0.2 MPa, with the normalized pull-off curves compared.","section":"§ Pull-off strength of fibrils with an initial adhesion defect, Eq. (10)"},{"comment":"The full pull-off curves are assembled from equilibrium states obtained on separate meshes with different pre-crack sizes c, unified by having equal δC. The manuscript states that this is 'a convenient one for obtaining accurate and reliable results for the unstable pull-off process' but does not validate that the quasi-static assembly reproduces the actual unstable detachment trajectory. Since the pull-off strength is read as the maximum of this assembled curve, this assumption is load-bearing. I recommend verifying the procedure for at least one representative case by a direct transient or path-following calculation, or by demonstrating that the ascending and descending branches form a single equilibrium manifold with a limit point at the claimed pull-off state.","section":"§ Detachment process and Fig. 2"},{"comment":"The mesh resolution criterion in Appendix A uses the strip-yield estimate (A3) and requires 'at least 5 finite elements' along the Dugdale zone, but this rule is heuristic, and the reported refinement study is localized to 'a representative case'. All quantitative claims, including the extended plateau in Fig. 3 and the LEFM comparisons in Figs. 7–9, depend on the accuracy of the cohesive-zone resolution. I ask for a quantitative convergence test for at least one edge-detachment case and one central-detachment case at high χ (e.g., χ ≈ 103), comparing pull-off strength between the present mesh and a mesh with half the interface element size.","section":"Appendix A, mesh resolution and convergence"}],"minor_comments":[{"comment":"The caption for Fig. 2 does not explain the meaning of the black triangle, the colored triangles, or the lines of circles and dashes; a concise legend or a sentence identifying these elements would improve readability.","section":"§ Results for fibril detachment, Fig. 2 caption"},{"comment":"The text contains a typographical error: 'Following Khaderi et al al [37]' should read 'Following Khaderi et al. [37]'.","section":"§ Discussion, paragraph on Khaderi et al."},{"comment":"The abstract contains the LaTeX artifact '\\c{hi}'; this should be typeset as the Greek letter χ.","section":"§ Abstract"},{"comment":"The reference to Aksak et al. as [21] for the optimal mushroom shape and as [24] for the fitting formula should be checked; the text appears to attribute the fit to a different study than the one cited, and the numbering should be made consistent.","section":"§ Introduction and references"},{"comment":"The phrase 'at the bottom of the line of circles and dashes having the same color' is unclear; specifying that the color-coded circles and dashes correspond to the post-peak equilibrium branch for a given χ would help the reader.","section":"§ Results for fibril detachment, text near Fig. 2(a)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a mostly sound numerical study, but it relies heavily on prior work by the same group (Tang et al., Balijepalli et al., Fleck et al., Khaderi et al.) without shipping code or data, and the mesh convergence is presented in only a localized form. For a journal that publishes computational mechanics, making the simulation data available and adding the requested verification runs would substantially increase confidence. The lack of a σ0-variation check for the defect cases is the main correctness risk; it is fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid, careful computational paper that delivers what it promises—complete pull-off curves for mushroom fibrils, a clean type 1/type 2 classification, and a closed-form defect-sensitivity relation. The qualitative superiority of wide, thin caps was already known, but the full unstable branches and the central-defect LEFM scaling are genuinely new.\n\nWhat it does well: the FE protocol is described in unusual detail, including element counts, interface mesh size, and a Dugdale-zone resolution criterion in Appendix A. The authors check sigma0-independence for two defect-free geometries and explicitly decline to fit strength curves, which is the right call. The comparison against Del Campo's experiment is honest: they predict a factor of 6, report the experimental factor of 20, and explain the discrepancy through conformal contact on a hemispherical glass surface. That is the kind of reporting I want to see.\n\nThe soft spots are mostly minor. First, Appendix Eq. (A2) appears to miss the beta-squared factor that appears in Eq. (8); as written it is inconsistent with the main-text LEFM derivation. Either the notation sigma_I is meant to be the stalk stress, or there is a typo. The plotted curves use Eq. (8), so it does not affect results, but it needs fixing.\n\nSecond, the stress-test concern is legitimate: the sigma0-independence check is only done for c=0 envelopes, while all central-defect strengths in Figs. 7–9 come from c>0 states. The dimensional reduction to chi almost certainly carries over—the problem is linear elastic with a step cohesive law, so normalized strength should depend only on chi and geometry, not on sigma0 and delta_C separately—but the authors never demonstrate it for the defect cases. A skeptic could ask for one extra set of runs at a second sigma0 with a central defect. This is a request for an additional check, not evidence of a wrong result.\n\nThird, no code or data are shipped, and global mesh convergence is not performed. The localized refinement study is reasonable, but for a paper whose value is in design curves, reproducibility would be better served by archiving the input files.\n\nWho is this for: anyone doing computational design of fibrillar adhesives, and experimentalists who want defect-sensitivity maps. It is not a breakthrough, but it is a useful and honest piece of engineering mechanics. With the appendix typo fixed and the defect-case dimensional check added, I would be happy to see it in print.\n\nRecommendation: send to peer review. The central claims hold up; the issues are addressable and none are load-bearing.","headline":"A careful and useful extension of mushroom-fibril cohesive-zone modeling, with a closed-form defect scaling that deserves referee time despite a few addressable gaps.","tokens_in":24541,"tokens_out":7154,"would_cite":true,"duration_ms":72163,"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 mushroom-shaped fibril's pull-off strength is set by one adhesion parameter and cap geometry; a wide, thin cap keeps it near the theoretical maximum over a huge size range, until a central defect takes over.","keywords":["mushroom-shaped fibrils","bio-inspired adhesion","pull-off strength","Dugdale cohesive zone","dimensionless adhesion parameter chi","adhesion defects","finite element simulation","design maps"],"falsifier":"Measure the pull-off strength of a $\\beta = 1.25$, $\\xi = 0.01$ PDMS fibril with tip diameter 6.52 mm and a central defect 3.26 $\\mu$m in diameter: the paper predicts a drop from 0.78 MPa to 0.27 MPa, and a substantial deviation would overturn the LEFM-based defect predictions.","tokens_in":23373,"feed_emoji":"🍄","tokens_out":10219,"duration_ms":97087,"temperature":0.7,"pith_summary":"The paper uses finite element simulations with a Dugdale cohesive zone to work out when and how mushroom-shaped adhesive fibrils detach from a rigid substrate. It claims that, for defect-free fibrils, the pull-off strength depends only on a single dimensionless adhesion parameter $\\chi$ together with the cap geometry, and that detachment is always unstable under load control once the interface separation reaches a critical value, whether that happens at the edge or at the center. A wide, thin cap (flange-to-stalk ratio $\\beta = 1.25$, thickness ratio $\\xi = 0.01$) keeps the pull-off strength at or above $0.78\\,\\sigma_0$ for $\\chi$ up to about 8000, a regime no cylindrical punch reaches. The paper also claims that a central adhesion defect can severely reduce that strength, especially at large $\\chi$, and that this reduction follows linear elastic fracture mechanics predictions once $\\chi$ is high.","feed_headline":"Mushroom fibrils keep 80% of ideal adhesion over a vast size range","feed_subtitle":"Wide, thin caps delay edge detachment, but a tiny central defect can cut that strength by two-thirds.","key_machinery":"The load-bearing machinery is the Dugdale cohesive zone model, in which the interface carries the constant stress $\\sigma_0$ until its separation reaches the critical value $\\delta_C$, combined with the dimensionless parameter $\\chi = \\sigma_0^2 D_f (1-\\nu^2)/(2\\pi E W_{\\mathrm{adh}})$ that measures the fibril diameter against the cohesive zone length. To capture unstable pull-off, the authors assemble each load-extension curve from hundreds of separate finite element solutions, each with a fixed partially detached region (a central circle or an edge annulus of radius $c$) and with $\\delta_C$ evaluated at the edge of that region. For central defects, a small-crack LEFM formula supplies the asymptotic strength at high $\\chi$ and provides the comparison that validates the numerical results.","core_discovery":"The central discovery is that the detachment behavior of these fibrils collapses onto a family of design curves: the normalized pull-off strength $\\sigma_C/\\sigma_0$ is a function only of $\\chi$, $h/D_f$, $D_f/D$, $R/D_f$ and, for a central defect, $2c/D_f$. In every geometry studied, pull-off occurs exactly when the maximum interface separation first reaches $\\delta_C$, and the ensuing separation is unstable under load control regardless of whether it initiates at the edge or at the center. Mushroom geometries with wide, thin caps suppress the edge stress concentration and move the initiation site to the center, which keeps the strength high over an enormous range of $\\chi$; however, edge initiation remains possible in every geometry under some conditions, while central initiation is not. When a central adhesion defect is present, it acts as a penny-shaped crack, and at high $\\chi$ the numerical strength curves converge to the paper's closed-form linear-elastic-fracture-mechanics expression (Eq. 8), validating the crack analogy.","pith_inferences":["Editorial inference: the LEFM asymptote implies that at large $\\chi$, for a fixed relative defect size $2c/D_f$, the pull-off strength falls roughly as $\\chi^{-1/2}$, so increasing the allowed fibril size for a fixed absolute defect costs more than a linear penalty in strength.","Editorial inference: the strength-versus-defect-size curves in Figures 7-9 are exactly the input a statistical array model needs, so combining them with a measured defect-size distribution would predict array-level adhesion and its scatter without new mechanics.","Editorial inference: the no-slip, rigid-substrate idealization likely makes these predictions upper bounds; a compliant substrate or sliding interface would relieve edge stress concentrations and could shift the edge-versus-center competition.","Editorial inference: because the paper verifies the single-parameter $\\chi$ reduction for only two defect-free geometries, extending the central-defect curves to other $\\beta$ and $\\xi$ would first require checking the same reduction in the defected regime."],"forward_implications":["If the central claim is right, the strength-versus-$\\chi$ curves constitute a design map: for a given material, the fibril tip diameter can be chosen so that $\\chi$ stays below $\\chi_{c1}$ (flaw-insensitive, strength equal to $\\sigma_0$) or below $\\chi_{c2}$ (central detachment, strength at least $0.78\\,\\sigma_0$).","The wide-thin-cap geometry $\\beta = 1.25$, $\\xi = 0.01$ permits tip diameters orders of magnitude larger than a punch for the same retained strength: about 6.5 mm for the PDMS example and about 1.7 mm for keratin.","Because detachment is unstable under load control, arrays of these fibrils should show abrupt, defect-triggered failure events rather than gradual softening.","Central defects of even a few micrometers can cut the pull-off strength of an otherwise optimal large fibril by 65% or more at high $\\chi$, so defect size, not just geometry, must enter sizing rules."],"supporting_citations":[{"why":"Supplies the dimensional-reduction argument that pull-off strength depends only on $\\chi$ and geometry, which the paper verifies and extends.","marker":"[22]"},{"why":"Defines the mushroom geometry with fillet radius and the thin-large-cap design whose central-detachment behavior is analyzed here.","marker":"[28]"},{"why":"Provides the straight-punch pull-off curve used as the baseline in Figure 3.","marker":"[26]"},{"why":"Supplies the penny-shaped interface crack stress-intensity factor used to derive the LEFM prediction in Eq. (8).","marker":"[33]"},{"why":"Provides punch-fibril detachment results with perimeter defects used to compare defect sensitivity with the mushroom.","marker":"[37]"},{"why":"Provides experimental pull-off measurements on mushroom and punch fibrils used to benchmark the model.","marker":"[14]"},{"why":"Provides the statistical array-strength framework that motivates needing strength as a function of defect size.","marker":"[31]"},{"why":"Supplies the Dugdale strip-yield model used to set mesh resolution and to estimate cohesive-zone size in Appendix A.","marker":"[27]"},{"why":"Provides experimental evidence of local detachment at the center of mushroom fibrils, justifying the central-defect focus.","marker":"[30]"}],"fun_headline_variants":["Mushroom fibrils obey design curves, never stable under load control","Wide thin caps shift mushroom fibril detachment to the center","Central defects slash mushroom fibril pull-off strength at high χ","Mushroom fibril pull-off collapses onto universal design curves","Edge detachment always possible, but wide caps delay it"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that varying $\\sigma_0$ and $\\delta_C$ separately is unnecessary and that only the combined parameter $\\chi$ plus shape ratios control pull-off strength; this was checked for only two defect-free geometries and never for the central-defect calculations, so if that dimensional reduction fails in the defective regime the design curves lose generality.","fun_headline_variants_meta":{"raw":{"variants":["Mushroom fibrils obey design curves, never stable under load control","Wide thin caps shift mushroom fibril detachment to the center","Central defects slash mushroom fibril pull-off strength at high χ","Mushroom fibril pull-off collapses onto universal design curves","Edge detachment always possible, but wide caps delay it"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000816,"raw_usage":{"total_tokens":3581,"prompt_tokens":956,"completion_tokens":2625,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":2539}},"tokens_in":572,"tokens_out":2625,"duration_ms":21041,"temperature":1.0,"reasoning_tokens":2539,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:44:00.206886+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the pull-off strength of a $\\beta = 1.25$, $\\xi = 0.01$ PDMS fibril with tip diameter 6.52 mm and a central defect 3.26 $\\mu$m in diameter: the paper predicts a drop from 0.78 MPa to 0.27 MPa, and a substantial deviation would overturn the LEFM-based defect predictions.","supporting_citations":[{"cited_title":"Hui, C","cited_arxiv_id":null,"evidence_quote":"Supplies the dimensional-reduction argument that pull-off strength depends only on $\\chi$ and geometry, which the paper verifies and extends."},{"cited_title":"G.; Begley, M","cited_arxiv_id":null,"evidence_quote":"Defines the mushroom geometry with fillet radius and the thin-large-cap design whose central-detachment behavior is analyzed here."},{"cited_title":"A.; Khaderi, S","cited_arxiv_id":null,"evidence_quote":"Provides the straight-punch pull-off curve used as the baseline in Figure 3."},{"cited_title":"K.; Bregman, A.M","cited_arxiv_id":null,"evidence_quote":"Supplies the penny-shaped interface crack stress-intensity factor used to derive the LEFM prediction in Eq. (8)."},{"cited_title":", McMeeking, R.M","cited_arxiv_id":null,"evidence_quote":"Provides punch-fibril detachment results with perimeter defects used to compare defect sensitivity with the mushroom."},{"cited_title":"Greiner, C.; Arzt, E","cited_arxiv_id":null,"evidence_quote":"Provides experimental pull-off measurements on mushroom and punch fibrils used to benchmark the model."},{"cited_title":"McMeeking R.M.; Foster K.; Tinnemann V.; Hensel R.; Arzt E","cited_arxiv_id":null,"evidence_quote":"Provides the statistical array-strength framework that motivates needing strength as a function of defect size."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Dugdale strip-yield model used to set mesh resolution and to estimate cohesive-zone size in Appendix A."},{"cited_title":"Hernández L.; Fischer S.C.L.; Arzt E.; Bennewitz R.; Hensel R","cited_arxiv_id":null,"evidence_quote":"Provides experimental evidence of local detachment at the center of mushroom fibrils, justifying the central-defect focus."}],"review_version":1}