Tổng quan luận án

Luận án tiến sĩ "Tương tác tán sắc gần các vật thể vĩ mô" do nghiên cứu sinh Nguyễn Dũng Chỉnh thực hiện tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Thành phố Hồ Chí Minh, thuộc chuyên ngành Vật lý lý thuyết và vật lý toán (mã số chuyên ngành: 62440103), dưới sự hướng dẫn khoa học của PGS. Hồ Trung Dũng, bảo vệ năm 2023.

Tương tác tán sắc là các tương tác điện từ có bản chất lượng tử xuất hiện giữa các vật thể trung hòa điện thông qua trạng thái chân không của trường điện từ. Tùy thuộc vào bản chất vi mô hay vĩ mô của hệ, tương tác tán sắc được phân loại thành: tương tác van der Waals (vdW) giữa các nguyên tử hoặc phân tử; tương tác Casimir-Polder (CP) giữa nguyên tử và vật thể vĩ mô; và tương tác Casimir giữa các vật thể vĩ mô với nhau. Mặc dù là tương tác yếu nhất trong thang đo điện từ, tương tác tán sắc đóng vai trò nền tảng trong vật lý phân tử, vật lý bề mặt, sự bám dính sinh học, bẫy nguyên tử và công nghệ thông tin lượng tử.

Khoảng trống nghiên cứu được tác giả xác định xuất phát từ giới hạn của các công trình đi trước:

  • Cấu trúc trụ là hình học cơ bản của sợi quang truyền thông và các buồng vi cộng hưởng hiện đại, nhưng hàm Green trong hệ tọa độ trụ có cấu trúc toán học phức tạp hơn nhiều so với hệ phẳng hay hệ cầu, dẫn đến việc nghiên cứu còn chưa đầy đủ.
  • Các công trình trước đây về hệ trụ phần lớn chỉ giới hạn ở mô hình hai hoặc ba lớp, hoặc áp dụng mô hình nguồn dòng dài vô hạn hai chiều, chưa phản ánh đầy đủ mô hình điện động lực học lượng tử ba chiều.
  • Cấu trúc trụ phản xạ Bragg nhiều lớp (Bragg reflector) với khả năng giam nhốt quang học và tạo các mode cộng hưởng sắc nét chưa từng được khảo sát toàn diện đối với tương tác CP và tốc độ rã tự phát.
  • Trong bài toán tương tác vdW giữa hai nguyên tử, các nghiên cứu trước đây chỉ khảo sát vị trí nguyên tử đặt trên trục hoặc song song với trục, bỏ qua cấu hình hai nguyên tử đặt đối xứng trên cùng một đường kính gần biên. Đồng thời, việc tính toán số hàm Green trên trục tần số ảo gặp khó khăn lớn do sự phân kỳ tại cận vô cùng.
  • Khi chuyển sang hệ trụ có chiều dài hữu hạn, các nghiên cứu gần đúng Born trước đó chỉ khảo sát moment lưỡng cực định hướng dọc theo trục ($z$), chưa phân tích đầy đủ các phương phân cực khác và chưa tính toán thế CP cho nguyên tử ở trạng thái cơ bản cũng như kích thích gần đầu mút khối trụ.

Luận án đặt ra các mục tiêu và nhiệm vụ nghiên cứu cụ thể, được tổ chức thành ba bài toán lớn:

  1. Bài toán 1: Khảo sát tốc độ rã tự phát và tương tác Casimir-Polder của một nguyên tử hai mức kích thích đặt ở lớp trong cùng của hệ trụ dài vô hạn có tường phản xạ Bragg nhiều lớp; xác định cấu trúc mode cộng hưởng và so sánh định lượng với mô hình tường phản xạ toàn phần (Perfectly Reflecting Wall - PRW).
  2. Bài toán 2: Nghiên cứu tương tác van der Waals giữa hai nguyên tử ở trạng thái cơ bản bên trong hệ trụ dài vô hạn cho hai cấu hình không gian: hai nguyên tử trên trục của hệ trụ Bragg và hai nguyên tử đặt đối diện nhau trên cùng một đường kính gần biên của hệ trụ hai lớp; phân tích ảnh hưởng của vật liệu thực (Silicon, Vàng) và phát triển kỹ thuật tính số trên trục tần số ảo.
  3. Bài toán 3: Khảo sát ảnh hưởng của biên hai đầu khối trụ lên tốc độ rã tự phát và tương tác Casimir-Polder đối với một nguyên tử đặt gần một đầu của hệ trụ có chiều dài hữu hạn, mở rộng cho mọi hướng phân cực của moment lưỡng cực điện và đánh giá khả năng tạo bẫy nguyên tử ba chiều.

Phạm vi và đối tượng nghiên cứu của luận án:

  • Đối tượng nghiên cứu: Nguyên tử hai mức (ở trạng thái kích thích hoặc trạng thái cơ bản), cặp nguyên tử trung hòa điện không phân cực vĩnh cửu, trường điện từ lượng tử hóa, và các cấu trúc vật thể vĩ mô dạng trụ.
  • Phạm vi không gian và hình học: Hệ trụ đa lớp dài vô hạn có cấu trúc gương phản xạ Bragg (số lớp $N$ từ 2 đến 25), hệ trụ phản xạ lý tưởng hai lớp, và hệ trụ điện môi có chiều dài hữu hạn ($H/\lambda_A = 10, R/\lambda_A = 1$); vị trí nguyên tử được đặt trên trục đối xứng, lệch trục, đối diện trên đường kính hoặc đặt gần mặt phẳng đầu mút khối trụ ($z_A/\lambda_A$ từ $4.96$ đến $5.0$).
  • Phạm vi vật liệu: Môi trường chân không, vật liệu phản xạ toàn phần lý tưởng ($\epsilon \to -\infty$), vật liệu điện môi lý tưởng không suy hao ($\epsilon_L = 1, \epsilon_H = 9$ hoặc $25$), và các vật liệu thực có tán sắc và hấp thụ gồm Silicon ($\mathrm{Si}$), Vàng ($\mathrm{Au}$), Silica ($\mathrm{SiO_2}$).
  • Phạm vi vật lý: Quá trình diễn ra tại nhiệt độ không tuyệt đối ($T = 0$), tương tác được mô tả thông qua gần đúng lưỡng cực điện và gần đúng sóng quay (đối với tốc độ rã tự phát).

Tổng quan tài liệu và vị trí của luận án

Tương tác tán sắc bắt đầu từ các khảo sát cơ học lượng tử của London (1930) [5] cho tương tác vdW giữa hai nguyên tử trung hòa điện, không phân cực trong không gian tự do ở khoảng cách nhỏ so với bước sóng chuyển mức nguyên tử, dẫn xuất quy luật lực hút tỷ lệ nghịch với lũy thừa bậc sáu của khoảng cách ($r^{-6}$). Casimir và Polder (1948) [6] phát triển lý thuyết lượng tử tổng quát tính đến hiệu ứng trễ do tốc độ ánh sáng hữu hạn, chỉ ra rằng ở khoảng cách lớn hơn nhiều bước sóng chuyển mức, thế tương tác chuyển sang quy luật tiệm cận $r^{-7}$, đồng thời công bố biểu thức cho tương tác giữa một nguyên tử và một tấm dẫn lý tưởng (tương tác CP). Lý thuyết vdW sau đó được mở rộng cho hai nguyên tử ở trạng thái kích thích (công trình [7-9]), hệ một nguyên tử kích thích và một nguyên tử ở mức cơ bản ([7, 10]), nguyên tử phân cực từ ([11, 12]), hệ có mặt photon nhiệt ([13]), và hai nguyên tử trong điện trường tĩnh ([14]).

Về ảnh hưởng của vật thể vĩ mô, các cấu hình phẳng ([21-23, 26]), cầu ([24-26]) và trụ ([27-30]) lần lượt được nghiên cứu. Điển hình, nghiên cứu [21] chỉ ra tương tác vdW giữa hai nguyên tử ở nhiệt độ khác không trong buồng cộng hưởng phẳng độc lập với hiệu ứng trễ. Đối với hệ trụ, công trình [27] chứng minh tương tác vdW giữa hai nguyên tử đặt giữa hai tường phản xạ toàn phần của hệ trụ ba lớp tăng đáng kể so với không gian tự do. Công trình [28] khảo sát hai nguyên tử trên trục của hệ trụ hai lớp phản xạ toàn phần, chỉ ra tương tác giảm theo hàm mũ theo khoảng cách dọc trục nhưng tăng nhiều bậc khi moment lưỡng cực định hướng dọc trục $z$ trong hệ trụ silicon. Khi hai nguyên tử lệch trục, công trình [29] phát hiện lực tăng nhiều bậc đối với lưỡng cực định hướng theo phương bán kính ở bên ngoài hệ trụ hai lớp. Đặc biệt, nghiên cứu [30] chứng minh tương tác vdW biến mất trong chân không khi hai moment lưỡng cực vuông góc nhau nhưng lại tiếp tục tồn tại khi có mặt hệ trụ hai lớp.

Trong lĩnh vực điện động lực học lượng tử buồng cộng hưởng (cavity QED), tốc độ rã tự phát và dịch chuyển Lamb (tương đương thế CP) được khảo sát bởi nhiều nhóm tác giả: Wylie và Sipe (1984) [51] dùng lý thuyết phản hồi tuyến tính cho nguyên tử gần tấm phẳng; Rippin và Knight (1996) [52] dùng phương pháp khai triển mode cho nguyên tử hai mức trong hệ trụ hai lớp vô hạn và gương Bragg dọc trục; Henkel và Sandoghdar (1998) [53] phân tích ảnh hưởng của độ nhám bề mặt điện môi; Klimov và Ducloy (2002, 2004) [54, 55] nghiên cứu dịch chuyển lưỡng cực và tứ cực gần hệ trụ kim loại và điện môi; Guzatov và Klimov (2006) [56] khảo sát trụ dẫn điện tiết diện elip; Bennett (2009) [57] dùng hàm Green và khai triển Born cho bề mặt quang khắc; Kien và Rauschenbeutel (2014) [58] khảo sát lưỡng cực quay gần bề mặt phẳng; Buhmann và Scheel (2008) [59] thảo luận điều kiện áp dụng lý thuyết Lifshitz cho lực CP nhiệt; công trình [60] dùng khai triển Born xét biên tấm phẳng; Ellingsen, Buhmann và Scheel (2010) [61] nghiên cứu buồng cộng hưởng trụ kim loại ở nhiệt độ khác không; Bender và cộng sự (2014) [62] đo đạc thực nghiệm tương tác CP; công trình [63] nghiên cứu tốc độ rã bên trong trụ hai lớp hữu hạn và loại bỏ kỳ dị ở tâm; Stourm và cộng sự (2019) [64] tính tốc độ rã và dịch chuyển Lamb của nguyên tử $^{87}\mathrm{Rb}$ gần sợi quang $\mathrm{SiO_2}$.

Đối với hệ trụ đa lớp, các nghiên cứu trước như [65, 66] xem xét hệ $N$ lớp nhưng giới hạn ở mô hình hai chiều do giả định nguồn là dòng dài vô tận; các công trình [67, 68] tập trung vào tiết diện tán xạ để ứng dụng tàng hình; công trình [69] khảo sát Tamm plasmon polaritons. Đối với hệ trụ hữu hạn, các mô hình thực nghiệm như nguồn phát đơn photon [71-73] và bẫy ion tích hợp sợi quang [74] đòi hỏi tính toán lý thuyết chính xác. Gần đúng Born từng được áp dụng để tính tốc độ rã của nguyên tử bên ngoài khối trụ hữu hạn [78] và bên trong khối trụ [63], nhưng cả hai đều chỉ giới hạn ở trường hợp moment lưỡng cực định hướng dọc theo trục $z$.

Khoảng trống mà luận án lựa chọn giải quyết bao gồm:

  1. Thiết lập bài toán ba chiều giải tích và giải số đầy đủ cho tốc độ rã tự phát và tương tác CP bên trong hệ trụ phản xạ Bragg $N$ lớp, không dùng mô hình nguồn dòng hai chiều.
  2. Khảo sát cấu hình hai nguyên tử đặt đối xứng trên cùng một đường kính trong ống dẫn sóng trụ đa lớp, đồng thời đề xuất giải pháp toán học xử lý tích phân hàm Green tán xạ trên trục tần số ảo.
  3. Mở rộng lý thuyết khai triển Born cho hệ trụ có chiều dài hữu hạn sang mọi phương định hướng lưỡng cực (phương $z$, phương bán kính $\rho$, phương góc $\phi$ và trung bình định hướng ngẫu nhiên), phân tích thế CP ba chiều ở vùng lân cận đầu mút khối trụ.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án xây dựng trên nền tảng điện động lực học lượng tử vĩ mô (macroscopic QED) trong môi trường vật chất có tán sắc và hấp thụ thỏa mãn hệ thức Kramers-Kronig, theo khuôn khổ lý thuyết của Dung, Ho Trung Dung và các cộng sự ([86, 88]).

Trường điện từ lượng tử được biểu diễn thông qua các toán tử vector thế và cường độ điện trường. Để mô tả tính hấp thụ của môi trường, các toán tử mật độ điện tích nhiễu $\rho(\mathbf{r}, \omega)$ và mật độ dòng nhiễu $\mathbf{j}(\mathbf{r}, \omega)$ được đưa vào phương trình Maxwell lượng tử, liên hệ với toán tử phân cực nhiễu $\mathbf{P}(\mathbf{r}, \omega)$ qua hệ thức $\rho(\mathbf{r}, \omega) = -\nabla \cdot \mathbf{P}(\mathbf{r}, \omega)$. Toán tử cường độ điện trường $\hat{\mathbf{E}}(\mathbf{r}, \omega)$ được biểu diễn qua tích phân không gian của các trường toán tử sinh/hủy hạt boson liên tục $\mathbf{f}(\mathbf{r}, \omega)$ và $\mathbf{f}^\dagger(\mathbf{r}, \omega)$, với nhân tích phân là hàm tensor Green $\boldsymbol{G}(\mathbf{r}, \mathbf{r}', \omega)$:

$$\hat{\mathbf{E}}(\mathbf{r}, \omega) = i \sqrt{\frac{\hbar}{\pi \epsilon_0}} \frac{\omega^2}{c^2} \int d^3\mathbf{r}' , \sqrt{\text{Im},\epsilon(\mathbf{r}', \omega)} , \boldsymbol{G}(\mathbf{r}, \mathbf{r}', \omega) \cdot \mathbf{f}(\mathbf{r}', \omega)$$

Toán tử Hamilton của hệ trường điện từ kết hợp môi trường hấp thụ và hạt tích điện được xây dựng dưới hai sơ đồ tương đương:

  • Sơ đồ tương tác tối thiểu (minimal-coupling scheme): Được áp dụng trong Chương 1 và Chương 2 kết hợp lý thuyết Wigner-Weisskopf để xác định tốc độ rã tự phát $\Gamma$ và độ dịch chuyển mức năng lượng của nguyên tử hai mức. Sau khi áp dụng gần đúng lưỡng cực điện và gần đúng sóng quay (rotating wave approximation), tốc độ rã tự phát $\Gamma$ và thế Casimir-Polder $U$ liên hệ trực tiếp với phần ảo và phần thực của hàm Green tán xạ tại vị trí nguyên tử $\mathbf{r}_A$:

$$\Gamma = \frac{2\omega_A^2}{\hbar \epsilon_0 c^2} \mathbf{d}_A \cdot \text{Im},\boldsymbol{G}(\mathbf{r}_A, \mathbf{r}_A, \omega_A) \cdot \mathbf{d}_A$$

$$U = -\frac{\omega_A^2}{\epsilon_0 c^2} \mathbf{d}_A \cdot \text{Re},\boldsymbol{G}^{(sc)}(\mathbf{r}_A, \mathbf{r}_A, \omega_A) \cdot \mathbf{d}_A$$

  • Sơ đồ tương tác đa cực (multipolar-coupling scheme): Được áp dụng trong Chương 1 và Chương 3 dựa trên biến đổi Power-Zienau-Woolley theo tiếp cận của Safari và cộng sự (2006) [23]. Sử dụng lý thuyết nhiễu loạn lượng tử bậc bốn đối với toán tử tương tác $\hat{H}_{int} = -\mathbf{d}_A \cdot \hat{\mathbf{E}}(\mathbf{r}_A) - \mathbf{d}_B \cdot \hat{\mathbf{E}}(\mathbf{r}_B)$, kết hợp với kỹ thuật tích phân vòng Cauchy trên nửa trên mặt phẳng phức tần số để chuyển tích phân tần số thực sang tần số ảo $\omega = iu$, thu được biểu thức tổng quát cho thế vdW giữa hai nguyên tử mức cơ bản:

$$U_{AB}(\mathbf{r}_A, \mathbf{r}_B) = -\frac{\hbar \mu_0^2}{\pi} \int_0^\infty du , u^4 \alpha_A(iu) \alpha_B(iu) , \text{Tr}\left[ \boldsymbol{G}(\mathbf{r}_A, \mathbf{r}_B, iu) \cdot \boldsymbol{G}(\mathbf{r}_B, \mathbf{r}_A, iu) \right]$$

Các phương pháp tính toán và xử lý số liệu thực tế bao gồm:

  1. Phương pháp giải tích hàm tensor Green trong hệ trụ: Tác giả sử dụng dạng tường minh hàm Green của Li và cộng sự (2000) [82] cho hệ trụ đa lớp vô hạn. Trường điện từ được phân tích theo hệ hàm vector sóng trụ $\mathbf{M}{e/o,n\eta}(h)$ và $\mathbf{N}{e/o,n\eta}(h)$. Các hệ số tán xạ $C_j^{(H,V)s}$ cho các mode sóng điện ngang TE ($H$) và từ ngang TM ($V$) được xác định thông qua ma trận truyền $2 \times 2$ giữa các mặt phân cách trụ liên tiếp.
  2. Kỹ thuật giải số trên trục tần số ảo: Trong bài toán tương tác vdW tại tần số ảo $\omega = iu$, hàm dưới tích phân có cận trên phân kỳ khi $u \to \infty$. Tác giả xử lý bằng cách bóc tách phần tiệm cận giải tích của hàm Green chân không và hàm Green tán xạ tương ứng với phản xạ cục bộ, tính tích phân giải tích phần tiệm cận và chỉ áp dụng phương pháp số cho phần hiệu số hội tụ nhanh.
  3. Phương pháp khai triển Born bậc nhất (First-order Born approximation): Áp dụng cho hệ trụ có chiều dài hữu hạn ở Chương 4, nơi điều kiện biên không cho phép tìm nghiệm giải tích chính xác của hàm Green. Khai triển chuỗi Born đối với độ cảm điện của khối trụ hữu hạn và ngắt chuỗi ở số hạng bậc nhất giúp đưa bài toán về dạng tích phân khối trên thể tích của vật thể hữu hạn.

Nội dung chính theo từng chương

Chương 1: Lý thuyết tổng quan

Chương 1 hệ thống hóa cơ sở lý thuyết về lượng tử hóa trường điện từ và tương tác tán sắc trong môi trường vật chất. Mở đầu bằng mô tả định tính về tương tác Casimir giữa hai tấm phẳng kim loại lý tưởng thông qua sự gián đoạn của mật độ trạng thái photon ảo trong buồng cộng hưởng, và tương tác van der Waals giữa hai nguyên tử trung hòa điện do sự thăng giáng lưỡng cực tức thời.

Tác giả trình bày chi tiết quy trình lượng tử hóa trường điện từ trong môi trường có tán sắc và hấp thụ tổng quát bằng phương pháp hàm tensor Green. Chương này dẫn xuất đầy đủ các biểu thức nền tảng:

  • Biểu thức tốc độ rã tự phát $\Gamma$ và độ dịch chuyển mức năng lượng của nguyên tử hai mức tương tác với trường theo lý thuyết Wigner-Weisskopf.
  • Biểu thức tương tác Casimir-Polder cho nguyên tử nhiều mức ở trạng thái kích thích hoặc cơ bản bằng lý thuyết nhiễu loạn bậc hai, phân tách tường minh thế CP thành phần cộng hưởng $U^{res}$ và phần phi cộng hưởng $U^{nres}$.
  • Biểu thức thế tương tác van der Waals giữa hai nguyên tử ở trạng thái cơ bản thông qua lý thuyết nhiễu loạn bậc bốn trong sơ đồ tương tác đa cực, chứng minh sự chuyển đổi tích phân từ trục tần số thực sang trục tần số ảo $\omega = iu$ nhờ định lý tích phân Cauchy.

Chương 2: Nguyên tử kích thích bên trong hệ trụ dài vô hạn

Chương 2 giải quyết bài toán tính toán tốc độ rã tự phát $\Gamma/\Gamma_0$ và thế Casimir-Polder $U$ của một nguyên tử hai mức kích thích đặt tại lớp trong cùng của một hệ trụ nhiều lớp có chiều dài vô hạn. Hệ trụ được thiết kế thỏa mãn điều kiện của gương phản xạ Bragg hình trụ:

$$R_m - R_{m-1} = \frac{\lambda_0}{4 \text{Re}\sqrt{\epsilon_{m+1}(\omega)}}$$

trong đó $\lambda_0$ là bước sóng trung tâm của vùng cấm quang học (band gap), $\epsilon_L = 1$ và $\epsilon_H = 9$ (hoặc thay đổi) là hằng số điện môi của các lớp kế tiếp nhau.

Tác giả phân tích hàm Green tán xạ thành ba thành phần tương ứng với các hướng phân cực của moment lưỡng cực nguyên tử:

  • Khi lưỡng cực định hướng dọc theo trục $z$, chỉ có các mode TM đóng góp vào hàm Green tán xạ.
  • Khi lưỡng cực định hướng theo phương bán kính $\rho$ hoặc phương góc $\phi$, cả mode TE và TM đều tham gia đóng góp.
Thông số khảo sát Cấu hình hệ trụ Bragg Cấu hình phản xạ toàn phần (PRW)
Số lớp ($N$) $N = 5, 15, 25$ $N = 2$ ($\epsilon_1 \to -\infty$)
Hằng số điện môi $\epsilon_L = 1, \epsilon_H = 9$ (hoặc thay đổi) $\epsilon = 1$ bên trong, thành dẫn lý tưởng
Điều kiện lớp Bề dày quang học mỗi lớp là $\lambda_0/4$ Thành phản xạ tại bán kính $R$
Mode đóng góp khi lưỡng cực hướng $z$ Chỉ có mode TM ($n = 0$ khi trên trục) Chỉ có mode TM ($p_{nm}$)
Mode đóng góp khi lưỡng cực hướng $\rho, \phi$ Cả mode TM và TE ($n = 1$ khi trên trục) Cả mode TM ($p_{nm}$) và TE ($q_{nm}$)
Vị trí mode cộng hưởng đầu tiên ($p_{01}$) $R = 0.3822591759\lambda_0$ Phân kỳ vô hạn tại nghiệm $J_0(p_{01}) = 0$

Các phát hiện chính của Chương 2:

  • Sự hình thành vùng cấm (Band gap): Khi nguyên tử đặt trên trục có lưỡng cực hướng theo trục $z$, vùng cấm bức xạ ($\Gamma/\Gamma_0 \ll 1$) xuất hiện rõ rệt ngay từ $N = 5$ lớp và biên dạng vùng cấm đạt trạng thái ổn định khi tăng lên $N = 15$ và $N = 25$. Khi tăng giá trị hằng số điện môi $\epsilon_H$, vùng cấm mở rộng cả về bề rộng tần số và độ suy giảm sâu của tốc độ rã.
  • Cấu trúc mode cộng hưởng: Khi cố định tần số nguyên tử $\omega_A = \omega_0$ và thay đổi bán kính lớp trong cùng $R$, các mode cộng hưởng xuất hiện với số lượng tăng dần khi bán kính tăng. Mode cộng hưởng đầu tiên $p_{01}$ xuất hiện tại giá trị bán kính chuẩn hóa $R = 0.3822591759\lambda_0$, tương ứng bề dày quang học của lớp trong cùng lớn hơn $\lambda_0/2$. Vị trí các đỉnh cộng hưởng của hệ Bragg hoàn toàn tương đồng với các vị trí phân kỳ của hệ phản xạ toàn phần (PRW), được xác định bởi điều kiện mẫu số hàm Green bằng không ($k_V = 0$ và $k_H = 0$).
  • Thế giam nhốt Casimir-Polder: Tại các tần số ứng với mode cộng hưởng, thế CP tạo thành giếng thế năng giam nhốt nguyên tử xung quanh trục của khối trụ. Luận án cũng mở rộng tính toán với vật liệu thực có tán sắc là Silicon ($\epsilon_{\mathrm{Si}} = 12.25$) và $\mathrm{SiO_2}$ ($\epsilon_{\mathrm{SiO_2}} = 2.25$).

Chương 3: Tương tác van der Waals giữa hai nguyên tử ở mức cơ bản bên trong hệ trụ dài vô hạn

Chương 3 khảo sát tương tác van der Waals giữa hai nguyên tử ở mức cơ bản bên trong hệ trụ dài vô hạn, sử dụng biểu thức tích phân hàm Green trên trục tần số ảo $\omega = iu$.

Nội dung tập trung vào hai cấu hình hình học cụ thể:

  1. Cấu hình 1: Hai nguyên tử đặt dọc theo trục đối xứng $z$ của hệ trụ Bragg ($N = 9, \epsilon_L = 1, \epsilon_H = \epsilon$) và hệ trụ hai lớp. Khảo sát được thực hiện với các giá trị $\epsilon = 9$ và $\epsilon = 25$, so sánh trực tiếp với trường hợp hệ phản xạ lý tưởng ($\epsilon \to -\infty$). Kết quả chỉ ra rằng thế vdW giảm theo quy luật hàm mũ theo khoảng cách giữa hai nguyên tử dọc theo trục.
  2. Cấu hình 2: Hai nguyên tử có bước sóng chuyển mức $\lambda_A = 780,\text{nm}$ ($\Omega/c = 2\pi/780,\text{nm}$) đặt đối diện nhau trên cùng một đường kính bên trong hệ trụ hai lớp, cách đều tâm và tiến sát biên bề mặt với các khoảng cách $\Delta$ khác nhau.

Các phát hiện chính của Chương 3:

  • Tác giả phát hiện sự tăng cường và ức chế mạnh mẽ của lực tán sắc vdW. Trong cấu hình hai nguyên tử đặt đối diện trên đường kính gần biên, độ lớn của lực tương tác vdW có thể tăng lên đến bốn bậc độ lớn ($10^4$ lần) so với tương tác trong không gian tự do.
  • Ảnh hưởng của vật liệu thực: Khảo sát đối với thành trụ bằng Vàng ($\mathrm{Au}$) và Silicon ($\mathrm{Si}$) cho thấy sự sai khác đáng kể so với mô hình phản xạ toàn phần. Khi bề rộng vạch hấp thụ (linewidth) của vật liệu thực tăng lên, biên độ tương tác vdW giảm do tiêu tán năng lượng trong môi trường. Đối với hệ trụ Bragg, khi hằng số điện môi tĩnh $\epsilon_0$ của Silicon tăng từ $11.87$ lên $16$ và $25$, thế tương tác vdW tiệm cận dần về kết quả của hệ phản xạ lý tưởng.

Chương 4: Tương tác Casimir-Polder của nguyên tử gần một đầu của hệ trụ dài hữu hạn

Chương 4 khảo sát ảnh hưởng của tính hữu hạn về chiều dài của khối trụ lên tốc độ rã tự phát và tương tác Casimir-Polder, giải quyết giới hạn của các mô hình trụ dài vô hạn. Hệ vật lý được khảo sát có bán kính chuẩn hóa $R/\lambda_A = 1$, chiều dài hữu hạn $H/\lambda_A = 10$, và hằng số điện môi $\epsilon(\omega_A) = 1$.

Do hàm Green cho hệ trụ hữu hạn không có nghiệm giải tích dạng đóng, tác giả áp dụng phương pháp khai triển Born bậc nhất. Luận án mở rộng việc thiết lập công thức khai triển Born cho cả hai trường hợp:

  • Nguyên tử đặt bên ngoài khối trụ hữu hạn.
  • Nguyên tử đặt bên trong khối trụ hữu hạn.

Các kết quả khảo sát chi tiết:

  • Tốc độ rã tự phát: Được khảo sát tại vị trí rất gần một đầu của khối trụ ($z_A/\lambda_A = 4.96$, trong khi đầu trụ tại $z = H/2\lambda_A = 5.0$). So sánh giữa khối trụ hữu hạn và khối trụ vô hạn cùng bán kính cho thấy hiệu ứng biên làm biến đổi đáng kể tốc độ rã tự phát khi nguyên tử dịch chuyển qua mặt phân cách bán kính $\rho_A/\lambda_A = 1$.
  • Thế Casimir-Polder lên nguyên tử kích thích: Thế chuẩn hóa $\bar{U} = U \cdot 8\pi\epsilon_0 / (\mu_0 p^2 \omega_A^3)$ được tính toán theo tọa độ bán kính $\rho_A$ và tọa độ trục $z_A$. Kết quả chỉ ra sự biến thiên sâu sắc của thế năng tại vùng không gian chuyển tiếp giữa bên trong và bên ngoài đầu mút khối trụ.
  • Thế Casimir-Polder lên nguyên tử ở mức cơ bản: Luận án phân tích chi tiết thế CP của nguyên tử mức cơ bản trên trục $z$ trong vùng lân cận đầu mút khối trụ. Tác giả khảo sát độc lập các trường hợp moment lưỡng cực định hướng vuông góc với bề mặt đầu trụ (phương $z$), song song với bề mặt (phương $\rho$ và $\phi$), và giá trị trung bình cho định hướng ngẫu nhiên.
  • Sự hình thành bẫy nguyên tử ba chiều: Kết quả tính toán thế CP chuẩn hóa theo các tọa độ không gian ($x_A, y_A, z_A$) chỉ ra rằng tại vùng không gian bên ngoài nằm ngay trên trục đối xứng gần đầu mút của khối trụ, xuất hiện một cấu trúc giếng thế thế năng Casimir-Polder ba chiều. Cấu trúc thế này có khả năng bẫy và giam giữ nguyên tử trong không gian mà không cần trường ngoài bổ trợ.

Kết quả và những đóng góp mới

Luận án tự xác định các đóng góp mới về mặt khoa học lý thuyết và giá trị ứng dụng thực tiễn:

  1. Về mặt lý luận và phương pháp học thuật:

    • Xây dựng hoàn chỉnh mô hình lý thuyết ba chiều mô tả tốc độ rã tự phát và tương tác Casimir-Polder của nguyên tử hai mức bên trong hệ trụ phản xạ Bragg đa lớp dựa trên kỹ thuật hàm tensor Green, khắc phục hoàn toàn hạn chế của các mô hình dòng hai chiều trước đó.
    • Phát triển giải pháp kết hợp tính toán giải tích và giải số để xử lý hàm Green tán xạ trên trục tần số ảo $\omega = iu$, giải quyết triệt để sự phân kỳ toán học tại cận vô cùng của tích phân, tạo ra công cụ tính toán tin cậy cho bài toán tương tác van der Waals giữa các nguyên tử trong ống dẫn sóng.
    • Mở rộng phương pháp khai triển Born bậc nhất cho hệ trụ có chiều dài hữu hạn để tính toán cho mọi hướng phân cực của moment lưỡng cực nguyên tử (phương trục $z$, phương bán kính $\rho$, phương góc $\phi$ và trung bình đẳng hướng), vượt qua giới hạn của các công trình tiền nhiệm vốn chỉ khảo sát được phương dọc trục.
  2. Về mặt thực tiễn và phát hiện vật lý:

    • Xác lập định lượng điều kiện hình thành và độ ổn định của vùng cấm quang học (band gap) đối với tốc độ rã tự phát trong hệ trụ Bragg ($N \ge 5$). Chứng minh sự tương đồng về vị trí mode cộng hưởng giữa hệ trụ Bragg và hệ phản xạ toàn phần (PRW), xác định mode cộng hưởng cơ bản $p_{01}$ xuất hiện tại bán kính $R = 0.3822591759\lambda_0$.
    • Chứng minh tương tác van der Waals giữa hai nguyên tử đặt đối diện trên đường kính gần biên trong hệ trụ hai lớp có thể được tăng cường vượt bậc lên tới bốn bậc độ lớn ($10^4$ lần) so với không gian tự do, hoặc bị ức chế tùy thuộc vào cấu hình hình học và hướng phân cực.
    • Phát hiện sự tồn tại của hệ thế năng Casimir-Polder ba chiều xuất hiện tại vùng không gian trên trục bên ngoài một đầu của khối trụ hữu hạn, cung cấp luận cứ vật lý cho việc thiết kế bẫy nguyên tử quang học.
  3. Kiến nghị và đề xuất giải pháp:

    • Luận án khuyến nghị sử dụng mô hình hệ trụ phản xạ toàn phần (PRW) như một công cụ giải tích thuận tiện để định vị nhanh các mode cộng hưởng của cấu trúc Bragg đa lớp, do tính chất trùng khớp của vị trí các đỉnh phổ.
    • Đề xuất ứng dụng cấu trúc trụ Bragg đa lớp trong việc kiểm soát bức xạ tự phát để chế tạo nguồn phát đơn photon (single-photon sources) và ứng dụng cấu trúc đầu mút sợi quang nano trong kỹ thuật bẫy nguyên tử trung hòa cho công nghệ thông tin lượng tử.

Hạn chế và hướng nghiên cứu tiếp

Các hạn chế được văn bản luận án thừa nhận:

  • Các mô hình hình học trong luận án chủ yếu tập trung vào hệ trụ tròn đối xứng lý tưởng, chưa bao quát các dạng tiết diện phức tạp khác hoặc bề mặt có độ nhám thực tế.
  • Đối với hệ trụ có chiều dài hữu hạn ở Chương 4, phương pháp khai triển Born mới chỉ được áp dụng ở phép gần đúng bậc nhất (first-order Born approximation) và ngắt chuỗi số hạng đầu; điều này giới hạn độ chính xác khi độ chênh lệch hằng số điện môi của khối trụ so với môi trường xung quanh trở nên lớn.
  • Mặc dù khung lý thuyết tổng quát ở Chương 1 thiết lập cho môi trường có cả tính chất từ ($\mu \neq 1$), toàn bộ các tính toán giải số và ứng dụng cụ thể trong hệ trụ ở các chương sau đều bỏ qua tính chất từ ($\mu = 1$).
  • Các tính toán số chi tiết đều giả định nhiệt độ bằng không tuyệt đối ($T = 0$), chưa tính đến sự đóng góp của các photon nhiệt ở nhiệt độ phòng.

Các hướng nghiên cứu tiếp theo được tác giả đề xuất:

  • Mở rộng phương pháp khai triển Born lên các số hạng bậc cao hơn nhằm nâng cao độ chính xác khi mô tả các hệ trụ hữu hạn có hằng số điện môi lớn.
  • Phát triển bài toán tương tác tán sắc và bức xạ tự phát trong hệ trụ ở điều kiện nhiệt độ hữu hạn ($T > 0$) để đánh giá đầy đủ tác động của hiệu ứng nhiệt Lifshitz.
  • Nghiên cứu tương tác tán sắc trong các môi trường vật liệu dị hướng, vật liệu biến hóa (metamaterials) hoặc vật liệu thuận tay trái (left-handed materials) có độ từ thẩm phức tạp.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu có giá trị đối với:

  • Nghiên cứu sinh và học viên cao học ngành Vật lý: Cung cấp tài liệu hoàn chỉnh về kỹ thuật lượng tử hóa trường điện từ trong môi trường điện môi vĩ mô có tán sắc và hấp thụ, cùng phương pháp toán học xử lý hàm tensor Green trong hệ tọa độ trụ.
  • Giảng viên và nhà nghiên cứu Quang học lượng tử, Cavity QED: Cung cấp cơ sở lý thuyết chuẩn xác về mối liên hệ giữa hàm Green tán xạ với tốc độ rã tự phát, độ dịch chuyển Lamb và thế tương tác tán sắc.
  • Các nhóm nghiên cứu công nghệ nano và quang lượng tử thực nghiệm: Cung cấp các thông số định lượng (vị trí mode cộng hưởng, kích thước bước sóng, khoảng cách biên) giúp định hướng thiết kế buồng cộng hưởng sợi quang, linh kiện quang tử nano, nguồn phát photon đơn và cấu hình bẫy nguyên tử gần đầu sợi quang.

Các phần nội dung đáng tham khảo nhất trong văn bản:

  • Mục 1.2 đến 1.5: Hệ thống công thức giải tích chuẩn xác từ toán tử Hamilton đến các đại lượng vật lý quan sát được.
  • Mục 2.2: Kỹ thuật xây dựng ma trận truyền $2 \times 2$ để xác định các hệ số tán xạ của sóng TE và TM trong cấu trúc trụ đa lớp.
  • Mục 3.2: Quy trình toán học kết hợp giải tích và số học để xử lý tích phân hàm Green trên trục tần số ảo.
  • Mục 4.1 và 4.2: Phương pháp triển khai gần đúng Born bậc nhất cho cấu trúc trụ hữu hạn áp dụng cho các hướng lưỡng cực khác nhau.

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

1. Vùng cấm quang học (band gap) đối với tốc độ rã tự phát trong hệ trụ Bragg xuất hiện từ số lớp bao nhiêu và biến đổi như thế nào khi thay đổi thông số cấu trúc? Theo kết quả giải số tại Hình 2.2 của luận án, vùng cấm quang học (nơi tốc độ rã tự phát chuẩn hóa $\Gamma/\Gamma_0 \ll 1$) xuất hiện rõ rệt ngay từ khi hệ có $N = 5$ lớp. Hình dáng của vùng cấm trở nên ổn định khi tăng số lớp lên $N = 15$ và $N = 25$. Khi giữ nguyên số lớp $N = 15$ và tăng giá trị hằng số điện môi cao $\epsilon_H$ (với $\epsilon_L = 1$), vùng cấm mở rộng ra về dải tần số đồng thời độ suy giảm của tốc độ rã tự phát trở nên sâu hơn.

2. Mode cộng hưởng đầu tiên trong hệ trụ Bragg xuất hiện tại giá trị bán kính nào và có đặc điểm gì so với hệ phản xạ toàn phần? Khi khảo sát nguyên tử trên trục có tần số chuyển mức trùng với tần số trung tâm $\omega_A = \omega_0$ trong hệ trụ 15 lớp ($\epsilon_L = 1, \epsilon_H = 9$), mode cộng hưởng cơ bản đầu tiên $p_{01}$ xuất hiện tại bán kính lớp trong cùng $R = 0.3822591759\lambda_0$. Giá trị này ứng với bề dày quang học của lớp trong cùng lớn hơn $\lambda_0/2$. Vị trí của đỉnh cộng hưởng này gần như trùng khớp hoàn toàn với vị trí phân kỳ của hệ phản xạ toàn phần (PRW), cho phép dùng hệ PRW để định vị mode cộng hưởng của hệ Bragg.

3. Sự định hướng của moment lưỡng cực nguyên tử quyết định những mode sóng nào đóng góp vào hàm Green tán xạ trong hệ trụ? Theo các biểu thức giải tích từ (2.21) đến (2.23) trong luận án:

  • Khi moment lưỡng cực của nguyên tử định hướng dọc theo trục đối xứng $z$, chỉ duy nhất các mode sóng từ ngang TM đóng góp vào hàm Green tán xạ.
  • Khi moment lưỡng cực định hướng theo phương bán kính $\rho$ hoặc phương góc $\phi$, cả hai loại mode sóng điện ngang TE và từ ngang TM đều đồng thời tham gia đóng góp.

4. Tương tác van der Waals giữa hai nguyên tử bên trong hệ trụ có thể tăng cường tối đa bao nhiêu lần so với không gian tự do? Trong cấu hình hai nguyên tử mức cơ bản đặt đối diện nhau trên cùng một đường kính gần biên của hệ trụ hai lớp, kết quả tính toán tại Chương 3 chỉ ra rằng thế tương tác và độ lớn của lực van der Waals có thể được tăng cường vượt bậc lên đến bốn bậc độ lớn ($10^4$ lần) so với tương tác giữa hai nguyên tử tương tự trong không gian tự do.

5. Cơ chế nào cho phép tạo bẫy nguyên tử ba chiều tại đầu mút của khối trụ hữu hạn? Tại Chương 4, việc khảo sát thế Casimir-Polder chuẩn hóa của nguyên tử đặt trên trục đối xứng $z$ ở vùng không gian bên ngoài gần sát đầu mút khối trụ hữu hạn ($H/\lambda_A = 10, R/\lambda_A = 1$) cho thấy năng lượng tương tác tạo thành một cực tiểu thế năng theo cả ba chiều không gian ($x_A, y_A, z_A$). Giếng thế năng cục bộ này xuất hiện do sự gián đoạn hình học tại mép biên khối trụ, tạo ra lực hút/đẩy cân bằng có khả năng giam giữ nguyên tử.

Kết luận

Luận án của nghiên cứu sinh Nguyễn Dũng Chỉnh đã nghiên cứu toàn diện tương tác tán sắc (Casimir-Polder và van der Waals) cùng tốc độ rã tự phát của nguyên tử gần các cấu trúc trụ vĩ mô bằng phương pháp hàm tensor Green trong điện động lực học lượng tử ba chiều. Công trình đã làm sáng tỏ cấu trúc mode cộng hưởng và vùng cấm quang học trong hệ trụ phản xạ Bragg đa lớp, chứng minh khả năng tăng cường tương tác van der Waals lên đến bốn bậc độ lớn trong cấu hình nguyên tử đối diện gần biên, đồng thời phát hiện cấu trúc giếng thế Casimir-Polder ba chiều tại đầu mút khối trụ hữu hạn bằng gần đúng Born. Các kết quả giải tích và giải số của luận án cung cấp cơ sở vật lý chuẩn xác cho việc kiểm soát tương tác nguyên tử - trường điện từ, định hướng ứng dụng trong chế tạo nguồn phát đơn photon và bẫy nguyên tử quang học.