Reflection resonances in surface-disordered waveguides: strong higher-order effects of the disorder

We study coherent wave scattering through waveguides with a step-like surface disorder and find distinct enhancements in the reflection coefficients at well-defined resonance values. Based on detailed numerical and analytical calculations, we can unambiguously identify the origin of these reflection resonances to be higher-order correlations in the surface disorder profile which are typically neglected in similar studies of the same system. A remarkable feature of this new effect is that it relies on the longitudinal correlations in the step profile, although individual step heights are random and thus completely uncorrelated. The corresponding resonances are very pronounced and robust with respect to ensemble averaging, and lead to an enhancement of wave reflection by more than one order of magnitude.

hand, is due to the correlations between the squares of the slopes (squares of the derivatives) of the profiles at the points x and x′. In most theoretical studies, this -correlator is, however, neglected since it constitutes a higher-order term in the weak disorder expansion where the disorder amplitude is the relevant expansion parameter. In our article we will provide conclusive evidence that this term, although being of higher order, can dominate the transmission through a surface-disordered waveguide and that it needs to be taken into account in a comprehensive description.
The first numerical and experimental indications that the -correlator plays, indeed, an important role have been put forward in a recent study on a specific waveguide geometry which was designed to highlight the presence of this new term [33]. Here we go an important step further by demonstrating that the influence of this correlator shows up not just for carefully chosen waveguide geometries, but in a quite general class of waveguides. In particular, we will show that waveguides with a step-like surface disorder which have been well studied by the community yield unambiguous and very pronounced signatures for the influence of the -correlator which, to the best of our knowledge, have so far been overlooked. In these waveguides (see figure 1), the surface disorder features steps of random height (in the direction transverse to propagation) and of constant width (in longitudinal direction). Such waveguides have been considered in quite a few recent studies [20,26,[36][37][38][39], as they are attractive model systems both for an experimental implementation as well as for a numerical computation. This is because waveguides with the above specifics can be easily built up by combining a series of rectangular waveguide stubs, each of which has no surface disorder but a randomly chosen height. Our analysis will show that in these concatenated systems the -correlator gives rise to well-defined resonances in the reflection coefficients, which are perfectly reproduced in a corresponding numerical study. At these resonant values the -correlator may strongly dominate over the lower-order -correlator such that a conventional description breaks down here. To show this not only by numerical evidence, but also by the corresponding analytical expressions, we significantly expand the existing theoretical framework presented in [34,35]. This is mainly because the surface derivatives which enter the -correlator diverge at the steps in the surface profile and thus require a special treatment. Another important extension of the theory which we take into account is due to multiple scattering events between the propagating modes in the waveguide which yield a significant contribution beyond the single-scattering terms that have been considered so far. In this sense, our combined analytical-numerical study Figure 1. Illustration of the considered surface-disordered waveguide of length L attached to semi-infinite collinear leads of width d. The step-like surface disorder is characterized by a constant step width Δ and a maximum disorder-strength δ, respectively (see text for details). Flux is injected from the left. An example of an antisymmetric geometry is shown in which the upper and the lower surface disorder are identical. not only reveals a new effect, but also contributes to an extension of the underlying theory to the point where the analytical formulas that we derive provide predictions which quantitatively match with the numerical calculations that we perform independently.

Model
We consider a simple, non-trivial model consisting of a quasi-1D corrugated waveguide (or conducting wire) with discrete steps in the surface profile. This rough waveguide of length L and average width d L ≪ is attached to infinite leads of width d on the left and right (see Altogether three different cases will be considered in terms of the symmetries of the boundary profiles with respect to the horizontal center axis at y = 0: Following the assumptions adopted in a few recent papers [20,26,[36][37][38][39], the functions x ( ) i σξ are chosen as sequences of horizontal steps of constant width Δ and random heights, around the upper (lower) boundary of the attached leads. In our numerical analysis we set d = 1 and 0.04 δ = , resulting in a variance of the disorder, 12 , which is small compared to the width of the waveguide, d σ ≪ . Note that we have realized, in the above way, a scattering system which is truly random yet features very strong spatial correlations in its surface disorder since the waveguide exhibits a potential step at each integer multiple of the step-width Δ.

Analytical method
According to the theory developed in [34,35], the correlations in the surface disorder enter the scattering properties of the system through two independent correlators. The first one is the binary correlator of the surface profile, which contains contributions only from the amplitude x ( ) ξ and the derivative of the surface profile x ( ) ξ′ . Correspondingly, the scattering mechanism that this correlator gives rise to is referred to as the amplitude-gradient-scattering (AGS) mechanism.
The other correlator contains scattering contributions which are independent of those in equation (4) and which are related to the square of the profileʼs derivative, x ( ) 2 ξ′ , in an effective potential description (see details in [34,35] . The corresponding scattering process is thus referred to as the square-gradient-scattering (SGS) mechanism. We emphasize here that the validity of the identity used in different contexts (see, e.g., [33][34][35]) is restricted to Gaussian random processes and cannot be applied for the present step-like surface profiles. Indeed, as we will see below, this simplification would lead to a severe underestimation of the SGS mechanism in the present context.
In our further analysis it will not be binary correlators themselves which will be the key quantities, but rather their Fourier transforms ( ) W k x and ( ) x ik x For the scattering system in figure 1 the -correlator can be obtained analytically, which is strongly peaked for surface points x and x′ which are closer to each other than the step width in the disorder, . For completeness and since it is a key parameter in [34,35], we want to stress the fact that, if defining the correlation length R as the variance of the binary correlator ( ) , the step width Δ and correlation length R are directly linked with each other, R 6 Δ = .
For the sake of simplicity we will use the step width Δ in all expressions in the following since it represents the quantity which we tune in our simulations and which is therefore the more natural parameter in our system.
The Fourier transform of ( ) then yields the analytical expression for the roughness-height power spectrum ( )  (8) and (9) implicitly assume that the rectangular steps in the profile boundary can be perfectly resolved by the scattering wave. However, due to the finite wavelength at which the scattering process takes place, the resolution of the surface profile also will always be finite. To accommodate this limited resolution, we introduce an effective smearing of the step profiles based on a Fermi-function ( ) ⎤ ⎦ (see the appendix for more details). The smoothness of this function is governed by the parameter ρ which leads to a smearing of a step profile over a region x 12 Δ ρ ≈ (see figure A1 in the appendix for a corresponding illustration). If we now estimate that a scattering wave with a wavelength λ is associated with a resolution of x 2 Δ λ ≈ , we obtain for the smearing parameter 0.03 Employing this value for all further calculations, a comparison with the numerical data suggests that this simple estimate already captures our simulations remarkably well. An exception is that in symmetric waveguides a reduced value of 0.01 ρ = yields better agreement. Based on these observations, we expect our theoretical framework to be applicable also to waveguides with a smoothened step boundary (provided that the smoothing is on or below the same scale as the wavelength).
When incorporating the smoothness of the steps into the roughness-height power spectrum given by equation (9), we can again obtain a simple analytical expression, which takes the following form (see the appendix for details), π ρ π ρ ≈ , yielding the result already obtained for infinitely sharp steps, equation (9).
The above approach involving a smearing of the step disorder turns out to be essential when considering the roughness-square-gradient power spectrum ( ) S k x . This is because, without the smearing, the corresponding expressions would diverge, as can easily be understood from the fact that the gradient turns into a delta function at the position of a step when an infinite resolution is assumed. This divergence is, however, conveniently tamed through the above procedure involving the Fermi-function, yielding the following analytical expression for In addition to the smearing parameter ρ, the above expression also contains the integer number N eff , which determines the number of steps N 2 1 eff + that are effectively involved in the scattering process. The notion of an effective number has been introduced here to take into account that the total number of steps in the waveguide, N L 2 1 Δ + = , is typically significantly larger than the value N 2 1 eff + which we find to reproduce our data. This difference, N N eff ≪ , can be attributed to the finite penetration depth of the propagating wave as a result of which the effective longitudinal dimension of the waveguide is greatly reduced. We shall thus determine the quantity N eff through a direct comparison with the numerical data to be presented below. Also note that we have used ensemble averaging for the derivation of the above equations (10) and (11) (to ensure convergence of equation (A.3) in the appendix) [40]. Recent work demonstrates, however, that an application of the predictions following from the two different correlators above also yields good quantitative agreement for individual disorder realizations as in single disordered waveguides [33].
A direct comparison of the expressions for the two correlators in equations (10) and (11) provides the insight that the SGS term ( ) S k x becomes large at exactly the same points at which the AGS term ( ) with M integer, the SGS term will thus dominate over the AGS term. As we will demonstrate below, this fact provides the key element for the occurrence of the pronounced resonances in reflection that we observe, and we will discuss how this resonance condition is realized for different symmetry classes. Note that these dominant SGS contributions in would be suppressed if we applied the customary approximation (used, e.g., in [33][34][35]) that the defining expression for the SGS term, , can be replaced by the simplified term With the above expressions equations (10) and (11) we now have the relevant quantities at hand for setting up the perturbation theory analysis of scattering in surface-disordered waveguides. For this analysis to be applicable, the perturbation induced by the surface disorder has to be weak, resulting in the following independent requirements, Here, L n is the partial attenuation length of the nth incoming mode (from the left) which takes into account both the scattering in forward direction (to the right) and in backward direction (to the left). The cycle length n Λ is the distance between two successive reflections of the nth mode from the unperturbed surfaces. Under the conditions in equation (12) the waves are weakly attenuated over the correlation length R, the step width Δ, and over the cycle length n Λ . Clearly, the correlation length must be smaller than the waveguide length, R L ≪ . When applying, in the above limit of weak disorder, the perturbative treatment following [34,35], we obtain the mode-specific inverse attenuation lengths for scattering from any incoming mode n into any mode n′ [35], All L nn′ can be decomposed into backward (b) and forward (f) scattering contributions as well as into terms which are associated with the AGS and SGS mechanism of surface scattering. In their full, detailed form we thus obtain for the terms in equation (13), Here the factors A nn′ and B nn′ depend on the symmetry between the two profiles x ( ) ξ ↑ and x ( ) ξ ↓ (see table 1), and the terms depending on k k n n + ′ contribute to backward scattering whereas those depending on k k n n − ′ result in forward scattering. The overall attenuation length of mode n can be obtained by means of the sum over all corresponding partial inverse mode-specific As one can see from equations (14) and (15)  S k x depend differently on the external parameters, in particular, on the wavenumber k x and on the module width Δ. We thus may arrive at the situation in which at specific values of the wavenumber the SGS-term in equation (15) ( 4 σ ∝ ) can be comparable to (or even larger than) the AGS-term in equation (14) ( 2 σ ∝ ). In particular, the points discussed above, where a peak value in ( ) coincides with a zero of ( ) W k x , can be expected to lead to interesting transmission characteristics.
To test this scenario explicitly, we performed extensive numerical simulations on transport through surface-disordered waveguides of all three symmetry classes.

Numerical method
For these numerical simulations we employ the efficient 'modular recursive Greenʼs function method' (MRGM) [17,41] to solve the Schrödinger equation for the Hamiltonian (in atomic units), on a discretized grid. The potential term V defines the surface potential, which is infinite outside the waveguide and flat (V = 0) inside, corresponding to hard-wall boundary conditions. Since the scattering problem in equation (16) is equivalent to the Helmholtz equation, our approach is not only suitable for electronic systems but can, e.g., also be applied to microwave systems as in [33,42], or quite generally to systems which satisfy a Helmholtz-like equation. The MRGM is particularly advantageous for the present setup since the vertical steps in the disorder profile allow us to assemble the waveguide by connecting a large number of rectangular elements, which will be referred to as 'modules'. These modules are chosen to have equal width Δ, but different heights. The computation is based on a finite-difference approximation of the Laplacian and proceeds such that we first calculate the Greenʼs functions for a number of modules with different heights. These Greenʼs functions are then connected to each other by way of a matrix Dyson equation [41]. It is the different heights of the modules and additionally introduced random vertical shifts between them that give rise to the desired random sequence of vertical steps in the surface profile. To satisfy the additional symmetry imposed on the waveguide, we arrange the modules in order to respect this specific symmetry.
The key element of our numerical approach is an 'exponentiation' algorithm [20] which allows us to simulate transport through extremely long waveguides at moderate numerical costs. Rather than connecting individual modules with each other until the length of the waveguide is reached, we first connect several sequences of randomly assembled modules. In a subsequent step these 'supermodules' are then randomly permuted and connected to each other to form a next generation of supermodules. Continuing this iterative procedure allows us to obtain the Greenʼs functions of waveguides with a length that increases exponentially with the number of generations. For waveguides of moderate lengths, we tested this supermodule technique against the conventional approach where the modules are assembled one after the other. We found that the disorder-averaged Greenʼs functions obtained in these two ways do not show any noticeable difference from each other [20].
To calculate the desired transmission ( ) t nn′ and reflection amplitudes ( ) r nn′ for incoming flux from the left lead, we project the Greenʼs function at the scattering wavenumber k onto the fluxcarrying lead modes in the left and right lead, respectively.

Comparison between analytical and numerical results
In order to compare the analytical predictions of equations (14) and (15)  and N 10 6 = (nonsymmetric waveguide), respectively. To suppress effects which are due to individual disorder realizations, we also average the transmission over 100 (symmetric and  Our fits are then performed with the disorder-averaged transmission curves (details are provided below). To keep the system at a manageable degree of complexity and to perform a direct comparison with equations (14) and (15) we restrict ourselves to the regime of two open waveguide modes, N 2 d = , by choosing the wavenumber k to be fixed at the value k d 2.55 π = . By varying the step width Δ in the surface disorder incrementally, we numerically scan through the module width dependence of the transmission (at each value of Δ an ensemble average over 50-100 waveguide realizations is performed).
We will now discuss the disordered waveguides with different symmetry separately, as both the predictions from equations (14) and (15) as well as the procedure to extract the attenuation lengths are specific for each symmetry.

Symmetric profiles
In symmetric waveguides the up-down symmetry of the entire scattering structure, , results in the fact that modes of different symmetry cannot scatter into each other. For the two-mode waveguide considered here, this means that the two modes n 1, 2 = scatter fully independently of each other with only intra-mode scattering (with n n = ′) being relevant and inter-mode scattering (with n n ≠ ′) being entirely absent. Correspondingly, the only scattering mechanism that attenuates an incoming wave in mode n is back-scattering into the same mode (forward-scattering in the same mode does not attenuate the mode and intermode scattering is forbidden). For our analysis we therefore need to consider only the intramode back-scattering (b) length L ( ) nn b which follows from equations (14) and (15)   where W ( ) · and S ( ) · are defined by equations (10) and (11), respectively. Due to the decoupling of the two modes, we are here in the 1D limit of single-channel scattering where all modes are localized and diffusion is absent (as in 1D bulk scattering systems [30]), resulting in an exponential decrease of the transmission T(L) with waveguide length L, For 1D scattering the localization length ξ is related to the mean free path as follows l 2 ξ = [30]. Identifying the mean free path for each mode with the specific backward scattering length L ( )  (17) and (18). The corresponding results for L 1 ( ) nn b as a function of Δ are shown in figure 3(a). We also plot the theoretical predictions given by equations (17) and (18), respectively, as well as the AGS terms L 1 ( alone. The agreement we find between the AGS terms and the numerical calculations is already remarkably good for most of the chosen parameters, such that the SGS contributions can be easily identified to be dominant at those specific parameter values where deviations from the AGS predictions occur (see vertical arrows in figure 3(a)). In full agreement with our theoretical analysis, we find that the values of the step width Δ where this happens are determined by the resonance condition k M 2 2 n Δ π = (with M an integer), which we already identified earlier as those points where the contribution of the AGS terms vanishes while the SGS terms are maximal. Note that this condition leads to different resonance values for each of the two modes with n = 1, 2, At these well-defined values we not only find that the theory solely based on the AGS terms deviates from the numerics (see figure 3(b)), but that the additional SGS terms fill the missing gaps in the theory very well in terms of resonant contributions to the inverse attenuation lengths L 1 nn (see figure 3(a)). Since maxima in the inverse attenuation length correspond to maxima in the reflection (i.e., minima in the transmission) we may thus conclude that the SGS mechanism leads to reflection resonances in the systems under study. While these resonances are already ). An automated fitting procedure yields mode-specific attenuation lengths L n which show excellent agreement with our analytical estimates (see below). In the above figures, 0.01 σ ≈ and k 2.55 π = have been employed, respectively.
clearly discernible in the symmetric waveguides, we will find that they are even more pronounced in the antisymmetric waveguides that we investigate in the next section.

Antisymmetric profiles
In the case of antisymmetric waveguide profiles, we have , i.e., the waveguide width is constant throughout the waveguide (see figure 1). The situation is more complicated than in symmetric waveguides as inter-mode scattering is allowed here. A proper description of transmission thus has to incorporate both the intra-and the inter-mode scattering contributions. For mode-specific values of the transmission, we again have to ask which scattering events contribute: consider, e.g., the transmission of a mode into itself, T nn . In a perturbative treatment this quantity is determined by all scattering mechanisms that scatter the incoming mode n into the other mode or reverse its direction of propagation. This happens both through forwardscattering from mode n into the second available mode n n ≠ ′ as well as through backwardscattering into any of the two modes n 1, 2 = ′ . The attenuation length extracted from the transmission T nn thus has to be compared to the predictions for the attenuation length L nn , which is given as The remaining question at this point is how to extract the attenuation lengths L nn from the numerical data for T nn when modes do not just localize as in the symmetric case. In the presence of inter-mode scattering, the wave injected into a disordered waveguide first propagates  (17) and (18) (both ) including both the AGS and the SGS terms. Note the very good agreement which we find between the numerical data and the analytical theory, in particular also for those resonant values k M 2 2 n Δ π = where the SGS contributions dominate (marked by arrows). Panel (b) shows the AGS predictions alone, i.e., when the SGS mechanism is omitted (both ). For all data shown the following parameter values were used: ballistically, then scatters diffusively, and eventually localizes at very long waveguide lengths. For the two-mode waveguide considered here, the diffusive regime is, however, not well pronounced such that the crossover region between ballistic scattering and localization is comparatively narrow. Since, additionally, in the localized regime the mode with the higher localization length ξ always dominates [20], extracting mode-specific attenuation lengths is best achieved in the ballistic regime where the transmission decreases linearly with the system length L, T L L 1 nn nn ≈ − . We will use this relation to extract the attenuation lengths L nn from the disorder-averaged numerical transmission values T nn in the ballistic regime. In practice, we use the criterion [ ] T 0.9, 1 nn ∈ to ensure that the requirement of ballistic transport is satisfied (see figure 2(b)). Figure 4 shows the numerically obtained results for L 1 nn , including a comparison with the predictions from equations (20) and (21). In panel (a), both modes are displayed, with yellow full circles corresponding to n = 1 and blue diamonds to n = 2, respectively. In the case of antisymmetric waveguides, a direct comparison of the numerical results for the two different modes reveals immediately where the SGS mechanism is at work (see figure 4(a)): since the terms in equations (20) and (21) associated with the AGS mechanism are identical for L 1 11 and L 1 22 , any difference between the two attenuation lengths can be expected to be due to the SGS mechanism. The numerical results reveal that around 2.5 Δ ≈ an extended region opens up in which the two modes decouple and their attenuation lengths are significantly different. To clarify whether this decoupling is, indeed, due to the SGS mechanism, we compare the numerical results with the corresponding analytical predictions in figures 4(b), (c). The agreement we obtain is, again, excellent, allowing us to identify the contributions of the SGS mechanism in detail. First of all, we find that the decoupling of modes is, indeed, due to the SGS mechanism as it is accurately reproduced when the SGS terms are included. Second, the theoretical analysis also predicts that the SGS terms should give rise to small resonant enhancements of the inverse attentuation length at the resonant values k M 2 2 n Δ π = (see arrows in figures 4(b), (c)). Also, these predictions are very well reproduced by the numerical data.
To corroborate the consistency of our above arguments on forward-and backwardscattering contributions, we also investigated the total mode transmissions T t n n N nn 2 d = ∑ ′ ′ , which are now different from the mode-to-mode transmissions T nm due to inter-mode scattering. In the ballistic regime the T n should be determined by backward-scattering alone, since forwardscattering just redistributes the flux which is incoming in one mode over all available rightmoving modes. Since, however, the right-moving modes are summed over in the expression for T n , any influence of forward-scattering drops out in our perturbative treatment. Only when taking into account higher-order forward/backward-scattering events (as in the diffusive or localized regime), the influence of forward-scattering should also be noticeable on the T n . In the ballistic regime, however, we should have T L L 1 . The agreement which we find between the predictions for L 1 n and our numerical results is, in part, remarkably good (see figure 5). A comparison with the expression for L L L 1 1 1 figure 5(d) reveals an excellent agreement between theory and simulation. With the help of figure 5(c), it can be understood that the SGS mechanism contributes by way of two distinct effects: most obviously, we obtain peaks indicating enhanced resonant back-scattering in our system for k , with M integer. Note that these peaks in L 1 2 lead to back-scattering lengths which are about one order of magnitude larger than the (conventional) AGS background. The SGS mechanism can, however, also be identified as a finite contribution to the inverse scattering length at values ( ) ′, exactly where the AGS term in equation (23) vanishes. It is therefore the SGS mechanism which prevents a perfect transparency of the waveguide.
When comparing the first mode data to equation (22), we find that the numerical curve cannot be fully reproduced by the corresponding analytical expression for L 1 1 (see figure 5(b)).
While the AGS contribution L 1 ( in antisymmetric waveguides, 1 nn versus the step width Δ, as obtained numerically for two-mode antisymmetric waveguides, with L 1 11 (yellow ) and L 1 22 (blue ). We find very good agreement with the analytical expressions in equations (20) and (21)  the SGS contribution L 1 ( is a factor 16 smaller (see figure 5(a)). How can this be reconciled with the numerical finding that L 1 1 and L 1 2 are mostly equal? We suspect higher-order terms in scattering to be responsible for these deviations which go beyond the first-order nature of the underlying theory, where the incident wave is assumed to scatter only once before leaving the scattering region. Our aim in the following will be to include such higher-order contributions based on the knowledge of the first-order scattering lengths. Consider here, e.g., the scattering length of the first mode, L 1 , which, as we have assumed so far, is attenuated by back-scattering from the first mode into the first (L ( ) b 11 ) and into the second mode (L ( ) b 12 ), respectively. The next higher-order contribution would be given by forward-scattering into the second mode (governed by L ( ) ). Based on the magnitude of the involved scattering lengths, the forward-scattering 1 2 ↔ occurs much more frequently than a back-scattering event, i.e., the propagating wave undergoes forward-scattering multiple times before it is back-scattered (see figure 6). Consequently, the modes can be assumed to be almost equally distributed between mode 1 or 2 before back-scattering occurs. As a result, the backscattering contribution should also be composed of both modes in equal shares. Since the ). We find a quantitative agreement with the predictions for L 1 2 (panel (d)), while a sizeable discrepancy is observed between numerics and analytical curves for the first mode in (b).
forward-scattering occurs in series and the back-scattering in parallel, this translates into an additional effective second-order term for the inverse scattering length, Equation (24) represents a simple qualitative estimate of second-order contributions to the inverse scattering lengths, and we expect that this expression can be made more quantitative by employing a full-fledged diagrammatic theory. Note that this correction term does not feature an explicit mode dependence since the redistribution of the flux is the same for both propagating modes.
Based on the above, the total inverse scattering lengths can be written as the sum of the following contributions, where the superscripts i ( ) denote the order of the contribution. A comparison of this result with the numerical data is shown in figure 7, yielding much better agreement than without the second-order contributions. In particular, we find (see figure 7(a)) that incorporating the new effective scattering length L 1 ( ) 2,eff resolves the discrepancy we found earlier for the inverse attenuation length of the first mode, L 1 1 . This result also allows us to understand the similarity between the numerical data for L 1 1 and L 1 2 while the first-order SGS contributions are very different for these two modes: the reason is apparently the strong intermode coupling induced by efficient forward-scattering L ( ) f which lets the inverse attenuation length of the first mode L 1 1 inherit the behaviour of the second mode L 1 2 . Correspondingly, the reason for the decoupling between the modes in figure 7 at around 2.5 Δ ≈ can also now be identified: in this parameter window the intermode scattering strength is strongly reduced, allowing the different attenuation lengths to maintain their mode-specific values. Everywhere else (outside this parameter window) the behaviour of L 1 1 is governed by L 1 2 . Figure 6. Illustration of scattering processes of different order. Two processes are shown that attenuate the forward-moving first mode, n = 1. One process is of first order and just consists of a single back-scattering event from the first mode into any of the two backward-moving modes. However, since in the case we consider, forward-scattering 1 2 ↔ is dominant as compared to back-scattering, it is much more likely that the first mode undergoes multiple scattering events in a forward direction before back-scattering occurs. Such terms can thus yield a sizeable contribution, although they formally are of higher order in the number of scattering events which they undergo.

Nonsymmetric profiles
Nonsymmetric waveguides represent the most general case for waveguide symmetries since, in contrast to the previous sections, both boundaries are not restricted by any symmetry requirement, i.e., we have . Turning at first to the partial attenuation lengths, table 1 and equations (14) and (15) We can also immediately write down the corresponding total scattering lengths   ) and including ( ) second-order corrections of equation (24). Note the quantitative agreement that is achieved through the inclusion of higherorder scattering terms. π σ π σ σ π π π σ = + + + + + + + As can be seen from these equations for L 1 n and L 1 nn , intra-and inter-mode as well as the AGS and the SGS scattering lengths now all contribute to the scattering process, in contrast to symmetric or antisymmetric waveguides where the coefficient matrices A nn′ and B nn′ from table 1 feature zeros at symmetry-specific entries. This fact underlines the role of nonsymmetric waveguides as the most general case to study in surface-corrugated systems.
For comparison with the numerical data, we determine the quantities L 1 nn and L 1 n by means of a fit to the transmission in the ballistic regime, T L L 1 nn nn ≈ − and T L L 1 n n ≈ − , in complete analogy to antisymmetric waveguides in the preceding section. Figure 8 shows a comparison between numerics and theory for nonsymmetric two-mode waveguides. Concentrating at first on L 1 nn (first row in figure 8), we note that we find very good agreement between the theoretical and the numerical curves for both the first and the second mode. Similar to the corresponding results in antisymmetric waveguides, a peak at 2.5 Δ ≈ is clearly visible for L 1 22 and, more hidden, also in L 1 11 . The other peaks emerging in the analytical expression for L 1 22 can also be found in our numerical data ( figure 8(b)), albeit slightly more concealed than in the previous section. The positions of these resonances can again be determined from the resonance condition in equation (19).
Turning to the assessment of our results for L 1 n , we can now, with the knowledge from the last section, also take into account the higher-order scattering contributions given by equation (24). Figures 8(c), (d) shows a comparison of the analytical expressions for the attenuation lengths L 1 n with the numerical data. Note that here we also include the first-order predictions as dot-dashed lines. As found before in antisymmetric two-mode waveguides, L 1 2 is already captured very well by equation (29). The arc-structure driven by the AG scattering mechanism again shows a remarkable agreement with the numerical curve; the same is true at resonant points where we find dominating contributions of the SGS mechanism. As before, a more elaborate argument incorporating second-order terms in the scattering is needed for explaining the behaviour of L 1 1 (figure 8(c)). Taking only first-order expressions from equation (28) into account results in theoretical predictions which deviate from our numerical data by about one order of magnitude. Moreover, the period of the oscillations in L 1 1 does not seem to coincide with the analytical predictions. Only after allowing for higher-order terms in L 1 n , i.e., where forward-scattering followed by back-scattering is taken into account by employing equation (24), agreement can be reestablished (figure 8(c) solid line).
Note that the reflection resonances which we observe for nonsymmetric and symmetric waveguides are not as pronounced as in the case of antisymmetric waveguides. This can be understood by the fact that the cross-section of an antisymmetric waveguide remains constant throughout the entire waveguide length such that the wavenumber k x also does not change in the course of propagation. As a consequence, the resonance condition k , with M integer, can be fulfilled very accurately in antisymmetric waveguides, while for waveguides with different symmetries the resonance condition is fulfilled only on average.

Summary
In summary, we have investigated waveguides with a step-like surface disorder supporting two propagating modes. Our study reveals a resonant enhancement of wave reflection in these systems, an effect which has, to the best of our knowledge, not yet been observed earlier, despite the popularity of the employed waveguide model. To manifest this effect, we performed extensive numerical calculations using a waveguide model with symmetric, antisymmetric, and nonsymmetric random profiles, respectively. We compare our numerical findings to a recently proposed surface scattering theory [34,35], which we extend to include higher-order scattering processes as well as to account for the limited resolution with which a scattering wave is sensitive to the surface disorder. We find very good agreement with this new theoretical framework and can thereby associate the origin of the reflection resonances with a higher-order term in the weak disorder expansion of the attenuation lengths. A detailed derivation of this socalled 'square gradient scattering' term is put forward, which, for the systems we consider, results in a fully analytical expression. We show that this previously neglected contribution is , with M integer, we find up to an order-of-magnitude enhancements of the reflection. Not only do our results constitute the first evidence of these resonances in waveguides, but they also provide the first unambiguous signatures of the square-gradient scattering mechanism in waveguides with arbitrary symmetries. The very good agreement which we find between numerical and analytical results provides a solid basis for a general understanding of wave transmission through waveguides with surface roughness. This knowledge may be particularly important in view of experimental possibilities to engineer the transmission characteristics of waveguides through their surface profiles [43,44].
In order to describe an effective smoothing of a step-like waveguide boundary due to a finite resolution capacity of the propagating wave, we consider a profile x ( ) ξ which consists of N 2 1 + steps of width Δ and random heights n α that feature zero mean and unit variance, The smoothing of the steps is modelled by assuming x ( ) Π ρ to be the sum of two Fermi- with the parameter ρ controlling the smearing of the steps, corresponding to the finite resolution of the propagating wave. In the limit of 0 ρ → , i.e., if we assume perfect resolution, the unit box function A1).

A.2. Roughness-height power spectrum W(k)
To calculate , we employ the Wiener-Khinchin theorem [40], Here, we approximate the truncated Fourier transform in equation (A.7) with k ( ) Π ρ , such that it is independent of the summation index n and thus can be pulled in front of the summation. For the parameters employed in the present paper, this step is very well justified and only leads to a vanishingly small error.