{"id":"049b8430-1d13-4a28-af74-c96828e021ff","arxiv_id":"1908.07663","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Transient ring-up of ultrasound contrast microbubbles lowers the predicted attenuation peak near resonance and changes the shell viscosity fitted from short-pulse data by a factor of about three.","lead":"This paper derives a new formula for how much sound a coated microbubble absorbs while it is still starting to ring up, before steady vibration begins. It shows that very short ultrasound pulses, common in contrast-agent measurements, can change the inferred shell viscosity by about a factor of three.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (15) uses an averaging window N = 2ln10/(πδ) that is tied to the fitted damping, not to the measured pulse duration; for broadband pulse-echo data the steady-state transfer function governs attenuation, so the 2/3 peak reduction may not apply to the experiments used.","rationale":"The reader's conditional verdict is appropriate, and the weakest assumption identified by the reader is the same one that this stress test finds most load-bearing: the measurement model behind Eq. (15). The paper explicitly applies its transient attenuation formula to broadband short-pulse attenuation data, but such data are conventionally analyzed in the frequency domain, where the steady-state transfer function already includes all transient effects through the causal impulse response. The internal inconsistency in the definition of N — derived from the fitted damping constant rather than from the actual transducer pulse duration — strengthens this concern and provides a concrete handle for testing it. Within the idealized monochromatic-tone-burst model, the oscillator derivation is plausible, and the paper's energy-balance logic is internally coherent. However, the magnitude of the claimed correction and the factor-of-three change in inferred shell viscosity depend on averaging over a window that is not the measured pulse length. The proposed simulation is a single, decisive check: it directly compares the spectral-ratio attenuation coefficient of a dilute UCA suspension with Medwin's formula and with Eq. (15). If the simulation confirms Medwin's steady-state result, the central correction does not apply to broadband attenuation measurements, and the paper would need to be reframed as a model of energy loss from a monochromatic tone burst, not as a correction to standard UCA characterization. The lack of comparison with the previously published Clarke-Leighton transient cross-section formula is an additional weakness but is secondary to the measurement-model issue. I therefore retain the reader's CONDITIONAL verdict rather than moving to REJECT, because the tone-burst version of the formula may still be correct and testable.","tokens_in":6273,"tokens_out":8056,"duration_ms":154833,"concrete_test":"Simulate a 2.25 MHz one-cycle Gaussian-enveloped pulse propagating through a dilute UCA suspension governed by the linearized RPE coupled to the wave equation, using the fitted parameters of Table 1, and estimate α(f) by the standard spectral-ratio method. If the resulting α(f) matches Medwin's steady-state Eq. (13) rather than Eq. (15), the transient suppression is an artifact of the non-physical averaging window and the Table 1 refit is invalid. As a secondary check, re-fit the same data with Eq. (15) but replace N by the actual measured pulse duration in cycles; if the fitted κ_s no longer drops by about a factor of three, the headline result depends on the arbitrary choice of N.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III defines N in Eq. (15) as the number of periods for the transient to decay to 1% (N = 2ln10/(πδ)) and then states that \"the period should be the pulse duration in a pulse-echo system.\" These two prescriptions are incompatible. The Xia et al. 2015 data to which the paper fits were acquired with broadband transducers; a 2.25 MHz transducer has a pulse duration of roughly 0.44 µs, about one cycle, whereas the decay-based N is many cycles for the fitted damping. Eq. (15) therefore time-averages over forced-oscillation cycles that occur after the driving pulse has ended, and the integration window itself depends on the shell dilatational viscosity κ_s being fitted, through δ. More fundamentally, broadband attenuation measurements are obtained from the spectral ratio of received pulses. For a linear time-invariant suspension, the frequency-dependent attenuation is determined by the steady-state transfer function (Medwin or Commander-Prosperetti); the ring-up transient is the time-domain manifestation of that same transfer function and does not supply an additional 1/3 reduction at resonance. If the paper's claim were restricted to true tone-burst energy measurements, Eq. (15) could be a legitimate extension, but the central application to broadband attenuation data and the factor-of-three change in dilatational viscosity in Table 1 depend on this unvalidated measurement model. No comparison with the Clarke-Leighton transient cross-section formula, no error analysis, and no independent validation are provided.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript revisits the standard practice of estimating the attenuation of ultrasound contrast agent (UCA) suspensions from the linearized Rayleigh-Plesset equation, which is usually treated as a driven harmonic oscillator in steady state (Medwin's formula). The author argues that for short-pulse interrogation, the transient (ring-up) part of the oscillator solution is not negligible and that its inclusion lowers the predicted attenuation peak to about two-thirds of the steady-state value and shifts the peak to the damped resonance frequency. A new formula, Eq. (15), is proposed for attenuation including transient contributions, and it is fitted to the broadband attenuation data of Xia et al. (2015). This fitting yields substantially different shell parameters, in particular a dilatational viscosity that is about three times smaller than that obtained with Medwin's formula. The paper concludes that steady-state attenuation formulas can mis-estimate UCA shell properties when short pulses are used.","tokens_in":6523,"tokens_out":5524,"duration_ms":263063,"significance":"If the central claim were fully validated, the paper would have practical value for UCA material characterization: it predicts that shell dilatational viscosities inferred from short-pulse attenuation data could be overestimated by a factor of about three, which would affect subsequent modeling of UCA dynamics. The analytical treatment of the damped oscillator is standard and the paper is clearly written, but the practical significance rests entirely on the validity of Eq. (15) as a model for the measured attenuation. That link is not established in the present manuscript, and the comparison with broadband pulse-echo data is problematic. The paper also does not compare against the existing Clarke-Leighton time-dependent cross-section formulation, which is the established treatment of the same transient effect. As a result, the contribution is presently more of an interesting research note than a fully supported method for shell parameter estimation.","major_comments":[{"comment":"The derivation of Eq. (15) is not verifiable as presented. The text says that substituting R(t)=R0+X(t) into Eq. (12) and neglecting 'all the second order terms' yields Eq. (15), but the displayed equation appears to contain a term proportional to the integral of X(t) cos(ωt) dt, which is itself second order in the acoustic amplitude. Please give a transparent step-by-step derivation from Eq. (12) to Eq. (15), explicitly state which terms are retained and which are neglected, and define every symbol and integration limit in the displayed formula.","section":"Section III, Eq. (15)"},{"comment":"The definition of N in Eq. (15) is self-referential and inconsistent with the stated pulse-echo application. N is defined as the number of periods for the transient to decay to 1%, i.e., N=2ln(10)/(πδ), where δ is the oscillator damping constant that includes the shell dilatational viscosity κ_s being fitted. The immediately following sentence, 'Note that the period should be the pulse duration in a pulse-echo system,' prescribes an independent, measured pulse duration. These two prescriptions are incompatible. For the 2.25 MHz transducer cited in the paper, the pulse duration is roughly 0.44 µs (about one acoustic cycle), whereas the decay-based N for the fitted damping is many cycles. Please specify which N is used in the fits of Table 1 and analyze the sensitivity of the fitted parameters to this choice.","section":"Section III, definition of N"},{"comment":"The central application to the broadband attenuation data of Xia et al. (2015) is not justified. Those data were acquired with broadband transducers and are typically obtained from spectral ratios of received pulses, so for a linear time-invariant suspension the frequency-dependent attenuation is governed by the steady-state transfer function of the bubbly medium. The ring-up transient is the time-domain manifestation of that same transfer function and does not provide an additional one-third reduction in the measured attenuation at resonance. The author must either derive how spectral-ratio attenuation measurements relate to the time-averaged energy absorption of Eq. (15), or restrict the claim to true tone-burst energy measurements and remove or reinterpret the fit to Xia et al. data. As written, Table 1 and Figure 4 do not support the conclusion that the shell dilatational viscosity is a factor of three lower than previously estimated.","section":"Section III and Figure 4"},{"comment":"The factor-of-three change in dilatational viscosity is the main quantitative claim, yet Table 1 reports no uncertainty estimates, no goodness-of-fit measures, and no residuals for either model. Please report confidence intervals (e.g., from a covariance analysis or bootstrap) and show the fit residuals over the frequency range used. In addition, the paper does not compare Eq. (15) with the Clarke-Leighton (2000) time-dependent cross-section formula, which already addresses the same transient effect; such a comparison would clarify what the new formula adds and would provide a sanity check on the magnitude of the predicted reduction.","section":"Table 1 and Section III"}],"minor_comments":[{"comment":"The abstract uses 'shot-pulse ultrasound'; this should read 'short-pulse ultrasound.'","section":"Abstract"},{"comment":"Eq. (1) is written for a free bubble, while the shell enters only later through the viscoelastic model. Please state explicitly which shell rheological equation is linearized to obtain Eq. (2), so that the damping and stiffness parameters in Eqs. (3)-(4) are defined unambiguously.","section":"Section I"},{"comment":"The caption distinguishes '(a) from steady state oscillation' and '(b) from transient oscillation,' but both panels compare the full solution with Medwin's formula; please clarify the intended distinction between the two panels.","section":"Figure 2"},{"comment":"The statement that the difference between the full solution and Eq. (15) is 'less than 8%' should be quantified with respect to the frequency range and the quantity being compared (e.g., peak attenuation or frequency-averaged difference).","section":"Section III, after Figure 3"},{"comment":"In the conclusion, 'alternations of the estimated dilatational viscosities' should read 'alterations of the estimated dilatational viscosities.'","section":"Section IV"},{"comment":"The experimental conditions of the Xia et al. (2015) data are not described here; please state the bubble size distribution, concentration, temperature, and transmit/receive configuration, or explicitly reference the methods section of the prior paper.","section":"Figure 4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reads like an early-stage research note on an interesting idea, but the connection between the proposed time-averaged energy formula and the actual broadband attenuation measurements is the decisive gap. If the author can either provide a convincing derivation of the measurement-model link or reframe the paper as a theoretical prediction for tone-burst energy measurements with new experimental validation, the contribution could become publishable. Without that, the factor-of-three change in Table 1 is not supported. The lack of comparison with Clarke-Leighton and the absence of uncertainty analysis are additional concerns that should be addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Lang Xia's short paper makes a claim that could matter: when you measure attenuation of UCA suspensions with short pulses, the bubble may not have rung up to steady state, so the standard Medwin formula overestimates attenuation near resonance. The derived Eq. (15) is a compact approximation for attenuation including transient effects, and the paper is honest that the underlying oscillator physics is old. Its own Figure 2 shows the peak dropping to about 2/3 and shifting to the damped resonance frequency. The math within the oscillator model is standard, and the approximation error relative to the full solution is less than 8%. That part is fine.\n\nThe issue is the bridge between Eq. (15) and the attenuation data. Eq. (15) integrates over N periods, where N is defined by the decay of the transient to 1% (N = 2ln10/(πδ)). Then the paper says \"the period should be the pulse duration in a pulse-echo system.\" Those two statements are incompatible. N defined via the damping constant is typically many cycles for the fitted parameters, while the 2.25 MHz transducer used in the fitted data has a pulse duration of about one cycle. The integration window itself depends on the very shell viscosity being fitted, which creates a self-consistency problem.\n\nMore fundamentally, broadband attenuation measurements are obtained from the spectral ratio of received pulses. For a linear time-invariant suspension, the frequency-dependent attenuation is governed by the steady-state transfer function (Medwin or Commander-Prosperetti). The ring-up transient is the time-domain manifestation of that same transfer function; it does not supply an additional reduction at resonance once you do the spectral analysis. So the factor-of-three change in dilatational viscosity in Table 1 depends on a measurement model the paper does not defend. If the claim were restricted to true tone-burst energy measurements, Eq. (15) would be a legitimate extension. But the paper applies it to broadband experiments, and that is where it falls down.\n\nThe paper also does not compare with the Clarke-Leighton transient cross-section formula, even though it cites it; it provides no error bars on the fitted parameters and no independent validation. Those are real but softer gaps.\n\nWho should read it: people who characterize UCA shells from attenuation data, and anyone thinking about pulse-duration effects in linear oscillator models. It is a short, clear paper that raises a legitimate concern, but the central application needs rethinking. I would not cite it as it stands, but I would ask the authors to resubmit after clarifying the measurement model, restricting claims to tone-burst measurements or deriving the spectral version, and validating on at least two datasets. The idea deserves serious referee attention.","headline":"A plausible transient-correction formula for UCA attenuation, but the application to broadband data conflates time-domain ring-up with the steady-state transfer function and likely overstates the shell-parameter correction.","tokens_in":7103,"tokens_out":2919,"would_cite":false,"duration_ms":550741,"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":"Short-pulse ultrasound overestimates microbubble attenuation because the bubbles never ring up to steady state before the pulse ends.","keywords":["ultrasound contrast agents","microbubbles","linearized Rayleigh-Plesset equation","transient oscillation","acoustic attenuation","extinction cross-section","dilatational viscosity","pulse duration"],"falsifier":"Use a monodisperse UCA suspension and measure attenuation at a fixed frequency near resonance with identical peak pressure but different pulse durations (for example, 0.2, 0.5, 1 and 5 µs). Eq. (15) predicts the attenuation coefficient falls with pulse duration and the peak is about two-thirds of the steady-state value for very short pulses; the steady-state formula predicts no dependence on pulse length. If the measured attenuation shows no systematic pulse-length dependence, the central correction is not supported.","tokens_in":5970,"feed_emoji":"🫧","tokens_out":9635,"duration_ms":86089,"temperature":0.7,"pith_summary":"Ultrasound contrast agents are coated microbubbles whose shell properties are usually extracted by comparing measured attenuation of a short ultrasound pulse with a theoretical curve that assumes the bubble has reached steady-state oscillation. This paper argues that clinical and laboratory pulses are often too short for the bubble to ring up, so the standard curve overestimates attenuation near resonance. Including the transient part of the linearized Rayleigh-Plesset solution lowers the predicted attenuation peak to about two-thirds of the steady-state value and shifts it to the damped resonance frequency. When the new formula is fitted to published attenuation data, the inferred shell dilatational viscosity drops by roughly a factor of three and the elasticity rises slightly. The practical claim is that shell characterization with steady-state formulas misreads shell physics when the interrogation pulse is short.","feed_headline":"Short pulses shrink microbubble attenuation by a third near resonance","feed_subtitle":"New transient formula cuts inferred shell viscosity about threefold for short-pulse attenuation data.","key_machinery":"The central object is the linearized Rayleigh-Plesset equation written as a forced, damped harmonic oscillator, $\\ddot{X}+2\\delta\\omega_0\\dot{X}+\\omega_0^2X=F\\cos(\\omega t)$, whose complete solution is the sum of a steady-state response and a decaying transient at the damped frequency $\\omega_d=\\omega_0\\sqrt{1-\\delta^2}$. The argument is carried by Eq. (15), the time-averaged extinction cross-section of the full solution over $N$ drive cycles, with $N=2\\ln(10)/(\\pi\\delta)$ chosen as the number of periods required for the transient to decay to 1% of its initial amplitude; the paper notes that in a pulse-echo measurement the integration period should be the pulse duration. Substituting only the steady-state solution into the same energy integral yields the standard attenuation formula as a special case, so Eq. (15) is the general statement and the steady-state formula is its long-pulse limit.","core_discovery":"The paper claims that for a UCA driven by a short tone burst, the acoustic attenuation should be computed from the full solution of the damped driven oscillator—steady-state part plus decaying transient—rather than from the steady-state response alone. The full solution behaves differently only while the transient is alive; for typical lipid-shell parameters and a transducer pulse around 0.4 µs, the bubble has not yet rung up before the pressure ends. Within this transient window, the extinction cross-section—the energy a bubble removes from the sound field per unit incident intensity—is reduced, with the resonance peak falling to roughly 2/3 of the steady-state peak and moving to the damped resonance frequency. Fitting the transient-corrected formula to measured attenuation gives a dilatational viscosity of $0.62 \\times 10^{-8}\\,\\mathrm{N\\,s/m}$ and an elasticity of $0.45\\,\\mathrm{N/m}$, versus $1.97 \\times 10^{-8}\\,\\mathrm{N\\,s/m}$ and $0.40\\,\\mathrm{N/m}$ from the steady-state formula; the viscosity change is large enough that predictions of nonlinear UCA dynamics would differ.","pith_inferences":["If attenuation data are reduced by spectral analysis of a broadband pulse, the frequency-dependent attenuation may be governed by the steady-state transfer function rather than by the single-frequency ring-up transient; Eq. (15)'s correction would then not apply to those experiments.","The factor-of-three viscosity shift implies that published shell parameters obtained with steady-state formulas from broadband data may be systematically high; re-fitting existing datasets with Eq. (15) could revise them.","The same transient-averaging idea could be extended to scattering cross-sections and to nonlinear diagnostics by convolving the forced-oscillator response with the actual pulse envelope, giving a testable prediction that attenuation near resonance should decrease monotonically with pulse shortening.","A decisive check would be to measure attenuation of a monodisperse UCA suspension at fixed frequency and pressure while varying only pulse length; Eq. (15) predicts a measurable drop, whereas the steady-state formula predicts none."],"forward_implications":["Short-pulse attenuation data for UCAs should be fitted with Eq. (15); otherwise the fitted dilatational viscosity is roughly three times too large for typical lipid-shell bubbles.","The observed resonance peak in short-pulse attenuation will sit below the undamped resonance frequency; treating that downshift as a change in shell elasticity conflates a ring-up artifact with a material property.","Shell parameters inferred with the steady-state formula may mispredict nonlinear bubble behavior, such as subharmonic and ultraharmonic emission amplitudes.","For long pulses or continuous-wave excitation the transient correction vanishes and the two formulas agree, so the correction matters only when the pulse duration is comparable to or shorter than the ring-up time.","Attenuation coefficients become pulse-duration dependent: the same suspension can appear to attenuate differently with different transducers unless the transient contribution is modeled."],"supporting_citations":[{"why":"Supplies the steady-state attenuation formula that the paper identifies as the standard and corrects.","marker":"Medwin 1977"},{"why":"Gives the energy-based extinction-cross-section derivation and the equivalence of oscillator and effective-medium approaches from which Eq. (15) is built.","marker":"Xia 2018"},{"why":"Provides the viscoelastic shell model used for the UCA's linear dynamics.","marker":"Sarkar, Shi et al. 2005"},{"why":"Contributes the measured attenuation data and previously fitted shell parameters used in the comparison.","marker":"Xia, Porter et al. 2015"},{"why":"Establishes that bubble attenuation before steady state depends on cycle length and drive amplitude, motivating the transient treatment.","marker":"Clarke and Leighton 2000"},{"why":"First reported the reduction in dissipation attributed to ring-up time effects.","marker":"Akulichev, Bulanov et al. 1986"},{"why":"Provides the effective-medium linear attenuation derivation that is equivalent to the oscillator energy method at low volume fractions.","marker":"Commander and Prosperetti 1989"},{"why":"Documents broadband-transducer attenuation measurements used in shell characterization, the setting where short pulses make the correction relevant.","marker":"Frinking and de Jong 1998"}],"fun_headline_variants":["Transient microbubble dynamics alter attenuation estimates","Short-pulse ultrasound reveals transient bubble attenuation dip","Transient Rayleigh-Plesset solution lowers attenuation, viscosity","Bubble attenuation under short pulses needs transient correction","Transient term slices microbubble attenuation, recalibrates shell"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's correction assumes that a short-pulse attenuation measurement corresponds to the time-averaged energy absorption of a single-frequency forced oscillator over the pulse duration, with the number of effective cycles set by the damping constant; if attenuation is instead extracted from the frequency spectrum of a broadband pulse, the steady-state response may already govern and the correction would not apply.","fun_headline_variants_meta":{"raw":{"variants":["Transient microbubble dynamics alter attenuation estimates","Short-pulse ultrasound reveals transient bubble attenuation dip","Transient Rayleigh-Plesset solution lowers attenuation, viscosity","Bubble attenuation under short pulses needs transient correction","Transient term slices microbubble attenuation, recalibrates shell"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000879,"raw_usage":{"total_tokens":3788,"prompt_tokens":918,"completion_tokens":2870,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":2805}},"tokens_in":534,"tokens_out":2870,"duration_ms":18535,"temperature":1.0,"reasoning_tokens":2805,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:01:01.788141+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use a monodisperse UCA suspension and measure attenuation at a fixed frequency near resonance with identical peak pressure but different pulse durations (for example, 0.2, 0.5, 1 and 5 µs). Eq. (15) predicts the attenuation coefficient falls with pulse duration and the peak is about two-thirds of the steady-state value for very short pulses; the steady-state formula predicts no dependence on pulse length. If the measured attenuation shows no systematic pulse-length dependence, the central correction is not supported.","supporting_citations":[{"cited_title":"Counting bubbles acoustically: a review","cited_arxiv_id":null,"evidence_quote":"Supplies the steady-state attenuation formula that the paper identifies as the standard and corrects."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the energy-based extinction-cross-section derivation and the equivalence of oscillator and effective-medium approaches from which Eq. (15) is built."},{"cited_title":"Characterization of ultrasound contrast microbubbles using in vitro experiments and viscous and viscoelastic interface models for encapsulation","cited_arxiv_id":null,"evidence_quote":"Provides the viscoelastic shell model used for the UCA's linear dynamics."},{"cited_title":"Interpreting attenuation at different excitation amplitudes to estimate strain- dependent interfacial rheological properties of lipid -coated monodisperse microbubbles","cited_arxiv_id":null,"evidence_quote":"Contributes the measured attenuation data and previously fitted shell parameters used in the comparison."},{"cited_title":"A method for estimating time-dependent acoustic cross - sections of bubbles and bubble clouds prior to the steady state","cited_arxiv_id":null,"evidence_quote":"Establishes that bubble attenuation before steady state depends on cycle length and drive amplitude, motivating the transient treatment."},{"cited_title":"Acoustic sounding of gas bubbles in sea water","cited_arxiv_id":null,"evidence_quote":"First reported the reduction in dissipation attributed to ring-up time effects."},{"cited_title":"Linear pressure waves in bubbly liquids: Comparison between theory and experiments","cited_arxiv_id":null,"evidence_quote":"Provides the effective-medium linear attenuation derivation that is equivalent to the oscillator energy method at low volume fractions."},{"cited_title":"Acoustic modeling of shell -encapsulated gas bubbles","cited_arxiv_id":null,"evidence_quote":"Documents broadband-transducer attenuation measurements used in shell characterization, the setting where short pulses make the correction relevant."}],"review_version":1}