Interference and Diffraction Laboratory Study Guide

Experimental Objectives and Overview

• Primary Objectives:

  • Examine wave phenomena including single-slit diffraction and double-slit interference.

  • Measure and analyze the locations of intensity maxima and minima produced by monochromatic light passing through narrow slits.

  • Calculate slit width ($a$) and slit separation ($d$) along with their corresponding uncertainties, comparing experimental findings to manufacturer-specified values.

  • Derive the full theoretical intensity equations for single-slit diffraction and double-slit interference from first principles.

  • Model the complete double-slit diffraction-interference intensity distribution using computational tools and compare theoretical predictions directly against measured experimental data.

  • Master formal laboratory report structure and numerical analysis protocols.

• Key Apparatus & Parameters:

  • Monochromatic Diode Laser source with wavelength \n\lambda = 650\,\text{nm}\n

  • Diffraction slit plates (single slits and double slits of varying width $a$ and separation $d$)

  • Viewing screen and paper mounting surface

  • Optical rail setup equipped with a combined light intensity sensor and linear travel sensor

  • Data collection interface utilizing Logger Pro software

Part I: Visualizing Diffraction Patterns and Measuring Extremes

Single Slit Diffraction Procedure & Data Analysis

• Data Collection Steps:

  • Attach a blank sheet of paper securely to the viewing screen.

  • Direct the diode laser beam through the single slit onto the screen.

  • Mark the exact locations of the dark fringes (diffraction minima) on the attached paper.

  • Mark the central maximum location using an arrow indicator.

  • Measure and record the perpendicular distance ($L$) from the slit to the screen using a precision ruler.

  • Remove the paper from the screen to commence quantitative data processing.

• Data Analysis Steps:

  • Measure and record the distance between symmetric pairs of minima across the central maximum in a spreadsheet.

  • Divide the measured separation between each pair by $2$ to compute the distance ($y_{nn}$) from the central maximum to the $n^{\text{th}}$ minimum ($y_{1n}, y_{2n}, y_{3n}, \dots$).

• Calculate the diffraction angle ($\theta_n$) corresponding to each minimum using standard trigonometry:     \n\tan(\theta_n) = \frac{y_{nn}}{L} \implies \theta_n = \arctan\left(\frac{y_{nn}}{L}\right)\n

• Apply the single-slit diffraction minima condition to calculate the slit width ($a$):     \n\sin(\theta_n) = \frac{n\lambda}{a} \implies a = \frac{n\lambda}{\sin(\theta_n)}\n     where $n = 1, 2, 3, \dots$ represents the integer order of the minimum.

• Estimate the measurement uncertainty for the calculated slit width ($a$).

Two Slit Interference Procedure & Data Analysis

• Data Collection Steps:

  • Attach a new sheet of paper to the viewing screen.

  • Direct the diode laser beam through the double slit set.

  • Mark the center positions of the bright interference spots (interference maxima).

  • Mark the exact position of the central bright maximum ($m = 0$) with an arrow.

  • Place distinct markings ("x") at the locations of the primary single-slit diffraction envelope minima modulating the pattern.

  • Measure and record the slit-to-screen distance ($L$).

• Data Analysis Steps:

  • Measure and record distances between symmetric pairs of single-slit envelope minima to calculate the single-slit width ($a$) using single-slit diffraction equations.

  • Measure and record the distances between symmetric pairs of fine double-slit interference maxima moving outward from the central maximum.

  • Divide each pair distance by $2$ to obtain the position $y_{ms}$ from the central maximum to the $m^{\text{th}}$ maximum ($y_{1s}, y_{2s}, y_{3s}, \dots$).

  • Calculate the angular displacement ($\theta_m$) for each interference maximum:     \n\theta_m = \arctan\left(\frac{y_{ms}}{L}\right)\n

  • Determine the slit separation ($d$) using the double-slit interference maximum condition:     \n\sin(\theta_m) = \frac{m\lambda}{d} \implies d = \frac{m\lambda}{\sin(\theta_m)}\n     where $m = 0, 1, 2, 3, \dots$ represents the integer order of the interference maximum.

  • Quantify the experimental uncertainty associated with the calculated separation ($d$).

Part II: Measuring Intensity Distribution and Mathematical Derivations

Fundamental Relationship Between Intensity and Electric Field

• Wave intensity ($I$) is directly proportional to the time average of the magnitude of the Poynting vector ($\vec{S}$), which in turn scales with the square of the electric field magnitude ($E$):   \nI = \langle S \rangle \propto \langle E^2 \rangle\n

• Time-average of a time-dependent function $f(t)$ over a period $T$ is defined as:   \n\langle f(t) \rangle = \frac{1}{T} \int_0^T f(t)\,dt\n   where period $T$ for an angular frequency $\omega$ is given by:   \nT = \frac{2\pi}{\omega}\n

• Time-average of a squared harmonic function over a full period satisfies:   \n\sin^2(\omega t) = \frac{1 - \cos(2\omega t)}{2}\n   \n\left\langle \sin^2(\omega t + \varphi_0) \right\rangle = \frac{1}{T} \int_0^T \sin^2(\omega t + \varphi_0)\,dt = \frac{1}{T} \left[ \frac{T}{2} - \frac{\sin(2\omega T + 2\varphi_0) - \sin(2\varphi_0)}{4\omega} \right] = \frac{1}{2}\n   since $\cos(2\omega t)$ oscillates symmetrically about zero over any full period $T$.

Complete Derivation of Single Slit Diffraction Intensity Distribution

• Huygens' Principle Sub-zone Division:

  • Consider a single slit of width $a$ illuminated by uniform, coherent plane waves.

  • Divide the slit width $a$ into $N$ equal continuous sub-zones, each having width:     \n\Delta y = \frac{a}{N}\n

  • The physical path difference ($\delta$) between rays emitted from adjacent sub-zones traveling toward point $P$ at angle $\theta$ is:     \n\delta = \Delta y \sin(\theta)\n

  • Corresponding phase difference ($\Delta \beta$) between adjacent sub-zones:     \n\frac{\Delta \beta}{2\pi} = \frac{\delta}{\lambda} = \frac{\Delta y \sin(\theta)}{\lambda} \implies \Delta \beta = \frac{2\pi}{\lambda} \Delta y \sin(\theta)\n

  • Total phase difference ($\beta$) across the entire slit width $a$ from sub-zone $1$ to sub-zone $N$:     \n\beta = N \Delta \beta = \frac{2\pi}{\lambda} N \Delta y \sin(\theta) = \frac{2\pi}{\lambda} a \sin(\theta)\n

• Superposition of Electric Fields from $N$ Zones:

  • Let the electric field contribution from the first zone at point $P$ be:     \nE_1 = E_{10} \sin(\omega t)\n

  • Subsequent zones contribute fields with progressive phase increments of $\Delta \beta$:     \nE_2 = E_{10} \sin(\omega t + \Delta \beta)\n     \nE_3 = E_{10} \sin(\omega t + 2\Delta \beta)\n     \nE_N = E_{10} \sin\left(\omega t + (N-1)\Delta \beta\right)\n

  • Total field $E$ at screen point $P$ is the vector sum:     \nE = \sum_{j=1}^N E_j = E_{10} \left[ \sin(\omega t) + \sin(\omega t + \Delta \beta) + \dots + \sin\left(\omega t + (N-1)\Delta \beta\right) \right]\n

• Trigonometric Summation Identity Method:

  • Utilize the trigonometric product-to-sum identity:     \n\cos(\alpha) - \cos(\beta) = -2 \sin\left(\frac{\alpha + \beta}{2}\right) \sin\left(\frac{\alpha - \beta}{2}\right)\n

  • Evaluate successive cosine differences:     \n\cos\left(\omega t - \frac{\Delta \beta}{2}\right) - \cos\left(\omega t + \frac{\Delta \beta}{2}\right) = 2 \sin(\omega t) \sin\left(\frac{\Delta \beta}{2}\right)\n     \n\cos\left(\omega t + \frac{\Delta \beta}{2}\right) - \cos\left(\omega t + \frac{3\Delta \beta}{2}\right) = 2 \sin(\omega t + \Delta \beta) \sin\left(\frac{\Delta \beta}{2}\right)\n     \n\cos\left(\omega t + \frac{3\Delta \beta}{2}\right) - \cos\left(\omega t + \frac{5\Delta \beta}{2}\right) = 2 \sin(\omega t + 2\Delta \beta) \sin\left(\frac{\Delta \beta}{2}\right)\n     \n\cos\left(\omega t + \left(N - \frac{1}{2}\right)\Delta \beta\right) - \cos\left(\omega t + \left(N - \frac{3}{2}\right)\Delta \beta\right) = 2 \sin\left[\omega t + (N-1)\Delta \beta\right] \sin\left(\frac{\Delta \beta}{2}\right)\n

  • Summing both sides creates a telescoping series where internal terms cancel out entirely:     \n\cos\left(\omega t - \frac{\Delta \beta}{2}\right) - \cos\left[\omega t + \left(N - \frac{1}{2}\right)\Delta \beta\right] = 2 \sin\left(\frac{\Delta \beta}{2}\right) \sum_{j=0}^{N-1} \sin(\omega t + j\Delta \beta)\n

  • Applying the sum-to-product identity to the remaining boundary terms on the left:     \n\cos(\alpha) - \cos(\beta) = -2 \sin\left(\frac{\alpha + \beta}{2}\right) \sin\left(\frac{\alpha - \beta}{2}\right)\n     \n\cos\left(\omega t - \frac{\Delta \beta}{2}\right) - \cos\left[\omega t + \left(N - \frac{1}{2}\right)\Delta \beta\right] = 2 \sin\left(\omega t + (N-1)\frac{\Delta \beta}{2}\right) \sin\left(\frac{N \Delta \beta}{2}\right)\n

  • Equating both expressions yields the closed-form sum of sine terms:     \n\sum_{j=0}^{N-1} \sin(\omega t + j\Delta \beta) = \frac{\sin\left[\omega t + (N-1)\frac{\Delta \beta}{2}\right] \sin\left(\frac{\beta}{2}\right)}{\sin\left(\frac{\Delta \beta}{2}\right)}\n

  • Resulting total electric field equation:     \nE = E_{10} \left[ \frac{\sin(\beta / 2)}{\sin(\Delta \beta / 2)} \right] \sin\left(\omega t + (N-1)\frac{\Delta \beta}{2}\right)\n

• Calculating Time-Averaged Intensity and Continuum Limit:

  • Take the time-average of $E^2$ using $\left\langle \sin^2(\dots) \right\rangle = \frac{1}{2}$:     \nI \propto \langle E^2 \rangle = \frac{1}{2} E_{10}^2 \left[ \frac{\sin(\beta / 2)}{\sin(\Delta \beta / 2)} \right]^2\n

  • Define peak intensity $I_0$ as the maximum central intensity ($\beta = 0$) where $I_0 = \frac{1}{2} E_{10}^2 N^2$:     \nI = \frac{I_0}{N^2} \left[ \frac{\sin(\beta / 2)}{\sin(\Delta \beta / 2)} \right]^2\n

  • Take the continuous limit ($N \to \infty$, $\Delta \beta \to 0$), applying the small-angle approximation $\sin(x) \approx x$:     \nN \sin\left(\frac{\Delta \beta}{2}\right) \approx N \left(\frac{\Delta \beta}{2}\right) = \frac{N \Delta \beta}{2} = \frac{\beta}{2}\n

  • Substitute back into the expression to obtain the definitive Single-Slit Intensity Distribution:     \nI = I_0 \left[ \frac{\sin(\beta / 2)}{\beta / 2} \right]^2 = I_0 \left[ \frac{\sin\left(\frac{\pi a \sin(\theta)}{\lambda}\right)}{\frac{\pi a \sin(\theta)}{\lambda}} \right]^2\n     where the ratio $\frac{\sin(x)}{x}$ is defined as the sinc function.

Complete Derivation of Double Slit Interference Intensity Distribution

• Superposition Principle:

  • Total field $\vec{E}$ at screen point $P$ is the vector sum of light fields from slit 1 and slit 2:     \n\vec{E} = \vec{E}_1 + \vec{E}_2\n

  • Instantaneous Poynting magnitude scales as:     \nS \propto E^2 = (\vec{E}_1 + \vec{E}_2)^2 = E_1^2 + E_2^2 + 2 \vec{E}_1 \cdot \vec{E}_2\n

  • Time-averaged intensity:     \nI = \langle S \rangle \propto \langle E_1^2 \rangle + \langle E_2^2 \rangle + 2 \langle \vec{E}_1 \cdot \vec{E}_2 \rangle\n

  • The cross-term $2 \langle \vec{E}_1 \cdot \vec{E}_2 \rangle$ represents wave correlation (interference):

  • Incoherent Sources: Phase difference fluctuates randomly over time, driving the correlation term to zero:       \nI_{\text{incoherent}} = I_1 + I_2\n

  • Coherent Constructive Interference ($\vec{E}_1 = \vec{E}_2$):       \nI_{\text{constructive}} = I_1 + I_1 + 2I_1 = 4I_1\n

  • Coherent Destructive Interference ($\vec{E}_1 = -\vec{E}_2$):       \nI_{\text{destructive}} = I_1 + I_1 - 2I_1 = 0\n

• Derivation for Coherent Harmonic Waves:

  • Define wave fields arriving at point $P$ from slits 1 and 2 with equal amplitude $E_0$:     \nE_1 = E_0 \sin(\omega t)\n     \nE_2 = E_0 \sin(\omega t + \phi)\n     where $\phi$ represents phase shift from path difference $\delta = d \sin(\theta)$.

  • Sum fields via trigonometric identity:     \n\sin(\alpha) + \sin(\beta) = 2 \sin\left(\frac{\alpha + \beta}{2}\right) \cos\left(\frac{\alpha - \beta}{2}\right)\n     \nE = E_1 + E_2 = E_0 \left[ \sin(\omega t) + \sin(\omega t + \phi) \right] = 2 E_0 \cos\left(\frac{\phi}{2}\right) \sin\left(\omega t + \frac{\phi}{2}\right)\n

  • Compute time-averaged intensity:     \nI \propto \langle E^2 \rangle = 4 E_0^2 \cos^2\left(\frac{\phi}{2}\right) \left\langle \sin^2\left(\omega t + \frac{\phi}{2}\right) \right\rangle = 2 E_0^2 \cos^2\left(\frac{\phi}{2}\right)\n

  • Setting peak intensity $I_0 = 2 E_0^2$ when $\phi = 0$:     \nI = I_0 \cos^2\left(\frac{\phi}{2}\right)\n

  • Express phase difference $\phi$ in terms of slit separation $d$:     \n\frac{\delta}{\lambda} = \frac{\phi}{2\pi} \implies \phi = \frac{2\pi}{\lambda} \delta = \frac{2\pi}{\lambda} d \sin(\theta)\n

  • Substitute $\phi$ to obtain Double-Slit Interference Intensity Distribution:     \nI = I_0 \cos^2\left(\frac{\pi d \sin(\theta)}{\lambda}\right)\n

  • Under small-angle approximation ($\sin(\theta) \approx \tan(\theta) = \frac{y}{L}$):     \nI = I_0 \cos^2\left(\frac{\pi d y}{\lambda L}\right)\n

Composite Intensity Distribution for Double Slit Diffraction

• Combining finite slit width diffraction ($a$) and double-slit separation interference ($d$):

  • Single-slit diffraction factor:     \nI_{\text{diffraction}} = \left[ \frac{\sin\left(\frac{\pi a \sin(\theta)}{\lambda}\right)}{\frac{\pi a \sin(\theta)}{\lambda}} \right]^2\n

  • Double-slit interference factor:     \nI_{\text{interference}} = \cos^2\left(\frac{\pi d \sin(\theta)}{\lambda}\right)\n

• The combined intensity function is the product of both terms:   \nI = I_0 \cos^2\left(\frac{\pi d \sin(\theta)}{\lambda}\right) \left[ \frac{\sin\left(\frac{\pi a \sin(\theta)}{\lambda}\right)}{\frac{\pi a \sin(\theta)}{\lambda}} \right]^2\n

• Physical Interpretation of Combined Terms:

  • Interference Factor: Provides high-frequency interference substructure fringes determined by slit separation $d$.

  • Diffraction Factor: Acts as an overarching intensity envelope determined by slit width $a$, modulating and limiting the amplitude of the inner interference peaks.

Data Processing, Analysis, and Report Requirements

• Experimental Intensity Normalization & Positioning:

  • Normalize measured raw intensity values by dividing all readings by peak central maximum intensity $I_0$, yielding relative intensity formatting ($\frac{I}{I_0}$).

  • Align the central maximum peak precisely at spatial position $x = 0$ by subtracting the peak offset coordinate from all position measurements.

• Theoretical Comparison Spreadsheet Construction:

  • Column 1: Position values incremented in small spatial steps covering the experimental range centered at zero.

  • Column 2: Theoretical single-slit diffraction intensity calculated using:     \nI_{\text{single}} = \left[ \frac{\sin\left(\frac{\pi a \sin(\theta)}{\lambda}\right)}{\frac{\pi a \sin(\theta)}{\lambda}} \right]^2\n

  • Column 3: Theoretical double-slit interference intensity calculated using:     \nI_{\text{double}} = \cos^2\left(\frac{\pi d \sin(\theta)}{\lambda}\right)\n

  • Column 4: Total double-slit diffraction-interference intensity calculated by multiplying Column 2 by Column 3.

  • Plotting: Graph single-slit, double-slit, and composite total curves alongside experimentally collected sensor curves for direct quantitative evaluation.

• Mandatory Lab Report Elements:

  • Tabulate experimental slit width ($a$) and slit separation ($d$) against manufacturer specifications, including explicit uncertainty bounds.

  • Provide a concise summary of the theoretical intensity derivations.

  • Demonstrate analytically that zeros of single-slit intensity correspond to single-slit diffraction minima:     \n\frac{\pi a \sin(\theta)}{\lambda} = n\pi \implies a \sin(\theta) = n\lambda\n

  • Demonstrate analytically that maxima of double-slit intensity correspond to interference maxima:     \n\frac{\pi d \sin(\theta)}{\lambda} = m\pi \implies d \sin(\theta) = m\lambda\n

  • Discuss potential systematic deviations between modeled theoretical predictions and experimental data (sensor saturation, beam alignment, aperture tolerances).