Kỹ thuật Nhiễu Loạn (Perturbation Theory) - Giải Gần Đúng Các Bài Toán Lượng Tử

1. Giới thiệu

Trong Vật lý Lượng tử, chúng ta thường phải đối mặt với các bài toán không có nghiệm giải tích chính xác. Một kỹ thuật mạnh mẽ để giải quyết các bài toán này là Kỹ thuật Nhiễu Loạn (Perturbation Theory). Ý tưởng cơ bản là coi hệ phức tạp như một sự "nhiễu loạn" nhỏ của một hệ đơn giản hơn mà chúng ta đã biết nghiệm.

2. Cơ sở lý thuyết

2.1. Hamiltonian và Phương trình Schrödinger

Hệ cơ lượng tử được mô tả bởi Hamiltonian $\hat{H}$, một toán tử năng lượng. Phương trình Schrödinger mô tả sự tiến triển của trạng thái hệ theo thời gian:

$$i\hbar \frac{\partial}{\partial t} |\Psi(t)\rangle = \hat{H} |\Psi(t)\rangle$$

Trong đó:

• $|\Psi(t)\rangle$ là hàm sóng mô tả trạng thái của hệ tại thời điểm $t$.
• $\hbar$ là hằng số Planck rút gọn.

Đối với các trạng thái dừng (trạng thái không thay đổi theo thời gian), phương trình Schrödinger trở thành:

$$\hat{H} |\psi_n\rangle = E_n |\psi_n\rangle$$

Trong đó:

• $|\psi_n\rangle$ là hàm sóng của trạng thái dừng thứ $n$.
• $E_n$ là năng lượng của trạng thái dừng thứ $n$.

2.2. Hệ Không Nhiễu Loạn và Hệ Nhiễu Loạn

Giả sử chúng ta có một hệ không nhiễu loạn (unperturbed system) mà chúng ta đã biết nghiệm, tức là chúng ta biết Hamiltonian $\hat{H}_0$, các trạng thái riêng $|\psi_n^{(0)}\rangle$ và các mức năng lượng tương ứng $E_n^{(0)}$:

$$\hat{H}_0 |\psi_n^{(0)}\rangle = E_n^{(0)} |\psi_n^{(0)}\rangle$$

Bây giờ, xét một hệ nhiễu loạn (perturbed system) với Hamiltonian $\hat{H}$ có thể viết dưới dạng:

$$\hat{H} = \hat{H}_0 + \lambda \hat{H}'$$

Trong đó:

• $\hat{H}'$ là toán tử nhiễu loạn.
• $\lambda$ là tham số nhiễu loạn, là một số nhỏ ( $|\lambda| \ll 1$ ).

Mục tiêu của Kỹ thuật Nhiễu Loạn là tìm nghiệm gần đúng (các mức năng lượng $E_n$ và các trạng thái riêng $|\psi_n\rangle$ ) của hệ nhiễu loạn dựa trên nghiệm đã biết của hệ không nhiễu loạn.

3. Kỹ thuật Nhiễu Loạn Độc Lập Thời Gian

Chúng ta sẽ xét trường hợp nhiễu loạn không phụ thuộc vào thời gian.

3.1. Khai triển Chuỗi

Chúng ta giả sử rằng các mức năng lượng và trạng thái riêng của hệ nhiễu loạn có thể được biểu diễn dưới dạng chuỗi lũy thừa của tham số nhiễu loạn $\lambda$:

$$\begin{aligned} E_n &= E_n^{(0)} + \lambda E_n^{(1)} + \lambda^2 E_n^{(2)} + ... \\ |\psi_n\rangle &= |\psi_n^{(0)}\rangle + \lambda |\psi_n^{(1)}\rangle + \lambda^2 |\psi_n^{(2)}\rangle + ... \end{aligned}$$

Trong đó:

• $E_n^{(1)}$ là hiệu chỉnh năng lượng bậc nhất.
• $E_n^{(2)}$ là hiệu chỉnh năng lượng bậc hai.
• $|\psi_n^{(1)}\rangle$ là hiệu chỉnh trạng thái bậc nhất.
• $|\psi_n^{(2)}\rangle$ là hiệu chỉnh trạng thái bậc hai.

3.2. Phương pháp

1. Thay thế các khai triển chuỗi vào phương trình Schrödinger:

   $$(\hat{H}_0 + \lambda \hat{H}') (|\psi_n^{(0)}\rangle + \lambda |\psi_n^{(1)}\rangle + \lambda^2 |\psi_n^{(2)}\rangle + ...) = (E_n^{(0)} + \lambda E_n^{(1)} + \lambda^2 E_n^{(2)} + ...) (|\psi_n^{(0)}\rangle + \lambda |\psi_n^{(1)}\rangle + \lambda^2 |\psi_n^{(2)}\rangle + ...)$$

2. Nhóm các số hạng theo lũy thừa của $\lambda$:

   • Bậc 0 ($\lambda^0$): $\hat{H}_0 |\psi_n^{(0)}\rangle = E_n^{(0)} |\psi_n^{(0)}\rangle$ (phương trình Schrödinger cho hệ không nhiễu loạn - đã biết nghiệm).
   • Bậc 1 ($\lambda^1$): $\hat{H}_0 |\psi_n^{(1)}\rangle + \hat{H}' |\psi_n^{(0)}\rangle = E_n^{(0)} |\psi_n^{(1)}\rangle + E_n^{(1)} |\psi_n^{(0)}\rangle$
   • Bậc 2 ($\lambda^2$): $\hat{H}_0 |\psi_n^{(2)}\rangle + \hat{H}' |\psi_n^{(1)}\rangle = E_n^{(0)} |\psi_n^{(2)}\rangle + E_n^{(1)} |\psi_n^{(1)}\rangle + E_n^{(2)} |\psi_n^{(0)}\rangle$
   • Và tiếp tục cho các bậc cao hơn.

3. Giải các phương trình theo từng bậc:

   • Hiệu chỉnh năng lượng bậc nhất: Nhân phương trình bậc 1 với $\langle\psi_n^{(0)}|$ từ bên trái:

     $$\langle\psi_n^{(0)}|\hat{H}_0 |\psi_n^{(1)}\rangle + \langle\psi_n^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle = E_n^{(0)} \langle\psi_n^{(0)}|\psi_n^{(1)}\rangle + E_n^{(1)} \langle\psi_n^{(0)}|\psi_n^{(0)}\rangle$$

     Sử dụng tính chất Hermitian của $\hat{H}_0$ và $\langle\psi_n^{(0)}|\hat{H}_0 = E_n^{(0)} \langle\psi_n^{(0)}|$, và điều kiện chuẩn hóa $\langle\psi_n^{(0)}|\psi_n^{(0)}\rangle = 1$, ta có:

     $$E_n^{(0)} \langle\psi_n^{(0)}|\psi_n^{(1)}\rangle + \langle\psi_n^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle = E_n^{(0)} \langle\psi_n^{(0)}|\psi_n^{(1)}\rangle + E_n^{(1)}$$

     Vậy:

     $$E_n^{(1)} = \langle\psi_n^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle$$

     Đây là giá trị trung bình của toán tử nhiễu loạn trong trạng thái không nhiễu loạn.

   • Hiệu chỉnh trạng thái bậc nhất: Khai triển $|\psi_n^{(1)}\rangle$ theo cơ sở các trạng thái riêng của hệ không nhiễu loạn:

     $$|\psi_n^{(1)}\rangle = \sum_{k \ne n} c_{kn} |\psi_k^{(0)}\rangle$$

     Lưu ý rằng chúng ta không bao gồm $k=n$ trong tổng vì thành phần đó có thể được điều chỉnh để đảm bảo chuẩn hóa.

     Thay vào phương trình bậc 1 và nhân với $\langle\psi_m^{(0)}|$ (với $m\ne n$) từ bên trái, ta được:

     $$\sum_{k \ne n} c_{kn} \langle\psi_m^{(0)}|\hat{H}_0 |\psi_k^{(0)}\rangle + \langle\psi_m^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle = \sum_{k \ne n} c_{kn} E_n^{(0)} \langle\psi_m^{(0)}|\psi_k^{(0)}\rangle + E_n^{(1)} \langle\psi_m^{(0)}|\psi_n^{(0)}\rangle$$

     Sử dụng tính trực giao $\langle\psi_m^{(0)}|\psi_k^{(0)}\rangle = \delta_{mk}$ và $\hat{H}_0 |\psi_k^{(0)}\rangle = E_k^{(0)} |\psi_k^{(0)}\rangle$, ta có:

     $$c_{mn} E_m^{(0)} + \langle\psi_m^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle = c_{mn} E_n^{(0)}$$

     Vậy:

     $$c_{mn} = \frac{\langle\psi_m^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle}{E_n^{(0)} - E_m^{(0)}} \quad (m \ne n)$$

     Do đó:

     $$|\psi_n^{(1)}\rangle = \sum_{m \ne n} \frac{\langle\psi_m^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle}{E_n^{(0)} - E_m^{(0)}} |\psi_m^{(0)}\rangle$$

   • Hiệu chỉnh năng lượng bậc hai: Nhân phương trình bậc 2 với $\langle\psi_n^{(0)}|$ từ bên trái và sử dụng các kết quả đã tìm được:

     $$E_n^{(2)} = \sum_{m \ne n} \frac{|\langle\psi_m^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle|^2}{E_n^{(0)} - E_m^{(0)}}$$

3.3. Tóm tắt kết quả

• Hiệu chỉnh năng lượng bậc nhất:

  $$E_n^{(1)} = \langle\psi_n^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle$$

• Hiệu chỉnh trạng thái bậc nhất:

  $$|\psi_n^{(1)}\rangle = \sum_{m \ne n} \frac{\langle\psi_m^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle}{E_n^{(0)} - E_m^{(0)}} |\psi_m^{(0)}\rangle$$

• Hiệu chỉnh năng lượng bậc hai:

  $$E_n^{(2)} = \sum_{m \ne n} \frac{|\langle\psi_m^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle|^2}{E_n^{(0)} - E_m^{(0)}}$$

3.4. Lưu ý quan trọng

• Điều kiện áp dụng: Kỹ thuật Nhiễu Loạn chỉ hiệu quả khi nhiễu loạn là nhỏ, tức là $|\lambda \hat{H}'| \ll |\hat{H}_0|$.
• Trạng thái suy biến: Khi có các trạng thái suy biến (có cùng năng lượng), cần có cách xử lý đặc biệt (xem mục 4).
• Sự hội tụ: Chuỗi nhiễu loạn không phải lúc nào cũng hội tụ.

4. Trường Hợp Trạng Thái Suy Biến

Khi các trạng thái của hệ không nhiễu loạn suy biến, tức là có nhiều trạng thái có cùng năng lượng, Kỹ thuật Nhiễu Loạn cần được điều chỉnh.

4.1. Bài toán

Giả sử có $g$ trạng thái $|\psi_{n,i}^{(0)}\rangle$ ($i = 1, 2, ..., g$) có cùng năng lượng $E_n^{(0)}$. Khi có nhiễu loạn, sự suy biến này có thể được phá vỡ, và chúng ta cần tìm các mức năng lượng mới và các trạng thái riêng tương ứng.

4.2. Phương pháp

1. Chọn cơ sở thích hợp: Thay vì sử dụng các trạng thái suy biến ban đầu $|\psi_{n,i}^{(0)}\rangle$, chúng ta cần chọn một cơ sở mới $\{|\phi_{n,i}\rangle\}$ trong không gian con suy biến sao cho toán tử nhiễu loạn $\hat{H}'$ là đường chéo trong cơ sở này. Tức là, chúng ta cần tìm các tổ hợp tuyến tính của $|\psi_{n,i}^{(0)}\rangle$:

   $$|\phi_{n,i}\rangle = \sum_{j=1}^g a_{ij} |\psi_{n,j}^{(0)}\rangle$$

   sao cho:

   $$\langle\phi_{n,i}|\hat{H}'|\phi_{n,j}\rangle = \delta_{ij} E_{n,i}^{(1)}$$

   trong đó $E_{n,i}^{(1)}$ là hiệu chỉnh năng lượng bậc nhất cho trạng thái $|\phi_{n,i}\rangle$.

2. Tìm các hệ số $a_{ij}$: Để tìm các hệ số $a_{ij}$, chúng ta xét phương trình bậc nhất trong chuỗi nhiễu loạn:

   $$\hat{H}_0 |\psi_n^{(1)}\rangle + \hat{H}' |\psi_n^{(0)}\rangle = E_n^{(0)} |\psi_n^{(1)}\rangle + E_n^{(1)} |\psi_n^{(0)}\rangle$$

   Nhân phương trình này với $\langle\psi_{n,k}^{(0)}|$ từ bên trái (với $k$ thuộc tập các trạng thái suy biến), ta được:

   $$\langle\psi_{n,k}^{(0)}|\hat{H}_0 |\psi_n^{(1)}\rangle + \langle\psi_{n,k}^{(0)}|\hat{H}' |\psi_n^{(0)}\rangle = E_n^{(0)} \langle\psi_{n,k}^{(0)}|\psi_n^{(1)}\rangle + E_n^{(1)} \langle\psi_{n,k}^{(0)}|\psi_n^{(0)}\rangle$$

   Do $|\psi_n^{(0)}\rangle$ là một tổ hợp tuyến tính của các trạng thái suy biến, chúng ta viết:

   $$|\psi_n^{(0)}\rangle = \sum_{j=1}^g a_{nj} |\psi_{n,j}^{(0)}\rangle$$

   Thay vào phương trình trên, ta có:

   $$\sum_{j=1}^g a_{nj} \langle\psi_{n,k}^{(0)}|\hat{H}' |\psi_{n,j}^{(0)}\rangle = E_n^{(1)} \sum_{j=1}^g a_{nj} \langle\psi_{n,k}^{(0)}|\psi_{n,j}^{(0)}\rangle = E_n^{(1)} a_{nk}$$

   Đặt $H'_{kj} = \langle\psi_{n,k}^{(0)}|\hat{H}' |\psi_{n,j}^{(0)}\rangle$, ta có hệ phương trình:

   $$\sum_{j=1}^g H'_{kj} a_{nj} = E_n^{(1)} a_{nk} \quad (k = 1, 2, ..., g)$$

   Đây là một bài toán trị riêng, chúng ta cần tìm các giá trị $E_n^{(1)}$ và các vector riêng tương ứng $a_n = (a_{n1}, a_{n2}, ..., a_{ng})$.

3. Giải bài toán trị riêng: Chúng ta viết hệ phương trình trên dưới dạng ma trận:

   $$\mathbf{H'} \mathbf{a} = E_n^{(1)} \mathbf{a}$$

   Trong đó $\mathbf{H'}$ là ma trận $g \times g$ với các phần tử $H'_{kj}$. Các giá trị $E_n^{(1)}$ là các trị riêng của ma trận $\mathbf{H'}$, và các vector $\mathbf{a}$ là các vector riêng tương ứng.

   Các trị riêng $E_{n,i}^{(1)}$ là các hiệu chỉnh năng lượng bậc nhất, và các vector riêng $\mathbf{a}_i$ cho chúng ta các hệ số $a_{ij}$ để xây dựng cơ sở mới $|\phi_{n,i}\rangle$.

4. Tính các hiệu chỉnh cao hơn (nếu cần): Sau khi tìm được cơ sở thích hợp và các hiệu chỉnh năng lượng bậc nhất, chúng ta có thể tính các hiệu chỉnh bậc cao hơn theo cách tương tự như trường hợp không suy biến, nhưng sử dụng cơ sở $|\phi_{n,i}\rangle$ thay vì $|\psi_{n,i}^{(0)}\rangle$.

4.3. Tóm tắt

Trong trường hợp suy biến, chúng ta cần:

1. Tìm ma trận nhiễu loạn $\mathbf{H'}$ trong không gian con suy biến.
2. Giải bài toán trị riêng để tìm các hiệu chỉnh năng lượng bậc nhất và các vector riêng.
3. Sử dụng các vector riêng để xây dựng cơ sở mới trong không gian con suy biến.
4. Tính các hiệu chỉnh cao hơn (nếu cần) trong cơ sở mới.

5. Ứng dụng

Kỹ thuật Nhiễu Loạn có rất nhiều ứng dụng trong Vật lý Lượng tử, bao gồm:

• Nguyên tử Hydro: Tính hiệu ứng Stark (sự dịch chuyển mức năng lượng do tác dụng của điện trường ngoài) và hiệu ứng Zeeman (sự dịch chuyển mức năng lượng do tác dụng của từ trường ngoài).
• Phân tử: Tính năng lượng liên kết và cấu trúc điện tử của phân tử.
• Vật lý chất rắn: Tính các tính chất điện và quang của chất rắn.
• Vật lý hạt nhân: Tính các mức năng lượng của hạt nhân.

6. Kết luận

Kỹ thuật Nhiễu Loạn là một công cụ mạnh mẽ để giải gần đúng các bài toán lượng tử phức tạp. Nó cho phép chúng ta hiểu sâu hơn về các hệ vật lý mà không có lời giải chính xác. Việc nắm vững lý thuyết và kỹ năng áp dụng Kỹ thuật Nhiễu Loạn là rất quan trọng trong Vật lý Lượng tử.