Phương Pháp Xấp Xỉ Pha Tĩnh (Stationary Phase Approximation)

1. Giới thiệu

Phương pháp xấp xỉ pha tĩnh là một kỹ thuật mạnh mẽ được sử dụng để ước lượng giá trị của các tích phân dao động, đặc biệt khi tần số dao động lớn. Trong Vật lý, phương pháp này thường được áp dụng trong các bài toán liên quan đến lý thuyết sóng, cơ học lượng tử, và quang học. Ý tưởng cơ bản của phương pháp là khi pha của hàm số trong tích phân biến đổi nhanh, sự đóng góp chính vào tích phân đến từ các điểm mà tại đó pha biến đổi chậm nhất, tức là "pha tĩnh".

2. Phát biểu toán học

Xét tích phân có dạng tổng quát:

$$I = \int_{a}^{b} A(x) e^{i k f(x)} dx$$

trong đó:

• $A(x)$ là một hàm biên độ (amplitude) biến đổi chậm.
• $f(x)$ là hàm pha, quyết định sự dao động của hàm số dưới dấu tích phân.
• $k$ là một tham số lớn, đại diện cho tần số dao động.
• $a$ và $b$ là các giới hạn tích phân.

Mục tiêu là ước lượng giá trị của $I$ khi $k \rightarrow \infty$.

3. Điểm pha tĩnh

Điểm pha tĩnh $x_0$ là điểm mà tại đó đạo hàm của hàm pha bằng không:

$$f'(x_0) = \frac{df(x)}{dx}\Big|_{x=x_0} = 0$$

Đây là các điểm mà pha $f(x)$ thay đổi chậm nhất.

4. Xấp xỉ Taylor xung quanh điểm pha tĩnh

Khai triển Taylor hàm pha $f(x)$ xung quanh điểm pha tĩnh $x_0$:

$$f(x) \approx f(x_0) + f'(x_0)(x - x_0) + \frac{1}{2}f''(x_0)(x - x_0)^2 + ...$$

Do $f'(x_0) = 0$, ta có:

$$f(x) \approx f(x_0) + \frac{1}{2}f''(x_0)(x - x_0)^2$$

Thay vào tích phân ban đầu, ta có:

$$I \approx \int_{a}^{b} A(x) e^{i k [f(x_0) + \frac{1}{2}f''(x_0)(x - x_0)^2]} dx$$

5. Đánh giá tích phân

Để đơn giản hóa, ta giả sử $A(x)$ biến đổi chậm và có thể xấp xỉ bởi giá trị của nó tại điểm pha tĩnh $A(x_0)$. Đồng thời, ta mở rộng giới hạn tích phân về $\pm \infty$ vì hàm số dao động sẽ nhanh chóng triệt tiêu bên ngoài một vùng nhỏ xung quanh $x_0$.

$$I \approx A(x_0) e^{i k f(x_0)} \int_{-\infty}^{\infty} e^{i \frac{k}{2}f''(x_0)(x - x_0)^2} dx$$

Đặt $u = \sqrt{\frac{k|f''(x_0)|}{2}}(x - x_0)$, ta có $dx = \sqrt{\frac{2}{k|f''(x_0)|}} du$.

Khi đó, tích phân trở thành:

$$I \approx A(x_0) e^{i k f(x_0)} \sqrt{\frac{2}{k|f''(x_0)|}} \int_{-\infty}^{\infty} e^{\pm i u^2} du$$

Tích phân Fresnel được định nghĩa là:

$$\int_{-\infty}^{\infty} e^{i u^2} du = \sqrt{\pi}e^{i\pi/4}$$ $$\int_{-\infty}^{\infty} e^{-i u^2} du = \sqrt{\pi}e^{-i\pi/4}$$

Do đó:

$$I \approx A(x_0) e^{i k f(x_0)} \sqrt{\frac{2\pi}{k|f''(x_0)|}} e^{i \frac{\pi}{4}sgn(f''(x_0))}$$

trong đó $sgn(f''(x_0))$ là hàm dấu của $f''(x_0)$.

6. Tổng quát hóa cho nhiều điểm pha tĩnh

Nếu có nhiều điểm pha tĩnh $x_{0,j}$ trong khoảng $(a, b)$, tích phân được xấp xỉ bằng tổng đóng góp từ mỗi điểm:

$$I \approx \sum_{j} A(x_{0,j}) e^{i k f(x_{0,j})} \sqrt{\frac{2\pi}{k|f''(x_{0,j})|}} e^{i \frac{\pi}{4}sgn(f''(x_{0,j}))}$$

7. Điều kiện áp dụng

Phương pháp xấp xỉ pha tĩnh có hiệu quả khi:

• $k$ đủ lớn, đảm bảo dao động nhanh.
• $A(x)$ biến đổi chậm so với $e^{i k f(x)}$.
• Các điểm pha tĩnh nằm cách xa nhau để các xấp xỉ Gauss xung quanh mỗi điểm không giao nhau đáng kể.

8. Ứng dụng

• Cơ học lượng tử: Tính toán biên độ tán xạ và các tích phân đường.
• Quang học: Nghiên cứu nhiễu xạ và giao thoa sóng ánh sáng.
• Lý thuyết sóng: Phân tích sự truyền sóng trong các môi trường khác nhau.

9. Ví dụ

Xét tích phân:

$$I = \int_{-\infty}^{\infty} e^{i k (x^4 - x^2)} dx$$

Tìm điểm pha tĩnh:

$$f(x) = x^4 - x^2$$ $$f'(x) = 4x^3 - 2x = 0 \Rightarrow x = 0, \pm \frac{1}{\sqrt{2}}$$

Tính $f''(x)$:

$$f''(x) = 12x^2 - 2$$

• Tại $x_0 = 0$, $f''(0) = -2$
• Tại $x_0 = \pm \frac{1}{\sqrt{2}}$, $f''(\pm \frac{1}{\sqrt{2}}) = 4$

Áp dụng công thức xấp xỉ pha tĩnh:

$$I \approx e^{i k f(0)} \sqrt{\frac{2\pi}{k|f''(0)|}} e^{i \frac{\pi}{4}sgn(f''(0))} + 2e^{i k f(\frac{1}{\sqrt{2}})} \sqrt{\frac{2\pi}{k|f''(\frac{1}{\sqrt{2}})|}} e^{i \frac{\pi}{4}sgn(f''(\frac{1}{\sqrt{2}}))}$$ $$I \approx \sqrt{\frac{\pi}{k}} e^{-i \frac{\pi}{4}} + 2\sqrt{\frac{\pi}{2k}} e^{-i k/4} e^{i \frac{\pi}{4}}$$

10. Kết luận

Phương pháp xấp xỉ pha tĩnh là một công cụ hữu ích để ước lượng các tích phân dao động trong nhiều lĩnh vực của Vật lý. Việc hiểu rõ các điều kiện áp dụng và cách tính toán các điểm pha tĩnh là chìa khóa để sử dụng phương pháp này một cách hiệu quả.