{"id":"da66556f-399f-4554-8031-2987716bd098","arxiv_id":"2601.22743","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A two-element fractional Kelvin-Voigt model calibrated to quasi-static rheometer and high-frequency MRE data describes porcine brain viscoelasticity from 0 to 2100 Hz, with corona radiata stiffest and viscosity dominating at high frequencies.","lead":"This study measured how three regions of pig brain respond to slow compression, tension, and shear as well as to fast vibrations up to 2100 Hz, and combined both datasets in one fractional viscoelastic model. The unified model gives brain simulators a single set of material parameters across a wide range of loading speeds.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed 0–2100 Hz unification rests on a single 0 Hz anchor extrapolated from nonlinear rheometer fits; no data constrain 0.053–300 Hz, so wideband validity is an interpolation across an unmeasured gap.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing point: the 0 Hz anchor is a model-extrapolated bridge between two different experiments, and the gap from 0.053 to 300 Hz is unmeasured. I agree with that assessment. My stress-test adds two specifics: (1) μ_e is structurally pinned by the anchor because G*(0)=μ_e, so the low-frequency agreement is partly by construction; (2) the low-frequency loss modulus is not constrained at all, since only G′(0) enters the calibration and G″ from the rheometer at 0.027/0.053 Hz is not used. The paper is otherwise careful: it discloses exclusions, reports individual fits, uses nonparametric statistics, and selects model forms with clear physical motivation. No fraud or p-hacking concerns. The conditional verdict is appropriate because a targeted intermediate-frequency measurement (or an MRE-only refit) could confirm or refute the interpolation. I recommend no change to the reader's verdict.","tokens_in":17510,"tokens_out":7930,"duration_ms":84445,"concrete_test":"Perform small-strain oscillatory shear rheometry (or DMA) on matched corona radiata, putamen, and thalamus samples at intermediate frequencies spanning 0.1–100 Hz, i.e., inside the unmeasured gap, and compare the measured G′(f) and G″(f) with the calibrated fractional Kelvin-Voigt model from Table 4. If the model systematically deviates—especially in G″ at low frequencies, which was not included in the calibration—the unified description fails because the 0 Hz anchor and MRE data alone do not determine the intervening behavior. As an auxiliary sensitivity check, refit the model using only MRE data (300–2100 Hz) and test whether its extrapolation to 0.027–0.053 Hz reproduces the rheometer-derived moduli; if not, the apparent unification is anchor-driven rather than a property of the model form.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Abstract, Sec. 5) is that a 5-parameter fractional Kelvin-Voigt model 'unifies' quasi-static rheometer and 300–2100 Hz MRE data into one comprehensive description. The only low-frequency constraint entering the calibration is G′(0), computed in Table 3 from Prony parameters identified by inverse-fitting a nonlinear Ogden-Prony model to large-strain compression/tension/shear data (Sec. 2.2.2). The paper itself states (Sec. 2.3.2) that the measured MRE response contains no information about the lower frequency regime; it then appends exactly one 0 Hz point. Because the fKV model has G*(0)=μ_e, μ_e is essentially fixed by this anchor (Table 4 matches Table 3 except CR 0.053 vs 0.057). No data exist between 0.053 Hz and 300 Hz, so the model is an interpolation through an unmeasured gap. The 0 Hz anchor itself depends on the assumptions that (i) the linear-equivalent storage modulus extracted from large-strain nonlinear tests is the correct zero-frequency limit of the linear MRE response, and (ii) samples measured at 1–8 h post mortem (rheometer) and 2–3.5 h post mortem (MRE) represent the same material. Also, the low-frequency loss modulus is entirely unconstrained: the 0 Hz anchor carries no G″ information, and rheometer-derived G″ at 0.027/0.053 Hz (Fig. 6) is not included in the calibration. A model can fit a 0 Hz anchor and high-frequency MRE while misrepresenting the intermediate band; the conclusion 'comprehensive characterization' is therefore stronger than the data support.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper characterizes three porcine brain regions (corona radiata, putamen, thalamus) with two techniques: multi-modal large-strain rheometer tests (compression/tension/shear) analyzed with a nonlinear Ogden-Prony model, and tabletop MRE at 15 frequencies between 300 and 2100 Hz analyzed with a Bessel-function wave solution. The authors then calibrate a two-element fractional Kelvin-Voigt (fKV) model, G*(ω)=μ_e + Σ c_i(iω)^β_i, using the MRE dynamic moduli together with a single 0 Hz storage-modulus value computed from the Prony fit to the quasi-static rheometer data. They report consistent regional differences across both domains, observe a transition from elastic-dominated to viscous-dominated behavior with increasing frequency, and claim that the calibrated fKV model unifies the two responses into one 'comprehensive characterization' over 0–2100 Hz.","tokens_in":17948,"tokens_out":4383,"duration_ms":45941,"significance":"If the unified-model claim were fully supported, the paper would provide a practical bridge between quasi-static and high-frequency measurements, with region-specific parameters for a single fractional constitutive law. The experimental work is careful and the regional comparison is a strength: the corona radiata is consistently stiffest, and the identified μ_e values reproduce the rheometer-derived 0 Hz anchor. The paper is also transparent about the two-stage calibration and the fact that the MRE data contain no low-frequency information. However, the central claim is currently stronger than the evidence: the 0 Hz anchor is model-derived rather than measured, there are no experimental data between 0.053 Hz and 300 Hz, and the low-frequency loss modulus is not used in the calibration. The model therefore interpolates across the entire intermediate band, and the 'comprehensive characterization' conclusion should be substantially tempered.","major_comments":[{"comment":"The 0 Hz anchor used in the fKV calibration is not a measured quantity. It is computed from the Ogden-Prony parameters identified in §2.2.2 via Eq. (7), i.e., G'(0)=μ_∞ = μ_0(1-Σg_k). Since the fKV model has G*(0)=μ_e, the resulting μ_e values in Table 4 are essentially identical to the Table 3 anchors (0.053 vs 0.057, 0.063 vs 0.063, 0.052 vs 0.052). The low-frequency branch of the model is therefore enforced by construction, not validated by data. The paper itself notes in §3.3 that μ_e correlates with the rheometer storage modulus 'since those values were added to the data set for the model calibration.' The abstract and conclusion claim a 'unified' and 'comprehensive' characterization; this overstates what is demonstrated. At most, the model provides a constrained interpolation anchored at one model-derived point at 0 Hz and MRE data from 300 Hz upward.","section":"§2.3.2, Tables 3 and 4"},{"comment":"There are no experimental data between 0.053 Hz and 300 Hz. The model is fit to one 0 Hz point and to 300–2100 Hz data; the excellent visual agreement in Fig. 8 is an interpolation across an unmeasured band that spans more than three orders of magnitude in frequency. Moreover, the low-frequency loss modulus is entirely unconstrained: the rheometer-derived G'' values at 0.027 and 0.053 Hz (shown in Fig. 6 and computed from Eq. 7) are not included in the calibration, and the 0 Hz anchor carries no G'' information. Consequently, the fKV model can match the endpoint storage modulus and the high-frequency response while misrepresenting the intermediate regime. The claim that the model 'very well captures the response over the whole frequency range' (Sec. 3.3) should be restricted to the measured frequencies, or the authors should add the low-frequency G'' points and/or a sensitivity analysis","section":"§2.3.2, Fig. 8; §3.3"},{"comment":"The fitted μ_2 values are 2.23e-16 kPa (CR), 1.21e-13 kPa (P), and 2.13e-10 kPa (T). These are numerically zero and span six orders of magnitude. With β_2 ≈ 0.93–0.96, the contribution c_2 = μ_2^(1-β_2) remains finite, but the parameter μ_2 itself is nonphysical and the model appears overparameterized for the available data. This is a direct symptom of the missing intermediate-frequency constraints: the second fractional element is tasked to bridge a 0 Hz anchor and high-frequency MRE without any data between. The authors should discuss identifiability, fix or constrain μ_2, or report the effective c_2 instead of μ_2.","section":"Table 4 and Fig. 10"},{"comment":"Rheometer measurements were conducted 1–8 h post mortem, while MRE measurements were performed 2–3.5 h post mortem. The discussion cites Weickenmeier et al. (2018) for post-mortem stiffening, so the two datasets may represent different tissue states. The 0 Hz anchor depends on the rheometer data, and the high-frequency MRE data on the MRE samples; if post-mortem time affects stiffness, the 'unified' model merges responses from different material states. The authors should either quantify the effect of post-mortem time or explicitly list this as a limitation of the combined calibration.","section":"§2.1"}],"minor_comments":[{"comment":"The discussion states that the storage modulus is 'related to the initial shear modulus μ_0'. This is not correct: the storage modulus at 0 Hz from Eq. (7) is the long-term shear modulus μ_∞ = μ_0(1-Σg_k), not μ_0. The initial shear modulus is the zero-frequency limit of the *relaxed* modulus, not the storage modulus. Please correct this wording.","section":"§4"},{"comment":"The phrase 'we included the storage modulus ... as an additional parameter at 0 Hz' is confusing: the storage modulus is a data point, not a parameter. Rephrase, e.g., 'we added a single data point at 0 Hz'.","section":"§2.3.2"},{"comment":"The upper panels show the model response from 0 to 100 Hz, but there are no data in this range. Plotting the model as a solid line across the unmeasured gap visually implies validation. Consider shading the interpolated region or plotting only measured frequencies.","section":"Fig. 9"},{"comment":"The entry for α_1, P vs T, reads '01'; this is likely a typo for '0.1' or '1'. Please check.","section":"Table A.7"},{"comment":"The p-values in the appendix are given without indicating which comparisons remain significant after Holm-Bonferroni correction. The text mentions the correction, but the tables show raw p-values. Please mark corrected significance levels.","section":"§2.4 / Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically sound in its two experimental components, but the central 'unified characterization' claim is not supported by the data as it stands. The most serious issue is not the absence of intermediate data per se — that is unavoidable in such multi-scale work — but the overstatement. The authors already acknowledge in Sec. 2.3.2 that the MRE data contain no low-frequency information; they should carry this acknowledgment into the abstract and conclusion. If they reframe the contribution as a framework for combining scales, add the low-frequency G'' points or a sensitivity analysis, and discuss the near-zero μ_2, the paper would be acceptable for publication. The current overclaim is fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one for the data, not for the master curve. The genuinely new thing is a single group measuring three porcine brain regions both quasi-statically (large-strain multi-modal rheometry with an Ogden–Prony inverse fit) and at 300–2100 Hz with tabletop MRE, then putting both into one two-element fractional Kelvin–Voigt model. The regional ordering (corona radiata stiffest, elastic-to-viscous crossover with frequency) is credible and useful, and the in-sample fits are good. That is a real contribution — most prior comparisons stop at lower frequencies or use phantoms or liver.\n\nWhat the paper does well: careful sample handling, documented post-mortem windows, disclosed exclusions (the CR 300 Hz points and gray-matter masking), region-wise statistics, and transparency that the MRE measurement contains no information below 300 Hz. The parameter tables are useful enough to reproduce and build on.\n\nSoft spots, in order. First, the 0 Hz anchor is not a measurement; it is G'(0) computed from Prony parameters that were themselves identified from large-strain nonlinear tests. The paper says this in Sec. 2.3.2, so the low-frequency end of the \"unified model\" is model-extrapolated. Second, there is no data between 0.053 Hz and 300 Hz. The fKV curve is an interpolation across that gap, not a validated prediction; \"comprehensive\" is stronger than the evidence supports. Third, the rheometer loss modulus at 0.027/0.053 Hz appears in Fig. 6 but is not included in the calibration, so low-frequency loss behavior is essentially unconstrained. Fourth, no uncertainty is reported on the final fKV parameters, and the fitted μ2 values (1e-16 to 1e-10 kPa) are so small they look like a numerical placeholder rather than a physically meaningful second element. None of these are fatal; they are calibration choices that should be stated as limitations.\n\nWho should read it: anyone building porcine brain models for surgical or TBI simulation, and anyone comparing MRE with rheometry. I would cite it for the regional data, not for the claim that a single model covers 0–2100 Hz without further validation.\n\nYes, send it to peer review. Ask for out-of-sample validation in the 0.1–300 Hz window, parameter uncertainties, and a softer abstract.","headline":"Useful regional porcine brain data and a sensible fKV fit, but the 0–2100 Hz \"unification\" is an interpolation across an unmeasured gap anchored by a model-derived 0 Hz point.","tokens_in":18499,"tokens_out":3085,"would_cite":true,"duration_ms":36692,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92C10","74D10"],"pacs":["87.19.R-","83.60.Bc"],"model":"deepseek-v4-flash","headline":"The paper claims that a five-parameter fractional Kelvin–Voigt model unifies quasi-static rheometer data and 300–2100 Hz MRE data for three porcine brain regions into one viscoelastic response curve.","keywords":["Brain tissue","Viscoelasticity","Fractional Kelvin-Voigt model","Magnetic resonance elastography","Rheometry","Porcine brain","Regional brain mechanics","Multi-scale characterization"],"falsifier":"Run small-strain oscillatory shear tests on the same porcine brain regions at intermediate frequencies (roughly 0.1–300 Hz) and compare the measured G' and G'' with the fractional Kelvin–Voigt model's prediction. If the measurements deviate systematically from the model's smooth power-law bridge, or if the 0 Hz anchor shifts with strain amplitude or post-mortem time, the unified characterization fails. A quicker check: perturb the 0 Hz G'(0) value by its estimated uncertainty and see whether the 300–2100 Hz MRE fit still holds.","tokens_in":17327,"feed_emoji":"🧠","tokens_out":8601,"duration_ms":80021,"temperature":0.7,"pith_summary":"Brain-tissue mechanics has long suffered from contradictory results depending on whether the tissue is loaded slowly or by high-frequency waves, and most studies characterize only one scale. This paper measures the corona radiata, putamen, and thalamus of porcine brain in both regimes—multi-modal rheometry (compression, tension, torsion, and relaxation) plus tabletop magnetic resonance elastography at 300–2100 Hz—and calibrates a fractional Kelvin–Voigt model with two springpot elements in parallel to a spring. The central claim is that this five-parameter model, anchored at 0 Hz by a storage modulus computed from the quasi-static data, reproduces the entire frequency response for all three regions, including the switch from elastic-dominated behavior at low frequencies to viscous-dominated behavior at high frequencies. A sympathetic reader cares because a single mechanical law across time scales would let simulations of surgery, injury, and development use one set of parameters instead of choosing between inconsistent measurements.","feed_headline":"Five parameters unify brain tissue response from 0 to 2100 Hz","feed_subtitle":"Porcine brain rheometry and MRE data collapse into one fractional model, capturing the elastic-to-viscous shift.","key_machinery":"The central object is the fractional Kelvin–Voigt model with two fractional elements (springpots) in parallel with an elastic spring: G*(ω)=μ_e + c1(iω)^{β1} + c2(iω)^{β2}, with 0≤β_i≤1 and β1<β2. Each springpot constant is c_i = μ_i^{1-β_i} η_i^{β_i} with η_i fixed to 1 Pa s, so the five free parameters are μ_e, μ_1, β_1, μ_2, β_2. The model is calibrated by adding the rheometer-derived storage modulus at 0 Hz as an extra data point to the MRE frequency sweep, then minimizing a least-squares residuum. The lower-exponent element (β between 0.34 and 0.46) sets the storage modulus and low-frequency behavior; the higher-exponent element (β between 0.93 and 0.96) sets the loss modulus at high fr","core_discovery":"On its own terms, the paper establishes that a single fractional Kelvin–Voigt model, G*(ω)=μ_e + c1(iω)^{β1} + c2(iω)^{β2}, unifies the viscoelastic response of ex vivo porcine brain tissue. The spring fixes the zero-frequency stiffness; one springpot shapes storage modulus and low-frequency behavior, the other the loss modulus at high frequencies. The 0 Hz anchor comes from the storage modulus computed from the rheometer Prony fits; with it, the model reproduces the MRE data (300–2100 Hz) and the low-frequency points for all three regions. The same regional ordering appears in both domains (corona radiata stiffest), and the loss modulus rises faster than the storage modulus, so the loss tan","pith_inferences":["If the same anchoring procedure were applied to human brain data from indentation and in vivo MRE, it might reconcile many contradictory stiffness values reported in the literature; the authors stop short of making that claim.","The model's smooth power-law bridge across the unmeasured 0.053–300 Hz gap is a concrete prediction: small-strain oscillatory shear tests in that range would show whether the true response follows the fractional curve or has additional features.","The fitted second springpot shear moduli are extremely small (near 10^-16 kPa for corona radiata), so the high-frequency element is nearly a pure dashpot; whether that degeneracy persists in human tissue or at even higher frequencies is a testable question."],"forward_implications":["A single five-parameter set describes brain tissue from slow surgical loading to MRE wave frequencies, so finite-element simulations can cover both regimes without re-fitting the material law.","Regional ordering detectable in both domains—corona radiata stiffest, putamen and thalamus similar—means whole-brain models should use region-specific parameters at all scales.","The model quantitatively predicts the loss tangent crossing from below one to above one, locating where brain tissue turns from elastic-dominated to viscous-dominated, which is relevant for impact and injury simulations.","The two springpot exponents separate roles (one for storage/elastic response, one for loss/viscous response), giving a compact way to summarize regional mechanical state.","The same bridging strategy worked for three brain regions, suggesting it could be applied to other soft tissues with scale-dependent disagreements."],"fun_headline_variants":["One model unifies brain stiffness from 0 to 2100 Hz","Fractional model ties brain rheology across scales","Single viscoelastic law spans brain tissue from 0 to 2100 Hz","Rheometry and MRE converge in fractional Kelvin-Voigt fit"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The unified curve rests on treating the storage modulus computed from large-strain, nonlinear rheometer fits as the true zero-frequency limit of the small-amplitude linear MRE response, across an unmeasured 0.053–300 Hz gap and using samples with different post-mortem times.","fun_headline_variants_meta":{"raw":{"variants":["One model unifies brain stiffness from 0 to 2100 Hz","Fractional model ties brain rheology across scales","Single viscoelastic law spans brain tissue from 0 to 2100 Hz","Rheometry and MRE converge in fractional Kelvin-Voigt fit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000557,"raw_usage":{"total_tokens":2510,"prompt_tokens":787,"completion_tokens":1723,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":531,"completion_tokens_details":{"reasoning_tokens":1647}},"tokens_in":531,"tokens_out":1723,"duration_ms":13518,"temperature":1.0,"reasoning_tokens":1647,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T06:26:25.171902+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run small-strain oscillatory shear tests on the same porcine brain regions at intermediate frequencies (roughly 0.1–300 Hz) and compare the measured G' and G'' with the fractional Kelvin–Voigt model's prediction. If the measurements deviate systematically from the model's smooth power-law bridge, or if the 0 Hz anchor shifts with strain amplitude or post-mortem time, the unified characterization fails. A quicker check: perturb the 0 Hz G'(0) value by its estimated uncertainty and see whether the 300–2100 Hz MRE fit still holds.","supporting_citations":[],"review_version":1}