pith. sign in

arxiv: 2212.04466 · v24 · submitted 2022-12-08 · 🧮 math.NA · cs.NA

A Full-waveform Approximation of Finite-Sized Acoustic Apertures: Forward and Adjoint Wavefields

Pith reviewed 2026-05-24 10:05 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords acoustic wave equationfull-waveform approximationadjoint operatorKirchhoff-Helmholtz integralreception operatorsDirichlet boundary dataphotoacoustic tomographytherapeutic ultrasound
0
0 comments X

The pith

Under Dirichlet boundary data the adjoint of the acoustic forward operator coincides with the time-reversed dipole integral on receiver surfaces.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper first shows that the analytic monopole and dipole surface integrals are exactly equivalent to their full-waveform approximations for finite apertures. It then defines reception operators that extract boundary measurements from free-space pressure fields obtained by solving the wave equation. Using this trace mapping the adjoint of the forward operator is derived, and the central result follows: under the standard assumption of Dirichlet-type boundary data the adjoint equals a constant multiple of the interior time-reversed dipole integral evaluated on the receivers. This equivalence supplies a practical route to accurate forward modeling and adjoint-state inversion when amplitude fidelity matters.

Core claim

The paper establishes equivalence between the analytic expressions of the monopole and dipole integral formulations and their full-waveform approximations. It introduces reception operators that map free-space pressure wavefields onto measured fields restricted to the boundary. The adjoint of the forward operator is derived from this mapping. Under the common practical assumption of Dirichlet-type boundary data, the adjoint operator coincides up to a constant factor with the interior-field time-reversed form of the dipole integral formula evaluated on the receiver surfaces.

What carries the argument

Reception operators that map free-space pressure wavefields to boundary measurements, from which the adjoint is obtained and shown to match the time-reversed dipole integral.

If this is right

  • Forward modeling of finite acoustic apertures can use the full-waveform approximations in place of the classical integrals.
  • Adjoint-state methods for inverse problems become implementable by simple time reversal on the receiver surfaces.
  • Amplitude-preserving reconstructions are available for therapeutic ultrasound optimization and attenuation imaging.
  • The same adjoint construction applies directly to photoacoustic tomography under the stated boundary assumption.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The result may allow existing time-reversal codes to serve as adjoint operators without additional surface-integral machinery.
  • If the Dirichlet assumption is relaxed, a different scaling or correction term would be needed to recover the adjoint.
  • The same trace-mapping approach could be examined for other linear wave equations where surface sources appear.

Load-bearing premise

Boundary data are of Dirichlet type so that pressure is prescribed on the surface.

What would settle it

A direct numerical comparison, for a known Dirichlet boundary setup, of the adjoint operator computed by the derived formula versus the interior time-reversed dipole integral; any systematic mismatch falsifies the claimed coincidence.

Figures

Figures reproduced from arXiv: 2212.04466 by Ashkan Javaherian, Seyed Kamaledin Setarehdan.

Figure 1
Figure 1. Figure 1: (a) A single emitter point and 39 receiver points positioned on a hemi￾sphere with a radius of 5.6cm, (b) source pulse, s, in the time domain, (c) source pulse, s, decomposed into amplitude and phase in the frequency domain. 6.1.2. Result [PITH_FULL_IMAGE:figures/full_fig_p018_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Analytic action of the causal Green’s function on a source defined on a single point. The analytic action approximated on receiver 10 and at all chosen frequencies, a) amplitude, b) phase. The analytic action approximated at the single frequency 1MHz and on all receiver points, c) amplitude d) phase. approximated using the full-field algorithm 1 and using a mass source discretised in space and time by Eq. … view at source ↗
Figure 3
Figure 3. Figure 3: a) A disc-shaped emitter and 64 receiver points. The positions are ordered by changing θ, φ and r. b) source pulse u n (rigid-baffle condition), c) source pulse p (soft-baffle condition). In addition, the azimuthal angles are set θ ∈ {π/4, 3π/4, 5π/4, 7π/4}. The receiver positions are ordered based on changing θ, φ and r in order. Accordingly, for receiver points 1-4, r = 6.5cm and φ = 0 are set, and the p… view at source ↗
Figure 4
Figure 4. Figure 4: The wavefield approximated on the plane x 1 = 3.1cm and at a single time 45µs after excitation of the disc-shaped emitter by a source pulse u n, which is shown in figure 3(b). The emitter disc’s centre is placed on the origin of the Cartesian coordinates, as shown in figure 3(a) by the yellow colour. (Not shown here.) a) Field II, b) Full-field (RE= 7.10%). by a factor 4. Consider that for both analytic an… view at source ↗
Figure 5
Figure 5. Figure 5: The wavefield approximated and recorded in time on the receivers after an excitation of the disc-shaped emitter by a source pulse u n, which is shown in figure 3(b). Receiver points: a) 1, b) 5, c) 9, d) 13. The full-field approach and the Field II toolbox were used to approximate the monopole integral formula (32). the dipole integral formula (34), which is equivalent to a summation of the actions of the … view at source ↗
Figure 6
Figure 6. Figure 6: The relative error (RE) of the wavefield approximated by the full-field approach in time and on the chosen receiver points after an excitation of the emitter disc by a source pulse u n, which is shown in figure 3(b). The wavefield analytically calculated using the Field II toolbox is used as the benchmark. The receiver positions are ordered by changing θ, φ and r. On-grid sampling. As explained in section … view at source ↗
Figure 7
Figure 7. Figure 7: The wavefield approximated on the plane x 1 = 3.1cm and at a time 45µs after excitation of the disc-shaped emitter by a source pulse p, which is shown in figure 3(c). The centre of the emitter disc is placed on the origin of the Cartesian coordinates, as shown in figure 3(a). (Not shown here.) a) Field II (far-field term) + Analytic (near-field term) , b) Full-field using a mass source discretised using Eq… view at source ↗
Figure 8
Figure 8. Figure 8: The wavefield approximated on the receivers in time after an excitation of the disc-shaped emitter by a source pulse p, which is shown in figure 3(c). Receiver points: a) 1, b) 5, c) 9, d) 13. The full-field approach that uses algorithm 1 and a momentum source defined by Eq. (52) approximates directly the original dipole formula (34). The Field II toolbox is used to approximate the far-field term of the di… view at source ↗
read the original abstract

The acoustic wave equation governs wave propagation induced by either volumetric radiation sources, or by surface sources of monopole or dipole type. For surface sources, boundary value problems yield wavefield representations via the Kirchhoff-Helmholtz or Rayleigh-Sommerfeld integrals. This study begins by establishing an equivalence between the analytic expressions of the associated monopole and dipole integral formulations and their full-waveform approximations. Leveraging this equivalence, we introduce reception operators that map free space pressure wavefields-obtained by solving the wave equation-onto measured fields restricted to the boundary. Building on this trace mapping, we derive the adjoint of the forward operator. We show that, under the common practical assumption of Dirichlet-type boundary data, the adjoint operator coincides-up to a constant factor-with the interior-field time-reversed form of the dipole integral formula, evaluated on the receiver surfaces. This study aims to advance the approximation of forward problems and the solution of inverse problems in acoustics, with a particular focus on applications that require accurate amplitude modeling, including therapeutic ultrasound optimization, attenuation reconstruction, and photoacoustic tomography.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper establishes an equivalence between the analytic expressions of the monopole and dipole integral formulations (via Kirchhoff-Helmholtz and Rayleigh-Sommerfeld identities) and their full-waveform approximations for the acoustic wave equation with surface sources. It defines reception operators as trace maps from free-space pressure solutions of the wave equation onto boundary-restricted measured fields, then derives the adjoint of this forward operator. Under the explicitly stated assumption of Dirichlet-type boundary data, the adjoint is shown to coincide, up to a constant factor, with the interior-field time-reversed form of the dipole integral formula evaluated on the receiver surfaces. The work targets improved forward approximations and adjoint-based inverse problems in acoustics.

Significance. If the derivations hold, the conditional equivalence supplies a parameter-free link between boundary-integral representations and full-waveform modeling for finite apertures, preserving amplitude information needed for therapeutic ultrasound optimization, attenuation reconstruction, and photoacoustic tomography. The explicit adjoint construction under a standard practical assumption enables direct use of time-reversal techniques without hidden regularity requirements or unstated limit interchanges. The paper's strength is the upfront declaration of scope, rendering the central claim falsifiable by direct numerical comparison of the two operator forms.

minor comments (2)
  1. [Abstract] The abstract would be strengthened by a single sentence stating the precise scaling factor in the adjoint identity or the section where the equivalence is proved.
  2. Figure captions should explicitly label the Dirichlet boundary data case when comparing forward and adjoint fields.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. No specific major comments appear in the provided report, so we offer no point-by-point rebuttals. We will address any editorial or minor suggestions in the revised manuscript.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivation begins from the acoustic wave equation and standard Kirchhoff-Helmholtz / Rayleigh-Sommerfeld integral identities, constructs explicit reception (trace) operators, and obtains the adjoint via the usual duality pairing. The central claim is an explicit conditional statement: under the declared Dirichlet boundary assumption the adjoint equals (up to scaling) the time-reversed interior dipole integral. No parameter is fitted to data and then relabeled a prediction, no self-citation supplies a uniqueness theorem or ansatz, and the equivalence is not obtained by redefinition. The chain is therefore self-contained against external functional-analysis benchmarks for adjoints of boundary-trace operators.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Review performed on abstract only; no explicit free parameters, invented entities, or additional axioms are stated beyond the Dirichlet boundary assumption and the claimed equivalence.

axioms (2)
  • domain assumption Equivalence holds between analytic monopole/dipole integral expressions and their full-waveform approximations
    Stated as the starting point of the study in the abstract.
  • domain assumption Dirichlet-type boundary data is the relevant practical case
    Invoked to obtain the adjoint coincidence result.

pith-pipeline@v0.9.0 · 5724 in / 1287 out tokens · 23560 ms · 2026-05-24T10:05:28.754605+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

74 extracted references · 74 canonical work pages

  1. [1]

    A k-space method for large-scale models of wave propagation in tissue

    T. D. Mast, L. P. Souriau, D. -L. D. Liu, M. Tabei, A. I. Nachman and R. C. Waag, “A k-space method for large-scale models of wave propagation in tissue”,IEEE Trans. Ultrason. Ferroelectr. Freq. , vol. 48, no. 2, pp. 341-354, March 2001, doi: 10.1109/58.911717

  2. [2]

    A k-space method for coupled first-order acoustic propagation equations

    M. Tabei, T. D. Mast, and R. C. Waag, “A k-space method for coupled first-order acoustic propagation equations”, J. Acoust. Soc. Am. vol. 111, pp. 53–63, 2002

  3. [3]

    k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields

    B. E. Treeby and B. T. Cox, “k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields”,J. Biomed. Opt. vol. 15, no. 2, 021314, 2010. THIS ARTICLE HAS NOT BEEN REVIEWED YET. 27

  4. [4]

    A causal and fractional all-frequency wave equation for lossy media

    S. Holm and S.P. Näsholm, “A causal and fractional all-frequency wave equation for lossy media”,J. Acoust. Soc. Am., Vol. 130, no. 4, pp. 2195-2202, 2011

  5. [5]

    A Matlab toolbox for the time domain simulation of acoustic wave fields

    B. Treeby and B. Cox, k-Wave user manual, “A Matlab toolbox for the time domain simulation of acoustic wave fields”, Version 1.1, 27th August 2016 (the last version)

  6. [6]

    Directional Sources in Wave-Based Acoustic Simulation

    S. Bilbao and B. Hamilton, “Directional Sources in Wave-Based Acoustic Simulation”, in IEEE/ACM Transactions on Audio, Speech, and Language Processing, vol. 27, no. 2, pp. 415-428, Feb. 2019, doi: 10.1109/TASLP.2018.2881336

  7. [7]

    Discretization of multipole sources in a finite difference setting for wave propagation problems

    M. J. Bencomo and W. W. Symes, “Discretization of multipole sources in a finite difference setting for wave propagation problems”,J. Comput. Phys. , vol. 386, pp. 296-322, 2019

  8. [8]

    Stability for semilinear wave equation in an inhomogeneous medium with frictional localized damping and acoustic boundary conditions

    M. M. Cavalcanti, V. N. Domingos Cavalcanti, C. L. Frota and A. Vicente, “Stability for semilinear wave equation in an inhomogeneous medium with frictional localized damping and acoustic boundary conditions”,SIAM J Control Optim, Vol. 58, no. 4, pp. 2411–2445, 2020

  9. [9]

    Nonlinear wave equation with Dirichlet and Acoustic boundary conditions: theoretical analysis and numerical simulation

    A. A. Alcântara, B. A. Carmo, H. R. Clark, R. R. Guardia and M. A. Rincon, “Nonlinear wave equation with Dirichlet and Acoustic boundary conditions: theoretical analysis and numerical simulation”,Comp. Appl. Math. Vol. 41, no. 141, 2022,https://doi.org/10.1007/s40314-022-01822-5

  10. [10]

    The importance of transfer learning in seismic modeling and imaging

    A. Siahkoohi, M. Louboutin and F. J. Herrmann, “The importance of transfer learning in seismic modeling and imaging”,Geophysics, Vol. 84, no. 6, pp. A47-A52, 2019,https://doi.org/10.1190/geo2019-0056.1

  11. [11]

    Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast

    A. Javaherian, F. Lucka and B. Cox, “Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast”,Inverse Problems, vol. 36, no. 12, 125010, 2020

  12. [12]

    Wave-equation migration velocity analysis via the optimal-transport-based objective function

    Y. Tang, B. Sun and T. Alkhalifah, “Wave-equation migration velocity analysis via the optimal-transport-based objective function”,Geophysics, Vol. 87, no. 3, pp. U109–U120, 2022

  13. [13]

    Inverse medium scattering problems with Kalman filter techniques

    T. Furuya and R. Potthast, “Inverse medium scattering problems with Kalman filter techniques”,Inverse Prob- lems, Vol. 38, no. 9, 095003, 2022

  14. [14]

    Recovery of wave speeds and density of mass across a heterogeneous smooth interface from acoustic and elastic wave reflection operators

    S. Bhattacharyya, M. V. de Hoop, V. Katsnelson and G. Uhlmann, “Recovery of wave speeds and density of mass across a heterogeneous smooth interface from acoustic and elastic wave reflection operators”,GEM-International Journal on Geomathematics , Vol. 13, no. 1, 2022

  15. [15]

    Seismic differential semblance-oriented migration velocity analysis — Status and the way forward

    H. Chauris and M. Farshad, “Seismic differential semblance-oriented migration velocity analysis — Status and the way forward”,Geophysics, Vol. 88, no. 6, pp. U81-U100, 2023

  16. [16]

    Moment tensor inversion of perforation shots using distributed-acoustic sensing

    M. Bader, R. G. Clapp, K. T. Nihei, and B. Biondi, “Moment tensor inversion of perforation shots using distributed-acoustic sensing”, Geophysics,Vol. 88, no. 6, pp. 1-36, 2023. https://doi.org/10.1190/ geo2021-0661.1

  17. [17]

    High- resolution static corrections derived from surface-wave tomography: Application to mineral exploration

    M. Papadopoulou, B. Brodic, E. Koivisto, A. Kaleshova, M. Savolainen, P. Marsden, and L. V. Socco, “High- resolution static corrections derived from surface-wave tomography: Application to mineral exploration”,Geo- physics, Vol. 88, no. 6, B317-B328, 2023

  18. [18]

    Field assessment of elastic full- waveform inversion of combined accelerometer and distributed acoustic sensing data in a vertical seismic profile configuration

    M. V. Eaid, S. D. Keating, K. A. Innanen, M. Macquet and D. Lawton, “Field assessment of elastic full- waveform inversion of combined accelerometer and distributed acoustic sensing data in a vertical seismic profile configuration”,Geophysics, vol. 88, no. 6, pp. WC163-WC180, 2023

  19. [19]

    Figueroa

    B. Kaltenbacher and W. Rundell, “On the simultaneous reconstruction of the nonlinearity coefficient and the sound speed in the Westervelt equation”,Inverse Problems, Vol. 39, no. 10, p. 105001, 2023, DOI 10.1088/1361- 6420/aceef2

  20. [20]

    Reweighted variational full-waveform inversions

    W. Wang, G. A. McMechan and J. Ma,“Reweighted variational full-waveform inversions”,Geophysics, Vol. 88, no. 4, R-499-R512https://doi.org/10.1190/geo2021-0766.1

  21. [21]

    An inverse boundary value problem arising in nonlinear acoustics

    G. Uhlmann and Y. Zhang, “An inverse boundary value problem arising in nonlinear acoustics”,SIAM Journal on Mathematical Analysis, Vol. 55, no. 2, pp. 1364-1404, 2023

  22. [22]

    Kaltenbacher and V

    B. Kaltenbacher and V. Nikolić, The vanishing relaxation time behavior of multi-term nonlocal Jor- dan–Moore–Gibson–Thompson equations,Nonlinear Analysis: Real World Applications , Vol. 76, pp. 103991, 2024, https://doi.org/10.1016/j.nonrwa.2023.103991

  23. [23]

    An Efficient Neumann series–based algorithm for thermoa- coustic and photoacoustic tomography with variable sound speed

    J. Qian, P. Stefanov, G. Uhlmann and H. Zhao, “An Efficient Neumann series–based algorithm for thermoa- coustic and photoacoustic tomography with variable sound speed”,SIAM J. Imaging Sci. , Vol.4, no.3, 2011. doi: 10.1137/100817280

  24. [24]

    Reconstructing absorption and scattering distributions in quantitative photoacoustic tomography

    T. Tarvainen, B. T. Cox, J. Kaipio, and S.R. Arridge, “Reconstructing absorption and scattering distributions in quantitative photoacoustic tomography”,Inverse Problems, Vol. 28, 2012 ,p.084009

  25. [25]

    Attenuation Models in Photoacoustics

    R. Kowar and O. Scherzer, “Attenuation Models in Photoacoustics”. In: H. Ammari (eds) Mathematical Modeling in Biomedical Imaging II. Lecture Notes in Mathematics, Vol. 2035, 2012. Springer, Berlin, Heidelberg.https: //doi.org/10.1007/978-3-642-22990-94. 28 THIS ARTICLE HAS NOT BEEN REVIEWED YET

  26. [26]

    Accurate model-based reconstruction algorithm for three-dimensional optoacoustic tomography

    X. L. Dean-Ben, A. Buehler, V. Ntziachristos and D. Razansky, “Accurate model-based reconstruction algorithm for three-dimensional optoacoustic tomography”,IEEE T-MI, Vol.31, no.10, pp.19221928, 2012

  27. [27]

    Model-Based Learning for Accelerated, Limited-View 3-D Photoacoustic Tomography

    A. Hauptmann et al., “Model-Based Learning for Accelerated, Limited-View 3-D Photoacoustic Tomography”, IEEE T-MI, vol. 37, no. 6, pp. 1382-1393, June 2018, doi: 10.1109/TMI.2018.2820382

  28. [28]

    Analysis of iterative methods in photoacoustic tomography with variable sound speed

    M.Haltmeier and L.V.Nguyen, “Analysis of iterative methods in photoacoustic tomography with variable sound speed”,SIAM Journal on Imaging Sciences , Vol.10,no.2,2017, doi: 10.1137/16M1104822

  29. [29]

    A continuous adjoint for photo-acoustic tomography of the brain

    A. Javaherian and S. Holman, “A continuous adjoint for photo-acoustic tomography of the brain”,Inverse Problems, vol. 34, no. 8, p. 085003, 2018

  30. [30]

    Direct quantitative photoacoustic tomography for realistic acoustic media

    A.Javaherian and S.Holman, “Direct quantitative photoacoustic tomography for realistic acoustic media”,Inverse Problems,vol.35,no.8,084004,2019

  31. [31]

    Deep-learning for photoacoustic tomography from sparse data

    S.Antholzer, M.Haltmeier and J.Schwab, “Deep-learning for photoacoustic tomography from sparse data”,Inverse Problems in Science and Engineering , Vol.27, no.7,pp.987-1005, 2019. doi: 10.1080/17415977.2018.1518444

  32. [32]

    Fully dense UNet for 2-D sparse photoacoustic to- mography artifact removal

    S. Guan, A. A. Khan, S. Sikdar and P. V. Chitnis, “Fully dense UNet for 2-D sparse photoacoustic to- mography artifact removal”,IEEE J. Biomed. Health. Inform. , vol. 24, no. 2, pp. 568-576, Feb. 2020, doi: 10.1109/JBHI.2019.2912935

  33. [33]

    Photoacoustic computed tomography for functional human brain imaging

    S. Na and L. V. Wang, “Photoacoustic computed tomography for functional human brain imaging”,Biomed. Opt. Express, Vol. 12, pp. 4056-4083, 2021

  34. [34]

    Analysis for full-field photoacoustic tomography with variable sound speed

    L. Nguyen, M. Haltmeier, R Kowar, and N. Do, “Analysis for full-field photoacoustic tomography with variable sound speed”,SIAM J. Imaging Sci. , Vol. 15, no. 3, 2022, 10.1137/21M1463409

  35. [35]

    Single-stageapproachforestimatingopticalparametersinspectral quantitative photoacoustic tomography

    M.Suhonen, A.Pulkkinen, andT.Tarvainen, “Single-stageapproachforestimatingopticalparametersinspectral quantitative photoacoustic tomography”,J. Opt. Soc. Am. , Vol. 41, no. 3, pp. 527-542, 2024

  36. [36]

    A frequency-domain model-based reconstruction method for transcranial photoa- coustic imaging: A 2D numerical investigation

    H. Park, J. Yao and Y. Jing, “A frequency-domain model-based reconstruction method for transcranial photoa- coustic imaging: A 2D numerical investigation”,Photoacoustics, Vol. 33, 2023, p. 100561

  37. [37]

    Scalable Funding of Bitcoin Micropayment Channel Networks

    A. Hauptmann ans T. Tarvainen, (2024). “Model-Based Reconstructions for Quantitative Imaging in Photoacous- tic Tomography”. In: Xia, W. (eds) Biomedical Photoacoustics. Springer, Cham.https://doi.org/10.1007/978-3- 031-61411-8_4

  38. [38]

    Full-wave image reconstruction in transcranial photoacoustic computed tomography using a finite element method,

    Y. Luo et al., “Full-wave image reconstruction in transcranial photoacoustic computed tomography using a finite element method,” inIEEE T-MI, doi:10.1109/TMI.2024.3456595

  39. [39]

    A forward Model incorporating elevation-focused transducer properties for 3-D full-waveform inversion in ultrasound computed tomography

    F. Li, U. Villa, N. Duric and M. A. Anastasio, “A forward Model incorporating elevation-focused transducer properties for 3-D full-waveform inversion in ultrasound computed tomography”,IEEE T-UFFC, vol. 70, no. 10, pp. 1339-1354, Oct. 2023, doi: 10.1109/TUFFC.2023.3313549

  40. [40]

    Frequency domain ultrasound waveform tomography: breast imaging using a ring transducer

    G. Y. Sandhu, C. Li, O. Roy, S. Schmidt and N. Duric, “Frequency domain ultrasound waveform tomography: breast imaging using a ring transducer”,Phys. Med. Biol. , Vol. 60, 5381–5398, 2015

  41. [41]

    Iterative methods for solving coefficient inverse problems of wave tomography in models with attenuation

    A. V. Goncharsky and S. Y. Romanov, “Iterative methods for solving coefficient inverse problems of wave tomography in models with attenuation”,Inverse Problems, vol. 33, pp. 025003, 2017

  42. [42]

    3-D Nonlinear Acoustic Inverse Scattering: Algorithm and Quantitative Results

    J. W. Wiskin, D. T. Borup, E. Iuanow, J. Klock and M. W. Lenox, “3-D Nonlinear Acoustic Inverse Scattering: Algorithm and Quantitative Results”,IEEE T-UFFC, vol. 64, no. 3, 2017

  43. [43]

    Full-waveform inversion imaging of the human brain

    L. Guasch, O. Calderón Agudo, M. Tang, P. Nachev, and M. Warner, “Full-waveform inversion imaging of the human brain”.Nature Digital Medicine, vol. 3, 28, 2020

  44. [44]

    Adjoint-state method for Hybridizable Discontinuous Galerkin discretization

    F. Faucher and O. Scherzer,“Adjoint-state method for Hybridizable Discontinuous Galerkin discretization”, ap- plication to the inverse acoustic wave problem”,Comput. Methods Appl. Mech. Eng. , Vol. 372, pp. 113406, 2020

  45. [45]

    Full- waveform inversion with resolution proxies for in-vivo ultrasound computed tomography

    I. E. Ulrich, S. Noe, C. Boehm, N K Martiartu, B Lafci, X. L. Dean-Ben, D. Razansky and A.Fitchner, “Full- waveform inversion with resolution proxies for in-vivo ultrasound computed tomography”, 2023 IEEE International Ultrasonics Symposium (IUS), Montreal, QC, Canada, 2023, pp. 1-4, doi: 10.1109/IUS51837.2023.10308297

  46. [46]

    Robust Imaging of Speed of Sound Using Virtual Source Transmission

    D. Schweizer, R. Rau, C. D. Bezek, R. A. Kubik-Huch and O. Goksel, “Robust Imaging of Speed of Sound Using Virtual Source Transmission”, inIEEE T-UFFC, Vol. 70, no. 10, pp. 1308-1318, 2023

  47. [47]

    Neural Born series operator for biomedical ultrasound computed Tomography

    Z. Zeng, Y. Zheng, Y. Zheng, Y. Li, Z. Shi Aand H. Sun, “Neural Born series operator for biomedical ultrasound computed Tomography”, 2023,https://arxiv.org/abs/2312.15575

  48. [48]

    Extending the search space of full-waveform inversion beyond the single-scattering Born approxi- mation: A tutorial review

    S. Operto, A. Gholami, H. S. Aghamiry, G. Guo, S. Beller, K. Aghazade, F. Mamfoumbi, L. Combe and A. Ribodetti, “Extending the search space of full-waveform inversion beyond the single-scattering Born approxi- mation: A tutorial review”,Geophysics, Vol. 88, no. 6, pp. R671–R702, 2023. doi:https://doi.org/10.1190/ geo2022-0758.1

  49. [49]

    Ultrasound Brain Tomography: Comparison of Deep Learning and Deterministic Methods

    M. Soleimani, T. Rymarczyk and G. Kłosowski, “Ultrasound Brain Tomography: Comparison of Deep Learning and Deterministic Methods”, inIEEE Trans. Instrum. Meas. , Vol. 73, pp. 1-12, 2024, Art no. 4500812, doi: 10.1109/TIM.2023.3330229. THIS ARTICLE HAS NOT BEEN REVIEWED YET. 29

  50. [50]

    Waveform inversion with a data driven estimate of the internal wave

    L. Borcea, J. Garnier, A. V. Mamonov and J. Zimmerling, “Waveform inversion with a data driven estimate of the internal wave”,SIAM J. Imaging Sci. , Vol. 16, no. 1, pp. 280-312, 2023,https://doi.org/10.1137/22M1517342

  51. [51]

    Numerical simulations of clinical focused ultrasound functional neurosurgery

    A. Pulkkinen, B. Werner, E. Martin and K. Hynynen, “Numerical simulations of clinical focused ultrasound functional neurosurgery”,Physics in Medicine & Biology , vol. 59, no. 7, p. 1679, 2014

  52. [52]

    Full-wave acoustic and thermal modeling of transcranial ultrasound propagation and investigation of skull-induced aberration correction techniques: a feasibility study

    A. Kyriakou, E. Neufeld, and B. Werner, G. Székely and N. Kuster, “Full-wave acoustic and thermal modeling of transcranial ultrasound propagation and investigation of skull-induced aberration correction techniques: a feasibility study”.J Ther Ultrasound , Vol. 3, no. 11, 2015.https://doi.org/10.1186/s40349-015-0032-9

  53. [53]

    Numerical evaluation of the skull for human neuromodulation with transcranial focused ultrasound

    J. K. Mueller, L. Ai, P. Bansal and W. Legon, “Numerical evaluation of the skull for human neuromodulation with transcranial focused ultrasound”,J. Neural Eng. , Vol. 14, p.066012 (19pp), 2017

  54. [54]

    A viscoelastic model for the prediction of transcranial ultrasound propagation: Application for the estimation of shear acoustic properties in the human skull

    S. Pichardo, C. Moreno-Hernandez, R. Andrew Drainville, V. Sin, L. Curiel, and K. Hynynen, “A viscoelastic model for the prediction of transcranial ultrasound propagation: Application for the estimation of shear acoustic properties in the human skull”,Phys. Med. Biol. Vol. 62, no. 17, pp. 6938–6962, 2017

  55. [55]

    Safety of Transcranial focused ultrasound stimulation: A systematic review of the state of knowledge from both human and animal studies

    C. Pasquinelli, L.G. Hanson, H.R. Siebner, H.J. Lee and A. Thielscher, “Safety of Transcranial focused ultrasound stimulation: A systematic review of the state of knowledge from both human and animal studies”,Brain Stimul., Vol. 12, no. 6, pp. 1367-1380, 2019. doi: 10.1016/j.brs.2019.07.024. Epub 2019 Jul 31. PMID: 31401074

  56. [56]

    Histologic safety of transcranial focused ultrasound neuromodulation and magnetic resonance acous- tic radiation force imaging in rhesus macaques and sheep

    P. Gaur, K.M. Casey, J. Kubanek, N. Li, M. Mohammadjavadi, Y. Saenz, G.H. Glover, D.M. Bouley and K.B. Pauly. “Histologic safety of transcranial focused ultrasound neuromodulation and magnetic resonance acous- tic radiation force imaging in rhesus macaques and sheep”.Brain Stimul . 2020 May-Jun;13(3):804-814. doi: 10.1016/j.brs.2020.02.017. Epub 2020 Feb ...

  57. [57]

    Comparison Between Ray-Tracing and Full-Wave Simulation for Transcranial Ultrasound Focusing on a Clinical System Using the Transfer Matrix Formalism

    T. Bancel et al., “Comparison Between Ray-Tracing and Full-Wave Simulation for Transcranial Ultrasound Focusing on a Clinical System Using the Transfer Matrix Formalism”, inIEEE T-UFFC , vol. 68, no. 7, pp. 2554-2565, July 2021, doi: 10.1109/TUFFC.2021.3063055

  58. [58]

    Benchmark problems for transcranial ultrasound simulation: Inter- comparison of compressional wave models

    J-F Aubry, O. Bates, C. Boehm, K. B. Pauly, D. Christensen, C. Cueto, P. Gélat, L. Guasch, J. Jaros, Y. Jing, R. Jones, N. Li, P. Marty, H. Montanaro, E. Neufeld, S. Pichardo, G. Pinton, A. Pulkkinen, A. Stanziola, A. Thielscher, B Treeby and E. V. Wout, “Benchmark problems for transcranial ultrasound simulation: Inter- comparison of compressional wave mo...

  59. [59]

    ITRUSST Consensus on Biophysical Safety for Transcranial Ultrasonic Stimulation

    J-F Aubry, D. Attali, M. Schafer, E. Fouragnan, C. Caskey, R. Chen, G. Darmani, E. J. Bubrick, J. Sallet, C. Butler, C. Stagg, M. Klein-Flugge, S-S Yoo, B. Treeby, L. Verhagen and K. B. Pauly, “ITRUSST Consensus on Biophysical Safety for Transcranial Ultrasonic Stimulation”, 2023.https://arxiv.org/abs/2311.05359

  60. [60]

    A practical guide to transcranial ultrasonic stimulation from the IFCN-endorsed ITRUSST consortium

    K. R. Murphy, T. Nandi, B. Kop, T. Osada, M. Lueckel, W. A. N’Djin, K. A. Caulfield, A. Fomenko, H. R Siebner, Y. Ugawa, L. Verhagen, S. Bestmann, E. Martin, K. B. Pauly, E. Fouragnan and T. O. Bergmann, “A practical guide to transcranial ultrasonic stimulation from the IFCN-endorsed ITRUSST consortium”, 2024. https://arxiv.org/abs/2407.07646

  61. [61]

    E. V. Wout, S. R. Haqshenas, P. Gélat, T. Betcke and N. Saffari, Boundary integral formulations for acoustic modelling of high-contrast media,COMPUT. MATH. APPL., Vol. 105, pp. 136-149, 2022.https://doi.org/10. 1016/j.camwa.2021.11.021

  62. [62]

    Acoustic multipole source–simplified lattice Boltzmann method for simulating acoustic propagation problems

    R. Chen, S. Cao, W. Liu, Q. Song and Y. You, “Acoustic multipole source–simplified lattice Boltzmann method for simulating acoustic propagation problems”,Int J Numer Meth Fluids . 2023; 95(7): 1174–1196. doi: 10.1002/fld.5178

  63. [63]

    Integral formulations for predicting acoustic radiation

    S. F. Wu, “Integral formulations for predicting acoustic radiation”.J. Acoust. Soc. Am. 1 October 2024; 156 (4): R7–R9.https://doi.org/10.1121/10.0030390

  64. [64]

    S. F. Wu, P. Zhou and Y. Lu, Determining excitation forces acting on the interior surface of an enclosure. Part I: Theory,J. Theor. Comput. Acoust. , Vol. 30, no. 1, 2022, 2250001

  65. [65]

    Universal back-projection algorithm for photoacoustic computed tomography

    M. Xu and L. V. Wang, “Universal back-projection algorithm for photoacoustic computed tomography”,PHYS- ICAL REVIEW E , Vol. 71, pp. 016706 (2005)

  66. [66]

    Exact and approximative imaging methods for photoacoustic tomography using an arbitrary detection surface

    P. Burgholzer, G. J. Matt, M. Haltmeier, Markus and G. Paltauf, “Exact and approximative imaging methods for photoacoustic tomography using an arbitrary detection surface”,Phys. Rev. E, vol. 75, no. 4, pp. 046706, 2007, doi : 10.1103/PhysRevE.75.046706

  67. [67]

    Mathematical Foundations of Imaging, Tomography and Wavefield Inversion

    A.J. Devaney, “Mathematical Foundations of Imaging, Tomography and Wavefield Inversion”. Cambridge Uni- versity Press; 2012

  68. [68]

    Ray-based inversion accounting for scattering for biomedical ultrasound tomography

    A. Javaherian and B. Cox, “Ray-based inversion accounting for scattering for biomedical ultrasound tomography”, Inverse Problems, vol. 37, no.11, 115003, 2021

  69. [69]

    Hessian-free ray-born inversion for high-resolution quantitative ultrasound tomography

    A. Javaherian, “Hessian-free ray-born inversion for high-resolution quantitative ultrasound tomography”, 2023. https://arxiv.org/abs/2211.00316. 30 THIS ARTICLE HAS NOT BEEN REVIEWED YET

  70. [70]

    Representing arbitrary acoustic source and sensor distributions in Fourier collocation methods

    E. S. Wise, B. T. Cox, J. Jaros and B. E. Treeby, “Representing arbitrary acoustic source and sensor distributions in Fourier collocation methods”,J. Acoust. Soc. of Am. , vol. 146, no. 1, pp. 278-288, 2019

  71. [71]

    Field: A Program for Simulating Ultrasound Systems

    J.A. Jensen: “Field: A Program for Simulating Ultrasound Systems”, Paper presented at the 10th Nordic-Baltic Conference on Biomedical Imaging Published inMed. Biol. Eng. Comput. , pp. 351-353, Vol. 34, Supplement 1, Part 1, 1996

  72. [72]

    Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers

    J.A. Jensen and N. B. Svendsen: “Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers”,IEEE T-UFFC., Vol. 39, pp. 262-267, 1992

  73. [73]

    A Time-Space Decomposition Method for Calculating the Nearfield Pressure Generated by a Pulsed Circular Piston

    J. F. Kelly and R. J. McGough. “A Time-Space Decomposition Method for Calculating the Nearfield Pressure Generated by a Pulsed Circular Piston”,IEEE T-UFFC, Vol. l53, pp. 1150-1159, 2006

  74. [74]

    SIMUS3: An open-source simulator for 3-D ultrasound imaging

    D. Garcia and F. Varray, “SIMUS3: An open-source simulator for 3-D ultrasound imaging”,Comput. Methods. Programs. Biomed., Vol. 250, pp. 108169, 2024.https://doi.org/10.1016/j.cmpb.2024.108169