2. Determining Intensity Extrema and Their Loci in Near-Field Single-Slit Diffraction Using the Half-Wave Zone Model
The half-period zone model is a conceptual framework used to explain light diffraction without relying on the direct calculation of Kirchhoff’s diffraction integrals. Instead, it is based on dividing the aperture into specific zones.While applying Fresnel's circular half-period zones to a circular aperture has proven successful, describing diffraction through a slit requires a slightly different approach. Specifically, in the slit model, the areas of the zones do not remain constant, meaning the number of secondary sources within each zone varies. My method for constructing these zones differs in the design of the central zone and can be found in the applet.
In this applet, we identify specific observation points where the boundaries of the outer zones coincide with the physical edges of the slit. We propose that these points are "close" to the critical points of the intensity function J = J(x, y), as calculated via the Kirchhoff integral for a slit. Comparisons with numerous Heatmaps generated using this integral for different slits widths confirm the validity of these assumptions (1). Furthermore, numerical calculations of the extrema align with the model's predictions with a high degree of accuracy (2, 3).
The primary objective is to use geometric transformations to determine the coordinates (x, y) of the intensity extrema J(x, y) within a diffraction pattern. These coordinates are expressed as functions of the integers m1 and m2, which define how the half-wave zones occupy the slit. For integer values of these filling numbers, the physical edges of the slit coincide with the boundaries of the full zones; in fact, the outermost boundaries of the zones should ideally coincide with the slit edges.
It has been established that these coordinates can be mapped using the functions: fx and fy: A(x, y) = (fx(m1, m2), fy(m1, m2)) (see Fig. 1). The derivation of these functions is based on a geometric analysis of the aspect ratios and coordinates within triangle △ABC (refer to the figures below):
● The vertex A(x, y) is the observation point. The base a is equal to the slit width b. Sides b and c are defined as: b: r1= y + Δ0 + m1 * λ/2 and c: r2= y + Δ0 + m2 * λ/2, respectively.
● m1 and m2: These represent the number of full zones fitting into the right and left parts of the slit, respectively, relative to the middle of the central zone (Point A').
● We define Δ0 as the slit constant. To determine its value, we consider an observation point located on the axis at the position of the final intensity J=J(x) maximum, as calculated via the Fraunhofer diffraction integral for the slit. This specific point is designated as the first focus (F1). In this configuration, the entire width of the slit is treated as the central zone. Consequently, Δ0 is defined as the optical path difference at this focal point between the central (axial) ray and the ray from the central zone boundary at the slit edge.
The results demonstrate that the extrema of the function J(x, y) in this model belong to three distinct families of curves. The actual extremum points are located at the intersections of these curves -see Fig. 6.
The resulting point system reveals three distinct directional groups, each following a characteristic curve.
1. Evidently, the points in the field for which the path difference Δ=r1-r2 remains constant lie on a family of hyperbolas. Given that Δ = (m1 - m2)*λ/2, where (m1 - m2) is a constant integer, these curves represent lines of constructive or destructive interference. The foci of these hyperbolas are the "sources" located at the ends of the slit, and the quantity (m1-m2) determines the specific branch of the family.
![[size=85][b][b]Fig. 1.[/b][u] Hyperbolic [/u][/b][u][b]"[color=#ff0000]Antin[/color][/b][/u][u][b][color=#ff0000]odal[/color][/b][/u][u][b]" [/b][/u][u][b]and [/b][/u][u][b]"N[/b][/u][color=#333333][u][b]odal[/b][/u][/color][u][b]" [/b][/u][u][b]lines [/b][/u][b]of [color=#ff0000]constructive[/color] and destructive interference of the [/b][color=#004dbb][i][b]two point source interference pattern [/b][/i][/color][i][b][color=#004dbb]from the edges [/color][color=#ff0000]B[/color][color=#004dbb] and [/color][color=#ff0000]C[/color][color=#004dbb] of a slit [/color][/b][/i][color=#333333][b]in single-slit diffraction pattern [/b][b]itself.
[/b] Formulas for determining the coordinates [b]x[sub]k[/sub][/b] and [b]y[sub]k[/sub][/b] of the "extrema" of the diffraction field can be found at the bottom of the figure.[/color][/size]](https://www.geogebra.org/resource/khxfqjp4/7eNf6JEPwgr5lIWL/material-khxfqjp4.png)
2. For the 'extreme' points of the model lying on each parabolic line, one integer index remains constant, while the other changes by one.

3. Implicit equation of the focal curve.
For any pair of extrema located symmetrically on the focal curve (e.g., points 0 and 4, or 1 and 3 from figure a), the sum of the parameters m1k, m2k is invariably equal to u - 1. Leveraging this symmetry, we can establish a system of equations (equations from figure b) by introducing a parameter k, which, in this specific case, varies from 4 to 0. By eliminating k between these equations, we derive the desired implicit equation.
Here the focal curve is defined as an implicit function. By adjusting the curve order (slider u), the curve is made to pass through all extrema corresponding to the focal points on the slit axis. The $15 and $18 lines (equations from figure b) describe the expressions for the focal curve in the cases of vertical and horizontal slit axes, respectively.

![[b]Fig. 3a and Fig. 3b [/b]](https://www.geogebra.org/resource/ssbrxyt2/8iTpK8dmTVeiqkQW/material-ssbrxyt2.png)
4. Vector Zone Analysis and Characteristic Features

5. Formation of zones for the "critical points" in the proposed model of light diffraction behind a slit
●Central Zone: A central zone of action within the slit is present for all "extreme points."
● Max: A characteristic feature of zone formation is that all zones of opposing actions (excluding the central zone) are compensated: at the ends of the slit there are red zones.
●Min: A characteristic feature of zone formation is that all zones of opposing actions (excluding the central zone) on either side of the central zone are uncompensated: at the ends of the slit there are blue zones.
●Zone formation for "saddle points":
-There are two distinct types of saddle points based on their location and formation -See the note for Fig. 2.
-Common Features: For both types, there is a lack of compensation (excluding the central zone): at the ends of the slit there are multi-colored(red-blue) zones.
![[b]Fig. 5.[/b]](https://www.geogebra.org/resource/pzqckgbw/IrvUK0LEIvVc4UVZ/material-pzqckgbw.png)
6. A Heatmap of the diffraction field, calculated using the Kirchhoff diffraction integral. Superimposed on the field are the "extreme points" derived from the proposed model of hyperbolic-parabolic light diffraction curves behind a slit

An interactive exploration of the described model is available in the applet: "Visualizing Near-Field Single-Slit Diffraction: A Tutorial Applet Based on the Half-Wave Zone Model"
