{"id":"93a1a4c2-af46-4831-9e0c-89e8356b03eb","arxiv_id":"2507.10841","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Hydrated collagen fibrils dissipate energy through rate-independent hysteresis with return point memory, captured by a generic model called EPICAL that predicts force and dissipated energy for arbitrary indentation paths.","lead":"Atomic force microscope measurements on hydrated collagen fibrils show that the energy lost in each push-pull cycle does not depend on how fast the tip moves, below 1 micrometer per second. A simple hysteresis model built from one large indentation cycle then predicts the force and energy loss for other indentation paths, pointing to a rate-independent memory mechanism instead of viscosity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The arbitrary-trajectory claim rests on one amplitude-independent unloading stiffness κT (eqn 1b/5); the only out-of-sample test is a single complex trajectory with acknowledged deviations, so κT's depth/amplitude dependence must be checked before the prediction is taken as general.","rationale":"The core observation—rate-independent hysteresis with return point memory over 0.01–1 µm/s—is supported by the collapse of Edis vs FM and by the complex-trajectory data. The proposed operator is a legitimate rate-independent hysteresis model, and the closed-form Edis expression matching the data is a real success. I searched for a more fundamental objection (e.g., time-dependent surface leveling corrupting rate independence, or the model being circular because fT is fit to one retraction and then used to predict another) and did not find one. The load-bearing weak point is the extrapolation from one calibrated κT to all possible return points. The paper itself flags the two places where this extrapolation fails. A single additional experiment with varied inner-loop depths, or a reanalysis of the existing cyclic data with depth-dependent κT, would settle whether the 'arbitrary trajectories' claim is quantitatively justified or needs to be a conditional claim. Since the reader already conditionalized on exactly this point and on data availability, my verdict is unchanged; the paper should not be rejected, but should not be accepted as demonstrating arbitrary-trajectory prediction without the κT(FM) check.","tokens_in":13289,"tokens_out":6272,"duration_ms":79730,"concrete_test":"Using the cyclic data behind Fig. 2b, measure the unloading slope κT from the first ~20 nm of each retraction curve as a function of FM (or of return-point depth) and recompute eqn (6) with κT(FM) in place of the single calibrated value. If the predicted Edis(FM) moves by more than the stated ~20% scatter for FM>5 nN, the constant-κT assumption is load-bearing and the arbitrary-trajectory claim must be qualified. Independently, run one additional complex trajectory with inner retract-approach loops whose reversal depths differ from the calibration cycle, and compare the measured force during each inner-loop retraction with fT(δ,δ_i*) using the calibrated κT; systematic deviations at the new depths would show fT must depend on δ_i*.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.8 constructs Γ̂[δ]=max{fR, min{fA, fT}} with fT(δ,δ_i*) = f(δ_i*) − κT(δ−δ_i*), and eqn (6) derives Edis(FM) from the same linear unloading. The least secure premise is that a single κT, calibrated by eqn (5) from the retraction slope of one large indentation, holds for every return point and every indentation amplitude. The paper's own evidence is limited to 'approximately linear' intermediate excursions with stiffness 'similar to kT' (Section 2.6), and Section 2.9 admits visible deviations on the first small loop and the final retraction (Fig. 3c,d). If κT varies with depth or with the position of the return point, both the force prediction for complex trajectories and the closed-form Edis prediction inherit a systematic error; Fig. 2c even shows a slight increase of kT for FM>5 nN. Since the out-of-sample demonstration is one complex trajectory at one velocity on one fibril, the 'arbitrary indentation trajectories' claim is under-determined. This is an empirical generalization that needs a direct test, not an internal inconsistency. No part of the argument requires rejecting the core rate-independence observation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports AFM nanoindentation measurements on hydrated collagen fibrils showing that the force-distance hysteresis and dissipated energy are independent of tip velocity in the range 0.01–1 um/s, and that the response exhibits return point memory. The authors construct a phenomenological hysteresis operator (eqns 1–3) whose parameters are extracted from one measured approach-retract cycle, and show that it predicts the dissipated energy as a function of maximum force (eqn 6) and reproduces the main features of a complex indentation trajectory with multiple intermediate retract-approach cycles (Fig. 3). The model combines elastoplastic indentation, a linear elastic unloading branch with constant stiffness kT, and an attractive capillary force floor, with an additional elementary hysteresis in the capillary region.","tokens_in":13580,"tokens_out":14060,"duration_ms":144871,"significance":"If the rate-independence result is robust, it is a significant finding for AFM-based nanomechanics of collagen and soft materials, because low-speed hysteresis is usually attributed to viscoelasticity. The closed-form Edis(FM) relation and the out-of-sample prediction of a complex trajectory are clear strengths, and the model is parameterized from a single cycle, which aids reproducibility. The main weaknesses are the limited validation of the 'arbitrary trajectories' claim and the reliance on an amplitude-independent kT. The paper is otherwise careful in distinguishing rate-independent from rate-dependent regimes and provides direct stress-relaxation and creep controls.","major_comments":[{"comment":"The claim that the hysteresis operator predicts the force and dissipated energy for 'arbitrary indentation trajectories' rests on a single complex trajectory measured at one velocity (vtip = 0.04 um/s) on one isolated fibril. The authors themselves note deviations on the first small retract-approach loop and the last retraction (Section 2.9). Please add a quantitative error analysis (e.g., RMS force error, relative Edis error) for the Fig. 3 prediction and, ideally, test additional trajectories with different return-point depths and loop amplitudes. If such data are not available, the wording in the abstract and Section 2.8 should be tempered to 'the tested trajectory class.'","section":"§2.9, Fig. 3"},{"comment":"The model's unloading branch fT uses a single constant kT determined from the retraction slope of one large indentation, yet Fig. 2c shows that kT increases for FM > 5 nN and Section 2.6 reports only that intermediate excursions have stiffness 'similar to kT.' Because this parameter controls the entire unloading behavior, a depth- or amplitude-dependent kT would propagate directly into the predicted force and Edis for arbitrary trajectories. Please report kT for return points at different indentation depths and amplitudes (e.g., the intermediate loops of Fig. 3) and quantify how the observed variation affects the predictions, or state the range of return-point amplitudes over which kT is effectively constant.","section":"§2.4, eqn (5)"},{"comment":"The quantitative model specification is incomplete for the attractive region 0 < delta < delta_off, where the text invokes an 'elementary hysteresis model' (ref. 43) without giving its equations. Since the full-trajectory prediction in Fig. 3 includes the attractive parts of the force signal, the model as written is not fully reproducible from the equations provided. Please include the explicit capillary-bridge operator or clearly restrict the quantitative predictions to delta <= 0.","section":"§2.8, eqns (1)-(3)"},{"comment":"The closed-form Edis(fM) derivation assumes that the linear unloading from the maximum-force point reaches the attractive limit fc before the tip returns to delta_s (i.e., delta_c = delta_min + (fM - fc)/kT <= delta_s). Please state this condition explicitly and confirm that all FM data points in Fig. 2b satisfy it; for large fM or small kT the loop-closure geometry changes and eqn (6) would need modification.","section":"§2.10, eqn (6)"}],"minor_comments":[{"comment":"The heading 'Results and disscusion' contains a typo; it should read 'Results and discussion.'","section":"Section 2 heading"},{"comment":"The caption does not mention error bars or the scatter reported as 'typically below 20%.' Please add representative error bars to the main figure or explicitly refer to Fig. S3 in the caption.","section":"Fig. 2b"},{"comment":"State the number of measurement positions and cycles per velocity used for the Edis statistics in the main text, rather than only in Fig. S3.","section":"Section 2.3"},{"comment":"In the compiled manuscript the formula appears with split exponents; please ensure the typeset equation is unambiguous, for example (f_M - f_c)^(1+1/gamma) / ((1+gamma) kappa_A^(1/gamma)).","section":"Eqn (6)"},{"comment":"The ESI footnote uses a placeholder DOI (10.1039/cXsm00000x/); update it to the actual supplementary information DOI.","section":"ESI footnote"},{"comment":"The return-point-memory discussion would benefit from an explicit demonstration of wiping out (a larger loop erasing the memory of a smaller one), since that property is central to the operator's definition.","section":"Section 2.6"}],"recommendation":"major_revision","confidential_remarks":"This manuscript is a previously published Soft Matter article (2024, 20, 2831–2839); the review concerns the arXiv posting. The rate-independence observation is valuable, but the 'arbitrary trajectories' claim is broader than the evidence presented. The requested additions—systematic kT characterization and additional trajectory validation—appear feasible within the authors' existing dataset. I see no grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe bottom line: this paper reports a genuinely new observation — rate-independent hysteresis with return point memory in AFM nanoindentation of hydrated collagen fibrils — and supports it with clean data. I'd send it to a serious referee.\n\nWhat's new: the claim that the shape of force-distance curves is independent of tip velocity between 0.01 and 1 µm/s, and that dissipated energy depends only on the maximum force, not on rate. The data in Fig 2 show this convincingly: overlapping curves at three speeds and a collapse of Edis vs FM. The return point memory demonstration in Fig 3 is also convincing: multiple small retract-approach loops return to the original trajectory, which is exactly what a generalized stop operator predicts.\n\nWhat the paper does well: the model is refreshingly simple. They take one large approach-retract cycle, fit a few parameters (approach curve, contact stiffness, adhesion), and then predict dissipated energy for a whole family of cycles with varying maximum force. That prediction matches the measured data across the full range. That's a genuine out-of-sample test, not a fit to the same data. They also honestly note where the model deviates: the first small loop and the final retraction.\n\nSoft spots: the model's 'arbitrary trajectories' claim is only demonstrated on one complex trajectory, with acknowledged deviations. The likely source is the assumption that the unloading stiffness κT is a single constant measured from one large cycle. Fig 2c shows κT drifts upward for FM > 5 nN, so a depth-dependent stiffness could improve the prediction. That's a scope issue rather than a fatal flaw. Second, the data availability is poor: no raw FD curves, no code for the model, and Fig 2b has no error bars (the text gives scatter values). For a result that challenges the default viscoelastic interpretation in a broad literature, the community will want to re-fit the model on their own data.\n\nCitation pattern is fine; self-citations are to their own methods papers, and they properly credit the hysteresis operator literature (Prandtl, Visintin, Brokate and Sprekels). I don't see circular fitting.\n\nWho should read this: anyone doing AFM nanoindentation of soft hydrated materials. It's a solid candidate for peer review, with the request to supply data and temper the generality claim.","headline":"The core rate-independence observation is solid and the model does real out-of-sample work; the 'arbitrary trajectories' claim is broader than the single tested trajectory supports.","tokens_in":14110,"tokens_out":3439,"would_cite":true,"duration_ms":35859,"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":"Hydrated collagen fibrils show rate-independent hysteresis with return point memory, and a three-function hysteresis operator fitted to one indentation cycle predicts the force and dissipated energy for arbitrary tip trajectories.","keywords":["rate-independent hysteresis","return point memory","collagen fibrils","atomic force microscopy","nanoindentation","energy dissipation","hysteresis operator","elastoplastic indentation"],"falsifier":"Measure cyclic force-distance curves with small retract-approach loops starting from several different return-point depths on the same fibril and compare the measured unloading slopes to the κT value from one deep cycle; if the initial slope of a small loop at shallow depth is systematically different, the straight-line unloading assumption fails. Alternatively, ramp the tip velocity above 1 µm/s, where inner hysteresis loops appear due to viscous phase lag, to test the stated rate-independent regime boundary.","tokens_in":13091,"feed_emoji":"🧬","tokens_out":5252,"duration_ms":60228,"temperature":0.7,"pith_summary":"The paper reports that when an atomic force microscope tip indents a hydrated collagen fibril slowly, at tip speeds between 0.01 and 1 µm/s, the measured force-distance curves do not change shape with speed: the dissipated energy per cycle is rate-independent and the system shows return point memory. It proposes a generic hysteresis model in which the entire force response is determined by three time-independent functions—the approach curve, a straight unloading branch with a fixed contact stiffness, and a constant attractive force limit—combined through the operator Γ̂[δ] = max{fR, min{fA, fT}}. The claim is that this operator, calibrated on one large approach-retract cycle, predicts force and dissipated energy for arbitrary indentation trajectories, unifying elastoplastic indentation, capillary adhesion, and surface leveling at these slow rates. A sympathetic reader would care because it identifies a previously unrecognized dissipation mechanism in a load-bearing connective-tissue component and provides a parameter-light predictive tool for nanoindentation.","feed_headline":"Nanoindentation of collagen: no rate dependence below 1 µm/s","feed_subtitle":"Force and dissipated energy under any indentation path follow from one measured approach-retract cycle.","key_machinery":"The central object is the hysteresis operator Γ̂ defined in eqn (3): for indentation δ and last return point δ*i, the output force is max{fR(δ), min{fA(δ), fT(δ, δ*i)}}. Here fA is a power-law approach curve of contact-model form, fT is a straight elastic unloading line with slope κT, and fR is the constant attractive (capillary and adhesion) force. The operator is rate-independent by construction because it depends on the sign of δ − δ*i rather than on time, and it implements return point memory by storing only the last return point. The same form is a generalized stop operator, with fA and fR serving as the upper and lower output bounds; the predicted forces, dissipated energies, and complex-trajectory behaviors all follow from plugging calibrated parameters from one measured cycle into this operator.","core_discovery":"On its own terms, the central discovery is that hydrated type I collagen fibrils, when indented perpendicular to their axis at velocities at or below 1 µm/s, dissipate energy through a rate-independent hysteretic process with return point memory rather than through viscous friction. The shape of force-distance curves is independent of tip velocity in this range, and the dissipated energy Edis depends only on the maximal force setpoint FM, with an analytical expression provided by the model. Return point memory manifests in trajectories with intermediate retract-approach cycles: at the end of each small cycle the force returns to the value at the cycle's start, so the system state depends only on the last return point. The model's operator reproduces the measured force and dissipated energy for complex trajectories, with deviations only at the first small internal loop and the final retraction.","pith_inferences":["If κT turns out to vary with indentation depth or with the position of the return point, a natural extension would be to replace the straight unloading line fT with a curved family parameterized by the return point; the max/min operator form could survive but would then encode memory beyond the last return point.","The operator's resemblance to a generalized stop operator suggests that superposition-based hysteresis descriptions in the Preisach family could be adapted to nanoindentation data if nested-loop measurements at multiple amplitudes were available.","A testable extension is to fit the EPICAL operator on one fibril and then check whether it transfers to other positions along the same fibril or to other hydrated biopolymer fibrils; if the mechanism is generic, the calibrated parameters should predict new force curves without refitting."],"forward_implications":["Given one calibrated large indentation cycle, the operator predicts the force output for any other indentation trajectory without fitting time-dependent viscoelastic parameters.","The analytical expression for Edis(FM) gives a velocity-free prediction of how much energy each indentation cycle dissipates as a function of indentation amplitude.","The model unifies three previously distinct phenomena—elastoplastic indentation, capillary adhesion and bridge formation, and surface leveling—into one rate-independent description.","Because native hydrated collagen fibrils are in the glassy state, the authors anticipate that rate-independent hysteretic energy dissipation will also be found in other soft condensed matter and biological materials at slow deformation rates."],"supporting_citations":[{"why":"Defines return point memory and memory wiping, the memory mechanism the model exploits.","marker":"29,30"},{"why":"Prandtl's elasto-plasticity model, the conceptual ancestor of the hysteresis operator.","marker":"34"},{"why":"Cyclic nanoindentation of Fe-based amorphous alloys, which motivates the elastoplastic interpretation of collagen fibril indentation.","marker":"35"},{"why":"Oliver-Pharr contact stiffness method, used to extract κT from the retraction slope.","marker":"40"},{"why":"Contact-model functional form used for the approach curve fA.","marker":"41"},{"why":"Procedure for locating d0, the onset of attractive forces, needed for calibration.","marker":"42"},{"why":"Elementary hysteresis model for capillary bridge formation and collapse in the region 0 < δ < δoff.","marker":"43"},{"why":"Generalized stop operator, to which the hysteresis operator Γ̂ is compared.","marker":"44,45"}],"fun_headline_variants":["One indentation cycle predicts collagen's rate-independent energy loss","Collagen's energy loss: no speed dependence, only memory","Hysteresis without rate: collagen's return point memory","Collagen at low speeds: rate-independent hysteresis with memory"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes that unloading from any return point follows a straight line with a single, depth-independent contact stiffness κT taken from one large indentation cycle; if κT actually varies with indentation depth or with the position of the return point, the predicted force and dissipated energy for complex trajectories will drift away from the data.","fun_headline_variants_meta":{"raw":{"variants":["One indentation cycle predicts collagen's rate-independent energy loss","Collagen's energy loss: no speed dependence, only memory","Hysteresis without rate: collagen's return point memory","Collagen at low speeds: rate-independent hysteresis with memory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002003,"raw_usage":{"total_tokens":7745,"prompt_tokens":808,"completion_tokens":6937,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":424,"completion_tokens_details":{"reasoning_tokens":6869}},"tokens_in":424,"tokens_out":6937,"duration_ms":51922,"temperature":1.0,"reasoning_tokens":6869,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:24:12.568255+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure cyclic force-distance curves with small retract-approach loops starting from several different return-point depths on the same fibril and compare the measured unloading slopes to the κT value from one deep cycle; if the initial slope of a small loop at shallow depth is systematically different, the straight-line unloading assumption fails. Alternatively, ramp the tip velocity above 1 µm/s, where inner hysteresis loops appear due to viscous phase lag, to test the stated rate-independent regime boundary.","supporting_citations":[],"review_version":1}