{"id":"3f98834b-cafd-440f-9568-982635a61561","arxiv_id":"2608.12200","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"The paper claims a focused ultrasound model yields 90 C tumor heating and a 0.75 beam-radius lesion, but those values are imposed inputs, not predictions.","lead":"This paper presents a mathematical model claiming that a one-second focused ultrasound pulse can heat a tumor core to 90 C while leaving surrounding tissue unharmed. The key numbers, however, are built into the model's assumptions rather than derived from independent physics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Even granting v2=0, the post-pulse containment proof is internally inconsistent: Eq. (36)'s Gaussian initial condition contradicts Eq. (33)'s parabolic end-of-pulse profile, and the derivative test at r_b does not bound injury outside r_b.","rationale":"I read the manuscript in good faith and attempted to evaluate the strongest version of its central argument: that a structurally fixed tumor suppresses second-order streaming, that the absorbed acoustic energy maps directly to a localized parabolic temperature field, and that the Arrhenius integral plus a Green's-function cooling analysis yield a sharp, contained lesion at 0.75w0. The kinematic suppression of streaming is indeed imposed rather than derived, as the Reader noted. However, I find an even more immediately decisive internal inconsistency downstream: the cooling proof does not use the temperature field produced by the heating phase. Eq. (36) replaces the parabolic profile of Eq. (33) with a different Gaussian, so the post-pulse dynamics are not the dynamics of the model's own lesion. At r_b the Gaussian initial condition is 67.2°C rather than 60°C, so the statement that 'the instant the pulse terminates' the boundary is below 60°C is false. The monotonicity argument only concerns that single radius, and a standard diffusing Gaussian initially increases outside the waist; no radial scan of T or Ω_cooling is provided. This is not a matter of disagreement with clinical consensus; it is a quantitative claim that can be checked by a 2D heat-conduction run. Because this flaw is internal and directly attacks the headline containment/localization result, it reinforces the REJECT verdict. If the cooling test were to pass, the v2=0 assumption and the circular Arrhenius-boundary construction would still need separate justification, so I see no basis to upgrade the verdict.","tokens_in":12911,"tokens_out":12860,"duration_ms":121234,"concrete_test":"Solve Eq. (35) for t′>0 with the actual end-of-pulse condition T(r,τ)=37+53(1−r^2/w0^2) for 0≤r≤w0 and T=37°C for r>w0 (or the extension implied by Eq. 33). Record max_{t′>0} T(r,t′) at r=0.8w0, 0.9w0, 1.0w0, and 1.2w0 and evaluate Ω_cooling(r) from Eq. (40) using Arrhenius constants consistent with a 60°C/1s threshold. If any monitored radius crosses 60°C after pulse cutoff or yields Ω_cooling(r)≥1, the containment claim and the >90% localization interpretation fail.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Granting the v2=0 kinematic assumption for the sake of argument, the central containment claim still does not follow from the paper's own equations. The cooling phase is initialized in Eq. (36) with θ(r,0)=ΔT_max exp(−r^2/w0^2), but the end-of-pulse temperature from Eq. (33) is θ(r,τ)=ΔT_max(1−r^2/w0^2). These differ everywhere except at r=0. In particular, at the claimed lesion boundary r_b=0.75w0, Eq. (34) fixes T=60°C, while Eq. (36) gives T=37+53 exp(−0.5625)≈67.2°C. So the 'immediate monotonic decay below 60°C' starts from a temperature that is already 7°C above the cytotoxic threshold; the cooling calculation is not the continuation of the heating problem. Moreover, the derivative test in Eqs. (41)–(42) is evaluated only at r=r_b. For a diffusing Gaussian, ∂θ/∂t′>0 for r>w0, so tissue just outside the beam initially heats up during cooling. The paper never evaluates the cooling injury integral Ω_cooling(r) for r>r_b, and the parabolic initial condition would produce a different diffusive spreading. The conclusion of 'complete structural containment' and the associated >90% localization rate is therefore not established by this analysis.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents an analytical 2D hydrothermodynamic model of HIFU ablation. It argues that in a dense anchored tumor, second-order acoustic streaming is exactly zero, so absorbed acoustic momentum becomes static pressure gradients and heat; that with a spherical focusing geometry and alpha=1/(2x0), a 1 s top-hat pulse produces a parabolic temperature profile; that the Arrhenius injury integral gives a lesion radius r_b=0.75w0; and that post-pulse Green's-function cooling proves thermal containment, with a volume-averaged tumor temperature of 72.1 C and a focal peak of 90 C. The paper is transparent about many of its assumptions and provides closed-form expressions throughout.","tokens_in":13243,"tokens_out":5468,"duration_ms":49345,"significance":"If the zero-streaming premise and the thermal predictions were established, the model would offer a simple analytical design rule for HIFU frequency-versus-depth and a sharp lesion criterion, which could be useful for device optimization. The paper is clearly written, and the closed-form expressions are easy to check. However, the central quantitative claims are not derived from the model equations in the way claimed: the lesion radius, the average temperature, and the containment proof follow from imposed input values or from an initial condition inconsistent with the heating profile. The potential significance is therefore not realized by the present manuscript.","major_comments":[{"comment":"The zero-streaming constraint v2=0 is asserted from the statement that 'structural biological boundaries are fixed and anchored'; it is not derived from the tissue mechanics or from any equation. This is load-bearing because Eq. (19) is simplified by dropping the convective term rho0 Cp v2 dot grad T, and the entire localization argument relies on that drop. The manuscript provides no estimate of a residual second-order velocity or the timescale on which v2 relaxes; without that, the claim of complete suppression is an assumption rather than a result.","section":"II, Eqs. (16)-(18)"},{"comment":"The lesion radius r_b=0.75w0 is not obtained by solving the Arrhenius injury integral. The text sets Omega(r_b)=1, but then substitutes T(r_b,tau)=60 C and T(0,tau)=90 C into the parabolic profile Eq. (33) and solves algebraically. The Arrhenius expression (29)-(32) plays no role in determining r_b beyond supplying the fixed 60 C threshold, and the 90 C peak is an assumed input chosen to match clinical practice. Thus Eq. (34) is a rearrangement of assumed boundary values, not a predicted lesion boundary.","section":"III, Eq. (34)"},{"comment":"The post-pulse problem is initialized with the Gaussian theta(r,0)=DeltaT_max exp(-r^2/w0^2), but the end-of-pulse profile from Eq. (33) is the parabola theta(r,tau)=DeltaT_max(1-r^2/w0^2). These agree only at r=0; at r_b=0.75w0 the Gaussian gives 37+53 exp(-0.5625)=67.2 C, not the 60 C that defines r_b. Consequently Eqs. (41)-(42) describe cooling of a different initial state, and the monotone decay at r_b does not imply containment for the actual parabolic profile. Furthermore, the derivative test only at r=r_b does not bound Omega_cooling(r) for r>r_b, where a diffusing Gaussian can initially heat; the paper never evaluates the cooling injury integral outside r_b.","section":"III, Eqs. (33), (36), (41)-(42)"},{"comment":"The reported average 72.1 C is computed from the same assumed DeltaT_max=53 C and gamma=0.75, so it is not an independent output. It also uses a 3D spherical volume average in a model that is otherwise presented as 2D and uses a 2D Green's function (Eq. (37)). Had the average been taken over a 2D disk, the prefactor would be different and the value would change (roughly 75.1 C for the same parabolic profile). The paper should either treat the problem consistently in 2D or provide a 3D derivation of the averaged quantity.","section":"III, Eqs. (43)-(46)"},{"comment":"The cooling injury integral Omega_cooling(r) is not evaluated. The text asserts that it 'drops to near zero within milliseconds' based on the derivative sign at r_b. This is not a calculation; the Arrhenius integrand is exponentially sensitive, and the paper supplies no bound on the integral over t' for r>=r_b. Without such a bound, the central claim of complete structural containment is unproven.","section":"III, Eq. (40)"}],"minor_comments":[{"comment":"The normalization of the intensity I0 is not defined in relation to the transducer aperture radius a and focal radius R_f; please specify exactly how I0 in Eq. (15) and Eq. (26) is obtained from the source parameters in Eq. (25).","section":"II, Eqs. (15), (25)-(26)"},{"comment":"The angular brackets in the average temperature label are rendered as 'T V' rather than as an explicit averaged quantity; please fix the notation for clarity.","section":"III, Fig. 1 caption"},{"comment":"The statement says that 'data and numerical codes' are available, but no code or repository is described anywhere in the paper; please either provide the code or revise the statement to match the content.","section":"Data Availability Statement"},{"comment":"The 'multi-scale asymptotic expansion' leading to Eq. (32) is stated without an error estimate; since the claim is that the Arrhenius integral is solved analytically, please quantify the neglected terms or cite a standard asymptotic result.","section":"III, Eq. (32)"}],"recommendation":"reject","confidential_remarks":"The central quantitative outputs of the paper—r_b=0.75w0, 72.1 C, and the containment proof—are not derived from the governing equations in a self-consistent way; they follow from assumed input values or from a cooling initial condition inconsistent with the heating profile. These are load-bearing issues, not presentation problems, so I do not see a path to revision within the scope of the manuscript. The paper also relies on a web database entry (reference [2]) for clinical localization rates; if a revised version is considered, the editor may wish to verify that source and request quantitative clinical data."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper's quantitative results are imposed rather than derived, and the post-pulse containment proof is internally inconsistent with the paper's own heating equations.\n\nWhat is actually new: not much. The optimization criterion α=1/(2x0) comes from your earlier work, adapted to a stationary tumor. The plane-wave comparison is clear and correct: plane waves overheat the surface, spherical focusing concentrates energy. The model is organized well enough that the assumptions are visible.\n\nThe core problem is that the central numbers are chosen, not calculated. The lesion radius r_b=0.75w0 follows from assuming T(0)=90°C and T(r_b)=60°C in the parabolic profile; the Arrhenius integral is not used to find the boundary. The average 72.1°C is then just the integral of that assumed profile over the assumed boundary.\n\nThe containment proof is worse. The cooling phase starts with a Gaussian initial condition (Eq. 36) that does not match the parabolic profile at the end of the pulse (Eq. 33). At r_b, the Gaussian gives about 67°C, not 60°C, and the derivative test at r_b only shows local decay. For r>w0, the temperature initially rises during diffusion, so the claim that tissue outside r_b never exceeds threshold is not established. This is a load-bearing internal inconsistency.\n\nEven granting the v2=0 kinematic suppression, the thermal reasoning does not deliver what the abstract promises. The v2=0 assumption itself is asserted without a scale analysis. There is also a 2D/3D mismatch: the model is 2D but the volumetric average uses 4πr^2.\n\nThe citation pattern is fine; the paper leans on your own prior work, which is acceptable when the prior work is publicly available.\n\nWho is this for? Possibly a reader interested in analytical HIFU approximations, but as a substantive result it is not reliable. The flaws are specific and potentially fixable: add a parameter table, derive the lesion boundary from an actual Arrhenius integration, and use a consistent initial condition for the cooling phase.\n\nMy recommendation: this should go to peer review, because the issues are concrete and a referee could force the necessary revisions. But in its current form the central claims do not hold up.","headline":"The paper's central numbers are imposed rather than derived, and the post-pulse containment proof is internally inconsistent with the paper's own heating profile.","tokens_in":13786,"tokens_out":4363,"would_cite":false,"duration_ms":37749,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A stationary tumor matrix suppresses acoustic streaming, converting absorbed ultrasound momentum into localized heat with a sharp lesion boundary at 0.75 w0.","keywords":["high-intensity focused ultrasound","acoustic streaming suppression","Pennes bioheat equation","Arrhenius injury integral","thermal ablation","second-order perturbation","spherical focusing","tumor ablation"],"falsifier":"Measure the time-averaged second-order velocity in a dense anchored tumor or tissue-mimicking phantom during a 1 s HIFU pulse. A nonzero $v_2$ with $\\rho_0 C_p v_2\\cdot\\nabla T$ comparable to $2\\alpha I$ would break the pressure-gradient heat conversion. Alternatively, record the radial temperature at $t=1$ s and shortly after: the model's cooling law (Eq. 39) predicts a monotone drop at $r_b=0.75\\,w_0$, so any measured temperature rise outside $r_b$ after pulse cutoff would falsify the containment claim.","tokens_in":12633,"feed_emoji":"🔥","tokens_out":6473,"duration_ms":52349,"temperature":0.7,"pith_summary":"The paper claims that in a dense, structurally anchored tumor, acoustic streaming is fully suppressed, so the momentum lost by a focused ultrasound wave cannot drive fluid motion; it is converted into static pressure gradients and, through viscous absorption, into heat at the same spot. Building on that constraint, the paper solves a simplified bioheat equation for a one-second top-hat pulse in a spherically focusing geometry, and derives a closed-form lesion boundary at $r_b=0.75\\,w_0$, a volumetric mean temperature of about $72.1\\,^\\circ$C with a $90\\,^\\circ$C central peak, and monotone post-pulse cooling below the $60\\,^\\circ$C necrosis threshold. The reason this matters is that it turns the observed high localization of HIFU ablation into a mechanistic statement: anchored tissue cannot advect heat away, so the absorbed wave energy stays where it was absorbed.","feed_headline":"Stationary tumors turn ultrasound momentum into heat, not streaming","feed_subtitle":"A one-second focused pulse gives a lesion edge at 75% of the beam waist, 72.1 °C average, and cooling below 60 °C.","key_machinery":"The carrying mechanism is the second-order expansion of the compressible viscous Navier-Stokes equations with the zero-streaming constraint $v_2=0$, imposed because the tumor matrix is treated as fixed and anchored. This converts the acoustic force density $F_2=2\\alpha I(r)/c_0\\,\\hat{r}$ into a static pressure gradient and removes the convective term $\\rho_0 C_p v_2\\cdot\\nabla T$ from the energy balance. The derived optimum $\\alpha=1/(2x_0)$, the spherical focusing factor $(a/r)^2$, the parabolic focal temperature profile $T(r,\\tau)=T_0+\\Delta T_{\\max}(1-r^2/w_0^2)$, the Arrhenius injury integral, and the 2D free-space Green's function convolution then combine to produce the lesion radius and the containment statement.","core_discovery":"The central claim is that the kinematic constraint $v_2=0$ inside a stationary tumor makes the time-averaged acoustic momentum flux $F_2 = 2\\alpha I(r)/c_0\\,\\hat{r}$ transfer entirely into a static second-order pressure gradient $\\nabla\\langle p_2\\rangle$, so the absorbed energy appears as heat rather than as acoustic wind. With that constraint, the Pennes bioheat equation reduces to a local, non-diffusive source on the pulse timescale; a spherical focusing geometry with geometric convergence $\\propto 1/r^2$ overrides exponential attenuation within a critical radius $r_{\\mathrm{crit}}=2x_0$; and the non-isothermal Arrhenius integral gives $r_b=0.75\\,w_0$. The stated thermal outcome is a peak focal temperature of $90^\\circ$C, a volumetric average of $72.1^\\circ$C inside the lesion, and post-pulse Green's-function cooling that stays monotonically below $60^\\circ$C at the boundary.","pith_inferences":["If streaming suppression is real, tumors with stiffer, more collagenous extracellular matrices should exhibit sharper thermal lesions than soft or fluid-rich tumors; comparing lesion sharpness across tumor stiffness would test this directly.","The $\\alpha=1/(2x_0)$ scaling suggests a simple patient-specific tuning rule: choose the HIFU frequency from target depth alone, then check whether the measured lesion radius matches $0.75\\,w_0$ across a range of depths.","The cooling phase replaces the parabolic heating profile by a Gaussian initial condition; whether real focal profiles are closer to Eq. (33) or Eq. (36) will determine how literally the $0.75\\,w_0$ boundary and the monotone cooling statement should be read."],"forward_implications":["Focal lesions in anchored tumors should have a sharp necrotic boundary at $r_b=0.75\\,w_0$ rather than a gradual damage gradient.","Optimal operating frequency for a given patient scales inversely with target depth through $\\alpha = 1/(2x_0)$.","Plane-wave HIFU is inherently limited because upstream tissue receives about $e$ times more heat than the tumor core, so spherical focusing is not optional.","After pulse termination the boundary temperature falls monotonically below $60^\\circ$C, which is why residual conduction injury is negligible on this timescale."],"supporting_citations":[{"why":"Frames the clinical HIFU ablation setting whose high localization rates the model aims to explain.","marker":"[1]"},{"why":"Reports the device operational parameters and the >90% localization benchmark used as the clinical reference.","marker":"[2]"},{"why":"Supplies the clinical intensity and focal pressure ranges underlying the high-amplitude pulse parameters.","marker":"[8]"},{"why":"Supplies the earlier hydrodynamic optimization framework and criterion alpha=1/(2x0) adapted here to stationary tumor domains.","marker":"[11]"},{"why":"Provides the tissue density and bulk-modulus contrast ranges for tumor versus healthy tissue.","marker":"[12]"},{"why":"Gives the bulk-modulus relation for sound speed used to define the tumor's acoustic refraction regimes.","marker":"[13]"}],"fun_headline_variants":["Tumor's stillness turns ultrasound into heat, not wind","Focused ultrasound hits 90°C in tumor core, spares boundary","One-second pulse yields sharp lesion edge at 0.75 beam width","Acoustic streaming suppressed: all momentum becomes localized heat","Spherical focusing beats attenuation: 72°C average tumor temp"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument rests on assuming, rather than deriving, that the tumor matrix is rigid enough to make the second-order streaming velocity exactly zero; if even a small $v_2$ exists, convective heat transport reappears and the claimed sharp containment no longer follows.","fun_headline_variants_meta":{"raw":{"variants":["Tumor's stillness turns ultrasound into heat, not wind","Focused ultrasound hits 90°C in tumor core, spares boundary","One-second pulse yields sharp lesion edge at 0.75 beam width","Acoustic streaming suppressed: all momentum becomes localized heat","Spherical focusing beats attenuation: 72°C average tumor temp"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00101,"raw_usage":{"total_tokens":4357,"prompt_tokens":1121,"completion_tokens":3236,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":737,"completion_tokens_details":{"reasoning_tokens":3147}},"tokens_in":737,"tokens_out":3236,"duration_ms":20267,"temperature":1.0,"reasoning_tokens":3147,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:13:27.835695+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the time-averaged second-order velocity in a dense anchored tumor or tissue-mimicking phantom during a 1 s HIFU pulse. A nonzero $v_2$ with $\\rho_0 C_p v_2\\cdot\\nabla T$ comparable to $2\\alpha I$ would break the pressure-gradient heat conversion. Alternatively, record the radial temperature at $t=1$ s and shortly after: the model's cooling law (Eq. 39) predicts a monotone drop at $r_b=0.75\\,w_0$, so any measured temperature rise outside $r_b$ after pulse cutoff would falsify the containment claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the device operational parameters and the >90% localization benchmark used as the clinical reference."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the clinical intensity and focal pressure ranges underlying the high-amplitude pulse parameters."},{"cited_title":"Unpinning of trapped oil droplets via non-resonant acoustic streaming in capillary tubes","cited_arxiv_id":"2606.27331","evidence_quote":"Supplies the earlier hydrodynamic optimization framework and criterion alpha=1/(2x0) adapted here to stationary tumor domains."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the tissue density and bulk-modulus contrast ranges for tumor versus healthy tissue."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the bulk-modulus relation for sound speed used to define the tumor's acoustic refraction regimes."}],"review_version":1}