Bối cảnh và vấn đề nghiên cứu

Việc giam giữ và bẫy các nguyên tử siêu lạnh trong các vùng không gian giới hạn là một chủ đề vật lý cơ bản thu hút sự quan tâm lớn trong nghiên cứu các hiệu ứng lượng tử của hệ vi mô. Các ứng dụng thực nghiệm nổi bật bao gồm việc nghiên cứu đặc tính phổ nguyên tử, đếm số lượng đồng vị phóng xạ, nâng cao độ chính xác của phép đo quang phổ và đặc biệt là khảo sát các hiệu ứng thống kê lượng tử ở nhiệt độ siêu thấp. Kể từ khi kỹ thuật làm lạnh bằng laser phát triển vào cuối thế kỷ XX, hiện tượng ngưng tụ Bose - Einstein (BEC) lần đầu tiên được Anderson và cộng sự quan sát trên nguyên tử $^{87}\text{Rb}$ vào năm 1995 ở nhiệt độ dưới $170\text{ nK}$, mở đường cho hàng loạt quan sát thực nghiệm tiếp theo trên các nguyên tử kiềm và kim loại kiềm thổ như $\text{Na}$, $\text{Li}$, $\text{H}$, $\text{He}$, $\text{K}$, $\text{Cs}$, $\text{Yb}$. Đến năm 2004, Kinoshita và cộng sự đã công bố kết quả thực nghiệm về sự hình thành trạng thái khí Tonks-Girardeau (TG) từ các boson tương tác mạnh trong không gian một chiều ($1\text{D}$), và năm 2008 nhóm của Cheinet đã quan sát thành công pha siêu lỏng và pha cách điện Mott của hệ boson trong mạng tinh thể quang học.

Mặc dù có nhiều công trình nghiên cứu về hệ khí lượng tử ít hạt, phần lớn các tiếp cận trước đây chỉ khảo sát các trường hợp riêng lẻ mà chưa đưa ra một phương pháp tổng quát giải quyết triệt để bài toán tương quan mạnh với tài nguyên tính toán tối ưu. Nhằm hoàn thiện bức tranh vật lý về hệ lượng tử ít hạt, khóa luận tập trung giải quyết đề tài: "Khảo sát một số tính chất dừng của hệ khí lượng tử ít hạt có tương tác trong các bẫy thế năng khác nhau ở $T = 0\text{ K}$".

Mục tiêu và các nhiệm vụ nghiên cứu cụ thể được xác định gồm:

  1. Xây dựng chương trình tính toán số dựa trên phương pháp chéo hóa ma trận chính xác cải tiến để giải phương trình Schrödinger dừng của hệ boson ít hạt trong không gian một chiều.
  2. Khảo sát các tính chất dừng trạng thái cơ bản của hệ boson ít hạt có tương tác tiếp xúc trong bẫy thế điều hòa (Harmonic Oscillator - HO) khi thay đổi cường độ tương tác từ không tương tác ($g=0$) đến giới hạn tương tác rất mạnh ($g=20$, $g=200$) với số hạt $N = 2, 3, 4$.
  3. Khảo sát các tính chất dừng trạng thái cơ bản của hệ boson trong bẫy mạng tinh thể quang học (Optical Lattice - OL) 3 giếng thế ($p=3$) với số hạt $N = 2, 3$ khi thay đổi cường độ tương tác từ $g=0$ đến $g=1$.
  4. Đối chiếu các kết quả giải số thu được với các nghiệm giải tích chính xác của mô hình khí Tonks-Girardeau và tiêu chuẩn ngưng tụ Penrose - Onsager để kiểm chứng độ chính xác của thuật toán.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Trạng thái cơ bản của hệ khí lượng tử gồm $N$ hạt boson siêu lạnh có khối lượng $m$, tương tác tầm ngắn (tiếp xúc) dạng hàm delta $\delta(x_i - x_j)$ tại nhiệt độ không tuyệt đối $T = 0\text{ K}$.
  • Phạm vi không gian và thế năng: Không gian một chiều ($1\text{D}$) với chiều dài tọa độ $L = 20$; hai dạng thế năng giam giữ gồm bẫy điều hòa ($V_{\text{ext}}(x) = \frac{1}{2}m\omega^2 x^2$ với chuẩn hóa $m = \omega = 1$) và mạng tinh thể quang học ($V_{\text{ext}}(x) = V_0 \sin^2(\frac{\pi p x}{L})$ với biên độ $V_0 = 1$, số giếng $p = 3$). Số lượng hạt khảo sát giới hạn ở $N = 2, 3, 4$ đối với bẫy điều hòa và $N = 2, 3$ đối với bẫy mạng tinh thể quang học.

Cơ sở lý thuyết và phương pháp

Khóa luận xây dựng trên nền tảng cơ học lượng tử nhiều hạt, cơ học thống kê lượng tử và các mô hình giải tích chuẩn trong không gian một chiều:

Khung lý thuyết và mô hình vật lý

  • Mô hình Hamiltonian 1D tổng quát: Hệ $N$ boson tương tác tầm ngắn được mô tả bởi toán tử năng lượng: $$\hat{H} = \sum_{i=1}^{N} \left[ -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x_i^2} + V_{\text{ext}}(x_i) \right] + g \sum_{i<j}^{N} \delta(x_i - x_j)$$ Trong đó $g > 0$ đặc trưng cho tương tác đẩy. Cường độ tương tác này có thể điều khiển trong thực nghiệm thông qua cộng hưởng Feshbach hoặc hiệu chỉnh tần số bẫy theo phương ngang.
  • Bài toán Lieb - Liniger: Mô tả hệ $N$ boson tương tác đẩy trong hộp thế tuần hoàn ($V_{\text{ext}} = 0$) với điều kiện biên tuần hoàn, được giải chính xác bằng phép thế sóng thử Bethe (Bethe Ansatz).
  • Mô hình Tonks - Girardeau (TG) và định lý ánh xạ Bose - Fermi: Trong giới hạn tương tác đẩy vô cùng lớn ($g \to \infty$), các boson không thể chiếm cùng một vị trí trong không gian ($|\Psi(x_1, \dots, x_N)| = 0$ khi $x_i = x_j$), hành xử như các boson lõi cứng (hard-core bosons). Theo định lý Girardeau (1960), hàm sóng nhiều hạt boson trạng thái cơ bản $\Psi_B$ được ánh xạ trực tiếp từ định thức Slater của hàm sóng hệ $N$ fermion không tương tác spin phân cực $\Psi_F$ thông qua hàm phản đối xứng đơn vị: $\Psi_B = |\Psi_F| = \prod_{1 \le j < k \le N} \text{sgn}(x_k - x_j) \Psi_F$.
  • Gần đúng trường trung bình Gross - Pitaevskii (GP): Áp dụng cho hệ nhiều hạt tương tác yếu ($N \gg 1, g \to 0$). Tuy nhiên, đối với hệ ít hạt tương tác mạnh, gần đúng trường trung bình bị phá vỡ hoàn toàn do không mô tả được các hiệu ứng tương quan lượng tử trực tiếp giữa hai hạt, đòi hỏi phải sử dụng phương pháp giải số nhiều hạt vi mô.
  • Các đại lượng đặc trưng:
    • Ma trận mật độ đơn hạt rút gọn (Reduced Single-Particle Density Matrix - RSPDM): $\rho_1(x, x') = N \int \Psi^*(x, x_2, \dots, x_N) \Psi(x', x_2, \dots, x_N) dx_2 \dots dx_N$. Đường chéo $\rho(x) = \rho_1(x, x)$ là hàm phân bố mật độ đơn hạt. Tính chất ngoài đường chéo đặc trưng cho trật tự tầm xa ngoài đường chéo (ODLRO) và tính kết hợp lượng tử.
    • Số chiếm đóng tự nhiên (Natural Occupancy): Các trị riêng $\lambda_j$ nhận được từ phương trình tích phân $\int \rho_1(x, x') \psi_j^{(\text{NO})}(x') dx' = \lambda_j \psi_j^{(\text{NO})}(x)$. Theo tiêu chuẩn Penrose - Onsager, hệ ở pha BEC hoàn toàn khi trị riêng lớn nhất $\lambda_0/N = 1$; nếu có nhiều trị riêng $\lambda_j$ cùng bậc với $N$, hệ rơi vào trạng thái phân mảnh (fragmentation).
    • Hàm phân bố cặp (Pair Distribution Function - PDF): $\rho_2(x_1, x_2) = N(N-1) \int |\Psi(x_1, x_2, x_3, \dots, x_N)|^2 dx_3 \dots dx_N$, thể hiện xác suất tìm thấy đồng thời hai hạt tại tọa độ $x_1$ và $x_2$.

Phương pháp tính toán số

  • Phương pháp chéo hóa ma trận chính xác cải tiến (Improved Exact Diagonalization - ED): Biểu diễn Hamiltonian sang hình thức lượng tử hóa lần hai bằng các toán tử sinh $\hat{a}k^\dagger$ và hủy $\hat{a}k$ của $M$ trạng thái đơn hạt cơ sở: $$\hat{H} = \sum{k, j=0}^{M-1} H{kj} \hat{a}k^\dagger \hat{a}j + \frac{1}{2} \sum{k, j, s, t=0}^{M-1} V{kjst} \hat{a}_k^\dagger \hat{a}_j^\dagger \hat{a}_t \hat{a}_s$$ Không gian Fock truyền thống có số chiều $D = \frac{(N+M-1)!}{N!(M-1)!}$. Khóa luận áp dụng thuật toán cải tiến của Trần Dương Anh Tài và cộng sự (2023): dựa vào tính chất đối xứng chẵn - lẻ của các trạng thái đơn hạt trong thế đối xứng, hàm sóng nhiều boson chỉ nhận đóng góp từ các trạng thái Fock có tính đối xứng chẵn, các trạng thái Fock lẻ bị triệt tiêu hoàn toàn. Việc loại bỏ các vector cơ sở lẻ giúp giảm số chiều không gian Hilbert $D'$ xuống xấp xỉ một nửa ($D' \approx D/2$), tiết kiệm đáng kể bộ nhớ RAM và thời gian tính toán.
  • Biểu diễn biến số gián đoạn (Discrete Variable Representation - DVR): Sử dụng kỹ thuật DVR theo sơ đồ Colbert - Miller (1992) với $N_{\text{DVR}} = 2^9 + 1 = 513$ điểm lưới trên đoạn $[-L/2, L/2]$ để chéo hóa Hamiltonian đơn hạt, tìm chính xác tập hàm riêng $\phi_k(x)$ và trị riêng $\varepsilon_k$.
  • Công cụ triển khai: Chương trình tính toán được lập trình bằng ngôn ngữ MATLAB, sử dụng hàm eigs để giải bài toán trị riêng ma trận thưa quy mô lớn, kết hợp module xây dựng không gian Fock kế thừa từ phần mềm TDED-O.

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

Chương 1: Tổng quan

Chương 1 hệ thống hóa lịch sử phát triển của lĩnh vực khí lượng tử siêu lạnh từ lý thuyết BEC năm 1925 của Einstein đến các mốc thực nghiệm trên hệ boson ($^{87}\text{Rb}, \text{Na}, \text{Li}, \dots$) và fermion siêu lạnh. Tác giả trình bày chi tiết nguyên lý bẫy nguyên tử trong mạng tinh thể quang học (OL) hình thành từ sự giao thoa laser kết hợp, đồng thời phân tích cơ chế điều khiển tương tác qua cộng hưởng Feshbach. Về mặt lý thuyết, chương này thiết lập mô hình toán học của hệ $1\text{D}$, làm rõ phạm vi áp dụng và giới hạn của mô hình giải tích Lieb - Liniger, mô hình Tonks - Girardeau với phép biến đổi sóng đối xứng - phản đối xứng, và phương trình Gross - Pitaevskii. Phần cuối chương định nghĩa toán học chặt chẽ cho các đại lượng đo lường lượng tử: ma trận mật độ đơn hạt rút gọn $\rho_1(x, x')$, hàm phân bố mật độ $\rho(x)$, hàm phân bố cặp $\rho_2(x_1, x_2)$, và số chiếm đóng tự nhiên $\lambda_j$ làm tiêu chí định lượng chuyển pha lượng tử.

Chương 2: Phương pháp chéo hóa ma trận chính xác cải tiến

Chương 2 tập trung vào việc thiết lập thuật toán và quy trình giải số. Tác giả trình bày chi tiết việc chuyển đổi Hamiltonian liên tục sang biểu diễn lượng tử hóa lần hai trên cơ sở các trạng thái Fock $|F_i\rangle = |n_0, n_1, \dots, n_{M-1}\rangle$ thỏa mãn điều kiện bảo toàn số hạt $\sum_{k=0}^{M-1} n_k = N$. Điểm mấu chốt của chương là việc chứng minh sự phân tách không gian Hilbert thành hai không gian con độc lập dựa trên tính đối xứng chẵn - lẻ. Bằng việc loại bỏ toàn bộ các trạng thái Fock lẻ không đóng góp vào hàm sóng nhiều boson đối xứng, số chiều ma trận cần chéo hóa giảm xấp xỉ 50%. Chương này cũng mô tả chi tiết biểu thức ma trận của các tích phân đơn hạt $H_{kj}$ và tích phân tương tác hai hạt $V_{kjst}$, quy tắc tác dụng của toán tử sinh hủy $\hat{a}_k^\dagger, \hat{a}k$, và cách trích xuất hàm sóng trạng thái cơ bản $|\Psi{\text{ground}}\rangle$ cùng các đại lượng vật lý liên quan từ vector riêng thu được.

Chương 3: Kết quả nghiên cứu

Chương 3 trình bày các kết quả mô phỏng số chi tiết thu được từ chương trình MATLAB với $M = 24$ trạng thái đơn hạt cơ sở và $N_{\text{DVR}} = 513$ điểm lưới:

  1. Thông số không gian Hilbert: Xác lập bảng đối chiếu kích thước ma trận trước và sau cải tiến cho các trường hợp số hạt $N = 2, 3, 4$, khẳng định hiệu quả giảm tải tài nguyên tính toán 50%.
  2. Hệ khí trong bẫy thế điều hòa (HO):
    • Năng lượng trạng thái cơ bản: Khi $g$ tăng từ 0 đến 20, năng lượng $E_0$ tăng đơn điệu từ mức không tương tác ($E_0 = N \times 0.5$) và tiệm cận giá trị giải tích của khí TG ($E_{\text{TG}} = 2.0$ cho $N=2$; $4.5$ cho $N=3$; $8.0$ cho $N=4$).
    • RSPDM và hàm phân bố mật độ $\rho(x)$: Ở $g=0$, $\rho(x)$ có dạng Gauss định xứ tại tâm bẫy. Khi $g$ tăng, tương tác đẩy làm mật độ trải rộng; tại $g=20$, đồ thị $\rho(x)$ phân tách rõ rệt thành đúng $N$ đỉnh độc lập, trùng khớp hoàn hảo với đường nghiệm giải tích TG. Đồng thời, các thành phần ngoài đường chéo của $\rho_1(x, x')$ suy giảm mạnh, báo hiệu sự triệt tiêu trật tự ODLRO.
    • Hàm phân bố cặp PDF: Tại $g=0$, hai hạt có xác suất chồng lấn cực đại tại cùng tọa độ ($x_1 = x_2$). Khi $g \to 20$, giá trị PDF dọc theo đường chéo chính $x_1 = x_2$ hoàn toàn triệt tiêu về 0, chứng minh thực nghiệm số cho điều kiện biên của boson lõi cứng.
    • Số chiếm đóng tự nhiên $\lambda_j$: Tại $g=0$, $\lambda_0/N = 1$ (100% ngưng tụ BEC). Khi $g=20$, $\lambda_0/N$ giảm mạnh và các trạng thái kích thích $\lambda_1, \lambda_2, \dots$ bắt đầu bị chiếm đóng. Khảo sát tiệm cận đến $g=200$ cho thấy số mũ suy giảm $\lambda_0 \sim N^{0.52}$, bám sát giá trị giải tích lý thuyết $N^{0.5}$ của Girardeau.
  3. Hệ khí trong bẫy mạng tinh thể quang học (OL):
    • Giải bài toán đơn hạt bằng DVR để xác định các mức năng lượng và hàm sóng đơn hạt trong thế $V_0 = 1, p = 3$.
    • Khi $g=0$, năng lượng hệ bằng $N \varepsilon_0$. Khi $g=1$ (tương tác mạnh), năng lượng tiệm cận giới hạn TG.
    • Hiệu ứng lấp đầy (Filling factor $\nu = N/p$): Với $N=3$ ($\nu = 3/3 = 1$ là số nguyên), tại giới hạn tương tác mạnh $g=1$, mật độ hạt phân bố đều tại 3 giếng thế, các phần tử ngoài đường chéo của $\rho_1(x, x')$ hoàn toàn biến mất, hệ chuyển đổi từ pha siêu lỏng sang pha cách điện Mott cục bộ với tỷ lệ chiếm đóng trên 3 quỹ đạo tự nhiên đầu tiên tiệm cận đều về $1/3$. Ngược lại, với $N=2$ ($\nu = 2/3$ không nguyên), mật độ ở giếng trung tâm cao hơn hai giếng biên, thành phần ngoài đường chéo của $\rho_1(x, x')$ vẫn tồn tại một phần do tồn tại thành phần phi định xứ ngăn cản sự hình thành pha cách điện Mott thuần túy.

Kết quả và đóng góp

Bảng số liệu tổng hợp trong nghiên cứu

Bảng 1: Số chiều không gian Hilbert theo phương pháp truyền thống ($D$) và cải tiến ($D'$) với $M = 24$ trạng thái đơn hạt

Số hạt ($N$) Số trạng thái đơn hạt ($M$) Số chiều truyền thống ($D$) Số chiều cải tiến ($D'$) Tỷ lệ giảm kích thước
$N = 2$ 24 300 156 $\approx 48.0%$
$N = 3$ 24 2600 1300 $50.0%$
$N = 4$ 24 17550 8814 $\approx 49.8%$

Bảng 2: Năng lượng trạng thái cơ bản trong bẫy điều hòa theo cường độ tương tác $g$ (đơn vị $\hbar\omega$)

Số hạt ($N$) $g = 0$ (Không tương tác) $g \to \infty$ (Giới hạn giải tích TG) Công thức giải tích $E_{\text{TG}}$
$N = 2$ 1.0000 2.0000 $E_{\text{TG}} = 0.5 + 1.5 = 2.0$
$N = 3$ 1.5000 4.5000 $E_{\text{TG}} = 0.5 + 1.5 + 2.5 = 4.5$
$N = 4$ 2.0000 8.0000 $E_{\text{TG}} = 0.5 + 1.5 + 2.5 + 3.5 = 8.0$

Bảng 3: So sánh đặc tính hệ khí trong mạng tinh thể quang học ($p = 3$ giếng)

Trường hợp Hệ số lấp đầy ($\nu = N/p$) Phân bố mật độ $\rho(x)$ ở $g=1$ Trạng thái ngoài đường chéo $\rho_1(x, x')$ Trạng thái pha lượng tử
$N = 2$ $2/3$ (Không nguyên) Giếng giữa cao hơn 2 giếng bên Còn tồn tại thành phần kết hợp Ngăn cản pha Mott, duy trì một phần phi định xứ
$N = 3$ $3/3 = 1$ (Số nguyên) 3 đỉnh bằng nhau tại 3 giếng Hoàn toàn triệt tiêu Chuyển pha sang cách điện Mott, $\lambda_0 \approx \lambda_1 \approx \lambda_2 \to 1/3$

Đóng góp khoa học chính

  • Về mặt phương pháp: Chứng minh thực nghiệm tính hiệu quả và độ chính xác cao của phương pháp chéo hóa ma trận chính xác cải tiến (loại bỏ trạng thái Fock lẻ), giúp giảm chính xác $50%$ số chiều không gian Hilbert mà không làm mất độ chính xác của trị riêng và hàm sóng trạng thái cơ bản.
  • Về mặt vật lý:
    • Mô tả định lượng quá trình chuyển pha liên tục từ ngưng tụ Bose - Einstein thuần túy ($\lambda_0/N = 1$) sang khí Tonks - Girardeau (các boson bị fermion hóa, tách thành $N$ đỉnh định xứ riêng biệt và triệt tiêu PDF tại $x_1=x_2$) trong bẫy điều hòa.
    • Khảo sát sự phụ thuộc tiệm cận của số chiếm đóng tự nhiên theo số hạt ở giới hạn tương tác mạnh ($g=200$), thu được quy luật $\lambda_0 \sim N^{0.52}$, phù hợp chặt chẽ với dự đoán lý thuyết của Girardeau ($\lambda_0 \sim N^{0.5}$).
    • Làm sáng tỏ ảnh hưởng cốt lõi của hệ số lấp đầy nguyên ($\nu = 1$) và không nguyên ($\nu = 2/3$) đối với sự hình thành pha cách điện Mott và sự triệt tiêu trật tự tầm xa ngoài đường chéo trong bẫy mạng tinh thể quang học ít hạt.

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

Hạn chế của đề tài

  • Phạm vi số chiều không gian: Nghiên cứu chỉ giới hạn trong không gian một chiều ($1\text{D}$), chưa xét đến các hiệu ứng hình học phức tạp hơn trong không gian 2D hoặc 3D.
  • Quy mô số hạt và trạng thái cơ sở: Do giới hạn về bộ nhớ RAM và năng lực tính toán xử lý ma trận lớn, số lượng hạt khảo sát chỉ dừng lại ở $N \le 4$ đối với bẫy điều hòa và $N \le 3$ đối với mạng tinh thể quang học trên cơ sở $M = 24$ trạng thái đơn hạt.
  • Điều kiện nhiệt độ và trạng thái động lực: Nghiên cứu mới chỉ tập trung vào các tính chất dừng của trạng thái cơ bản tại $T = 0\text{ K}$, chưa khảo sát các trạng thái kích thích, hiệu ứng nhiệt độ hữu hạn ($T > 0\text{ K}$) hoặc động lực học phụ thuộc thời gian khi thay đổi đột ngột thế bẫy (quench dynamics).

Hướng phát triển tiếp theo

  • Mở rộng thuật toán và mã nguồn tính toán số để khảo sát hệ boson chuyển động trong không gian hai chiều ($2\text{D}$) và ba chiều ($3\text{D}$).
  • Ứng dụng nền tảng tính toán hiện tại để nghiên cứu các hệ khí lượng tử hỗn hợp (hỗn hợp Bose - Fermi, Bose - Bose) và khảo sát động lực học lượng tử phụ thuộc thời gian của hệ ít hạt trong các bẫy thế quang học động.

Giá trị tham khảo

Khóa luận là tài liệu học thuật và kỹ thuật tính toán có giá trị tham khảo thiết thực cho:

  • Sinh viên và học viên cao học ngành Vật lý: Đặc biệt là chuyên ngành Vật lý lý thuyết, Vật lý tính toán và Vật lý nguyên tử - phân tử - quang học (AMO), cần nắm bắt cách áp dụng hình thức lượng tử hóa lần hai và phương pháp chéo hóa ma trận chính xác vào các bài toán lượng tử nhiều hạt.
  • Người làm nghiên cứu về phương pháp số lượng tử: Có thể tham khảo trực tiếp Chương 2 và Chương 3 để triển khai thuật toán biểu diễn biến số gián đoạn (DVR Colbert - Miller) kết hợp với kỹ thuật tối ưu hóa không gian Fock theo đối xứng chẵn - lẻ nhằm giải quyết các bài toán Hamiltonian ma trận kích thước lớn trên MATLAB hoặc Python.
  • Các nhóm nghiên cứu về mô phỏng lượng tử: Cung cấp các bộ dữ liệu chuẩn (benchmark) về năng lượng trạng thái cơ bản, ma trận mật độ $\rho_1$, hàm phân bố cặp $\rho_2$, và số chiếm đóng tự nhiên $\lambda_j$ để đối chiếu khi phát triển các phương pháp gần đúng khác như Monte Carlo lượng tử (QMC) hay mạng tensor (DMRG/MPS).

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

1. Điểm cải tiến mấu chốt của phương pháp chéo hóa chính xác (ED) trong khóa luận này là gì?
Phương pháp cải tiến dựa trên công trình của Trần Dương Anh Tài và cộng sự (2023), tận dụng tính chất đối xứng không gian của thế năng để phân tách không gian Hilbert thành hai nhánh chẵn và lẻ. Do hàm sóng nhiều boson có tính đối xứng hoàn toàn theo tọa độ, các trạng thái cơ sở Fock có tính chất lẻ không đóng góp vào hàm sóng trạng thái cơ bản. Việc loại bỏ hoàn toàn các trạng thái Fock lẻ giúp giảm xấp xỉ $50%$ số chiều không gian Hilbert ($D'$ so với $D$), từ đó giảm tải một nửa nhu cầu bộ nhớ RAM và tăng tốc độ giải số ma trận Hamiltonian.

2. Tại sao phương trình Gross - Pitaevskii (GP) không được sử dụng để giải bài toán trong nghiên cứu này?
Phương trình GP là phương trình Schrödinger phi tuyến dựa trên gần đúng trường trung bình, chỉ phù hợp cho hệ khí nhiều hạt ($N \gg 1$) và tương tác yếu ($g \to 0$), khi tất cả các hạt cùng ngưng tụ trên một hàm sóng đơn hạt. Đối với hệ khí ít hạt ($N = 2, 3, 4$) và ở vùng tương tác trung bình đến rất mạnh ($g = 20, g = 200$), tương quan tiếp xúc trực tiếp giữa từng cặp hạt đóng vai trò chi phối và làm phân mảnh hệ (phá vỡ trạng thái ngưng tụ đơn hạt), khiến gần đúng trường trung bình của GP không còn chính xác.

3. Khí Tonks - Girardeau (TG) là gì và khóa luận kiểm chứng trạng thái này qua những tiêu chuẩn nào?
Khí TG là trạng thái của hệ boson $1\text{D}$ trong giới hạn tương tác đẩy cực mạnh ($g \to \infty$), khi đó các boson hành xử như "boson lõi cứng" và có tính chất tương tự hệ fermion không tương tác (định lý ánh xạ Bose - Fermi). Khóa luận kiểm chứng trạng thái này qua 3 tiêu chuẩn:

  • Năng lượng trạng thái cơ bản tiến sát giá trị giải tích $E_{\text{TG}} = \sum_{j=0}^{N-1} \varepsilon_j$.
  • Hàm phân bố mật độ $\rho(x)$ phân tách thành đúng $N$ đỉnh định xứ riêng biệt.
  • Hàm phân bố cặp $\rho_2(x_1, x_2)$ dọc theo đường chéo chính $x_1 = x_2$ hoàn toàn triệt tiêu về 0 (xác nhận hai hạt không thể chiếm cùng vị trí).

4. Tại sao có sự khác biệt về ma trận mật độ ngoài đường chéo $\rho_1(x, x')$ giữa trường hợp $N=2$ và $N=3$ trong mạng tinh thể quang học 3 giếng?
Sự khác biệt bắt nguồn từ hệ số lấp đầy $\nu = N/p$:

  • Với $N=3, p=3 \implies \nu = 1$ (hệ số lấp đầy nguyên), khi tương tác mạnh ($g=1$), mỗi giếng chứa đúng 1 hạt, hệ chuyển pha hoàn toàn sang cách điện Mott cục bộ, triệt tiêu tính kết hợp lượng tử giữa các giếng, khiến phần tử ngoài đường chéo của $\rho_1(x, x')$ biến mất.
  • Với $N=2, p=3 \implies \nu = 2/3$ (không nguyên), mật độ hạt không thể phân bố đều ở cả 3 giếng mà tập trung ở giếng giữa nhiều hơn. Sự không tương thích này tạo ra phần mật độ phi định xứ, ngăn cản sự hình thành pha Mott hoàn toàn và giữ lại các thành phần kết hợp ngoài đường chéo.

5. Các thông số kỹ thuật lưới và cơ sở đơn hạt được sử dụng trong chương trình tính toán số là bao nhiêu?
Toàn bộ tính toán số sử dụng số trạng thái đơn hạt cơ sở $M = 24$, số điểm lưới không gian gián đoạn $N_{\text{DVR}} = 2^9 + 1 = 513$ điểm trên đoạn tọa độ có chiều dài $L = 20$, bước lưới $\Delta x = L/(N_{\text{DVR}} - 1)$. Việc chéo hóa ma trận được thực hiện bằng thuật toán eigs trong môi trường MATLAB kết hợp thuật toán sinh cơ sở Fock từ phần mềm TDED-O.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Ngọc Giàu đã khảo sát thành công các tính chất dừng trạng thái cơ bản của hệ khí boson siêu lạnh ít hạt ($N = 2, 3, 4$) trong bẫy điều hòa và mạng tinh thể quang học một chiều tại $T = 0\text{ K}$. Việc ứng dụng phương pháp chéo hóa ma trận chính xác cải tiến kết hợp biểu diễn biến số gián đoạn DVR đã giúp tối ưu hóa $50%$ kích thước không gian Hilbert mà vẫn đảm bảo độ hội tụ và độ chính xác cao khi đối chiếu với nghiệm giải tích của mô hình Tonks - Girardeau. Nghiên cứu đã làm sáng tỏ cơ chế fermion hóa của boson khi tăng cường độ tương tác đẩy và ảnh hưởng mang tính quyết định của hệ số lấp đầy đối với quá trình định xứ Mott trong bẫy quang học. Kết quả và hệ chương trình tính toán phát triển trong khóa luận là tiền đề kỹ thuật vững chắc để mở rộng khảo sát động lực học và các mô hình khí lượng tử trong không gian nhiều chiều.