User Tools
app:diffraction_grating
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
|
app:diffraction_grating [2018/01/30 12:58] pklapetek |
app:diffraction_grating [2018/01/30 13:05] pklapetek |
||
|---|---|---|---|
| Line 21: | Line 21: | ||
| Near-to-far field transformation (NFFF) integrated from plane behind the aperture is used to calculate far field radiation pattern. Preserving both TSF source and NFFF integration plane close to metal plane helps to minimize influence of finite plane wave and metal plane size. | Near-to-far field transformation (NFFF) integrated from plane behind the aperture is used to calculate far field radiation pattern. Preserving both TSF source and NFFF integration plane close to metal plane helps to minimize influence of finite plane wave and metal plane size. | ||
| - | Far field radiation pattern is evaluated on a sphere (or its part) with prescribed radius, here approximately 3 cm. Spherical coordinates are used for this, so the output of calculation spans some values of θ (in range of 0-360) and φ (in range of 0-180), so the results can contain different parts of a sphere as shown for different settings on the left image. As a result we get a 2D diagram of scattered light intensity that can be visualised or processed in Gwyddion. | + | Far field radiation pattern is evaluated on a sphere (or its part) with prescribed radius, here approximately 3 cm. Spherical coordinates are used for this, so the output of calculation spans some values of θ (in range of 0-360) and φ (in range of 0-180), so the results can contain different parts of a sphere as shown for different settings on the left image. As a result we get a 2D diagram of scattered light intensity that can be visualised or processed in Gwyddion. |
| There are two possible ways how to create source. We can use a monochromatic plane wave (e.g. 633 nm) and evaluate the scattered light at this wavelength directly. Alternatively, we can use a broadband pulse, here a gaussian function, determine its spectral dependency by Fast Fourier Transform (FFT) and determine spectral dependency of the computation output (again by FFT). This gives us the spectral dependency of scattered light from a single computation. | There are two possible ways how to create source. We can use a monochromatic plane wave (e.g. 633 nm) and evaluate the scattered light at this wavelength directly. Alternatively, we can use a broadband pulse, here a gaussian function, determine its spectral dependency by Fast Fourier Transform (FFT) and determine spectral dependency of the computation output (again by FFT). This gives us the spectral dependency of scattered light from a single computation. | ||
| As a result we obtain a far field radiation pattern as shown below on the right image (analytical result is in the left image). Both images are in logarithmic scale for intensity false color plot. Only a single lobe is calculated (i.e. a quarter of the whole pattern that is four times repeated in 0-360 degrees due to the fact that aperture is square). We can clearly see the diffraction pattern and we can compare it with analytical result. | As a result we obtain a far field radiation pattern as shown below on the right image (analytical result is in the left image). Both images are in logarithmic scale for intensity false color plot. Only a single lobe is calculated (i.e. a quarter of the whole pattern that is four times repeated in 0-360 degrees due to the fact that aperture is square). We can clearly see the diffraction pattern and we can compare it with analytical result. | ||
| + | {{ :app:aperture_results.png?400 |}} | ||
| Comparison can be even better done on single profile across the direction of diffraction pattern as shown below. Note the effect of limited angular resolution of the calculation which can be one of reasons of small differences (besides limited NFFF integration area and other numerical errors). | Comparison can be even better done on single profile across the direction of diffraction pattern as shown below. Note the effect of limited angular resolution of the calculation which can be one of reasons of small differences (besides limited NFFF integration area and other numerical errors). | ||
| + | {{ :app:aperture_graphs.png?600 |}} | ||
| ==== Transmission through a grating ==== | ==== Transmission through a grating ==== | ||
| - | + | {{:app:a_grating_mrizka.png?200 |}} | |
| Grating is formed by equally spaced set of apertures. An image of ideal diffraction grating aperture locations is shown in the left figure (ideally the number of apertures in every direction would be infinite). Grating is again made from non-transparent material and as light goes through its slits (shown by black color) it is diffracted. Due to constructive and destructive interference the transmitted light has strong dependency of intensity on propagation angle, so if we place a screen behind the grating we will see equally spaced set of points - diffraction maxima - modulated by diffraction of single aperture forming the grating. | Grating is formed by equally spaced set of apertures. An image of ideal diffraction grating aperture locations is shown in the left figure (ideally the number of apertures in every direction would be infinite). Grating is again made from non-transparent material and as light goes through its slits (shown by black color) it is diffracted. Due to constructive and destructive interference the transmitted light has strong dependency of intensity on propagation angle, so if we place a screen behind the grating we will see equally spaced set of points - diffraction maxima - modulated by diffraction of single aperture forming the grating. | ||
| Line 40: | Line 41: | ||
| To evaluate far field scattering pattern, a periodic NFFF computational boundary is used, which means that local fields on integration boundary are copied as many times as requested, simulating much larger periodic structure. | To evaluate far field scattering pattern, a periodic NFFF computational boundary is used, which means that local fields on integration boundary are copied as many times as requested, simulating much larger periodic structure. | ||
| - | The whole pattern as simulated for grating spacing to wavelength ratio of approx 10:1 and finite size of 13x13 holes is shown below, both in normal and logarithmic scale (note that in contrast to aperture simulation results now the θ axis is the horizontal one). | + | The whole pattern as simulated for grating spacing to wavelength ratio of approx 10:1 and finite size of 13x13 holes is shown below. |
| + | {{ :app:simulated_grating_1313.png?400 |}} | ||
| Analytically this can be written for diffraction angles θ1 (in the x direction) and θ2 (in the y direction) as | Analytically this can be written for diffraction angles θ1 (in the x direction) and θ2 (in the y direction) as | ||
| + | |||
| + | {{ :app:eq_grating.png?600 |}} | ||
| where a1 and a2 are grating hole spacing in x,y direction and N1 and N2 are number of holes in these directions; the rest of symbols is same as in aperture equation. | where a1 and a2 are grating hole spacing in x,y direction and N1 and N2 are number of holes in these directions; the rest of symbols is same as in aperture equation. | ||
| Even if it would be in principle possible to calculate the same images also using FDTD, the number of far field points with the same resolution would be around million which is already significantly slow in present version of GSvit. Therefore we had compared result only on a single profile in x direction again, similarily to how the aperture intensity graphs were obtained. Results are show below (x axis is in degrees) | Even if it would be in principle possible to calculate the same images also using FDTD, the number of far field points with the same resolution would be around million which is already significantly slow in present version of GSvit. Therefore we had compared result only on a single profile in x direction again, similarily to how the aperture intensity graphs were obtained. Results are show below (x axis is in degrees) | ||
| + | |||
| + | {{ :app:a_grating_results.png?600 |}} | ||
app/diffraction_grating.txt · Last modified: 2018/08/29 12:25 by pklapetek
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Noncommercial 4.0 International
