{"id":"edce19a1-07ae-4445-a257-28e36b0642ee","arxiv_id":"2512.17575","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A multipole-expansion model predicts the frequency response and resonance shifts of a bubble near a spherical inclusion of arbitrary size and mechanical nature.","lead":"This paper derives a linear acoustics model for how a gas microbubble's resonance changes when it sits near a spherical particle that is rigid, fluid-like, or viscoelastic. It could enable a new form of micro-elastography where bubbles are used to probe the stiffness of cells.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model excludes direct insonification of the spherical inclusion; for viscoelastic cells this can alter the predicted resonance fingerprints and undermines the inverse-elastography claim.","rationale":"The reader's weakest-assumption identification is directly supported by the manuscript's own statement about the air sphere and by the absence of an incident-field term in the sphere boundary conditions. This is the most load-bearing concern because the paper's novel application claim is bubble-mediated mechanosensing and inversion, which requires the model to describe the actual experimental situation. If the incident wave directly excites the sphere, the predicted frequency response (especially the cell-resonance coupling in Fig. 10b and the fingerprints in Figs. 11–12) could be materially different. I agree with the reader's CONDITIONAL verdict: the analytical framework and its limiting-case validations are solid, but the experimental relevance of the excitation assumption must be demonstrated before the inverse-elastography claim is accepted. A computational test that adds the incident wave to the sphere boundary conditions, or a full-wave simulation, would settle the matter. I do not see a stronger internal inconsistency in the derivations, and the concern is not about disagreement with consensus but about a specific modeling assumption that may fail in the target application.","tokens_in":26934,"tokens_out":5310,"duration_ms":61944,"concrete_test":"Re-solve the linear system for the viscoelastic-sphere case with an incident plane-wave field included in the boundary conditions at the sphere (and consistently at the bubble), or equivalently run a full FEM/BEM simulation of a 10-μm radius bubble at h = 1.5R10 from CMM spheres of R20 = 4R10 and 20R10, driven by a plane wave sweeping f/f0 from 0.2 to 1.4. Compare the predicted radial-mode frequency response, resonance shift, and the secondary CMM resonance feature in Fig. 10(b) with the present model. If the shift changes by more than about 10%, the current fingerprints are configuration-specific and need revision before inverse use.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the analytical model predicts the frequency response of a bubble near a rigid/fluid/viscoelastic sphere and that scanning this response can recover inclusion properties. A necessary condition is that the model's excitation matches the measurement protocol. The model applies the incident pressure only at the bubble surface; the sphere is driven exclusively by the bubble's scattered field. This is explicit for the air sphere in Sec. III B 2 ('the air sphere is assumed to oscillate only in response to the acoustic field radiated by the bubble, without being influenced by the external driving wave') and is implicit in the sphere boundary conditions (Eqs. B25–B28, B49–B52), where no incident-field term appears. In a real elastography experiment, a traveling or standing ultrasound wave insonifies both objects. For a small viscoelastic cell (e.g., CMM, R20 = 4R10), the incident wave can directly drive internal compressional and shear resonances; the paper itself finds cell resonances near the bubble resonance (Fig. 10b). If the inclusion is directly driven, the bubble's response gains additional scattering contributions and the resonance shifts and quality-factor fingerprints in Figs. 11–12 will change. Thus the inverse-characterization claim depends on an untested assumption about the experimental configuration.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents an analytical model for the linear frequency response of a gas microbubble near a spherical inclusion (rigid, fluid, or viscoelastic) in a viscous compressible liquid. The model is derived from the linearized equations of motion via Helmholtz decomposition, with boundary conditions expressed through multipole expansions and coordinate transformations between the bubble and sphere centers. Both radial and shape modes are considered. The authors validate the model against the linearized Rayleigh–Plesset equation, the Lamb shape-mode frequencies, and the planar-wall limits of Strasberg and Hay. They then present forward simulations showing how resonance frequencies and quality factors change with sphere size, distance, and material, and propose scanning this response as a basis for inverse mechanical characterization of microscopic objects such as biological cells.","tokens_in":27271,"tokens_out":4213,"duration_ms":51384,"significance":"If the central claims hold, the paper provides a unified, parameter-free analytical framework for bubble–sphere interaction that goes beyond existing planar-wall models and handles curvature, viscosity, compressibility, and viscoelasticity consistently. The careful benchmarking against established limits is a genuine strength, and the observation that cell-mimicking inclusions can exhibit resonances near the bubble resonance (Fig. 10b) is physically interesting. However, the inverse-elastography claim is not yet demonstrated: the paper shows forward 'fingerprints' but no actual inversion, no noise analysis, and no sensitivity or identifiability study. The application-oriented part of the abstract and conclusion is therefore stronger than the evidence presented.","major_comments":[{"comment":"The central application claim is that scanning the bubble frequency response allows the inclusion's mechanical properties to be recovered 'through inverse modeling' (Abstract, Conclusion). The paper, however, only presents forward maps: resonance frequency and quality factor as functions of material parameters (Figs. 11–12). There is no inversion of synthetic or experimental data, no study of noise sensitivity, and no discussion of whether the mapping is one-to-one (e.g., glycerin vs PMMA appear nearly degenerate in Fig. 11, and the authors themselves note ambiguity that they propose to resolve with an additional dimension). As written, the inverse-elastography claim is speculative. I would ask the authors to either provide a proof-of-principle inversion (e.g., recover known material parameters from synthetic noisy response curves) or explicitly reframe the conclusion as 'forward modelin","section":"Section III C and Conclusion"},{"comment":"The model assumes the spherical inclusion is excited only by the wave scattered by the bubble, not by the external driving pressure. This is explicit for the air sphere in Sec. III B 2 ('the air sphere is assumed to oscillate only in response to the acoustic field radiated by the bubble, without being influenced by the external driving wave') and is implicit in the sphere boundary conditions, where no incident-field term appears. In a real elastography experiment the incident ultrasound wave insonifies both the bubble and the inclusion. For a viscoelastic cell (e.g., CMM, Fig. 10b) with resonances near the bubble resonance, direct driving of the cell would add scattering contributions and could shift or alter the very fingerprints used for characterization. This is not a flaw of the model as a mathematical derivation, but it is a load-bearing assumption for the proposed measurement proto","section":"Sec. III B 2; Eqs. (B25)–(B28), (B49)–(B52)"},{"comment":"Shape-mode resonance frequencies are not obtained as poles of the unforced system but by forcing the bubble with Pac Σ Pn(µ1) and reading the peak of the modal amplitude sn. The text explicitly acknowledges (Sec. II B 3) that 'the proposed method does not allow the characterization of the resonance frequency of these modes' before introducing the special forcing. For a damped, driven oscillator the peak response frequency differs from the undamped natural frequency, and in a coupled multiple-scattering system the peak of a particular sn can be affected by neighboring modes. Thus the quantities reported in Figs. 3, 5, 7(c), 9(c) as 'resonance frequencies' are operationally defined as forced-response maxima, not derived natural frequencies. This does not invalidate the forward predictions if that is the intended definition, but the distinction should be stated clearly and the term 'natural","section":"Sec. II B 3 and Sec. III (preamble)"}],"minor_comments":[{"comment":"The symbol δ is used both for the viscous penetration depth in Eq. (10) and for the total damping coefficient in the Rayleigh–Plesset comparison Eq. (57). These are different quantities; please rename one to avoid confusion.","section":"Eqs. (10) and (57)"},{"comment":"The statement that 'a threshold value for the sphere radius exists ... above which convergence is not reached' is too vague. Which physical and numerical parameters determine the threshold, and how does the truncation N_t affect the reported results? Please provide a short convergence analysis or at least report the N_t values used, so the reader can assess the reliability of the large-sphere results.","section":"Sec. III B (preamble)"},{"comment":"The word 'insert' in several captions (e.g., Fig. 5(a), Fig. 10) should be 'inset'. Also, in Fig. 10(b) the inset is described as 'the particle's response' but the units and normalization of u_0 are not given; please add.","section":"Figures 5, 7, 9, 12"},{"comment":"The sentence 'as the radius of the rigid sphere increases, the resonance frequency curve progressively tends to that of the infinite wall' is followed by the observation that even R20 = 40R10 differs noticeably from the plane-wall case. The wording is confusing; please distinguish 'tends in the limit' from 'is still measurably different at the largest computed radius.'","section":"Sec. III B 1"}],"recommendation":"major_revision","confidential_remarks":"I found no evidence that the derivation is circular or that the authors deliberately fit to the target results; the core linear model is carefully derived and well benchmarked. The main gap is that the headline application (inverse characterization) is not actually demonstrated and rests on an untested assumption about the experimental illumination. These are fixable within the manuscript's scope, so I do not recommend rejection, but the claims need to be scaled back or substantiated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is a unified analytical model for a bubble oscillating near a finite spherical inclusion—rigid, fluid, or viscoelastic—in a viscous compressible liquid, including both radial and shape modes. The authors extend the established multipole program (Strasberg, Doinikov, Hay) to curved interfaces of arbitrary size and material, which nobody had done in this generality. The derivation is careful, and the checks against Rayleigh–Plesset, Lamb, Strasberg, and Hay in the appropriate limits give real confidence that the core model is right. That is the paper’s strength, and it is solid.\n\nThe soft spot is the inverse-elastography claim. The stress-test flags it correctly: the sphere is driven only by the bubble’s scattered field, with the incident pressure applied only at the bubble surface. The paper states this explicitly for the air sphere (Sec. III B 2) and the boundary conditions contain no incident-field term at the sphere. In any real experiment the ultrasound pulse hits the cell directly, and for the CMM case the cell has resonances near the bubble resonance—so direct driving would alter the very fingerprints Figs. 11–12 are meant to characterize. That is not a minor detail for the application; it is a missing piece of the measurement model. Relatedly, the abstract says the approach “demonstrates” inverse modeling, but there is no inversion, no noise analysis, no sensitivity study—the conclusion more honestly says “offers a way” and “opening the way.” The word choice in the abstract oversells.\n\nMinor issues: the convergence analysis is mentioned but no numbers are given (truncation order, tolerance), and the claim that resonance frequencies are found from the maximum of the amplitude response is fine but should be stated with the caveat that this is an estimate for a damped system. Neither changes the core finding.\n\nThe model itself is a legitimate contribution to the bubble-acoustics literature and deserves a serious referee. The application section needs real revision: either include direct insonification in the model, or clearly scope the claim to configurations where the sphere is shielded from the incident wave and justify that experimentally. I would not desk-reject this; I would send it out with an expectation of major revision, mostly targeted at the interpretation and the missing inversion analysis.","headline":"Solid multipole model for bubble-sphere resonance, but the mechanosensing application is overclaimed: the inclusion is never directly insonified and no inversion is shown.","tokens_in":27674,"tokens_out":1558,"would_cite":true,"duration_ms":18325,"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":"This paper claims that a single analytical model, built from the linearized equations of a viscous compressible liquid, predicts how a gas microbubble's resonance frequencies and amplitudes shift when it sits near a rigid, fluid, or viscoel","keywords":["microbubble","resonance frequency","spherical inclusion","viscoelastic","multipole expansion","elastography","viscous compressible fluid"],"falsifier":"Measure the frequency response of a 10-micron bubble near a well-characterized viscoelastic microsphere whose stiffness and viscosity are known independently, while the incident ultrasound insonifies both objects, and check whether the bubble's resonance-frequency shift and quality factor match the model's predictions and whether the sphere's own resonances appear at the predicted peaks. A mismatch in the sphere-resonance peaks would indicate the passive-inclusion assumption fails.","tokens_in":26849,"feed_emoji":"🫧","tokens_out":3219,"duration_ms":35038,"temperature":0.7,"pith_summary":"The paper aims to establish that the linear frequency response of a gas microbubble near a spherical inclusion of arbitrary size and mechanical nature is governed by one analytical model based on the first-order equations of a viscous compressible liquid. It predicts how the resonance frequencies and amplitudes of radial and shape modes change with bubble size, sphere radius, material properties, and distance. The central result is that the bubble's spectral response acts as an acoustic fingerprint of the neighboring object, so scanning it offers a route to recover the inclusion's mechanical properties by inverse modeling. A sympathetic reader would care because this could turn a standard ultrasound contrast agent into a local probe for tissue and cell stiffness at the microscale.","feed_headline":"Bubble resonances fingerprint a nearby sphere's material","feed_subtitle":"A new analytical model links a microbubble's frequency shifts to a nearby sphere's material properties.","key_machinery":"The central object is the coupled multipole expansion of the velocity potentials in the two spherical coordinate systems, with the bubble's scattered wave re-expressed in the sphere's coordinates and vice versa using translation formulas involving Clebsch-Gordan coefficients. This lets the authors assemble a linear system of boundary conditions—normal velocity and tangential stress at the bubble surface, and velocity/stress continuity at the sphere—whose solution yields the scattering coefficients of every mode and hence the bubble's modal amplitudes as functions of frequency and distance.","core_discovery":"The paper claims to derive, from linearized mass and momentum conservation in a viscous compressible liquid, a complete analytical description of a small gas bubble oscillating near a sphere of arbitrary size and mechanical nature. The model accounts for both the bubble's radial breathing mode and its nonspherical surface modes, and it treats the sphere as rigid, as a viscous compressible fluid, or as a viscoelastic solid by matching velocities and stresses at the sphere surface. It reproduces known limits far from the sphere and near planar walls, and it predicts that the radial resonance frequency shifts down as a rigid sphere approaches, shifts up near a compliant air sphere, and shows a","pith_inferences":["If the model's scanning claim holds, a natural extension is to invert not just the resonance peak but the full measured frequency response to extract multiple mechanical parameters (Young modulus and viscosity) simultaneously, since each parameter shifts amplitude, peak position, and Q differently.","The model's assumption that the sphere is only driven by the bubble's scattered field suggests that in experiments with free-floating cells the incident ultrasound may directly excite the cell, so a testable extension would deliberately compare this passive-inclusion model with an insonified-inclusion version to bound the error.","The same multipole machinery could likely be extended to multiple neighboring spheres or to nonspherical inclusions by summing translation contributions, permitting bubble-based mapping of heterogeneous microstructures rather than a single local object."],"forward_implications":["A single frequency sweep of a bubble can in principle distinguish among air, rigid, and viscoelastic neighboring objects through changes in resonance frequency, amplitude, and quality factor.","Near a soft viscoelastic sphere of comparable size, resonance peaks of the sphere itself can appear in the bubble's response, so the bubble can detect internal acoustic resonances of a nearby object.","For sufficiently large spheres, the model's predictions reduce to known planar-boundary results, meaning curvature effects matter mainly for small inclusions comparable to the bubble size.","Shape-mode resonance shifts are generally small (a few percent and less) and appear only at very short distances, so radial-mode information is the more robust sensor channel.","The (resonance-frequency shift, quality factor, distance) space provides material-specific signatures that can support inversion of mechanical properties."],"fun_headline_variants":["Microbubble resonance reads sphere's viscoelasticity","Bubble frequency shifts reveal adjacent sphere's stiffness","Model maps bubble's resonance to spherical inclusion properties","Resonant bubble near a sphere probes material traits","Bubble resonances fingerprint nearby sphere material"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The nearby sphere is assumed to move only because of the wave the bubble radiates, not because of the external ultrasound wave that drives the bubble, which may not hold in real experiments where the incident wave hits the sphere directly.","fun_headline_variants_meta":{"raw":{"variants":["Microbubble resonance reads sphere's viscoelasticity","Bubble frequency shifts reveal adjacent sphere's stiffness","Model maps bubble's resonance to spherical inclusion properties","Resonant bubble near a sphere probes material traits","Bubble resonances fingerprint nearby sphere material"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000865,"raw_usage":{"total_tokens":3527,"prompt_tokens":628,"completion_tokens":2899,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":372,"completion_tokens_details":{"reasoning_tokens":2826}},"tokens_in":372,"tokens_out":2899,"duration_ms":24162,"temperature":1.0,"reasoning_tokens":2826,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T15:12:55.233244+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the frequency response of a 10-micron bubble near a well-characterized viscoelastic microsphere whose stiffness and viscosity are known independently, while the incident ultrasound insonifies both objects, and check whether the bubble's resonance-frequency shift and quality factor match the model's predictions and whether the sphere's own resonances appear at the predicted peaks. A mismatch in the sphere-resonance peaks would indicate the passive-inclusion assumption fails.","supporting_citations":[],"review_version":1}