The surface states of intrinsic higher order topological phases are protected by the spatial symmetries of a finite sample. This property makes the existing scattering theory of topological invariants inapplicable: the scattering geometry is either incompatible with the symmetry or does not probe the bulk topology. We resolve this obstacle by using a symmetric scattering geometry that probes transport from the inside to the outside of the sample. We demonstrate that the intrinsic higher order topology is captured by the flux dependence of the reflection matrix. Our finding follows from identifying the spectral flow of a flux line as a signature of higher order topology. We show how this scattering approach applies to several examples of higher order topological insulators and superconductors. Our theory provides an alternative approach for proving bulk–edge correspondence in intrinsic higher order topological phases, especially in the presence of disorder.

We demonstrate that the classical dynamics influence the localization behaviour of Majorana wavefunctions in Majorana billiards. By using a connection between Majorana wavefunctions and eigenfunctions of a normal state Hamiltonian, we show that Majorana wavefunctions in both p-wave and s-wave topological superconductors inherit the properties of the underlying normal state eigenfunctions. As an example, we demonstrate that Majorana wavefunctions in topological superconductors with chaotic shapes feature quantum scarring. Furthermore, we show a way to manipulate a localized Majorana wavefunction by altering the underlying classical dynamics using a local potential away from the localization region. Finally, in the presence of chiral symmetry breaking, we find that the Majorana wavefunction in convex-shaped Majorana billiards exhibits caustics formation, reminiscent of a normal state system with magnetic field.