Law of refraction from Scalar Diffraction Theory I am trying to figure out if Snell's Law for refraction can be derived from Scalar Diffraction Theory.
The setup is this: light  (plane wave, with wave vector $\vec k_i = (k_x, k_y, k_z)$ ) falls on a flat interface which is taken to be x-y plane; the incident side has refractive index $n_i$ . I want to figure out the refracted ray wave vector $\vec k$, on the otherside of the interface, which has refractive index $n$. Of course, I assume that the magnitutes of wave vectors are proportional to the refractive indices of respective media; i.e. $ |k|/|k_i| = n/n_i$.
I start with the Fresnel Diffraction Integral, with plane wave light $$U(\xi,\eta,z=0) = \exp( i \left[k_x \xi + k_y \eta + 0 \right] )$$ and  got to the point of showing that $$U(x,y,z) = \exp( i \left[ k_x x + k_y y + k_z^{\prime} z \right] ),$$ where $$ k_z^{\prime} = k - \frac{k_x^2 + k_y^2}{2k}. $$
I was happy to see that this field has the same $k_x$ and $k_y$ values, but no matter how I play with $k_z^{\prime}$, I am unable to show that it is consistent with Snell's Law.
Any pointers? 
 A: In scalar wave theory you're simply going to impose the continuity of the scalar field across the boundary. 
The reasoning is the same whether or not you consider reflexions, whose magnitude and form can be gotten from scalar theory by imposing a further boundary condition of continuous normal derivatives across the interface.
With or without reflexions, the scalar variation along the interface on the incidence side is $\exp(i\,k\,n_i\,\sin\theta_i\,x)$, where $n_i$ is the refractive index of the medium on the incoming side of the interface, $k$ the freespace wavenumber and $x$ is measured along the interface plane in the plane of incidence (plane containing wavevector and interface normal). The $\sin\theta_i$ factor arises simply because the wavevector is skewed relative to the interface normal: it is the component of the wavevector in the interface plane.
By the same reasoning, the scalar variation along the interface on the output (transmission) side is simply $\exp(i\,k\,n_t\,\sin\theta_t\,x)$, where $n_t$ is the refractive index on that side.
So, for continuity of the scalar field, we must have:
$$\exp(i\,k\,n_i\,\sin\theta_i\,x) = \exp(i\,k\,n_t\,\sin\theta_t\,x),\;\forall x\in\mathbb{R}$$
which can only be true if
$$n_i\,\sin\theta_i = n_t\,\sin\theta_t$$
which is of course Snell's Law.
The full, vector field picture is to be found in the derivation of the Fresnel equations as given, for example, in section 1.5 of the sixth edition of:
Born and Wolf, "Principles of Optics"
however the simpler derivation above is more general, since it works for scalar waves such as sound, too.
A: Well I'm not sure where in your math the refractive indices of the media are hidden, or even the propagation velocities, so it's not clear to me that Snell's law should fall out of your diffraction theory.
Snell's law can commonly be obtained from Fermat's principle, by simply solving for the minimum path length across a media boundary.
A: Answering my own question: George's comment and Rod's answer, while not the exact answer, directly lead to me figure out the missing was the "math" steps! (which is below)
The Fresnel Diffraction Formula basically (somewhere in its derivation) Taylor expands square root i.e. $r=\sqrt(x^2 + y^2 + z^2) = z (1 + \frac{x^2 + y^2}{2z^2})$, with an approximation $x^2 + y^2 << z^2$. So what I had for normal component of the wave-vector after refraction was $ k - \frac{k_x^2 + k_y^2}{2k}$ basically the same approximation for $\sqrt{k^2 - (k_x^2 + k_y^2)}$, showing that the the wave-vector after refraction is $(k_x, k_y, \sqrt{k^2 - (k_x^2+k_y^2})$, where the wave vector before refraction was $(k_x, k_y, k_z)$.  Given that $\sqrt{k_x^2+k_y^2+k_z^2}/k = n_i/n$, the wave vectors satisfy Snell Law for refraction!
