{"id":"cf49db91-b332-4726-a1b7-44966433ca86","arxiv_id":"2506.03816","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An acoustoelastic formula connects the difference of two perpendicular shear-wave speeds to the local stress difference, enabling stress imaging in hydrogels and muscle without constitutive parameters.","lead":"Researchers measured local stress inside soft materials by tracking shear waves sent in two perpendicular directions with an ultrasound transducer, without knowing the material's stiffness. The method could make it possible to image stresses in tissues and soft machines non-invasively, using equipment already found in clinics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (2) is only established for coaxial stress and anisotropy; the muscle test stretches along the fiber, so non-coaxial tissue orientations remain an untested failure mode for the claim of general applicability.","rationale":"The paper's central theoretical result is the parameter-free identity (2). The stress identity in SM Note 1 is general and correct, and the principal-axis speeds equal the appropriate moduli when the measurement axes are symmetry directions. The critical step is SM Note 2's coaxiality assumption. Our reading of the acoustic tensor algebra shows the stress identity alone is not enough: the modes must not mix. The muscle experiment deliberately aligns the stretch with the fibers, so it cannot detect failure when the fiber direction has an out-of-plane component. Since the abstract promises applications to tissues where stress and fiber directions may not coincide, and no analytical or experimental treatment is given for that case, the central claim's generality is conditional. The reader's weakest-assumption selection of coaxiality matches our concern, though we would phrase it as 'measurement plane must be a symmetry plane.' We therefore recommend keeping the CONDITIONAL verdict. The typo in Eq. (1) vs SM Eq. (13) is real but does not affect Eq. (2).","tokens_in":15538,"tokens_out":32801,"duration_ms":281434,"concrete_test":"Solve the Christoffel equation for the SM Note 5 transversely isotropic acoustoelastic model with fiber direction m=(cos45°,sin45°,0) and uniaxial stress along x. Determine the SV-mode speeds for n=(1,0,0) and n=(0,0,1) and test whether ρ(vx²−vz²)=σ1−σ3. If the equality fails, the coaxiality assumption in SM Note 2 is load-bearing and Eq. (2) requires the measurement plane to be a material symmetry plane. An experimental complement would use an anisotropic phantom with fibers at 45° out of the imaging plane and compare inferred and applied stresses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"SM Note 2 posits coaxiality of stress, deformation, and structural anisotropy so that the waves along x and z are pure shear modes with ρv_x²=A0_1313 and ρv_z²=A0_3131. The stress identity A0_1313−A0_3131=σ11−σ33 (SM Eq. 4) itself holds without that assumption, but the identification of measured speeds with those moduli does not. The muscle experiment stretches along the fiber direction, so it only validates the fully coaxial case. If fibers are not contained in the x–z measurement plane (e.g., m=(cosφ,sinφ,0)), the x–z plane is not a symmetry plane, the in-plane SV and out-of-plane SH modes couple, and the z-polarized particle velocity tracked by ultrasound does not propagate at the pure-mode speed. Eq. (2) then fails, so the abstract's broad claim to measure stress in soft tissues is not supported for realistic non-coaxial fiber architectures.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes and tests an acoustoelastic shear-wave method to measure the difference of principal stresses in soft materials without knowing their constitutive parameters. The theoretical basis is Eq. (2), derived from a constitutive-independent identity (SM Eq. 4) together with the assumption that all anisotropy is coaxial with the stress, so that waves along the two principal directions propagate as pure shear modes. The authors implement the method with a medical ultrasound transducer programmed to apply six successive acoustic radiation forces, producing both lateral and vertical shear waves; vertical speeds are extracted from z > 7 mm data with an interference correction. They validate the method on a PVA hydrogel under uniaxial and bending stress (errors ~5%) and on ex vivo porcine muscle under passive uniaxial stretch along the fiber direction (errors ~15%), and they analyze viscoelastic effects using a QLV model.","tokens_in":15707,"tokens_out":5719,"duration_ms":62448,"significance":"The stress identity and the basic experimental demonstration are valuable: if the coaxiality condition holds, the method offers a parameter-free route to local stress imaging, and the hydrogel and muscle experiments compare against externally applied loads rather than fitted stress values. The programmed-ARF excitation is a useful technical advance. However, the paper's broad claim of applicability to soft tissues is not yet supported for non-coaxial fiber/stress architectures, and the vertical-wave measurement relies on a correction/cutoff whose influence is not quantified. With those gaps addressed or the claims appropriately restricted, the work would be a solid contribution.","major_comments":[{"comment":"SM Note 2 assumes the coaxiality conditions Cτ = τC, CMMᵀ = MMᵀC, and τMMᵀ = MMᵀτ, and Eq. (12) restricts the moduli; without these assumptions the measured ρv_x² and ρv_z² cannot be identified with A⁰₁₃₁₃ and A⁰₃₁₃₁. The muscle experiment stretches along the fiber, so it only validates the coaxial case. The abstract states that the method will find broad applications for diagnosing diseases that alter stresses in soft tissues, but in general tissue fibers need not be coaxial with the in vivo stress; in that case the x–z plane is not a symmetry plane and the tracked SV mode is not a pure mode, so Eq. (2) fails. The authors should either extend the theory to non-coaxial situations, provide a validation with fibers not aligned with the principal stress, or clearly restrict the claimed domain of applicability.","section":"SM Note 2 and main text Eq. (2)"},{"comment":"The vertical speed is not measured as a direct plane-wave speed: the source interference makes the apparent speed approach vz only for z > 7 mm, and the paper relies on this cutoff without a sensitivity analysis. Since Eq. (2) uses v_x² − v_z², a systematic bias in vz maps directly into a bias in the inferred stress. The zero-stress isotropy check (vx ≈ vz) and the hydrogel ground-truth comparisons provide indirect support, but a quantitative study of how the z cutoff and the interference correction affect the inferred stress would be needed to support the claimed accuracy in less controlled, anisotropic samples.","section":"Generating shear waves / Fig. 2G"},{"comment":"The central result provides only the difference of two in-plane principal stresses, σ₁ − σ₃; for a general biaxial stress state, individual stress components require additional assumptions or measurements, and only for uniaxial stress (σ₃ = 0) does Eq. (2) directly give σ₁. The paper sometimes refers broadly to measuring 'mechanical stresses', and this limitation should be stated more prominently in the abstract and discussion to avoid overstating what the method delivers in realistic tissue loading conditions.","section":"Scope of Eq. (2)"}],"minor_comments":[{"comment":"The powers in Eq. (1) of the main text appear different from those in SM Eq. (13): the main text writes ρv² = α cos²θ + 2β cos²θ sin²θ + γ sin²θ while the SM has cos⁴θ and sin⁴θ. Please align the notation and clarify whether 'cos2θ' means cos²θ or cos(2θ).","section":"Eq. (1) vs SM Eq. (13)"},{"comment":"Equation (14) in SM Note 2 writes σ₁ − σ₂ on the left-hand side, but the main text and the preceding derivation use σ₁ − σ₃; this appears to be a typo.","section":"SM Note 2, Eq. (14)"},{"comment":"The phrase 'the so-call longitudinal shear waves' should be 'the so-called longitudinal shear waves'.","section":"Results, 'Generating shear waves...' paragraph"},{"comment":"The caption uses 'AFRs' where 'ARFs' (acoustic radiation forces) is meant.","section":"Fig. 4C caption"},{"comment":"The identified stresses in the muscle experiment are shown without error bars; reporting the standard deviation of the five measurements would help readers judge the statistical significance of the ~15% maximum error.","section":"Fig. 4E"}],"recommendation":"major_revision","confidential_remarks":"The claim of constitutive-parameter-free stress measurement is strong and will attract attention. The main risk is overstatement: the method is validated only in coaxial geometries, while the clinical motivation requires non-coaxial cases. The paper should be pushed to either demonstrate the non-coaxial case in a tissue-mimicking phantom with controlled fiber orientation and oblique loading, or temper the abstract and discussion. The vertical-wave measurement procedure also needs robustness analysis before the accuracy claims can be fully trusted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper has a clean theoretical core, Eq. (2), and it earns its claim in the isotropic hydrogel and in muscle stretched along its fibers. The big question is whether the coaxiality assumption will hold in real tissues, where fiber direction and stress direction don't have to align.\n\nWhat is actually new: the stress identity σ1−σ3 = ρ(vx²−vz²)/cos 2θ0 was already in Li et al. JASA 2020; the paper says so and generalizes it to anisotropic coaxial solids in SM Note 2. The genuinely new pieces are the programmed ARF excitation that produces measurable vertical shear waves, and the muscle demonstration showing the method works in an anisotropic material under the coaxial condition. The validation against externally applied loads is the right way to test the claim, and the reported errors are honest: about 5% in hydrogel, 15% in muscle. They also discuss viscoelasticity explicitly instead of burying it.\n\nThe soft spots are real but localized. The weakest link is coaxiality. SM Note 2 assumes C commutes with the initial stress τ and the structural tensor MMᵀ, and Eq. (2) relies on the measured vx and vz being pure-mode speeds. If the fiber direction is not in the x–z measurement plane, the modes couple and the formula fails. The muscle test stretches along the fiber, so it only validates the aligned case. Real tissues often have fiber and stress directions mismatched, so the abstract's broad claim outruns the evidence. That should be stated as a limitation or tested. Second, the vertical shear wave speed is measured after a source-interference correction and a z > 7 mm cutoff; if that correction is biased, the inferred stress is biased. It looks reasonable, but I'd want a sensitivity analysis. There is also a typo-level inconsistency between main-text Eq. (1), which has cos²θ and sin²θ, and SM Eq. (13), which has cos⁴θ and sin⁴θ. Not load-bearing, but needs fixing. No code or raw data is provided, though the experimental description is detailed enough that replication seems plausible.\n\nWho this is for: soft matter physicists and biomechanics people working on elastography, acoustoelasticity, or non-invasive stress measurement. It deserves a serious referee; the central idea is solid and the validation is meaningful. I'd accept with revisions: either demonstrate non-coaxial validity or narrow the claims, add the sensitivity analysis, and fix the equation inconsistency.","headline":"A clean constitutive-independent stress formula with real experimental validation in aligned cases, but the coaxiality assumption leaves a genuine untested gap for realistic tissue architectures.","tokens_in":16282,"tokens_out":1928,"would_cite":true,"duration_ms":21036,"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":"This paper derives and demonstrates a constitutive-parameter-free formula that turns two perpendicular shear wave speeds into the principal stress difference in an incompressible soft material.","keywords":["acoustoelasticity","shear wave elastography","stress measurement","soft materials","acoustic radiation force","incompressible elasticity","skeletal muscle","ultrafast ultrasound"],"falsifier":"Take a transversely isotropic phantom whose fibres are deliberately inclined at a known angle to an applied uniaxial stress, measure $v_x$ and $v_z$ along the transducer axes, and compare the value $\\rho(v_x^2-v_z^2)/\\cos 2\\theta_0$ with the known applied principal stress difference. If the two agree within measurement error even when fibres and stress are misaligned, the coaxiality premise is unnecessary; if they diverge systematically with misalignment, the premise is load-bearing and the method's range is limited to aligned tissues.","tokens_in":1610,"feed_emoji":"〰️","tokens_out":2212,"duration_ms":87521,"temperature":0.7,"pith_summary":"The paper sets out to show that the difference of the two in-plane principal stresses inside an incompressible soft material can be measured non-invasively from the speeds of two shear waves travelling in perpendicular directions, without knowing the material's constitutive parameters. The load-bearing identity is $\\sigma_1-\\sigma_3=\\rho(v_x^2-v_z^2)/\\cos 2\\theta_0$, which the paper derives from the acoustoelastic wave equation under a coaxiality assumption. The authors implement the measurement by programming six acoustic-radiation-force pushes so that lateral and vertical shear waves are generated simultaneously, then imaging them with ultrafast ultrasound. They demonstrate the method on a hydrogel under uniaxial and bending stress and on ex vivo skeletal muscle under passive tension, with maximum reported errors of about 5% and 15% respectively.","feed_headline":"Two shear wave speeds reveal stress with no material constants","feed_subtitle":"A pair of perpendicular shear wave speeds yields the stress difference directly, with no calibration of the material.","key_machinery":"The machinery is equation (2) plus a way to excite and time two perpendicular shear waves. A focused ultrasound beam applies acoustic radiation force; by stepping the focus through six points at $d=1$ mm spacing, the authors create a supershear moving load (Mach number about 10) whose Huygens-Fresnel interference amplifies the otherwise weak vertical shear wave. Speeds $v_x$ and $v_z$ are extracted from spatiotemporal maps by Radon transformation, and because phase and group speeds coincide along the principal axes, the measured group speeds feed directly into equation (2). The derivation itself rests on the stress identity $\\alpha-\\gamma=\\sigma_1-\\sigma_3$ for incompressible solids and on the coaxiality condition that structure tensors, initial stress, and deformation commute.","core_discovery":"The central discovery is that for an incompressible soft solid, the principal stress difference is encoded in two shear wave speeds exactly through $\\sigma_1-\\sigma_3=\\rho(v_x^2-v_z^2)/\\cos 2\\theta_0$, with no dependence on the strain-energy function or on third-order elastic constants. The identity comes from $\\alpha-\\gamma=\\sigma_1-\\sigma_3$, where $\\alpha=A^0_{1313}$ and $\\gamma=A^0_{3131}$ are Eulerian elastic moduli, combined with the wave-speed relations $\\rho v_x^2=\\alpha$ and $\\rho v_z^2=\\gamma$ for waves along the principal axes. It holds for arbitrary anisotropy, including fibre reinforcement and initial stress, as long as all sources of anisotropy are coaxial with the stress and deformation. The paper presents this as a parameter-free route to stress, validated by imaging uniaxial stress, bending stress, and passive muscle tension.","pith_inferences":["The identity gives only the difference of principal stresses; reconstructing the full in-plane stress tensor would require one more condition, such as a known stress-free boundary or an independent measurement of one principal stress, which the paper does not address.","Because the acoustic-radiation-force excitation is only one way to launch shear waves, the same formula should work with mechanical shakers or other sources, making the method portable to non-ultrasound settings.","A direct test of the coaxiality premise would be a phantom with fibres deliberately misaligned from the applied stress; the paper's muscle experiment only covers the aligned case.","For living tissues, passive tension along fibres is the favourable case; active contraction or residual stresses that are not fibre-aligned may violate the condition, so extending to those settings would need a multi-angle version of the measurement."],"forward_implications":["In hydrogels and other weakly viscoelastic soft solids, uniaxial and bending stress fields can be imaged remotely with errors around 5%, without calibrating the constitutive law.","The method works in anisotropic, fibre-reinforced tissues when fibres are aligned with the stress, as demonstrated by the passive muscle experiment, opening a route to in vivo stress estimation.","Because only the principal stress difference is obtained, quantitative use requires either a known zero-stress direction, such as a free surface, or a second independent condition; the paper uses $\\sigma_3=0$ for uniaxial cases.","Viscoelasticity biases the inferred stress downward, but the paper's analysis with a quasi-linear model indicates errors below about 10% for soft materials whose stress relaxation is less than about 50%.","The measured quantity is the group speed along the principal axes, where phase and group speeds coincide, so the method transfers directly to standard shear-wave-elastography systems."],"supporting_citations":[{"why":"establishes the isotropic version of the stress-wave-speed result that this paper generalizes to coaxial anisotropy","marker":"[25]"},{"why":"provides the theory of incremental waves in initially stressed incompressible solids","marker":"[28]"},{"why":"supplies the elliptical wavefront and group-versus-phase speed analysis used to justify measuring group speeds","marker":"[34]"},{"why":"provides the Radon-transform method used to extract shear wave speeds from spatiotemporal data","marker":"[41]"},{"why":"supplies the acousto-visco-elastic model used to quantify the bias introduced by tissue viscosity","marker":"[45]"},{"why":"gives the representation of Cauchy stress and Eulerian moduli from which the stress identity is derived","marker":"[49]"},{"why":"provides the equation of motion for plane shear waves in prestressed anisotropic materials","marker":"[50]"}],"fun_headline_variants":["Shear wave pair reads stress without material laws","Stress from two wave speeds, no constants needed","Two shear speeds map stress in soft matter","Parameter-free stress imaging via shear waves"],"cache_read_input_tokens":18432,"weakest_assumption_plain":"The whole method collapses if the material's structural anisotropy, its initial stress, and its deformation are not aligned with one another, because then the two measured waves no longer sit exactly on the principal stress axes and equation (2) no longer follows.","fun_headline_variants_meta":{"raw":{"variants":["Shear wave pair reads stress without material laws","Stress from two wave speeds, no constants needed","Two shear speeds map stress in soft matter","Parameter-free stress imaging via shear waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000724,"raw_usage":{"total_tokens":3226,"prompt_tokens":907,"completion_tokens":2319,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":523,"completion_tokens_details":{"reasoning_tokens":2263}},"tokens_in":523,"tokens_out":2319,"duration_ms":17087,"temperature":1.0,"reasoning_tokens":2263,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:54:42.976464+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a transversely isotropic phantom whose fibres are deliberately inclined at a known angle to an applied uniaxial stress, measure $v_x$ and $v_z$ along the transducer axes, and compare the value $\\rho(v_x^2-v_z^2)/\\cos 2\\theta_0$ with the known applied principal stress difference. If the two agree within measurement error even when fibres and stress are misaligned, the coaxiality premise is unnecessary; if they diverge systematically with misalignment, the premise is load-bearing and the method's range is limited to aligned tissues.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"establishes the isotropic version of the stress-wave-speed result that this paper generalizes to coaxial anisotropy"},{"cited_title":"Shams, M","cited_arxiv_id":null,"evidence_quote":"provides the theory of incremental waves in initially stressed incompressible solids"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the elliptical wavefront and group-versus-phase speed analysis used to justify measuring group speeds"},{"cited_title":"Rouze, M","cited_arxiv_id":null,"evidence_quote":"provides the Radon-transform method used to extract shear wave speeds from spatiotemporal data"},{"cited_title":"Berjamin, R","cited_arxiv_id":null,"evidence_quote":"supplies the acousto-visco-elastic model used to quantify the bias introduced by tissue viscosity"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the representation of Cauchy stress and Eulerian moduli from which the stress identity is derived"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the equation of motion for plane shear waves in prestressed anisotropic materials"}],"review_version":1}