Translation of transient acoustic fields

Michael J. Carley

REVIEW 3 major objections 5 minor 16 references

Translation of transient acoustic fields

T0 review · 3 major / 5 minor · reviewed 2026-08-02 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A point-and-shoot translation method uses the Kirchhoff–Helmholtz integral to shift spherical harmonic acoustic expansions between spheres, converging to six digits and scaling as N^2.

desk verdict A genuinely direct time-domain SH translation operator with honest error tests and a plausible core, but the paper's N^2 complexity claim is internally inconsistent with its own algorithm and scaling data. read the letter →

arxiv 2607.15138 v1 pith:IA4LMUYD submitted 2026-07-16 physics.comp-ph

Michael J. Carley This is my paper · ORCID

classification physics.comp-ph

keywords acousticfieldtranslationsphericalharmonicsfastmultipolemethodKirchhoff–Helmholtzintegraltime-domainacousticspoint-and-shoottransientwavesLebedevquadrature

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to establish a practical translation operator for time-domain fast multipole methods: given a spherical harmonic expansion of transient acoustic data on one sphere, it computes the expansion on another, arbitrarily oriented sphere. The method splits the shift into a rotation, an axial translation via precomputed ring impulse responses, and a rotation back. The author shows numerically that the error converges to six digits for suitable sphere radii and that per-time-step cost scales approximately as N^2, where N is the spherical harmonic order. If correct, this removes a known bottleneck in transient FMM and makes time-domain translations comparable in cost to frequency-domain ones.

What carries the argument

The key object is the set of precomputed weighted impulse responses I^(m)_i (r1,r2,z) for each source-target ring pair, derived from the Kirchhoff–Helmholtz integral and advanced-time Lagrange weights. Along with the rotation matrices R = A Y_R (evaluating the expansion at rotated nodes), these operators turn the shift into a sequence of matrix-vector products: rotate, apply axial translation ring by ring, rotate back. The axial translation is O(N^2) because the azimuthal mode order is preserved, and only N source modes are transferred to N target modes per ring.

What would settle it

Take a source sphere with expansion order N=40, translate by the (2,1,0) shift with a rotation angle that does not align with quadrature symmetries, and compare the translated expansion against direct evaluation from the source for a point source. If the coefficient error does not reach about 1e-6 for sphere radius a=1.2*sqrt(3)d while increasing N to 48, the rotation or translation assumption fails. Alternatively, rotate a known band-limited spherical harmonic expansion by an arbitrary angle using R=AY_R and check whether the coefficients are preserved to machine precision as the number of qu

Watch

Extended reading notes

Core claim

The central claim is that the Kirchhoff–Helmholtz integral, applied to spherical harmonic surface data, gives a complete shift operator for transient fields: after rotating the source expansion so the translation axis is the z-axis, the integral factorizes over coaxial rings, each modal interaction being a precomputable impulse response. Applying the advanced-time Lagrange interpolation, these responses update the target modal coefficients at each time step. Two rotations, implemented by resampling at rotated quadrature nodes, handle arbitrary orientations. The method is shown to converge to six digits for local, upward, and downward shifts when sphere radii are chosen away from the closest

Load-bearing premise

The rotation of spherical harmonic expansions is done by resampling the function at rotated quadrature nodes with the same integration weights, which is exact only if the quadrature rule stays exact after an arbitrary rotation—something the paper does not prove or bound.

Editorial extensions

If this is right

Time-domain FMM shift operations become practical at expansion order up to about 48, with per-step cost scaling as N^2.1 and a small leading constant.
Six-digit accuracy holds for local, upward, and downward shifts in a parameter window of sphere radius roughly 1.15 to 1.25 times the minimum enclosing radius.
The very closest local interaction, the (2,0,0) shift, is not computed accurately and would need direct source evaluation; the method is limited to non-overlapping spheres.
The precomputed impulse responses are reusable across all time steps, making the method well suited to long transient simulations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the rotation step error is bounded as N grows, the same point-and-shoot structure could extend to other wave equations, such as vector or elastic waves, by replacing the scalar Kirchhoff–Helmholtz kernels with their tensorial analogues.
The observed N^2 scaling suggests the axial translation is essentially a convolution in azimuth; one might replace the ring loop with a fast transform to push the exponent lower for very high N.
The failure of the (2,0,0) case hints that accuracy is controlled by the ratio of sphere radius to box size; a user could choose sphere radii adaptively per interaction rather than globally.
One could test the rotation component in isolation by rotating a known band-limited function and measuring coefficient error as a function of quadrature order, which would separate rotation error from translation error.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Read it. The one thing to know: the core translation idea is real and the numerical section is mostly honest, but the abstract oversells the scaling.

What's new: the paper translates time-domain spherical-harmonic expansions between arbitrarily oriented, non-intersecting spheres directly, without Laplace transforms or concentric geometry. The Kirchhoff-Helmholtz ring-to-ring impulse responses are precomputed; the usual point-and-shoot rotations are then matrix-vector products per time step. Error tests cover local, upward, and downward shifts, and show roughly six-digit errors inside a parameter window. The paper is also honest that the closest local interaction fails and would have to be handled by direct evaluation. That is a useful, credible numerical primitive.

What worries me: the abstract's 'approximately N^2' claim does not survive contact with the paper. Section 2.7 says the rotation matrices are N(N+1) by N(N+1), which is O(N^4), and Figure 12 fits rotation time as N^3.5. Figure 12 fits axial translation as N^2.1, but Algorithm 4 is nested over source rings, target rings, and azimuthal modes. If the number of rings grows with N, that loop is closer to O(N^3). The measured slope describes the tested range, not an asymptotic law. The rotation itself is approximate—evaluation at rotated quadrature nodes—and no error bound is supplied, even though every translation uses two rotations. Less serious: only one source configuration is tested, the printed axial expressions are worth checking for a surface-area normalization factor, and code is only 'available upon request.'

Bottom line: the method is a legitimate contribution for time-domain FMM/BEM implementers and deserves a serious referee. The referee should ask for a corrected complexity statement, some rotation error analysis or experiments, and a second source distribution. I would cite it if I worked in this area, and I would bring it to reading group mainly to discuss the scaling gap.

Referee Report

3 major / 5 minor

Summary. The paper develops a method to translate transient acoustic field data between two non-intersecting spheres with arbitrary relative orientation. The field on each sphere is represented as a truncated spherical harmonic expansion. The translation uses a 'point-and-shoot' strategy: rotate the source expansion into a frame aligned with the source-to-target axis, perform an axial translation using precomputed Kirchhoff–Helmholtz ring impulse responses, and rotate the resulting target expansion back. Numerical tests compare the translated coefficients against direct evaluation from a point source and report close to six-digit accuracy for a range of parameters, except for the very closest interaction. The paper claims computational effort scales approximately as N^2, where N is the spherical harmonic order.

Significance. If the method is both accurate and truly O(N^2) per time step, it would provide a practical translation operator for time-domain fast multipole methods, filling a known gap. The paper has several genuine strengths: the Kirchhoff–Helmholtz derivation is standard, the error is measured against a direct analytic point-source evaluation with no fitted constants, and the method is presented in algorithmic detail with precomputation steps. The author's honest reporting of the failure for the closest interaction is also a positive feature. However, the complexity claim is contradicted by the paper's own analysis and timing data, and the rotation step lacks a rigorous accuracy guarantee, so the central practical claim needs significant revision.

major comments (3)

[Abstract; Sec. 2.7; Fig. 12] The central claim 'computational effort scales approximately as N^2' is inconsistent with the paper's own complexity analysis and measurements. Section 2.7 states that rotation matrices are of size N(N+1) × N(N+1), which implies O(N^4) work per matrix-vector multiplication, and Fig. 12 fits the rotation time as N^3.5. Even if the constant is small, the asymptotic scaling is not N^2. The abstract and conclusions should be corrected to state the actual scaling, or the rotation must be replaced by an algorithm with provably lower complexity.
[Sec. 2.4; Algorithm 4] The claimed O(N^2) for the axial translation is not supported by the algorithm as written. Algorithm 4 contains three nested loops over source rings i, target rings j, and azimuthal modes m. If the number of rings on each sphere grows with N, as required for accurate quadrature, this loop structure is O(R_s R_t N) = O(N^3), not O(N^2). The text in Sec. 2.7 counts only the 'generation of azimuthal modes' and 'transfer of N source modes' but omits the double sum over rings. The measured N^2.1 slope may reflect a fixed or slowly growing ring count over the tested range, but it cannot be taken as an asymptotic scaling law.
[Sec. 2.5; Eq. (12)] The rotation method R = A Y_R, following Lessig et al., is exact only if the quadrature rule used to build A remains exact after the arbitrary rotation of the evaluation nodes. For Lebedev-type quadratures, exactness relies on the specific node set and weights; rotating the nodes does not in general preserve exactness. Since every translation applies both a forward and an inverse rotation, any rotation error propagates into the six-digit convergence claim. The paper offers no error bound or separate convergence test for the rotation. I ask that the authors quantify the rotation error (e.g., by rotating a known expansion through a set of angles and measuring the coefficient error) and either show it is negligible in the reported parameter window or replace the rotation with an exact Wigner-D based method.

minor comments (5)

[Sec. 2.4, Eq. (15)] The notation in Eq. (15) uses Q(1)_m and N(1)_m for source ring modes but the target coefficients are denoted Q(2)_m and N(2)_m. This is understandable but should be made explicit, especially in the Fourier coefficient recovery step that follows.
[Algorithm 4] Typo: 'Algorith 1' should be 'Algorithm 1'. Also, the ordering of loops over i and j is not described; clarifying whether the impulse responses are stored as a flat array would improve reproducibility.
[Sec. 2.3, Eq. (14)] The number N_φ of azimuthal quadrature points is never specified or analyzed. Since the ring impulse responses are computed with this discrete sum, its accuracy is part of the overall convergence. Please state the values used in the tests and provide guidance on how N_φ should scale with N.
[Sec. 3, Fig. 12] The figure reports computation time for rotation and translation separately, but the translation time appears to include only the axial step, not the two rotations. Since the abstract claims N^2 for the whole operation, the total time (rotations plus translation) should be plotted and fitted separately, or the wording clarified.
[Sec. 5] Data availability: 'Code implementing the method... is available upon request' is a weaker statement than providing a permanent repository. Given that the numerical claims are central, a versioned code release would substantially strengthen the paper.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the translation method is derived from an external integral theorem and validated against direct evaluation, with no fitted parameters or self-referential definitions.

full rationale

The paper's central derivation uses the Kirchhoff–Helmholtz integral (Eq. 1) as an external acoustic theorem, then discretizes it into ring impulse responses (Eqs. 13–14) and a spherical-harmonic shift operator. No parameter is fitted to the target data; the only inputs are the source expansion coefficients, quadrature rules, and geometry. The claimed accuracy is measured against a direct evaluation from the point source (Eq. 17), so the 'prediction' of the target expansion is not equivalent to any fitted input. The self-citations [8, 10] provide prior numerical quadrature and advanced-time interpolation machinery, but these are supporting tools, not premises that assume the target result. The rotation step cites Lessig et al. [13], an independent external method, and although its accuracy is not bounded, that is a correctness/error-analysis concern rather than circularity. Likewise, the abstract's O(N^2) scaling claim is contradicted by the paper's own Section 2.7 rotation complexity and Figure 12 timing, but that is an internal consistency issue, not a definitional or fitted-input circularity. No step in the derivation reduces to its own input by construction, and no load-bearing uniqueness theorem is imported from the authors' prior work.

Assumptions & free parameters 4 free parameters · 7 assumptions · 0 invented entities

The method itself introduces no free physical parameters; the listed quantities are numerical discretization choices that the user must set, and the six-digit accuracy claim is explicitly conditional on them. The axioms are standard wave-theorem results plus domain assumptions of the FMM setting. No invented entities are postulated.

free parameters (4)

Sphere radius ratio a/sqrt(3)d
Accuracy window found in tests: six-digit convergence for 1.15 <~ a/sqrt(3)d <~ 1.25 (Figs. 7-11); closest interactions and small radii lose accuracy. This is a user-chosen parameter, not fitted to data.
Spherical harmonic order N = 32-40
Tests use N=32,36,40; error decreases with N until a time-discretization plateau (Fig. 8). Choice is a discretization parameter, not fitted.
Time discretization (number of time steps n_t) = 128,256,512
Error plateaus near 1e-6 due to time discretization and Lagrange interpolation; reported accuracy is conditional on n_t.
Lagrange interpolation order K = 4
Fourth-order interpolation used in all tests; accuracy of the advanced-time method depends on K.

assumptions (7)

standard math Kirchhoff-Helmholtz integral theorem for the scalar wave equation (Eq. 1) with retarded time.
Foundation of the translation; gives exterior/interior pressure from surface pressure and normal derivative. Standard theorem (Pierce [9]).
domain assumption All sources lie strictly inside source sphere S1, so surface data on S1 completely determine the field at points on/near S2.
FMM configuration; required for the K-H integral to be valid. Violated if sources leak outside S1.
domain assumption Source and target spheres do not intersect.
Required to avoid singularities in the K-H kernel R; the method cannot handle S1 and S2 overlapping. Explained in Sec. 2 and Fig. 1.
domain assumption For field points inside S1, the K-H sign is reversed.
One-line note after Eq. 1; needed in the downward-pass tests but not explicitly carried through Eqs. 13-15.
domain assumption Spherical harmonic expansion truncation at order N is adequate for the source data.
The method operates on truncated SH coefficients; convergence requires smooth field and sufficiently large N (see error discussion around Fig. 8).
domain assumption Lebedev quadrature and the rotation method R=AY_R accurately represent integration and rotation of SH expansions.
Rotation is performed by evaluating the function at rotated nodes and re-expanding (Lessig et al. [13]); exactness under arbitrary rotation is not proven.
standard math Advanced-time Lagrange interpolation accurately approximates retarded-time values.
Used in Sec. 2.2 for impulse responses; standard interpolation, accuracy depends on order K and time step.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Translation of transient acoustic fields." pith.science (2026). https://pith.science/paper/IA4LMUYD

@misc{pith2026260715138,
  author       = {Pith},
  title        = {Pith review of: Translation of transient acoustic fields},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IA4LMUYD}},
  note         = {Machine review of arXiv:2607.15138}
}

abstract

A method is presented for the translation of acoustic field data from a source to a target region. Field data are represented as spherical harmonic expansions on spheres surrounding the source and target regions respectively and expansions are translated using a ``point and shoot'' method using the Kirchhoff-Helmholtz integral to carry out an axial translation from one sphere to the other. The principal motivation for the method is its use in a time-domain Fast Multipole Method, and test cases reflective of this application are presented. The method converges to six digits for appropriate values of parameters and computational effort scales approximately as $N^{2}$ where $N$ is the order of spherical harmonic expansion for the field data.

Figures

Figures reproduced from arXiv: 2607.15138 by the authors.

**Figure 2.** Interpolation of time signal where i accounts for the time delay, j is incremented between time steps and wk are the weights of a Lagrange interpolation rule [16, 8] of order K. In the remainder of the paper, we write the weight vector as w(R/c). We will also need to evaluate the first and second time derivatives of the pressure, which we do using the Lagrange differentiation weights [16] with the corresponding weig… view at source ↗

**Figure 3.** Radiation from source ring (1) to target ring (2) [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

Figures from the paper (9 more)

**Figure 4.** Figure 4: Translation from source to target sphere [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]

**Figure 5.** Figure 5: Indexing of boxes for local shift operations: source box is at bottom left; [PITH_FULL_IMAGE:figures/full_fig_p014_5.png]

**Figure 6.** Figure 6: Configuration for upward and downward translation operations. For the up [PITH_FULL_IMAGE:figures/full_fig_p014_6.png]

**Figure 7.** Figure 7: Error in local shift operations against sphere radius with [PITH_FULL_IMAGE:figures/full_fig_p016_7.png]

**Figure 8.** Figure 8: Error in (2, 1, 0) shift operations against sphere radius, nt = 512; bullet N = 32, box N = 36, cross N = 40. 1 1.1 1.2 1.3 -2 -3 -4 -5 -6 -7 a/√ 3d log10 [PITH_FULL_IMAGE:figures/full_fig_p017_8.png]

**Figure 9.** Figure 9: Error in upward shift operation against sphere radius with [PITH_FULL_IMAGE:figures/full_fig_p017_9.png]

**Figure 10.** Figure 10: Error in downward shift operation against sphere radius with [PITH_FULL_IMAGE:figures/full_fig_p017_10.png]

**Figure 11.** Figure 11: Error against sphere radius for (2, 1, 0) (bullet), upward (box) and downward (cross) shifts with N = 32, nt = 512. 16 24 32 40 48 10−2 10−1 100 101 N t/s [PITH_FULL_IMAGE:figures/full_fig_p018_11.png]

**Figure 12.** Figure 12: Computation time for translation (bullets) and rotation operations (boxes) [PITH_FULL_IMAGE:figures/full_fig_p018_12.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references

[1]

Fast evaluation of asymptotic waveforms from gravitational perturbations.Classical and Quantum Gravity, 30, 2013

Alex G Benedict, Scott E Field, and Stephen R Lau. Fast evaluation of asymptotic waveforms from gravitational perturbations.Classical and Quantum Gravity, 30, 2013

2013
[2]

Field and Stephen R

Scott E. Field and Stephen R. Lau. Fast evaluation of far-field signals for time- domain wave propagation.Journal of Scientific Computing, 64:647–669, 2015

2015
[3]

P. A. Martin. Acoustic scattering by a sphere in the time domain.Wave Motion, 67:68–80, 2016

2016
[4]

P. A. Martin. The pulsating orb: solving the wave equation outside a ball.Pro- ceedings of the Royal Society of London. A., 472:20160037, 2016

2016
[5]

The solution of the scalar wave equation in the exterior of a sphere.Journal of Computational Physics, 274:191–207, 2014

Leslie Greengard, Thomas Hagstrom, and Shidong Jiang. The solution of the scalar wave equation in the exterior of a sphere.Journal of Computational Physics, 274:191–207, 2014

2014
[6]

N. A. Gumerov and R. Duraiswami.Fast Multipole Methods for the Helmholtz equation in three dimensions. Elsevier, Oxford, 2004. 19

2004
[7]

An enhancement of the fast time-domain boundary element method for the three-dimensional wave equation.Computer Physics Communications, 271:108229, 2022

Toru Takahashi, Masaki Tanigawa, and Naoya Miyazawa. An enhancement of the fast time-domain boundary element method for the three-dimensional wave equation.Computer Physics Communications, 271:108229, 2022

2022
[8]

M. J. Carley. Numerical evaluation of the Kirchhoff–Helmholtz integral outside a sphere.Journal of the Acoustical Society of America, 158(3):2622–2630, 09 2025

2025

Show all 16 references

[9]

A. D. Pierce.Acoustics: An introduction to its physical principles and applica- tions. Acoustical Society of America, New York, 1989

1989
[10]

Fast computation of time-dependent acoustic fields.Journal of the Acoustical Society of America, 140(5):3963–3970, November 2016

Michael Carley and Ghader Ghorbaniasl. Fast computation of time-dependent acoustic fields.Journal of the Acoustical Society of America, 140(5):3963–3970, November 2016

2016
[11]

V . I. Lebedev. Spherical quadrature formulas exact to orders 25–29.Siberian Mathematical Journal, 18:99–107, 1977

1977
[12]

C. H. L. Beentjes. Quadrature on a spherical surface. Technical report, Mathemat- ical Institute, University of Oxford, 2015

2015
[13]

Lessig, T

C. Lessig, T. de Witt, and E. Fiume. Efficient and accurate rotation of finite spher- ical harmonics expansions.Journal of Computational Physics, 231(2):243–250, 2012

2012
[14]

Casalino

D. Casalino. An advanced time approach for acoustic analogy predictions.Journal of Sound and Vibration, 261(4):583–612, 2003

2003
[15]

Source-time dominant aeroacoustics.Com- puters and Fluids, 33:791–800, 2004

Manuel Kessler and Siegfried Wagner. Source-time dominant aeroacoustics.Com- puters and Fluids, 33:791–800, 2004

2004
[16]

Trefethen

Jean-Paul Berrut and Lloyd N. Trefethen. Barycentric Lagrange interpolation. SIAM Review, 46(3):501–517, 2004. 20

2004

Pith tools

Reviewed August 2, 2026 · model on record in the stance chip above.

[1] [1]

Fast evaluation of asymptotic waveforms from gravitational perturbations.Classical and Quantum Gravity, 30, 2013

Alex G Benedict, Scott E Field, and Stephen R Lau. Fast evaluation of asymptotic waveforms from gravitational perturbations.Classical and Quantum Gravity, 30, 2013

2013

[2] [2]

Field and Stephen R

Scott E. Field and Stephen R. Lau. Fast evaluation of far-field signals for time- domain wave propagation.Journal of Scientific Computing, 64:647–669, 2015

2015

[3] [3]

P. A. Martin. Acoustic scattering by a sphere in the time domain.Wave Motion, 67:68–80, 2016

2016

[4] [4]

P. A. Martin. The pulsating orb: solving the wave equation outside a ball.Pro- ceedings of the Royal Society of London. A., 472:20160037, 2016

2016

[5] [5]

The solution of the scalar wave equation in the exterior of a sphere.Journal of Computational Physics, 274:191–207, 2014

Leslie Greengard, Thomas Hagstrom, and Shidong Jiang. The solution of the scalar wave equation in the exterior of a sphere.Journal of Computational Physics, 274:191–207, 2014

2014

[6] [6]

N. A. Gumerov and R. Duraiswami.Fast Multipole Methods for the Helmholtz equation in three dimensions. Elsevier, Oxford, 2004. 19

2004

[7] [7]

An enhancement of the fast time-domain boundary element method for the three-dimensional wave equation.Computer Physics Communications, 271:108229, 2022

Toru Takahashi, Masaki Tanigawa, and Naoya Miyazawa. An enhancement of the fast time-domain boundary element method for the three-dimensional wave equation.Computer Physics Communications, 271:108229, 2022

2022

[8] [8]

M. J. Carley. Numerical evaluation of the Kirchhoff–Helmholtz integral outside a sphere.Journal of the Acoustical Society of America, 158(3):2622–2630, 09 2025

2025