pith. sign in

arxiv: 2604.17617 · v1 · submitted 2026-04-19 · ⚛️ physics.optics · cond-mat.dis-nn

Universal compression of wave fields in weakly scattering media

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

classification ⚛️ physics.optics cond-mat.dis-nn
keywords wave field compressionweak scatteringFourier dispersion shelllossy compressionintensity correlationsoptical sensitivityelectromagnetic simulationson-shell representation
0
0 comments X

The pith

Wave fields in weakly scattering media compress by hundreds of times by confining their Fourier content to a thin dispersion shell, allowing accurate second-order calculations directly on the reduced data.

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

The paper shows that wave fields propagating in media where scattering is weak become confined in Fourier space to a narrow shell around the dispersion relation. This universal property stems from the long scattering mean free path relative to wavelength and produces two scale separations that enable lossy compression. The resulting scheme preserves coherent interference, so quantities like intensity and correlations can be obtained by convolution without full field reconstruction. Full-wave simulations in large volumes create cubic scaling storage demands that block ensemble work; this approach cuts data size by factors up to 380 while keeping errors small enough for practical use in optics and acoustics.

Core claim

OSCAR exploits the universal confinement of the Fourier representation of wave fields to a thin dispersion shell, a direct consequence of wave propagation when the scattering mean free path significantly exceeds the wavelength. The resulting compression ratio reflects on-shell confinement and excess Fourier volume from sub-wavelength discretization. Second-order quantities such as intensity, correlations, and optical sensitivity remain accurate when computed via convolution entirely in compressed space because coherent interference between independently compressed fields is preserved.

What carries the argument

On-shell confinement of the wave-field Fourier representation to a thin dispersion shell, which encodes the propagation physics of weak scattering and permits both data reduction and direct operations on the compressed representation.

If this is right

  • Compression ratios reach approximately 380 times with sub-percent field reconstruction error in 2D and 3D electromagnetic simulations.
  • Intensity, field correlations, and optical sensitivity can be obtained by direct convolution in compressed space without reconstructing the full fields.
  • The method supports routine ensemble averaging at system sizes relevant to biomedical optics, seismology, and underwater acoustics.
  • The approach is universal for any weakly scattering medium and independent of specific scatterer realizations.

Where Pith is reading between the lines

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

  • The same shell confinement should apply to acoustic or elastic waves in weakly scattering solids or fluids, suggesting direct transfer of the compression scheme.
  • Time-dependent or pulsed fields could be handled by stacking thin shells at different frequencies, provided the mean-free-path condition holds across the bandwidth.
  • Integration with existing mesh or wavelet compressors could further reduce storage for problems that combine sub-wavelength features with large overall scale.

Load-bearing premise

The scattering mean free path is much longer than the wavelength, which forces the Fourier content of the wave field onto a narrow shell.

What would settle it

A numerical simulation in a medium with mean free path comparable to wavelength that shows the Fourier spectrum of the field spreading well beyond the expected thin shell and producing large errors in reconstructed intensity would falsify the claimed compression accuracy.

Figures

Figures reproduced from arXiv: 2604.17617 by Alexey Yamilov, Pablo Jara.

Figure 1
Figure 1. Figure 1: Illustration of on-shell structure of wave-field in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of the OSCAR method. Left: real-space field [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Power fraction P/P∞ retained within the charac￾teristic shell as a function of normalized width α = ∆k ℓs. Solid line: analytical prediction, Eq. (7). Circles: numerical results averaged over Nrlz = 1000 realizations. The agree￾ment validates the Lorentzian spectral model and the sample￾independence of the energy normalization factor. ally define the intensity error ϵI = [PITH_FULL_IMAGE:figures/full_fig_… view at source ↗
Figure 4
Figure 4. Figure 4: Reconstruction accuracy in 2D. (a) Phase-sensitive [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Second-order quantities from compressed 2D fields. (a) Schematic of the convolution geometry: optical sensitivity [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Three-dimensional rendering of beam propaga [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Azimuthally averaged spectral density log10⟨|E˜x(k ′ )| 2 ⟩φ in the (k ′ z, k′ ⊥) plane for a representa￾tive 3D disorder realization after anisotropic shell masking. Dashed line: dispersion shell |k ′ | = k. Anisotropic mask, see Eq. (17), with α=4, ffwd=9, σfwd=0.1. The enhanced mask width near k ′ z ≈k captures the forward-scattering peak, while the narrow baseline width at large angles retains only wea… view at source ↗
Figure 8
Figure 8. Figure 8: Compression of 3D vector fields. Original (top) and reconstructed (bottom) cross-sections of log [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Reconstruction error ϵϕ, Eq. (13), vs boundary trim distance ∆trim normalized to shell width in the forward direc￾tion ∆kfwdℓs =∆k(1+ffwd)ℓs =20, 30, 40, and 60 in 3D. Error decreases monotonically with trimming as boundary Gibbs artifacts are excluded, and saturates beyond ∆trim·∆kfwd ≳1. The required trim is set by the shell thickness and is indepen￾dent of λ and ∆x. V. CONCLUSION We have introduced OSCA… view at source ↗
read the original abstract

Advances in computational methods have made full-wave simulations in large disordered media increasingly feasible, but the resulting field data, scaling with the cube of the ratio of system size to wavelength, creates a severe storage and post-processing bottleneck. Generic compression methods are sample-specific and preclude operations on compressed data. We introduce OSCAR (On-Shell Compression And Reconstruction), a physics-based lossy compression scheme for weakly scattering media. OSCAR exploits the universal confinement of the Fourier representation of wave fields to a thin dispersion shell, a direct consequence of wave propagation when the scattering mean free path significantly exceeds the wavelength. The resulting compression ratio reflects two distinct scale separations: on-shell confinement due to weak scattering, and the excess Fourier-space volume introduced by sub-wavelength discretization of the scatterers. Crucially, second-order quantities such as intensity, correlations, and (optical) sensitivity can be computed via convolution entirely in compressed space and remain accurate even when individual field reconstruction incurs appreciable error, because coherent interference between independently compressed fields is preserved. Numerical simulations of electromagnetic waves in 2D and 3D confirm compression ratios up to ${\sim}380\times$ with sub-percent field error, enabling routine ensemble studies at scales relevant to biomedical optics, seismology, and underwater acoustics.

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

1 major / 3 minor

Summary. The manuscript introduces OSCAR, a physics-based lossy compression scheme for wave fields in weakly scattering media. It exploits the confinement of the Fourier representation to a thin dispersion shell around |k|=k0, a direct consequence of the wave equation under the condition that the scattering mean free path greatly exceeds the wavelength. This yields compression ratios up to ~380x, with the key feature that second-order quantities (intensity, correlations, optical sensitivity) can be obtained via convolutions performed entirely in the compressed representation, remaining accurate even when individual field reconstructions have appreciable error. The approach is validated through 2D and 3D electromagnetic simulations.

Significance. If the accuracy claims for second-order observables hold, the method would provide a substantial practical advance for large-scale wave simulations by reducing storage demands while preserving the ability to compute ensemble statistics without full decompression. The parameter-free, derivation-based nature (no fitting or self-referential definitions) and explicit preservation of coherent interference are notable strengths that distinguish it from generic compressors.

major comments (1)
  1. Abstract: the claim that second-order quantities remain accurate via convolution in compressed space, even with appreciable field reconstruction error, lacks an explicit bound or scaling analysis for off-shell Fourier components. The Helmholtz equation implies that the field transform has non-zero support outside the dispersion shell whose amplitude scales with the scattering source strength; without a derivation showing when the resulting cross terms in quadratic forms fall below the reported sub-percent threshold, the central accuracy assertion is incompletely supported.
minor comments (3)
  1. The reported compression ratios and sub-percent errors would be strengthened by explicit details on how shell thickness and discretization parameters were selected, together with direct comparisons against standard compressors (e.g., wavelet or DCT-based) on identical wave-field datasets.
  2. Clarify the precise error metric used for the 'sub-percent field error' (L2 norm, maximum absolute error, or ensemble-averaged) and whether it is reported for the reconstructed fields or for the second-order quantities themselves.
  3. The abstract would benefit from a brief statement of the achieved compression ratios separately for the 2D and 3D cases to illustrate scaling with dimensionality.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive evaluation of the significance of OSCAR and for the constructive comment on the abstract. We address the concern below and have revised the manuscript to strengthen the theoretical support for the accuracy of second-order observables.

read point-by-point responses
  1. Referee: Abstract: the claim that second-order quantities remain accurate via convolution in compressed space, even with appreciable field reconstruction error, lacks an explicit bound or scaling analysis for off-shell Fourier components. The Helmholtz equation implies that the field transform has non-zero support outside the dispersion shell whose amplitude scales with the scattering source strength; without a derivation showing when the resulting cross terms in quadratic forms fall below the reported sub-percent threshold, the central accuracy assertion is incompletely supported.

    Authors: We agree that the abstract claim would benefit from an explicit scaling argument. In the revised manuscript we have added a short derivation in Section II (Theory) showing that, under the weak-scattering condition ℓ_s ≫ λ, the off-shell Fourier amplitude is bounded by O((k_0 ℓ_s)^{-1}) relative to the on-shell component. Consequently, the cross terms that appear in quadratic forms (intensity, correlations, sensitivity) are suppressed by O((k_0 ℓ_s)^{-2}). For the regime explored in our simulations (k_0 ℓ_s ≳ 30), this places the relative error comfortably below the reported sub-percent level. We have also inserted a concise reference to this bound in the abstract and updated the caption of Figure 3 to note the scaling. The numerical results already demonstrate that the convolution-based observables remain accurate even when individual-field reconstruction error is several percent; the added analysis now supplies the missing analytic justification. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained from wave equation

full rationale

The paper's core mechanism—on-shell confinement of the Fourier representation of wave fields—follows directly from the Helmholtz/wave equation under the weak-scattering condition (mean free path significantly exceeds wavelength), as stated in the abstract. OSCAR compression and the claim that second-order quantities (intensity, correlations, sensitivity) remain accurate via convolution in compressed space are presented as consequences of this physics-based scale separation and preservation of coherent interference, without any fitted parameters, self-referential definitions, or load-bearing self-citations. No step reduces by construction to its own inputs; the derivation chain is independent and externally grounded in standard wave propagation principles.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that weak scattering (mean free path >> wavelength) produces strict on-shell confinement in Fourier space; no free parameters are explicitly fitted in the abstract, and no new physical entities are postulated beyond the standard wave-propagation framework.

axioms (1)
  • domain assumption Wave propagation in media where the scattering mean free path significantly exceeds the wavelength confines the Fourier representation of the field to a thin dispersion shell.
    Invoked in the abstract as the direct physical basis for the compression scheme.

pith-pipeline@v0.9.0 · 5518 in / 1346 out tokens · 40886 ms · 2026-05-10T05:04:17.679576+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

54 extracted references · 54 canonical work pages

  1. [1]

    (5): P P∞ = 2 π tan−1(α) +O (kℓs)−1 ,(7) whereα= ∆k ℓ s

    Power normalization Retaining spatial Fourier components within a shell of width ∆kcentered at|k ′|=kcaptures only a part of the total spectral weight in Eq. (5): P P∞ = 2 π tan−1(α) +O (kℓs)−1 ,(7) whereα= ∆k ℓ s. Forα= 1, 50% of the spectral weight is retained;α= 3 recovers≈80%. Reconstructed fields must be rescaled by p P∞/Pto preserve total energy. Fo...

  2. [2]

    For a uniform grid with spacing ∆x, the continu- ous relationk 2 =P i k′2 i becomes [17] ω c 2 ¯ε∆x2 = 4 dX i=1 sin2 k′ i∆x 2 ,(8) which reduces to the continuum limit as ∆x→0

    Discrete dispersion relation Finite-difference discretization modifies the dispersion re- lation. For a uniform grid with spacing ∆x, the continu- ous relationk 2 =P i k′2 i becomes [17] ω c 2 ¯ε∆x2 = 4 dX i=1 sin2 k′ i∆x 2 ,(8) which reduces to the continuum limit as ∆x→0. At finite resolution, the iso-surface deviates from a sphere (circle in 2D) toward...

  3. [3]

    A Gaussian beam impinging on a semi-infinite scattering medium is an example of such a situation

    Boundary/source-induced anisotropy In finite systems with directional and/or spatially con- fined sources, the azimuthal distribution of spectral weight along the shell might become nonuniform. A Gaussian beam impinging on a semi-infinite scattering medium is an example of such a situation. In this case, fields near the injection boundary and/or in the do...

  4. [4]

    Since each factor ˜E(k′) (spatial Fourier transform) is confined to a shell of width∼1/ℓ s centered at|k ′|=k, their convolution should only be appreciable in the|q|≲∆k region

    Intensity as convolution As an example, we consider intensity of a scalar fieldE(r) found in a specific realization of disorder (no ensemble averaging) I(r) =E(r)E ∗(r).(10) In Fourier space, this product becomes a convolution: ˜I(q) = Z ddk′ (2π)d ˜E(k′) ˜E∗(k′ −q),(11) whereqis the spatial frequency of the intensity pattern andd=2,3 is the dimensionalit...

  5. [5]

    Optical sensitivity Optical sensitivityS(r 0) quantifies the response of de- tected flux to a localized absorptive perturbation at po- sitionr 0 inside the medium [39–41]. For a point pertur- bationδε(r) =i δϵ δV δ(r−r 0), it takes the form S(r0) =−k 2δVRe [E out(r0)E0(r0)],(12) whereE 0(r0) is the field atr 0 due to illumination from the source andE out(...

  6. [6]

    Ensemble averaging without decompression For disorder-averaged quantities I(r) or correlation func- tions E(r)E ∗(r′), the sample-independence of shell ge- ometry becomes essential. Because all compressed quan- tities in all disorder realizations share the same shell structure, ensemble operations can be performed on com- pressed data by averaging the con...

  7. [7]

    B. M. Notaroˇ s, F. P. Andriulli, and H. Bagci, Guest edi- torial frontiers in computational electromagnetics, IEEE Transactions on Antennas and Propagation71, 9175 (2023)

  8. [8]

    A. P. Mosk, A. Lagendijk, G. Lerosey, and M. Fink, Con- trolling waves in space and time for imaging and focusing in complex media, Nature Photonics6, 283 (2012)

  9. [9]

    J.-H. Park, Z. Yu, K. Lee, P. Lai, and Y. Park, Perspec- tive: Wavefront shaping techniques for controlling multi- ple light scattering in biological tissues: Toward in vivo applications, APL Photonics3, 100901 (2018)

  10. [10]

    Gigan, O

    S. Gigan, O. Katz, H. B. de Aguiar, E. R. Andresen, A. Aubry, J. Bertolotti, E. Bossy, D. Bouchet, J. Brake, S. Brasselet, Y. Bromberg, H. Cao, T. Chaigne, Z. Cheng, W. Choi, T. ˇCiˇ zm´ ar, M. Cui, V. R. Curtis, H. De- fienne, M. Hofer, R. Horisaki, R. Horstmeyer, N. Ji, A. K. LaViolette, J. Mertz, C. Moser, A. P. Mosk, N. C. P´ egard, R. Piestun, S. Pop...

  11. [11]

    H. Cao, A. P. Mosk, and S. Rotter, Shaping the propa- gation of light in complex media, Nature Physics18, 994 (2022)

  12. [12]

    A. F. Fercher, W. Drexler, C. K. Hitzenberger, and T. Lasser, Optical coherence tomography—principles and applications, Reports on Progress in Physics66, 239 (2003)

  13. [13]

    M. K. Kim, Principles and techniques of digital holo- graphic microscopy, SPIE Reviews1, 018005 (2010)

  14. [14]

    Kraszewski and J

    M. Kraszewski and J. Pluci´ nski, Coherent-wave monte carlo method for simulating light propagation in tissue, inProceedings of SPIE, Vol. 10063 (2016) p. 100631T

  15. [15]

    J. P. McCoy, L. E. Matloff, P. J. Vorster, D. X. Hammer, and R. D. Ferguson, Finite-difference time- domain (FDTD) optical simulations: A primer for life sciences and biophotonics practitioners, Micron142, 102989 (2021)

  16. [16]

    Virieux and S

    J. Virieux and S. Operto, An overview of full-waveform inversion in exploration geophysics, Geophysics74, WCC1 (2009)

  17. [17]

    Liu and Y

    Q. Liu and Y. J. Gu, Seismic imaging: From classical to adjoint tomography, Tectonophysics566–567, 31 (2012)

  18. [18]

    F. B. Jensen, W. A. Kuperman, M. B. Porter, and H. Schmidt,Computational Ocean Acoustics, 2nd ed. (Springer, 2011)

  19. [19]

    Y.-T. Lin, J. M. Collis, and T. F. Duda, A three- dimensional parabolic equation model of sound propaga- tion using higher-order operator splitting and Pad´ e ap- proximants, Journal of the Acoustical Society of America 132, EL364 (2012)

  20. [20]

    Graham, T

    I. Graham, T. Hou, O. Lakkis, and R. Scheichl, eds.,Nu- merical analysis of multiscale problems, Lecture Notes in Computational Science and Engineering, Vol. 83 (Springer, Berlin, Heidelberg, 2012)

  21. [21]

    Mache and I

    S. Mache and I. M. Vellekoop, Domain decomposition of the modified Born series approach for large-scale wave propagation simulations, Journal of Computational Physics550, 114619 (2026)

  22. [22]

    K. S. Yee, Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic me- dia, IEEE Transactions on Antennas and Propagation 14, 302 (1966)

  23. [23]

    Taflove and S

    A. Taflove and S. C. Hagness,Computational Electro- dynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, Boston, MA, 2005)

  24. [24]

    J. M. Ball and W. Li, Using high-resolution microscopy data to generate realistic structures for electromagnetic FDTD simulations from complex biological models, Na- ture Protocols19, 1348 (2024)

  25. [25]

    Krischer, A

    L. Krischer, A. Fichtner, S. Zukauskaite, and H. Igel, Large-scale seismic inversion framework, Seismological Research Letters86, 1198 (2015)

  26. [26]

    X. Zhao, A. Curtis, and X. Zhang, Physically structured variational inference for Bayesian full waveform inver- sion, Journal of Geophysical Research: Solid Earth129, e2024JB029557 (2024)

  27. [27]

    T. C. A. Oliveira, Y.-T. Lin, and M. B. Porter, Under- water sound propagation modeling in a complex shallow water environment, Frontiers in Marine Science8, 751327 (2021)

  28. [28]

    Lee and M

    D. Lee and M. H. Schultz,Numerical Ocean Acoustic Propagation in Three Dimensions(World Scientific, Sin- gapore, 1995)

  29. [29]

    Akkermans and G

    E. Akkermans and G. Montambaux,Mesoscopic Physics of Electrons and Photons(Cambridge University Press, 2007)

  30. [30]

    W. Wang, W. Zhang, and T. Lei, Compression of seismic forward modeling wavefield using tuckermpi, Computers & Geosciences172, 105298 (2023)

  31. [31]

    Y. Ni, M. A. Denolle, Q. Shi, B. P. Lipovsky, S. Pan, and J. N. Kutz, Wavefield reconstruction of distributed acous- tic sensing: Lossy compression, wavefield separation, and edge computing, Journal of Geophysical Research: Machine Learning and Computation1, e2024JH000247 (2024)

  32. [32]

    J. B. Muir and Z. Zhan, Seismic wavefield reconstruc- tion using a pre-conditioned wavelet–curvelet compres- sive sensing approach, Geophysical Journal International 227, 303 (2021)

  33. [33]

    Boehm, M

    C. Boehm, M. Hanzich, J. de la Puente, and A. Ficht- ner, Wavefield compression for adjoint methods in full- waveform inversion, Geophysics81, R385 (2016)

  34. [34]

    Kukreja, J

    N. Kukreja, J. H¨ uckelheim, M. Louboutin, J. Wash- bourne, P. H. J. Kelly, and G. J. Gorman, Lossy check- point compression in full waveform inversion: a case study with ZFPv0.5.5 and the overthrust model, Geo- scientific Model Development15, 3815 (2022)

  35. [35]

    Ishimaru,Wave Propagation and Scattering in Ran- dom Media(Academic Press, New York, 1978)

    A. Ishimaru,Wave Propagation and Scattering in Ran- dom Media(Academic Press, New York, 1978)

  36. [36]

    Sheng,Introduction to Wave Scattering, Localization and Mesoscopic Phenomena, 2nd ed., Springer Series in Materials Science, Vol

    P. Sheng,Introduction to Wave Scattering, Localization and Mesoscopic Phenomena, 2nd ed., Springer Series in Materials Science, Vol. 88 (Springer, 2006)

  37. [37]

    Lagendijk and B

    A. Lagendijk and B. A. van Tiggelen, Resonant multiple scattering of light, Physics Reports270, 143 (1996)

  38. [38]

    S. L. Jacques, Optical properties of biological tissues: a review, Physics in Medicine and Biology58, R37 (2013)

  39. [39]

    M. C. W. van Rossum and T. M. Nieuwenhuizen, Multi- ple scattering of classical waves: microscopy, mesoscopy, and diffusion, Reviews of Modern Physics71, 313 (1999). 13

  40. [40]

    J. W. Goodman,Speckle Phenomena in Optics: Theory and Applications(Roberts and Company, Englewood, CO, 2007)

  41. [41]

    J. W. Goodman,Introduction to Fourier Optics, 3rd ed. (Roberts and Company, Englewood, CO, 2005)

  42. [42]

    P. Jara, H. C. Lin, C. W. Hsu, H. Cao, and A. Yamilov, Simulation of coherent remission in planar disordered medium, 2023 International Applied Computational Electromagnetics Society Symposium (ACES) , 1 (2023)

  43. [43]

    Bender, A

    N. Bender, A. Goetschy, C. W. Hsu, H. Yilmaz, P. Jara, A. Yamilov, and H. Cao, Coherent enhancement of opti- cal remission in diffusive media, Proceedings of the Na- tional Academy of Sciences119, 1 (2022)

  44. [44]

    P. Jara, A. Goetschy, H. Cao, and A. Yamilov, Harnessing coherent-wave control for sensing applications, Physical Review Applied24, 054027 (2025)

  45. [45]

    S. R. Arridge, Optical tomography in medical imaging, Inverse Problems15, R41 (1999)

  46. [46]

    S. R. Arridge and J. C. Schotland, Optical tomogra- phy: forward and inverse problems, Inverse Problems25, 123010 (2009)

  47. [47]

    Durduran, R

    T. Durduran, R. Choe, W. B. Baker, and A. G. Yodh, Diffuse optics for tissue monitoring and tomography, Re- ports on Progress in Physics73, 076701 (2010)

  48. [48]

    Sarma, A

    R. Sarma, A. G. Yamilov, S. Petrenko, Y. Bromberg, and H. Cao, Control of energy density inside a disordered medium by coupling to open or closed channels, Physical Review Letters117, 086803 (2016)

  49. [49]

    Bender, A

    N. Bender, A. Yamilov, H. Yılmaz, and H. Cao, Fluc- tuations and correlations of transmission eigenchannels in diffusive media, Physical Review Letters125, 165901 (2020)

  50. [50]

    H.-C. Lin, Z. Wang, and C. W. Hsu, Fast multi-source nanophotonic simulations using augmented partial fac- torization, Nature Computational Science2, 815 (2022)

  51. [51]

    The Mill HPC Cluster, Missouri Univer- sity of Science and Technology (2024), https://scholarsmine.mst.edu/the-mill/1/

  52. [52]

    Osnabrugge, S

    G. Osnabrugge, S. Leedumrongwatthanakun, and I. M. Vellekoop, A convergent Born series for solving the inho- mogeneous Helmholtz equation in arbitrarily large media, Journal of Computational Physics322, 113 (2016)

  53. [53]

    Lagendijk, B

    A. Lagendijk, B. van Tiggelen, and D. S. Wiersma, Fifty years of Anderson localization, Physics Today62, 24 (2009)

  54. [54]

    Yamilov, S

    A. Yamilov, S. E. Skipetrov, T. W. Hughes, M. Minkov, Z. Yu, and H. Cao, Anderson localization of electromag- netic waves in three dimensions, Nature Physics , 1308– (2023)