{"id":"50567e62-12c6-4195-b4b3-a93c45dadbe3","arxiv_id":"2507.20107","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A guided-wave optical coherence elastography technique quantifies layer-specific, nonlinear, anisotropic viscoelastic properties of arterial walls and shows that stretching reduces arterial viscosity.","lead":"Researchers measured how pig aorta walls damp vibrations using a high-resolution optical technique and found that stretching the artery makes it less viscous and more elastic. This could lead to a non-invasive way to profile arterial stiffness and disease risk, complementing existing ultrasound tests.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The stretch-driven decline in fitted δ may reflect parameter trade-off rather than a real viscosity decrease; a stretch-independent-viscosity refit is needed.","rationale":"The reader identified the shared-viscosity assumption in the two-layer model as the weakest point. I agree that this assumption needs scrutiny, but the more load-bearing concern is whether any stretch-dependent viscosity change is identifiable from the data at all. The single-layer model, which avoids the shared-viscosity constraint, also shows the same decreasing trend in δ, so relaxing that constraint may not change the conclusion. A stronger test is to compare against a null model with stretch-independent viscosity. The attenuation data are a helpful independent check, but they cover only part of the stretch range and are interpreted through the same viscoelastic framework. The proposed refit would directly test whether the dispersion data require the viscosity parameters to change with stretch, which is the crux of the central claim. This warrants conditional acceptance with a request for the control analysis, consistent with the reader's verdict.","tokens_in":35945,"tokens_out":6251,"duration_ms":65899,"concrete_test":"Refit the single-layer and two-layer viscoelastic models to the dispersion data at λ=1.2, 1.3, and 1.4 with η and δ fixed to their fitted values at λ=1.0 (or to a single shared constant), while allowing α, β, and γ to vary freely. Compare the residual sum of squares or AIC against the unconstrained fits. If the constrained model is not significantly worse (e.g., F-test p>0.05 or RMS increase <10%), the data do not require stretch-dependent viscosity, and the central claim should be weakened. If the constrained model clearly fails, the viscosity trend is confirmed.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that arterial viscoelasticity decreases with stretch is inferred from the KVFD model parameters η and δ, whose fitted values show a large drop in δ (e.g., axial δ from 0.45±0.06 at λ=1.0 to 0.09±0.02 at λ=1.4). However, the dispersion data are limited to 1–20 kHz, with incomplete S0 mode data at high stretch ratios, and the authors acknowledge that parameter uncertainties exceed 50% due to interdependence. The paper does not report a control fit in which viscosity parameters (η, δ) are held constant at their λ=1 values while elastic moduli are allowed to vary with stretch. Without such a restricted-model comparison, the apparent decrease in δ could be an artifact of parameter trade-off against the large stretch-dependent increases in α, β, and γ. The attenuation measurements provide some independent support for reduced dissipation, but they are shown only for λ=1 to 1.2, and the attenuation includes radiation losses into the fluid that depend on elasticity; the separation of intrinsic viscosity from leakage relies on the same viscoelastic model whose parameters are in question. Therefore, the quantitative claim of reduced viscosity at higher stretches (λ=1.3–1.4) is not securely established by the current analysis.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a guided-wave optical coherence elastography (OCE) method for characterizing ex vivo porcine thoracic aorta. Phase velocities of the A0 and S0 guided wave modes were measured over 1-20 kHz in the axial and circumferential directions at equibiaxial stretch ratios from 1.0 to 1.4. The data are fitted sequentially with a single-layer elastic model, a single-layer Kelvin-Voigt fractional derivative (KVFD) viscoelastic model, and a two-layer viscoelastic model in which the media and adventitia have distinct elastic moduli but shared viscosity parameters. The authors report stretch-dependent increases in shear and tensile moduli, greater adventitial stiffening under load, and decreases in the fitted fractional-order parameter, loss modulus, and wave attenuation with increasing stretch, which they interpret as a reduction in arterial viscosity under prestress. Collagen removal with CNBr reduces elastic moduli and the GOH stiffening exponent while leaving the fitted viscosity parameters roughly unchanged.","tokens_in":36258,"tokens_out":5925,"duration_ms":60506,"significance":"If the central claim holds, this is a useful advance: it demonstrates that a single OCE platform can constrain nonlinear, anisotropic, and layered viscoelastic properties of arteries, and it provides a quantitative data point on stretch-dependent dissipation. The strengths of the paper are the progressive model hierarchy (elastic to viscoelastic to two-layer), the internal consistency check via attenuation prediction, the comparison with literature constitutive parameters, and the detailed analytic derivations in the Supplementary Notes. The main risk is identifiability: the viscosity trend is inferred from a many-parameter fit with acknowledged uncertainties exceeding 50%, and the paper does not currently demonstrate that the trend survives a restricted-model comparison.","major_comments":[{"comment":"The central claim that arterial viscosity decreases with stretch is supported mainly by the fitted KVFD parameters in Table 3, where the fractional order delta decreases from 0.45 +/- 0.06 to 0.09 +/- 0.02 (axial) and from 0.46 +/- 0.06 to 0.05 +/- 0.03 (circumferential) as lambda goes from 1.0 to 1.4. However, no control fit is reported in which eta and delta are fixed at their lambda = 1 values while only the elastic parameters are allowed to vary with stretch. Because the elastic moduli increase substantially over the same range and the authors state in the Discussion that parameter uncertainties exceeded 50% due to interdependence, the apparent delta decline could reflect parameter trade-off rather than a real reduction in dissipation. Please add such a restricted-model comparison, or an equivalent identifiability analysis such as profile likelihood or bootstrap over the joint parameter distribution, and report the resulting fit quality for the higher-stretch cases.","section":"Viscoelastic two-layer wave model analysis; Table 3"},{"comment":"The two-layer model assumes identical eta and delta across media and adventitia, stated as a way 'to reduce the number of free parameters.' Yet the abstract and discussion describe layer-specific viscoelastic characterization. If the true layer-specific rheologies differ, the fitted shared viscosity parameters, and hence the stretch-dependent trend, may be biased averages. The assumption should be justified, relaxed in a sensitivity analysis, or the viscosity claim should be explicitly framed as an effective whole-wall property rather than a layer-specific intrinsic property.","section":"Viscoelastic two-layer wave model analysis"},{"comment":"The attenuation evidence in Fig. 6c is limited to lambda = 1 to 1.2, whereas the largest drop in the fitted delta occurs at lambda = 1.3-1.4. Moreover, the predicted attenuation is computed from the same viscoelastic model and includes radiation losses into the fluid that depend on the elastic parameters. As presented, this validation does not independently confirm the high-stretch viscosity decrease; please extend the attenuation measurements to higher stretch ratios or state this limitation explicitly.","section":"Stretch-dependent viscosity parameters of arterial tissues; Fig. 6c"},{"comment":"In Table 3, the circumferential adventitia shear modulus alpha is listed as 48 +/- 1 kPa at lambda = 1.4, after values of 32, 110, 220, and 325 kPa at lower stretches. This is inconsistent with the monotonic stiffening described in the text and shown in Fig. 5b. If this is a typographical error (e.g., 480 kPa), it must be corrected; otherwise it undercuts the layer-specific stiffening claim.","section":"Table 3"}],"minor_comments":[{"comment":"The phrase 'gold clinical gold standard' is a typo; revise to 'gold standard'.","section":"Introduction"},{"comment":"The exponential notation in the text accompanying Eq. (3) is inconsistent: the exponent should be sky + i(kx - omega t), not the form currently printed as e^{six} e^{i(kx - omega x)}.","section":"Equation (3)"},{"comment":"The sentence 'incompressibility yields lambda_x = (lambda_x lambda_z)^{-1}' should read lambda_y = (lambda_x lambda_z)^{-1}, or equivalently lambda_y = 1/(lambda_x lambda_z).","section":"Methods, Analytic modeling"},{"comment":"The statement that the S0 mode attenuation curves are similar is only demonstrated in Supplementary Fig. S2 at lambda = 1; please clarify whether the comparison was made at all stretch ratios.","section":"Fig. 6c and Supplementary Fig. S2"},{"comment":"The phrase 'A p-values less than 0.05' should be 'A p-value less than 0.05'.","section":"Statistical analysis"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies on the authors' own viscoelastic theory (reference 52), and the inversion code and data are not provided. Combined with the parameter-identifiability concern, this makes independent verification difficult. I do not regard this as a scientific flaw, but the journal may wish to encourage a data and code availability statement before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here's my read of arXiv:2507.20107.\n\nThe paper is a serious step forward in elastography: it adapts the authors' multi-wave OCE to the arterial wall, builds a two-layer guided-wave model with fractional viscoelasticity, and extracts layer-specific shear and tensile moduli under biaxial stretch. The experimental campaign is substantial — 10 porcine aortas, multiple stretch states, two directions, 1–20 kHz. The modeling is detailed, with explicit secular equations and a clear link to the GOH microstructure. The attenuation prediction from the fitted parameters (Fig 6c) is a legitimate internal check that the model is doing something right.\n\nThe headline claim — arterial viscosity drops as stretch increases — is plausible and consistent with the loss modulus and attenuation trends, but it is not nailed down. The fractional order δ falls from ~0.45 to ~0.1 between λ=1 and λ=1.4, while the elastic moduli more than double. With a KVFD model, δ is tradeable against the elastic parameters, and the paper does not report a refit where η and δ are held at their λ=1 values and only elastic parameters are adjusted. That control is necessary to rule out parameter migration. The two-layer model also forces identical η and δ across media and adventitia, which could bias the layer-specific elastic parameters and the viscosity trend. The authors honestly note that some parameter uncertainties exceed 50%, and the data/code are not public. The attenuation support is only shown up to λ=1.2 and includes radiation losses that share the same model.\n\nThat said, the elastic trends — circumferential stiffer than axial, adventitia taking over load at high stretch, collagen depletion reducing both stiffness and nonlinearity — are well supported and match the biaxial literature. The viscosity finding is likely correct in trend but needs a tighter identifiability argument; this is a fixable weakness, not a fatal one.\n\nThis paper deserves a serious referee. It is a real method with a real dataset. For the revision, I'd ask for the control fit, a sensitivity analysis on the equal-viscosity assumption, and open data or at least a simulated-data tutorial problem. The readership is biomechanics and elastography; anyone modeling arterial tissue would cite it, though I'd wait for the revision before citing the viscosity result quantitatively.","headline":"Strong OCE methodology with a real dataset, but the headline viscosity-with-stretch result needs a control fit before I'd trust the numbers.","tokens_in":36784,"tokens_out":2038,"would_cite":true,"duration_ms":21417,"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":"Arterial walls become markedly less viscous and more elastic as they are stretched, and the collagen-rich adventitia overtakes the media in stiffness under load.","keywords":["arterial biomechanics","optical coherence elastography","Lamb waves","nonlinear viscoelasticity","layer-specific mechanical properties","fractional viscoelastic model","collagen and elastin contributions","biaxial stretch"],"falsifier":"Perform dynamic mechanical analysis or slow stress-relaxation tests on the same porcine aorta tissue at the same equibiaxial stretch levels ($\\lambda = 1.0$ to $1.4$) and measure the loss tangent directly; if dissipation does not systematically fall as stretch increases, the OCE-inferred viscosity reduction is a modeling artifact. A cheaper check: refit the two-layer guided-wave model with independent viscosity parameters for media and adventitia and see whether the stretch-dependent decline in $\\eta$ and $\\delta$ survives the extra degrees of freedom.","tokens_in":35718,"feed_emoji":"🫀","tokens_out":12143,"duration_ms":108392,"temperature":0.7,"pith_summary":"Using broadband guided-wave optical coherence elastography on biaxially stretched porcine aortas, the paper claims that arterial walls become markedly less viscous as they are stretched: elastic shear and tensile moduli rise with prestretch while loss moduli, a fractional viscoelastic order ($\\delta$), and wave attenuation all fall. The claim matters because vascular stiffness and energy dissipation are central to diseases like hypertension and aneurysm, yet standard clinical measures such as pulse wave velocity cannot resolve layer-specific or viscoelastic behavior. The paper also argues that the collagen-rich adventitia is softer than the media at zero load but overtakes it in stiffness once stretch engages its collagen fibers, and that chemically removing collagen collapses the elastic moduli while leaving the viscoelastic 'spring-pot' nearly unchanged — evidence that elastin carries most low-strain viscous damping.","feed_headline":"Stretched arteries turn more elastic and less viscous","feed_subtitle":"Guided-wave optical elastography shows prestretch cuts viscous loss while collagen drives layer-specific stiffening.","key_machinery":"The load-bearing machinery is the Kelvin-Voigt fractional derivative (KVFD) model — a 'spring-pot' with complex modulus $\\Omega = \\eta(i\\omega)^\\delta$ coupled in parallel with an elastic spring — embedded in a generalized acousto-viscoelastic theory of small waves superimposed on a finitely pre-stressed solid. This replaces the purely elastic incremental stress tensor with a frequency-dependent one, and the dispersion of the two guided Lamb modes, A0 (quasi-antisymmetric, shear/bending dominated at low frequency) and S0 (quasi-symmetric, tensile/dilatational dominated), is solved through secular equations: 5×5 for a single layer and 9×9 for the two-layer media-adventitia model. The two-layer model gives each layer independent elastic moduli ($\\alpha_1, \\beta_1, \\gamma_1$ for the media; $\\alpha_2, \\beta_2, \\gamma_2$ for the adventitia) while sharing one pair of viscous parameters ($\\eta$, $\\delta$), and its high-frequency asymptotics divide into four regimes by the modulus ratio of adventitia to media, explaining why at zero stretch the modes asymptote to the adventitia's Scholte and shear velocities but under stretch approach the media's Rayleigh-wave limits. The Gasser-Ogden-Holzapfel constitutive model then converts the stretch-dependent moduli into fiber parameters: ground-matrix shear modulus $\\mu_0$, fiber stiffness $k_1$, nonlinear exponent $k_2$, fiber angle $\\varphi$, and dispersion $\\kappa$.","core_discovery":"The core discovery, stated on the paper's own terms, is that arterial viscoelasticity is strongly deformation-dependent: increasing prestress reduces tissue viscosity. Measuring the dispersion of the A0 (shear-dominated) and S0 (tensile-dominated) guided wave modes over 1-20 kHz under equibiaxial stretch from $\\lambda = 1.0$ to $1.4$, and fitting a two-layer viscoelastic guided-wave model, the authors find that shear and tensile moduli of both the media and adventitia rise with stretch, the adventitia stiffening faster than the media once $\\lambda$ exceeds about 1.1. Simultaneously, the imaginary (loss) parts of the complex shear and tensile moduli, the fractional order $\\delta$ of the Kelvin-Voigt fractional derivative model, and the directly measured wave attenuation all decrease with stretch, so the wall becomes more elastic and less energy-dissipating under load. Degrading collagen with cyanogen bromide reduces both shear and tensile moduli substantially but leaves the viscous parameters comparable to intact tissue, which the authors read as evidence that collagen governs nonlinear elastic stiffening while the elastin network dominates low-strain viscoelastic dissipation.","pith_inferences":["If the viscosity drop holds in vivo under pulsatile pressure, then damping of pulse waves should vary over the cardiac cycle as the wall loads and unloads; this predicts measurable cycle-dependent changes in wave attenuation that pulse-wave imaging could test — a consequence the paper does not state.","The paper leaves the KVFD extrapolation from 1-20 kHz to heartbeat frequencies near 1 Hz unvalidated; a direct test would compare OCE-derived loss moduli with slow stress-relaxation or dynamic mechanical analysis on the same samples to see whether a single power law spans both frequency ranges.","A testable prediction follows for aged or glycated arteries: if collagen crosslinking raises the nonlinear exponent $k_2$, the same protocol should show sharper adventitial stiffening while low-strain viscous damping stays roughly constant — a dissociation the authors do not draw.","Because the two layers share one viscosity pair, an identifiability analysis over the 1-20 kHz band would settle whether the apparent viscosity reduction reflects a tissue property or a constraint imposed by the shared-parameter model."],"forward_implications":["Under physiological prestretch, arterial walls dissipate less energy per pulse cycle than at low pressure, so vascular damping is pressure-dependent rather than a fixed tissue constant.","Constitutive models of arteries should treat viscosity as a stretch-dependent quantity: the KVFD parameters $\\eta$ and $\\delta$ and the loss tangent all change with prestress, so a single-valued arterial viscosity is incomplete.","The mechanical load shifts between layers: at $\\lambda = 1$ the media is slightly stiffer than the adventitia, but by $\\lambda = 1.4$ the adventitial tensile modulus is roughly 3.5 times the medial value, so collagen engagement is the dominant nonlinearity.","Collagen and elastin play separated roles: collagen drives the nonlinear rise of elastic moduli with stretch, while the elastin network accounts for most low-strain viscoelastic damping, since collagen digestion removes stiffness but not the spring-pot.","Guided-wave OCE can non-destructively map layer-specific, anisotropic, nonlinear viscoelasticity at roughly 10 µm resolution over 1-20 kHz, a capability that standard pulse-wave-velocity and ultrasound methods lack."],"supporting_citations":[{"why":"Supplies the generalized acousto-viscoelastic theory, including the KVFD-modified incremental stress tensor that the single- and two-layer viscoelastic wave models are built on.","marker":"[52]"},{"why":"Provides the incremental statics-and-dynamics formalism for pre-stressed elastic solids and the Eulerian elasticity tensor used in the elastic base model.","marker":"[53]"},{"why":"The prior multi-wave (S0/A0) OCE method this study extends, and the source of the mode interpretation linking A0 to shear and S0 to tensile deformation.","marker":"[51]"},{"why":"Establishes the fractional-order viscoelastic description of human arteries from stress relaxation, including the delta range (0.1-0.3) the extracted values are compared against.","marker":"[54]"},{"why":"The Gasser-Ogden-Holzapfel constitutive model used to translate stretch-dependent moduli into collagen fiber parameters.","marker":"[77]"},{"why":"The biaxial tensile study of human descending aorta whose adventitial stiffening trend is cited for agreement with the measured layer behavior.","marker":"[22]"},{"why":"Supplies the elastin-anisotropy framework used to interpret the collagen-degraded (elastin-dominated) samples.","marker":"[17]"}],"fun_headline_variants":["Stretch cuts arterial viscosity, boosts elasticity","Arteries under stretch become less viscous, more elastic","Collagen stiffens stretched arteries, cuts viscous loss","Prestretch reduces artery viscosity, stiffens layers","Stretch makes arteries more elastic, less lossy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that viscosity falls with stretch rests entirely on the assumption that one shared pair of Kelvin-Voigt fractional derivative parameters ($\\eta$, $\\delta$), applied identically to the media and adventitia, faithfully describes the wall's rheology — if the two layers dissipate energy differently, the apparent viscosity loss under stretch could be a modeling artifact rather than a real tissue property.","fun_headline_variants_meta":{"raw":{"variants":["Stretch cuts arterial viscosity, boosts elasticity","Arteries under stretch become less viscous, more elastic","Collagen stiffens stretched arteries, cuts viscous loss","Prestretch reduces artery viscosity, stiffens layers","Stretch makes arteries more elastic, less lossy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000174,"raw_usage":{"total_tokens":1285,"prompt_tokens":952,"completion_tokens":333,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":568,"completion_tokens_details":{"reasoning_tokens":258}},"tokens_in":568,"tokens_out":333,"duration_ms":3904,"temperature":1.0,"reasoning_tokens":258,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:49:01.436174+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform dynamic mechanical analysis or slow stress-relaxation tests on the same porcine aorta tissue at the same equibiaxial stretch levels ($\\lambda = 1.0$ to $1.4$) and measure the loss tangent directly; if dissipation does not systematically fall as stretch increases, the OCE-inferred viscosity reduction is a modeling artifact. A cheaper check: refit the two-layer guided-wave model with independent viscosity parameters for media and adventitia and see whether the stretch-dependent decline in $\\eta$ and $\\delta$ survives the extra degrees of freedom.","supporting_citations":[{"cited_title":"& Yun, S.-H","cited_arxiv_id":null,"evidence_quote":"Supplies the generalized acousto-viscoelastic theory, including the KVFD-modified incremental stress tensor that the single- and two-layer viscoelastic wave models are built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the incremental statics-and-dynamics formalism for pre-stressed elastic solids and the Eulerian elasticity tensor used in the elastic base model."},{"cited_title":"J., Atienza, J","cited_arxiv_id":null,"evidence_quote":"Establishes the fractional-order viscoelastic description of human arteries from stress relaxation, including the delta range (0.1-0.3) the extracted values are compared against."},{"cited_title":"C., Ogden, R","cited_arxiv_id":null,"evidence_quote":"The Gasser-Ogden-Holzapfel constitutive model used to translate stretch-dependent moduli into collagen fiber parameters."}],"review_version":2}