REVIEW 3 major objections 5 minor 44 references
Tailoring reflectionless complex media for non-Abelian braiding of acoustic modes
T0 review · 3 major / 5 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read The paper claims that reflectionless complex media can be given an arbitrary unitary transmission matrix, so cascading them multiplies their T-matrices and performs targeted unitary operations—demonstrated as non-Abelian braiding of four ac
desk verdict Inverse-designed acoustic T-matrices work in experiment, but the 'arbitrary' claim is overbroad and the loss normalization needs scrutiny. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing identity is the multiplication rule for reflectionless sections: if the reflection blocks of each section vanish, the total transmission matrix of stacked media is the ordered product of the individual T-matrices, turning propagation into a matrix product. To realize each T, the optimization uses the generalized Wigner-Smith operator (GWSO): for every cylinder, the derivative of the scattering matrix under a small displacement equals i S Q_r, and Q_r is a boundary integral of pressure and tangential gradient over the cylinder's rigid surface, proportional to the acoustic radiation force. The gradient-descent objective g(r) = 1 - |Tr(S_obj^dagger S)|^2/(4N^2) drives the actua
What would settle it
Measure one optimized section's raw scattered power with a calibrated absolute-energy probe at 3.2 kHz, then repeat the measurement after adding a known absorbing liner to the waveguide walls; if the normalized T-matrix changes substantially or the recovered zero reflection fails to hold under increased loss, the singular-vector-invariance assumption is disproved.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that the transmission matrix of a reflectionless complex medium is not a nuisance to be averaged over but a designable object: when R = 0, cascading K sections yields T_total = T_K ... T_2 T_1, so each section is a linear operator on the four waveguide modes. The authors implement this by minimizing the distance g(r) between the actual and target S-matrices, with gradients supplied by the generalized Wigner-Smith operator, a boundary integral involving the acoustic radiation force on each cylinder. The optimized sections realize the braid generators sigma1, sigma2, sigma3 with a complex phase eta = pi/6 in a 2x2 block Y(eta), experimentally
Load-bearing premise
The experimental demonstration assumes that weak dissipation only reduces the singular values of the scattering matrix while leaving its singular vectors unchanged, and that configurations with dwell time below 0.9 times the ensemble average are therefore representative; if that loss-compensation assumption fails, the reported zero reflection and exact T-matrix entries could be artifacts of the normalization rather than properties of the physical medium.
Editorial extensions
If this is right
- Because reflectionless sections multiply their transmission matrices, any sequence of unitary operations on the N modes can be composed by stacking sections, making a stack of such media a linear wave-based processor.
- Non-Abelian braiding of acoustic modes can be produced by static passive scatterers, without requiring the adiabatic evolution or parity-time-symmetric Hamiltonians used in earlier waveguide braiding demonstrations.
- The same optimization can encode discrete Fourier transforms and two-qubit-style gates such as CNOT into the transmission matrix, offering a route to multimode acoustic signal processing and wave-based logic.
- Since the scheme relies only on linear interference and rigid scatterers, it should transfer to microwaves, elastic waves, and photonics wherever the corresponding Wigner-Smith operator can be computed.
Reading between the lines
- I would read the paper's broader bet as: if the T-matrix is fully programmable, then a stack of such sections is a programmable multimode linear processor; a natural next test is whether reconfiguring cylinder positions in real time can switch the unitary operation while preserving zero reflection, which the static demonstrations do not prove.
- The loss-compensation step is the part most worth probing: the claimed exactness of the T-matrix entries after singular-value normalization assumes losses only scale the singular values. A stronger test would compare normalized S-matrices against full-wave simulations that include thermoviscous losses in the air and boundary layers rather than a uniform complex frequency shift.
- The optimization demonstrably works for four modes and eighty cylinders; whether the scheme scales to larger modal counts is open, since the S-matrix dimension grows and the gradient landscape may develop more local minima. One could test convergence statistics for N = 8 and N = 16 targets.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes an inverse-design method, based on generalized Wigner–Smith operators, for engineering reflectionless complex acoustic media whose transmission matrices are prescribed unitary operations. Concatenating such media multiplies their T-matrices, enabling wave-based unitary transformations. The authors optimize 80-cylinder sections in a four-mode waveguide to realize the braid-group generators σ1, σ2, σ3 with a nontrivial phase η=π/6, experimentally demonstrate the non-commutativity of σ1σ2 versus σ2σ1 and the commutativity of σ1σ3 versus σ3σ1, and implement BNOT and CNOT gates. A numerical 8×8 DFT example is also presented. The central experimental claims rely on measured S-matrices that are post-processed by singular-value normalization to remove loss.
Significance. If the experimental results are valid, the work would substantially extend the range of operations achievable with multiple-scattering media: instead of treating disorder as a scrambler, the authors show that tailored complex media can perform targeted unitary operations, including non-commuting braid-group operations, without relying on adiabatic evolution. The acoustic GWSO derivation in Supplementary Note 4 is a useful contribution in its own right. The method is, in principle, transferable to other wave platforms. However, the paper does not provide a rigorous argument for the 'arbitrary unitary T-matrix' claim, and the experimental demonstration depends on a loss-removal normalization whose mathematical justification is questionable. The paper ships no public code or data, so the reproducibility of the optimization and measurements rests on the 'available upon request' statement.
major comments (3)
- [Supplementary Note 5, Eqs. (S47)–(S52)] The loss-removal normalization is not justified. Eq. (S49) assumes that in the basis that diagonalizes the Wigner–Smith operator Qω, the symmetric unitary S(ω0) also admits a Takagi decomposition with the same basis. This is equivalent to [Qω,S]=0, which does not hold for a generic smooth family of symmetric unitary matrices. For example, S(ω)=R(ω) diag(e^{iφ1},e^{iφ2}) R(ω)^T with a real rotation R and φ1≠φ2 gives [Q,S]≠0. Then the weak-loss term (δ/2)Q does not simply rescale the singular values of S; it rotates its singular vectors. Replacing the measured singular values by unity then yields a unitary matrix that is not the physical lossless S. Since the experimental T-matrices in Figs. 2(b2–d2), 4(a3–d3), and 5(a,c) are all shown after this normalization, the reported R≈0 and exact T-matrix entries could be artifacts. The authors should present raw, unnormalized S-matrices or provide
- [Abstract; Results, 'Realization of braiding'; Supp. Note 3] The headline claim of arbitrary unitary transmission matrices is not demonstrated. The manuscript shows inverse-design optimization for three braid generators, BNOT, CNOT, and one numerical 8×8 DFT. There is no mathematical or scaling argument that any N×N unitary T can be realized by a sufficiently large complex medium, nor a study of how optimization success depends on N or target structure. The word 'arbitrary' in the abstract should be relaxed to 'a broad class of target unitary operations' unless a constructive proof or systematic scalability evidence is added.
- [Methods, 'The effect of dissipation'; Supp. Note 5] The experimental configurations are post-selected from 40 optimizations by the dwell-time criterion τ/⟨τ⟩≤0.9, but the manuscript does not report how many of the 40 realizations satisfied this criterion or the distribution of τ. This selection makes it impossible to assess the success rate of the design method, and it concentrates the demonstration on cases in which the loss correction is smallest. The criterion should be reported with the number of candidates, and ideally raw S-matrix data for unselected configurations should be shown to establish that the method is not cherry-picked.
minor comments (5)
- [Eq. (4) and Fig. 3] The term 'braiding' may be misleading because no physical exchange of scatterers occurs; the operations are static matrices satisfying braid relations. Consider clarifying in the introduction or discussion that this is braiding in the sense of matrix representations of the braid group, not topological braiding of scatterers.
- [Figs. 2(b2–d2), 4(a3–d3), 5(a,c)] The experimental S-matrix figures show magnitude/color only. Including numerical values or a quantitative error metric (e.g., deviation from the target T-matrix) would help support claims such as 'exactly σ1(π/6)'.
- [Supp. Note 2] The Yang–Baxter relation is only validated numerically, and the text states 'Experimental validation is currently beyond our capability'. This limitation should be clearly stated in the main text as well.
- [Methods, 'The effect of dissipation'] The global phases φ1, φ2, φ3 are stated but the extraction procedure is not described. Please indicate how these phases are determined from the measured S-matrices.
- [Data/Code availability] The 'available upon request' statements are weak. Depositing the optimization code and the raw experimental S-matrix data would materially improve reproducibility, especially given the loss-normalization step.
Circularity Check
No significant circularity: the design target is an input specification, and the experimental S-matrices are independent measurements.
full rationale
The paper's core derivation (Eq. (2) / Methods Eq. (10)) is a theorem of the Redheffer star product under the condition R=0; this condition is a design target, not an output fitted to the final claim. The optimization objective g (Eq. (5)) is an inverse-design metric that prescribes Sobj; achieving g<1e-4 is a fabrication/control demonstration, not a prediction extracted from the measured data. The experimental S-matrices are reconstructed from independent microphone/loudspeaker measurements and therefore provide a genuine check of the optimized configurations. No experiment-fitted parameter is renamed as a prediction: the reported global phases are measured phases, and the loss coefficient enters only to justify a normalization. The one substantive risk is in Methods/'The effect of dissipation' and Supplementary Note 5 (Eqs. S47-S52): recovering the lossless S(omega0) by singular-value normalization requires simultaneous diagonalization of S and the Wigner-Smith operator Q_omega, which is assumed ('In the same WS basis...') but not generally guaranteed. If this assumption fails, the normalized matrices would not equal the physical lossless S and the displayed R/T values would be artifacts. That is a correctness/validity concern, not circularity, because the normalization is not fitted to the target T-matrices and does not by construction impose R=0. Self-citations (refs. 22, 30) provide prior context and a standard S-matrix reconstruction recipe; neither is load-bearing for the central claim. Hence no circular step is identifiable.
Assumptions & free parameters
free parameters (4)
- Braiding phase eta =
pi/6
- Loss coefficient delta =
0.0032 omega_0
- Dwell-time selection threshold =
tau/<tau> <= 0.9
- Scatterer positions (160 per section) =
optimized
assumptions (6)
- domain assumption Lossless reciprocal scattering: S is unitary and symmetric, S = S^T.
- domain assumption Weak losses only scale the singular values of S, leaving singular vectors invariant.
- ad hoc to paper The S-matrix of an 80-cylinder section is a smooth function of positions and gradient descent reaches S_obj.
- domain assumption Two-line, 16-microphone sampling uniquely determines the four-mode S-matrix.
- domain assumption The 2D finite-element model with rigid cylinders represents the 3-cm-high waveguide.
- ad hoc to paper Any unitary T can be realized by a sufficiently large complex medium.
Cite this review
Pith. "Pith review of Tailoring reflectionless complex media for non-Abelian braiding of acoustic modes." pith.science (2026). https://pith.science/paper/JINMYMGR
@misc{pith2026260215366,
author = {Pith},
title = {Pith review of: Tailoring reflectionless complex media for non-Abelian braiding of acoustic modes},
year = {2026},
howpublished = {\url{https://pith.science/paper/JINMYMGR}},
note = {Machine review of arXiv:2602.15366}
}
read the original abstract
Multiple scattering of sound and light can be tailored for diverse applications. Despite the great progress enabled by technologies such as time-reversal propagation and wavefront shaping, the full control of the transmission matrix remains a significant challenge. In this work, we propose a multi-scattering-based approach to design reflectionless complex media with an arbitrary unitary transmission matrix. As such, the perfect transmission of waves through such a medium performs a unitary operation. Based on this principle, we experimentally demonstrated braiding of multiple waveguide modes in an acoustic waveguide via multiple scattering and showed non-Abelian characteristics arising from the concatenation of distinct complex media. Furthermore, we show that the principle can be extended for realizing arbitrary unitary operations beyond braiding. Our scheme uses generalized Wigner-Smith operators to design the optimal acoustic complex media with near-arbitrary targeted functionalities. The scheme is generally applicable beyond acoustics, with broad implications to other wave types. Our results demonstrate unprecedented control over multiple-scattering waves and establish complex media as a viable route to modal control and operations on compact platforms, providing a new design paradigm for multimode, reconfigurable wave-based devices with potential applications in multiplexed communication, imaging, and the manipulation of quantum waves.
Reference graph
Works this paper leans on
-
[1]
Introduction to Wave Scattering, Localization, and Mesoscopic Phenomena
Sheng, P. Introduction to Wave Scattering, Localization, and Mesoscopic Phenomena. (Springer, Berlin ; New York, 2006)
2006
-
[2]
& Gigan, S
Rotter, S. & Gigan, S. Light fields in complex media: Mesoscopic scattering meets wave control. Rev. Mod. Phys. 89, 015005 (2017)
2017
-
[3]
Imaging and computing with disorder
Gigan, S. Imaging and computing with disorder. Nat. Phys. 18, 980–985 (2022)
2022
-
[4]
Cao, H., Mosk, A. P. & Rotter, S. Shaping the propagation of light in complex media. Nat. Phys. 18, 994–1007 (2022)
2022
-
[5]
& Fink, M
Derode, A., Roux, P. & Fink, M. Robust acoustic time reversal with high-order multiple scattering. Phys. Rev. Lett. 75, 4206–4209 (1995)
1995
-
[6]
& Prada, C
Fink, M. & Prada, C. Acoustic time-reversal mirrors. Inverse Problems 17, R1–R38 (2001)
2001
-
[7]
& Fink, M
Draeger, C. & Fink, M. One-channel time reversal of elastic waves in a chaotic 2D-silicon cavity. Phys. Rev. Lett. 79, 407–410 (1997)
1997
-
[8]
Vellekoop, I. M. & Mosk, A. P. Focusing coherent light through opaque strongly scattering media. Opt. Lett. 32, 2309 (2007)
2007
Show all 44 references
-
[9]
P., Lagendijk, A., Lerosey, G
Mosk, A. P., Lagendijk, A., Lerosey, G. & Fink, M. Controlling waves in space and time for imaging and focusing in complex media. Nat. Photonics 6, 283–292 (2012)
2012
-
[10]
M., Lagendijk, A
Vellekoop, I. M., Lagendijk, A. & Mosk, A. P. Exploiting disorder for perfect focusing. Nat. Photonics 4, 320–322 (2010)
2010
-
[11]
G., Petrenko, S., Bromberg, Y
Sarma, R., Yamilov, A. G., Petrenko, S., Bromberg, Y. & Cao, H. Control of energy density inside a disordered medium by coupling to open or closed channels. Phys. Rev. Lett. 117, 086803 (2016)
2016
-
[12]
Bender, N. et al. Depth-targeted energy delivery deep inside scattering media. Nat. Phys. 18, 309–315 (2022). 20
2022
-
[13]
& Lerosey, G
Del Hougne, P., Fink, M. & Lerosey, G. Optimally diverse communication channels in disordered environments with tuned randomness. Nat. Electron. 2, 36–41 (2019)
2019
-
[14]
Zhang, H., Wang, Q., Fink, M. & Ma, G. Optimizing multi-user indoor sound communications with acoustic reconfigurable metasurfaces. Nat. Commun. 15, 1270 (2024)
2024
-
[15]
& Lerosey, G
del Hougne, P. & Lerosey, G. Leveraging chaos for wave-based analog computation: Demonstration with indoor wireless communication signals. Phys. Rev. X 8, 041037 (2018)
2018
-
[16]
W., del Hougne, P., de Rosny, J., Lerosey, G
Matthès, M. W., del Hougne, P., de Rosny, J., Lerosey, G. & Popoff, S. M. Optical complex media as universal reconfigurable linear operators. Optica 6, 465 (2019)
2019
-
[17]
& Gigan, S
Rafayelyan, M., Dong, J., Tan, Y., Krzakala, F. & Gigan, S. Large-scale optical reservoir computing for spatiotemporal chaotic systems prediction. Phys. Rev. X 10, 041037 (2020)
2020
-
[18]
Dorokhov, O. N. On the coexistence of localized and extended electronic states in the metallic phase. Solid State Commun. 51, 381–384 (1984)
1984
-
[19]
A., Pereyra, P
Mello, P. A., Pereyra, P. & Kumar, N. Macroscopic approach to multichannel disordered conductors. Ann. Phys. 181, 290–317 (1988)
1988
- [21]
-
[22]
Chen, Z.-G., Zhang, R.-Y., Chan, C. T. & Ma, G. Classical non-Abelian braiding of acoustic modes. Nat. Phys. 18, 179–184 (2022)
2022
-
[23]
Zhang, X.-L. et al. Non-Abelian braiding on photonic chips. Nat. Photon. 16, 390–395 (2022)
2022
-
[24]
Chen, Y. et al. High-dimensional non-Abelian holonomy in integrated photonics. Nat Commun 16, 3650 (2025). 21
2025
-
[25]
Wang, X.-M., Xu, J., Wang, X., Li, Z. & Ma, G. Topological Braiding of Bloch Eigenmodes Protected by Non-Abelian Quaternion Invariants. Preprint at (2025)
2025
-
[28]
& Fleury, R
Orazbayev, B., Malléjac, M., Bachelard, N., Rotter, S. & Fleury, R. Wave-momentum shaping for moving objects in heterogeneous and dynamic media. Nat. Phys. 20, 1441–1447 (2024)
2024
-
[29]
On the relation of transmission‐line theory to scattering and transfer
Redheffer, R. On the relation of transmission‐line theory to scattering and transfer. J. Math. Phys. 41, 1–41 (1962)
1962
-
[35]
& Kumar, S
Zaheer, M., Reddi, S., Sachan, D., Kale, S. & Kumar, S. Adaptive Methods for Nonconvex Optimization. in Proceedings of the 32nd International Conference on Neural Information Processing Systems 9815–9825 (Curran Associates Inc., Red Hook, NY, USA, 2018). 22
2018
-
[37]
virtual velocity
Gérardin, B., Laurent, J., Derode, A., Prada, C. & Aubry, A. Full transmission and reflection of waves propagating through a maze of disorder. Phys. Rev. Lett. 113, 173901 (2014). 1 SUPPLEMENTARY INFORMATION Tailoring multiple scattering acoustic media with perfect transmissio...
2014
-
[38]
Zhang, H., Delaporte, T. & Ma, G. Scattering matrices of two-dimensional complex acoustic media. npj Acoust. 1, 17 (2025)
2025
-
[39]
Ambichl, P. et al. Focusing inside disordered media with the generalized wigner-smith operator. Phys. Rev. Lett. 119, 033903 (2017)
2017
-
[40]
Horodynski, M. et al. Optimal wave fields for micromanipulation in complex scattering environments. Nat. Photonics 14, 149–153 (2020)
2020
-
[41]
& Davy, M
Horodynski, M., Kühmayer, M., Ferise, C., Rotter, S. & Davy, M. Anti-reflection structure for perfect transmission through complex media. Nature 607, 281–286 (2022)
2022
-
[42]
Faugno, W. N. & Ozawa, T. Non-abelian braiding of classical waves
-
[43]
Feijo, G. R. & Pinsky, P. M. Shape sensitivity calculations for exterior acoustics problems. Engineering Computations 18, 376–391 (2001)
2001
-
[44]
Shape sensitivity in 2D acoustic scattering
Sjöberg, M. Shape sensitivity in 2D acoustic scattering. (ETH Zürich, 2023)
2023
-
[45]
Hiptmair, R. & Li, J. Shape derivatives for scattering problems. Inverse Probl. 34, 105001 (2018)
2018
-
[46]
Andrade, M. A. B., Pérez, N. & Adamowski, J. C. Review of progress in acoustic levitation. Braz J Phys 48, 190–213 (2018). 16
2018
-
[47]
Acoustofluidics 7: The acoustic radiation force on small particles
Bruus, H. Acoustofluidics 7: The acoustic radiation force on small particles. Lab Chip 12, 1014 (2012)
2012
-
[48]
Wei, W., Thiessen, D. B. & Marston, P. L. Acoustic radiation force on a compressible cylinder in a standing wave. J. Acoust. Soc. Am. 116, 201–208 (2004)
2004
-
[49]
Beenakker, C. W. J. & Brouwer, P. W. Distribution of the reflection eigenvalues of a weakly absorbing chaotic cavity. Physica E: Low-dimensional Systems and Nanostructures 9, 463–466 (2001)
2001
-
[50]
Smith, F. T. Lifetime Matrix in Collision Theory. Phys. Rev. 118, 349–356 (1960)
1960
-
[51]
R., Mao, Y
Patel, U. R., Mao, Y. & Michielssen, E. Wigner–Smith time delay matrix for acoustic scattering: Theory and phenomenology. J. Acoust. Soc. Am. 153, 2769 (2023)
2023
-
[52]
Ikramov, Kh. D. Takagi’s decomposition of a symmetric unitary matrix as a finite algorithm. Comput. Math. Math. Phys. 52, 1–3 (2012)
2012
-
[53]
Patel, U. R. & Michielssen, E. Wigner-Smith Time-Delay Matrix for Electromagnetics: Theory and Phenomenology. IEEE Trans. Antennas Propagat. 69, 902–917 (2021)
2021
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.