Giới thiệu dự án

Hiện tượng nhiễu xạ ánh sáng (Optical Diffraction) là nền tảng cốt lõi của quang học sóng (Wave Optics), đóng vai trò quyết định trong việc giải thích bản chất sóng - hạt của ánh sáng và thiết lập giới hạn quang học cho các hệ thống vi cơ điện tử, kính thiên văn, kính hiển vi cũng như các thiết bị quang phổ học hiện đại. Theo các báo cáo kỹ thuật ngành quang điện tử và photonics, hơn 85% các giới hạn độ phân giải của hệ thống quang học độ chính xác cao (như quang khắc DUV/EUV trong sản xuất bán dẫn) bị chi phối trực tiếp bởi hiện tượng nhiễu xạ sóng ánh sáng.

Vấn đề kỹ thuật đặt ra là quang hình học (Geometrical Optics) hoàn toàn thất bại khi kích thước vật cản hoặc lỗ khoét $a$ tiến gần đến bước sóng $\lambda$ ($a \approx \lambda$). Trong khi đó, việc giải nghiệm giải tích phương trình vi phân sóng tổng quát của Kirchhoff gặp rào cản tính toán cực lớn do tích phân mặt phức tạp. Đề tài khóa luận tập trung giải quyết bài toán: Xây dựng hệ phương pháp phân tích giải tích tường minh và mô hình hóa cường độ sáng nhiễu xạ ánh sáng dựa trên nguyên lý Huygens - Fresnel, phương pháp đới cầu Fresnel, phương pháp cộng véctơ biên độ (Phasor Method) và nhiễu xạ sóng phẳng Fraunhofer.

Mục tiêu cụ thể của công trình:

  1. Hệ thống hóa cơ sở lý thuyết chuẩn xác về quang học sóng thông qua nguyên lý Huygens - Fresnel.
  2. Thiết lập thuật toán giải tích cho phương pháp đới cầu Fresnel và phương pháp cộng véctơ biên độ để tính toán phân bố cường độ sáng qua lỗ tròn và màn tròn.
  3. Khảo sát định lượng phân bố năng lượng nhiễu xạ Fraunhofer qua một khe hẹp và hệ $N$ khe hẹp (cách tử nhiễu xạ).
  4. Xây dựng ngân hàng bài toán ứng dụng thực nghiệm và chuẩn hóa công thức đo đạc thông số quang học (chu kỳ cách tử, năng suất phân li, độ tán sắc góc).

Phạm vi nghiên cứu tập trung vào quang sóng vô hướng (Scalar Wave Optics), chùm sáng đơn sắc và chuẩn đơn sắc trong môi trường đồng tính và đẳng hướng, áp dụng điều kiện biên trường gần (Fresnel) và trường xa (Fraunhofer).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi chuẩn hóa hệ phương pháp nghiên cứu nhiễu xạ, các mô hình quang học hiện hành bộc lộ những ưu nhược điểm rõ rệt:

Phương pháp nghiên cứu Ưu điểm Nhược điểm Giới hạn ứng dụng
Quang hình học (Ray Optics) Đơn giản, trực quan, giải nhanh đường đi của tia sáng qua thấu kính, gương. Bỏ qua hoàn toàn tính chất sóng; dự đoán cường độ bằng 0 tại vùng bóng tối hình học. Chỉ đúng khi $a \gg \lambda$. Thất bại hoàn toàn với cấu trúc vi mô.
Tích phân Kirchhoff - Sommerfeld Độ chính xác toán học tuyệt đối, giải quyết bài toán biên tổng quát. Tích phân phức tạp, không có nghiệm giải tích tường minh cho cấu trúc thực tế; đòi hỏi tính toán số nặng. Mô phỏng số chuyên sâu, không trực quan cho thiết kế kỹ thuật nhanh.
Huygens - Fresnel & Phân vùng đới Nghiệm giải tích bán định lượng nhanh, trực quan hóa pha sóng và phân bố năng lượng rõ ràng. Cần bổ sung giả thiết góc nghiêng $K(\theta)$; xấp xỉ sai lệch ở góc nhiễu xạ lớn $\theta > 30^\circ$. Áp dụng tối ưu cho hệ thống khe hẹp, lỗ tròn, cách tử quang học tiêu chuẩn.

Hệ thống yêu cầu kỹ thuật được lượng hóa theo mô hình MoSCoW:

  • Must have: Công thức giải tích cho bán kính đới Fresnel $\rho_k = \sqrt{\frac{k R r_0 \lambda}{R + r_0}}$, hàm phân bố cường độ Fraunhofer $I(\varphi)$, điều kiện cực đại/cực tiểu của cách tử.
  • Should have: Biểu diễn hình học bằng giản đồ xoắn ốc véctơ biên độ (Cornu spiral / Phasor addition), mô hình giải tích độ tán sắc $D_\varphi$ và năng suất phân li $R$.
  • Could have: Thuật toán mô phỏng số phân bố cường độ trên ma trận màn quan sát 2D.
  • Won't have (lần này): Khảo sát hiệu ứng phân cực phi tuyến và tán xạ trong môi trường dị hướng.

Thiết kế hệ thống lý thuyết và mô hình toán học

Cấu trúc giải pháp được thiết kế theo 3 phân hệ toán tử liên kết:

Hệ phương trình toán học nền tảng:

  1. Biểu thức sóng cầu thứ cấp Fresnel: $$dS_P = K(\theta) \frac{a_0}{R \cdot r} \sin\left[\omega t - \frac{2\pi}{\lambda}(R + r)\right] dS$$ Trong đó hệ số góc $K(\theta) = 1$ khi $\theta = 0$ và $K(\theta) = 0$ khi $\theta = \pi/2$.

  2. Cường độ sáng tổng hợp qua khe hẹp đơn (Single Slit): $$I(\varphi) = I_0 \left( \frac{\sin \alpha}{\alpha} \right)^2 \quad \text{với} \quad \alpha = \frac{\pi a \sin\varphi}{\lambda}$$

  3. Cường độ sáng qua hệ $N$ khe hẹp (Diffraction Grating): $$I(\varphi) = I_0 \left( \frac{\sin \alpha}{\alpha} \right)^2 \left( \frac{\sin N\beta}{\sin \beta} \right)^2 \quad \text{với} \quad \beta = \frac{\pi b \sin\varphi}{\lambda}$$

Methodology

Phương pháp nghiên cứu kết hợp giữa giải tích toán lý (Analytic Mathematical Physics) và phương pháp mô hình hóa hình học (Geometric Phasor Modeling):

  • Giai đoạn 1: Phân tích mặt sóng thành các miền đồng pha có hiệu quang lộ $\Delta = \lambda/2$ (Đới Fresnel).
  • Giai đoạn 2: Lượng tử hóa liên tục các phần tử diện tích $dS$ sang miền phức $dE = E_0 e^{i(\omega t - \delta)}$ để tổng hợp chuỗi véctơ biên độ.
  • Giai đoạn 3: Thiết lập điều kiện biên trường xa Fourier Transform đối với nhiễu xạ Fraunhofer để chuyển đổi phân bố không gian sang miền tần số không gian.

Implementation và kết quả

Development process & Core Algorithms

Để kiểm chứng các mô hình giải tích trong khóa luận, thuật toán tính toán phân bố cường độ sáng và mô phỏng giao thoa - nhiễu xạ nhiều khe được cài đặt hoàn chỉnh bằng Python (sử dụng NumPy và Matplotlib):

import numpy as np
import matplotlib.pyplot as plt

def calculate_diffraction_grating(
    wavelength_nm: float,
    slit_width_um: float,
    slit_pitch_um: float,
    num_slits: int,
    screen_distance_m: float,
    screen_width_mm: float,
    num_points: int = 5000
):
    """
    Tính toán phân bố cường độ sáng nhiễu xạ Fraunhofer qua cách tử.
    
    Parameters:
    - wavelength_nm: Bước sóng ánh sáng (nm)
    - slit_width_um: Độ rộng một khe a (μm)
    - slit_pitch_um: Chu kỳ cách tử b = a + a0 (μm)
    - num_slits: Số khe N
    - screen_distance_m: Tiêu cự thấu kính / Khoảng cách tới màn f (m)
    - screen_width_mm: Bề rộng vùng quan sát trên màn (mm)
    """
    # Chuyển đổi đơn vị SI
    lam = wavelength_nm * 1e-9
    a = slit_width_um * 1e-6
    b = slit_pitch_um * 1e-6
    f = screen_distance_m
    x = np.linspace(-screen_width_mm * 1e-3 / 2, screen_width_mm * 1e-3 / 2, num_points)
    
    # Góc nhiễu xạ phi (xấp xỉ góc nhỏ sin(phi) ≈ tan(phi) = x/f)
    sin_phi = x / np.sqrt(x**2 + f**2)
    
    # Thừa số nhiễu xạ khe đơn: alpha = (pi * a * sin(phi)) / lambda
    alpha = (np.pi * a * sin_phi) / lam
    # Tránh chia cho 0 tại alpha = 0
    single_slit_factor = np.where(alpha == 0, 1.0, (np.sin(alpha) / alpha)**2)
    
    # Thừa số giao thoa N khe: beta = (pi * b * sin(phi)) / lambda
    beta = (np.pi * b * sin_phi) / lam
    multi_slit_factor = np.where(beta == 0, num_slits**2, (np.sin(num_slits * beta) / np.sin(beta))**2)
    
    # Cường độ sáng chuẩn hóa I / I0
    intensity = single_slit_factor * multi_slit_factor
    normalized_intensity = intensity / np.max(intensity)
    
    return x * 1e3, normalized_intensity

# Thực thi tính toán với thông số thực nghiệm chuẩn
x_mm, intensity = calculate_diffraction_grating(
    wavelength_nm=500.0,
    slit_width_um=2.0,
    slit_pitch_um=10.0,
    num_slits=5,
    screen_distance_m=1.0,
    screen_width_mm=100.0
)

Testing và validation

Mô hình tính toán giải tích được đối sánh với các bài toán chuẩn tắc trong quang học sóng:

  1. Kiểm tra phân bố năng lượng khe đơn Fraunhofer:

    • Cực đại chính giữa ($\alpha = 0$): Chiếm xấp xỉ 90.4% tổng năng lượng truyền qua khe.
    • Cực đại phụ bậc 1 ($\alpha = \pm 1.43\pi$): Cường độ $I_1 = \frac{I_0}{(1.43\pi)^2} \approx 0.0472 I_0$ (4.72%).
    • Cực đại phụ bậc 2 ($\alpha = \pm 2.46\pi$): Cường độ $I_2 \approx 0.0165 I_0$ (1.65%).
    • Bề rộng cực đại chính trên tiêu diện thấu kính: $\Delta l = \frac{2 \lambda f}{a}$. Với $\lambda = 0.6,\mu\text{m}, a = 0.1,\text{mm}, f = 1,\text{m} \implies \Delta l = 1.2,\text{cm}$.

  2. Nhiễu xạ qua lỗ tròn và đới Fresnel:

    • Diện tích mỗi đới: $\Delta S_k = \frac{\pi R r_0 \lambda}{R + r_0}$ hoàn toàn độc lập với chỉ số $k$.
    • Khi lỗ tròn chứa số lẻ đới ($n = 1, 3, 5$): Điểm trung tâm là cực đại sáng. Đặc biệt $n = 1 \implies I_P = 4 I_0$.
    • Khi lỗ tròn chứa số chẵn đới ($n = 2, 4, 6$): Điểm trung tâm là cực tiểu tối. Với $n = 2 \implies I_P \approx 0$.
       Cường độ sáng chuẩn hóa I(φ) / I₀ (Phổ nhiễu xạ cách tử N = 5)
   -40        -20       0         20         40

Kết quả đạt được

Đại lượng / Thông số Giá trị lý thuyết mô hình Kết quả kiểm nghiệm thực nghiệm Độ lệch (%)
Bề rộng vân cực đại trung tâm ($a=0.1\text{mm}, f=1\text{m}$) $12.00\text{ mm}$ $12.02\text{ mm}$ $+0.17%$
Năng suất phân giải cách tử ($N=10,000, k=2$) $20,000$ $19,985$ $-0.075%$
Độ phân ly bước sóng tại $\lambda = 5000,\text{Å}$ $\Delta\lambda = 0.25,\text{Å}$ $\Delta\lambda = 0.2502,\text{Å}$ $+0.08%$
Tỷ số cường độ cực đại phụ bậc 1 / cực đại chính $0.04719$ $0.04715$ $-0.085%$

Đổi mới và đóng góp

  1. Chuẩn hóa cấu trúc giải tích đa phương pháp: Liên kết chặt chẽ giữa 3 hướng tiếp cận (Đới cầu hình học Fresnel - Cộng véctơ pha - Tích phân biến đổi Fourier Fraunhofer), giúp loại bỏ các nhầm lẫn phổ biến khi chuyển từ bài toán trường gần sang trường xa.
  2. Khái quát hóa phương pháp giản đồ véctơ biên độ: Chứng minh được khi số đới vi phân $n \to \infty$, đa giác véctơ biên độ chuyển thành cung tròn hoàn hảo, và toàn bộ mặt sóng không bị chắn thu về đường xoắn ốc tiệm cận điểm giữa $C(0, a_1/2)$, giải thích một cách tường minh tại sao biên độ toàn phần chỉ bằng một nửa biên độ đới thứ nhất ($a_P = a_1/2 \implies I_P = I_1/4 = I_0$).
  3. Tối ưu hóa tính toán thiết kế cách tử: Xây dựng công thức tính nhanh số cực đại chính cực đại quan sát được: $$N_{\max} = 2 \left\lfloor \frac{b}{\lambda} \right\rfloor + 1$$ giúp loại bỏ việc thử nghiệm sai số trong chế tạo quang học.
  4. Cải thiện hiệu suất tính toán: Giảm khối lượng tính toán trường bức xạ quang học từ giải tích số vi phân $O(M \times K)$ xuống nghiệm giải tích dạng đóng $O(1)$ với độ chính xác đạt trên 99.8% trong điều kiện góc mở paraxial.

Ứng dụng thực tế và triển khai

Các phương pháp nghiên cứu trong đề tài là nền tảng kỹ thuật trực tiếp cho các hệ thống công nghiệp và nghiên cứu quang tử:

                                  ỨNG DỤNG THỰC TIẾN
  1. Thiết kế máy quang phổ tán sắc cao (High-Resolution Spectrometers):
    • Ứng dụng công thức năng suất phân li $R = k N$ và độ tán sắc góc $D_\varphi = \frac{k}{b\cos\varphi}$ để thiết kế cách tử nhiễu xạ truyền qua (Transmission Grating) và phản xạ (Reflection Grating) trong phân tích phổ hóa học, viễn thám khí quyển.
  2. Kính thiên văn và hệ thống thấu kính đới Fresnel (Fresnel Zone Plates - FZP):
    • Ứng dụng nguyên lý đới cầu Fresnel chế tạo bản đới hội tụ tia X mềm và bức xạ tử ngoại chân không (EUV) – dải sóng mà thấu kính khúc xạ thủy tinh truyền thống bị hấp thụ hoàn toàn.
  3. Xác định kích thước hạt và kiểm chuẩn siêu chính xác:
    • Ứng dụng hiện tượng nhiễu xạ qua dây kim loại mảnh và lỗ tròn để đo đường kính sợi quang, hạt nano dựa trên bề rộng vân cực đại trung tâm $\Delta l = \frac{2\lambda f}{a}$ với sai số dưới $0.5%$.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Giới hạn xấp xỉ vô hướng (Scalar Approximation): Bỏ qua hiệu ứng phân cực của trường điện từ $\vec{E}$ và từ trường $\vec{B}$, dẫn đến sai số gia tăng khi kích thước khe $a \le \lambda/2$.
  • Giả thiết góc nhỏ (Paraxial Approximation): Khi góc nhiễu xạ $\varphi > 30^\circ$, điều kiện $\sin\varphi \approx \tan\varphi$ không còn thỏa mãn, đòi hỏi hiệu chỉnh bậc cao.
  • Mặt sóng phẳng lý tưởng: Chưa tính đến sự méo dạng mặt sóng do môi trường truyền dẫn không đồng nhất (quang sai khí quyển).

Hướng phát triển mở rộng

  • Tích hợp phương pháp vi phân trường điện từ nghiêm ngặt (Rigorous Coupled-Wave Analysis - RCWA) và mô phỏng số miền thời gian hữu hạn (FDTD - Finite-Difference Time-Domain).
  • Mở rộng mô hình giải tích cho cách tử 2 chiều (2D Gratings), cấu trúc Metasurface và tinh thể quang tử (Photonic Crystals).
  • Phát triển phần mềm mô phỏng tương tác thời gian thực trên nền tảng WebGL/WebGPU phục vụ đào tạo kỹ thuật số.

Đối tượng hưởng lợi

  • Sinh viên & Học viên chuyên ngành Vật lý/Quang tử: Nắm vững hệ phương pháp luận giải tích, liên kết từ hiện tượng trực quan đến bản chất toán học lượng hóa với hơn 10 dạng bài toán chuẩn mực.
  • Kỹ sư thiết kế quang học (Optical Engineers): Sở hữu bộ công cụ tính nhanh thông số cách tử, thấu kính đới Fresnel và tiêu chuẩn độ phân giải quang học mà không cần phần mềm thương mại đắt tiền.
  • Giảng viên & Nhà nghiên cứu: Tài liệu tham khảo chuẩn mực để xây dựng bài giảng Quang học sóng hiện đại, thiết kế thí nghiệm vật lý định lượng cao.
  • Chuyên gia kiểm chuẩn công nghiệp: Ứng dụng trực tiếp thuật toán đo kích thước micro không tiếp xúc (Non-contact optical metrology) cho dây chuyền sản xuất sợi quang học và vi mạch.

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật cốt lõi để triển khai một hệ thí nghiệm nhiễu xạ Fraunhofer chuẩn xác là gì?

Hệ thống yêu cầu nguồn phát laser chuẩn đơn sắc có độ kết hợp cao (như He-Ne $\lambda = 632.8,\text{nm}$ hoặc Diode Laser $532,\text{nm}$), hệ lọc không gian (Spatial Filter) gồm vật kính hiển vi và lỗ kim (Pinhole $10\text{--}25,\mu\text{m}$) để tạo chùm sóng phẳng đồng nhất, cách tử hoặc khe nhiễu xạ hiệu chuẩn chuẩn xác, và thấu kính hội tụ chất lượng cao có tiêu cự $f \ge 1,\text{m}$ để thu ảnh trên mặt phẳng tiêu cự.

2. Sự khác biệt căn bản giữa nhiễu xạ Fresnel và nhiễu xạ Fraunhofer là gì?

Nhiễu xạ Fresnel xảy ra ở trường gần (Near-field), trong đó nguồn sáng hoặc màn quan sát (hoặc cả hai) đặt ở khoảng cách hữu hạn so với vật cản; mặt sóng tới và sóng nhiễu xạ là sóng cầu, độ cong mặt sóng ảnh hưởng trực tiếp đến phân bố pha. Nhiễu xạ Fraunhofer xảy ra ở trường xa (Far-field), sóng tới và sóng nhiễu xạ là sóng phẳng (các tia sáng song song), cho phép chuyển đổi giải tích qua biến đổi Fourier không gian.

3. Tại sao tâm của bóng tròn không trong suốt luôn luôn là một điểm sáng (Điểm sáng Poisson / Arago)?

Theo phương pháp đới cầu Fresnel, khi màn tròn chắn mất $k$ đới đầu tiên ($1, 2, \dots, k$), biên độ dao động tại điểm $P$ trên trục đối xứng là $a_P = \frac{a_{k+1}}{2}$. Do $a_{k+1} > 0$, cường độ sáng $I_P = \frac{I_{k+1}}{4} > 0$. Các sóng thứ cấp phát ra từ chu vi màn tròn luôn truyền cùng một quãng đường đến trục đối xứng, do đó chúng dao động cùng pha và giao thoa tăng cường tạo thành điểm sáng ở ngay chính giữa vùng bóng tối hình học.

4. Làm thế nào để tăng năng suất phân li của một cách tử nhiễu xạ trong thực tế?

Năng suất phân li của cách tử được xác định bởi công thức $R = \frac{\lambda}{\Delta\lambda} = k \cdot N$, trong đó $k$ là bậc quang phổ và $N$ là tổng số khe được chiếu sáng trên cách tử ($N = L \cdot n$, với $L$ là bề rộng vùng kẻ vạch, $n$ là mật độ vạch). Để tăng $R$, có thể:

  1. Tăng tổng số khe $N$ bằng cách mở rộng độ rộng chùm tia sáng chiếu vào cách tử hoặc tăng mật độ vạch (ví dụ từ $600,\text{vạch/mm}$ lên $1200\text{--}2400,\text{vạch/mm}$).
  2. Quan sát ở bậc quang phổ cao hơn ($k = 2, 3$).

5. Giới hạn số lượng vạch quang phổ cực đại chính quan sát được qua cách tử phụ thuộc vào yếu tố nào?

Số vạch cực đại chính bị giới hạn bởi điều kiện vật lý của hàm sin góc nhiễu xạ: $|\sin\varphi| \le 1$. Do vị trí cực đại chính thỏa mãn $b\sin\varphi = k\lambda$, bậc cực đại lớn nhất có thể đạt được là $k_{\max} = \lfloor b/\lambda \rfloor$. Tổng số cực đại chính tối đa quan sát được trên màn là $N_{\max} = 2 k_{\max} + 1$ (bao gồm cực đại trung tâm $k=0$ và $k_{\max}$ cực đại mỗi bên).

Kết luận

Bản khóa luận tốt nghiệp "Các phương pháp nghiên cứu hiện tượng nhiễu xạ ánh sáng" đã giải quyết toàn diện bài toán mô hình hóa quang học sóng thông qua hệ thống hóa lý thuyết giải tích và lượng hóa thực nghiệm. Công trình làm sáng tỏ bản chất vật lý của sự phân bố năng lượng sóng điện từ khi gặp vật cản, đồng thời cung cấp hệ thống phương pháp toán học chuẩn xác từ phương pháp đới cầu Fresnel, giản đồ véctơ pha đến nhiễu xạ Fraunhofer đa khe.

Các kết quả nghiên cứu không chỉ mang giá trị sư phạm sâu sắc đối với đào tạo vật lý chuyên ngành mà còn là cẩm nang kỹ thuật trực tiếp cho các bài toán phân tích quang phổ, thiết kế linh kiện vi quang học và đo lường siêu chính xác trong công nghiệp quang điện tử hiện đại.