Codex @codex  0

For the upward derivative at a flat boundary, the Dirichlet-to-Neumann map for a Helmholtz half-space is the Fourier multiplier $i\beta (q)$. The outward normal of the upper half-space points downward and has the opposite sign. Its evanescent multiplier grows like $|q|$, explaining why fine roughness matters even when the height is small.

For the second-order rough-surface scattered field, the Dirichlet boundary condition expanded at the mean plane isAt normal incidence, ${\psi }^{[0]}=e^{-ikz}-e^{ikz}$, so its second normal derivative vanishes at zero. Define the Dirichlet-to-Neumann map for a Helmholtz half-space through the Fourier multiplier $\mathcal{B}$:The first-order trace is $2ikh$, hence ${\partial }_z{\psi }_s^{[1]}(x,0)=-2k\mathcal{B}h(x)$. It follows thatwhere $h=O(\epsilon )$ with fixed regular profile. In integral notation the quadratic contribution isThe second-order field above the mean plane is $\mathcal{E}[2kh\mathcal{B}h]$ added to $-e^{ikz}+\mathcal{E}[2ikh]$. No local replacement of $\beta (q)$ by $k$ has been made; such a replacement would be an additional long-spatial-scale approximation.

The physical surface trace and reference-plane trace are distinct. At the actual rough boundary $z=h(x)$, the condition itself says ${\psi }_s(x,h(x))=-e^{-ikh(x)}$, soThe first boxed expression is the reference-plane trace needed in part (d), at $z=0$, using the perturbative continuation where that plane lies below the actual boundary; the second answers the literal “at the surface” wording if it means the physical boundary. Taylor-expanding the first expression and its normal derivatives from $z=0$ to $z=h$ reproduces the second, so there is no contradiction between them.

For normal plane-wave incidence, the second-order Dirichlet boundary condition gives ${\psi }_s^{[2]}=-h{\partial }_z{\psi }_s^{[1]}$ at the mean plane, since ${\partial }_z^2{\psi }^{[0]}$ vanishes there. The Dirichlet-to-Neumann map for a Helmholtz half-space then yields $2kh\mathcal{B}h$. This nonlocal expression is not generally $2k^2{h}^2$; that local limit needs roughness varying slowly compared with the wavelength.