Difference between revisions of "Huygens Principle for a Planar Source"
| Line 82: | Line 82: | ||
We invoke the principle that the integral of a derivative is the function evaluated at the end points to give us the equation<br><br><math>\mathcal{J}(r')=2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}\left[aJ_1(ka\sin{\theta}')-0J_1(0k\sin{\theta}')\right]=2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}aJ_1(ka\sin{\theta}')</math><br><br> and <math>A(r')=\frac{z'\dot{A_0}}{2\pi c}e^{-ikct'} 2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}aJ_1(ka\sin{\theta}')=\frac{z'\dot{A_0}a}{c}\frac{e^{ikr'-ikct'}}{k\sin{\theta}'r'^2}J_1(ka\sin{\theta}')</math><br><br> | We invoke the principle that the integral of a derivative is the function evaluated at the end points to give us the equation<br><br><math>\mathcal{J}(r')=2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}\left[aJ_1(ka\sin{\theta}')-0J_1(0k\sin{\theta}')\right]=2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}aJ_1(ka\sin{\theta}')</math><br><br> and <math>A(r')=\frac{z'\dot{A_0}}{2\pi c}e^{-ikct'} 2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}aJ_1(ka\sin{\theta}')=\frac{z'\dot{A_0}a}{c}\frac{e^{ikr'-ikct'}}{k\sin{\theta}'r'^2}J_1(ka\sin{\theta}')</math><br><br> | ||
To find the angle to the first diffraction minimum, we must find the first zero of this amplitude function. This will occur when <math>J_1(ka\sin{\theta}')=0 \quad</math> | To find the angle to the first diffraction minimum, we must find the first zero of this amplitude function. This will occur when <math>J_1(ka\sin{\theta}')=0 \quad</math> | ||
| − | [[Image: | + | [[Image:Bessels_J0.svg|thumb|300px|right|Plot of Bessel function of the first kind, J<sub>α</sub>(x), for integer orders α=0,1,2.]] |
Revision as of 21:52, 8 July 2009
We start off with Maxwell's Equation in the Lorentz gauge:
where we use the metric signature (+,+,+,-) and
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A^\mu = (\mathbf{A},\frac{\Phi} {c})}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \square^2=\part_\mu \part^\mu = \nabla^2 - \frac{1}{c^2} \frac{\part^2}{\part t^2}}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle j^\mu = (\mathbf{j},c\rho), \part_\mu= (\mathbf{\nabla}, \frac{1}{c} \frac{\part}{\part t})}
The gauge condition for the Lorentz gauge is
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \part_\mu A^\mu = 0 \rArr \mathbf{\nabla} \cdot \mathbf{A}-\frac{1}{c^2} \frac{\part\Phi}{\part t}=0}
Introduce the Green's function at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r=(\mathbf{r},t)} from some impulse source at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r'=(\mathbf{r}',t')}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \square^2_rG(r,r')=\delta^4(r-r')}
and its Fourier transform
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \tilde{G} (q) = \frac{1}{(2\pi)^2} \int d^4r\, e^{-iq\cdot r} G(r,0)}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r,0)=\frac{1}{(2\pi)^2} \int d^4q\, e^{iq\cdot r} \tilde{G}(q)}
Translational symmetry implies
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r-r',0)=G(r,r')\quad}
so that
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r,r')=\frac{1}{(2\pi)^2}\int d^4q\, e^{iq\cdot (r-r')} \tilde{G} (q)}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \square^2_rG(r,r')=\frac{1}{(2\pi)^2}\int d^4q\,(-q^2)e^{iq\cdot (r-r')}\tilde{G}(q)}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \square^2_rG(r,r')=\frac{1}{(2\pi)^2}\int d^4q\, e^{iq\cdot (r-r')}(-k^2+\frac{\omega^2}{c^2})}
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle q=(\mathbf{k},\frac{\omega}{c})} . But
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \square^2_rG(r,r')=\delta^4(r-r')=\frac{1}{(2\pi)^4}\int d^4q\, e^{iq\cdot (r-r')}}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \tilde{G}(q)=\frac{(2\pi)^2}{(2\pi)^4}\frac{1}{-q^2}= \frac{-1}{(2\pi)^2q^2}}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r,r')=\frac{-1}{(2\pi)^4} \int d^4q\, e^{iq\cdot (r-r')} \frac{1}{(k+\frac{\omega}{c})(k-\frac{\omega}{c})}}
Chose the "retarded" solution, such that the function is zero unless t>t'.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r,r')=\frac{1}{(2\pi)^4}\int d^3k\, e^{i\mathbf{k}\cdot (r-r')}\int d(\frac{\omega}{c}) \frac{e^{-i\omega(t-t')}}{(\frac{\omega}{c}-k)(\frac{\omega}{c}+k)}\Theta(t-t')}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle =\frac{-2\pi}{(2\pi)^4}\int_0 \frac{k^2dk}{k} \sin\left({ck(t-t')}\right) 2\pi\int_{-1}^1 dze^{ik|\mathbf{r}-\mathbf{r'}|z}\Theta}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle =\frac{-1}{(2\pi)^2}\left(\frac{1}{|\mathbf{r}-\mathbf{r'|}}\right)2\int_0 dk \sin(ck(t-t')) \sin(k|\mathbf{r}-\mathbf{r'}|)\Theta}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle =\frac{1}{(2\pi)^2}\frac{2}{|\mathbf{r}-\mathbf{r}'|}\frac{2\pi}{4} \left[\delta(|\mathbf{r}-\mathbf{r}'|+c(t-t'))-\delta(|\mathbf{r}-\mathbf{r}'|-c(t-t'))\right]\Theta}
But the term Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta(|\mathbf{r}-\mathbf{r}'|+c(t-t'))\rightarrow 0 \quad\forall\quad t>t'}
so that
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G(r,r')=\frac{-1}{4\pi}\quad \frac{\delta(|\mathbf{r}-\mathbf{r}'|-c(t-t'))}{|\mathbf{r}-\mathbf{r}'|}}
Now to get the in the half-space with z>0 with the boundary condition Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G_1\quad }
atFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_3=z=0 \quad}
we take the difference:
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G_1(r,r')=\frac{-1}{4\pi}\left(\frac{\delta(|\mathbf{r}-\mathbf{r}'|-c(t-t'))}{|\mathbf{r}-\mathbf{r}'|}-\frac{\delta(|\mathbf{r}-\mathbf{r}'+2z'\hat{e_3}|-c(t-t'))}{|\mathbf{r}-\mathbf{r}'+2z'\hat{e_3}|}\right)}
Now use Green's theorem, with the generating function
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F^\mu=A(r)\part_\mu G_1(r,r')-G_1(r,r')\part_\mu A(r)}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \part_\mu \part^\mu G_1(r,r')=\delta^4(r-r')}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \part_\mu \part^\mu A(r)= \mu j(r)}
, let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle j(r)=0 \quad}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int \part_\mu F_\mu d^4r=A(r')}
Now invoke the divergence theorem on the half space :
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A(r')=-\int d^2r\int cdt\left[A(r)\frac{\part}{\part z}G_1(r,r')-G_1(r,r')\frac{\part}{\part z}A(r)\right]}
, where the last term is zero by the condition ofFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G_1(z=0,r')=0 \quad}
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A(r')=-c\int dt\int d^2rA(r)\frac{\part}{\part z}G_1(r,r')}
To do the t integral, I need to bring out the z derivative. To do this, I first turn it into a z' derivative, using the relation
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G_1(r,r')=\frac{-1}{4\pi}\left(\frac{\delta(|\mathbf{r}-\mathbf{r}'|-c(t-t'))}{|\mathbf{r}-\mathbf{r}'|}-\frac{\delta(|\mathbf{r}-\mathbf{r}''|-c(t-t'))}{|\mathbf{r}-\mathbf{r}''|}\right)}
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}''=\mathbf{r}'-2z'\hat{e_3}}
∴ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A(r')=\frac{-1}{4\pi}\frac{\part}{\part z'}\int_{z=0} d^2r\left(2\frac{A(\mathbf{r},t'-\frac{|\mathbf{r}-\mathbf{r}'|}{c})}{|\mathbf{r}-\mathbf{r}'|}\right)}
At Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z=0 \quad }
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\mathbf{r}-\mathbf{r}'|=\sqrt{r^2+z'^2}=S, dS=\frac{rdr}{\sqrt{r^2+z'^2}}}
If is independent of position, as in a plane wave propagating along the z axis, then:
This gives us uniform translation of waves at velocity c. More generally:
In our case, we consider only those waves which drop off as , so
In cylindrical coordinates, . Without loss of generality, we consider a harmonic solution with a particular frequency ω = kc.
Special Case
Picture an opaque screen with a circular aperture of radius a.
Let
Then
But
so that and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{1}{|\mathbf{r}-\mathbf{r}'|^2} \approx \frac{1}{r'^2}\left(1+\frac{2r\sin\theta'\cos\phi}{r'}\right)}
In this particular case, we are dealing with far-field effects only, so Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{2r\sin\theta'\cos\phi}{r'}\rightarrow 0 }
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{1}{|\mathbf{r}-\mathbf{r}'|^2} \approx \frac{1}{r'^2}}
So, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{J}(r')=\int_0^a rdr\int_0^{2\pi} d\phi\, \frac{e^{ik|\mathbf{r}-\mathbf{r}'|}}{|\mathbf{r}-\mathbf{r}'|}=\frac{e^{ikr'}}{r'^2}\int_0^a rdr\int_0^{2\pi} d\phi\, e^{-ikr\sin{\theta}'\cos{\phi}}}
The integral Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int_0^{2\pi} d\phi\, e^{-ikr\sin{\theta}'\cos{\phi}}}
is the integral representation of the zero order Bessel function of the first kind with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle kr\sin{\theta}' \quad}
as the argument. This gives us the equation:
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{J}(r')=\frac{e^{ikr'}}{r'^2}\int_0^a rdr 2\pi J_0(kr \sin{\theta}') }
To simplify the math, we make use of the fact that we can represent this Bessel functions as the derivative of a Bessel function of a different order. In general, the formula to compute this derivative is
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z^{v-k}J_{v-k}(z)=\left(\frac{1}{z}\frac{\part}{\part z}\right)^kz^vJ_v(z)}
In this case, we take Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v=k=1 \quad}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z=kr\sin{\theta}' \quad}
. So
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle J_0(kr\sin{\theta}')=\left(\frac{1}{kr\sin{\theta}'}\frac{\part}{\part (kr\sin{\theta}')}\right)(kr\sin{\theta}')J_1(kr\sin{\theta}')=\frac{d}{k\sin{\theta}'dr} J_1(kr\sin{\theta}')}
This gives us the equationFailed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{J}(r')=2\pi\frac{e^{ikr'}}{r'^2}\int_0^a rdr \frac{d}{k\sin{\theta}'dr} J_1(kr\sin{\theta}')}
We invoke the principle that the integral of a derivative is the function evaluated at the end points to give us the equation
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{J}(r')=2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}\left[aJ_1(ka\sin{\theta}')-0J_1(0k\sin{\theta}')\right]=2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}aJ_1(ka\sin{\theta}')}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A(r')=\frac{z'\dot{A_0}}{2\pi c}e^{-ikct'} 2\pi\frac{e^{ikr'}}{k\sin{\theta}'r'^2}aJ_1(ka\sin{\theta}')=\frac{z'\dot{A_0}a}{c}\frac{e^{ikr'-ikct'}}{k\sin{\theta}'r'^2}J_1(ka\sin{\theta}')}
To find the angle to the first diffraction minimum, we must find the first zero of this amplitude function. This will occur when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle J_1(ka\sin{\theta}')=0 \quad}