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

Kỹ thuật lượng tử (quantum engineering) và tính toán lượng tử (quantum computing) là các lĩnh vực nghiên cứu trọng điểm của vật lý hiện đại với tiềm năng ứng dụng trong khoa học thông tin, học máy và trí tuệ nhân tạo. Lĩnh vực này gắn kết chặt chẽ với quang tử học và vật lý các nguyên tử siêu lạnh (ultracold atoms). Nghiên cứu khí lượng tử siêu lạnh bắt đầu phát triển mạnh mẽ từ năm 1995 khi Eric Cornell, Carl Wieman và cộng sự thực hiện thành công thí nghiệm ngưng tụ Bose - Einstein (Bose-Einstein Condensation - BEC) trên khí $^{87}\text{Rb}$ loãng ở nhiệt độ $170\text{ nK}$. Sau đó, hiện tượng ngưng tụ tiếp tục được quan sát trên các nguyên tố kiềm khác như $^{23}\text{Na}$, $^{7}\text{Li}$, $^{1}\text{H}$, $^{85}\text{Rb}$ và $^{133}\text{Cs}$.

Sự tiến bộ của các kỹ thuật thực nghiệm như bẫy quang từ (MOT), vi mạch nguyên tử (atom chip) và cộng hưởng Feshbach cho phép giới hạn chuyển động của nguyên tử về không gian một chiều (1D), đồng thời kiểm soát chính xác số lượng hạt (từ hệ ít hạt đến nhiều hạt) và cường độ tương tác. Trong không gian 1D, tương quan lượng tử gia tăng rõ rệt, dẫn đến nhiều hiện tượng vật lý phi cổ điển như khí tương tác mạnh Tonks - Girardeau (TG), fermion hóa tổng hợp (Composite Fermionization - CF), fermion hóa hoàn toàn (Full Fermionization - FF) và hiện tượng tách pha (Phase Separation - PS).

Các nghiên cứu trước đây về hệ boson siêu lạnh ít hạt hai thành phần chủ yếu khảo sát trong bẫy thế điều hòa. Khóa luận này được thực hiện nhằm mở rộng việc nghiên cứu các tính chất dừng của hệ ít hạt hai thành phần sang trường thế năng của mạng tinh thể quang học một chiều.

Mục tiêu và nhiệm vụ nghiên cứu của đề tài bao gồm:

  1. Xây dựng mô hình lý thuyết mô tả hệ khí boson siêu lạnh hai thành phần trong mạng tinh thể quang học một chiều bằng phương pháp lượng tử hóa lần hai.
  2. Áp dụng phương pháp giải số chéo hóa ma trận chính xác (Exact Diagonalization - ED) kết hợp biểu diễn biến số rời rạc (Discrete Variable Representation - DVR) và thuật toán cải tiến để giải phương trình Schrödinger nhiều hạt.
  3. Kiểm chứng chương trình tính toán thông qua bài toán hệ một thành phần ở giới hạn tương tác mạnh Tonks - Girardeau giải tích.
  4. Khảo sát các đại lượng đặc trưng (ma trận mật độ đơn hạt rút gọn, hàm phân bố cặp, số chiếm đóng tự nhiên) của hệ boson hai thành phần ($2+2$ hạt) tại trạng thái cơ bản trong mạng tinh thể quang học 3 giếng thế ($p=3$) dưới các điều kiện tương tác và chiều cao thế khác nhau.
  5. Khảo sát và phân tích các đại lượng đặc trưng tương tự đối với mạng tinh thể quang học 2 giếng thế ($p=2$), làm rõ sự khác biệt về trạng thái cách điện Mott và cấu trúc tách pha do tính đối xứng và hệ số lấp đầy.

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

  • Đối tượng nghiên cứu: Hệ khí boson siêu lạnh ít hạt gồm hai thành phần $A$ và $B$ ($N_A = 2, N_B = 2$, tổng cộng 4 hạt).
  • Phạm vi không gian: Mạng tinh thể quang học một chiều có bề rộng $L = 40$ ($x \in [-20; 20]$), số lượng giếng thế khảo sát là $p = 2$ và $p = 3$.
  • Phạm vi tương tác và thế năng: Độ cao giếng thế $V_0$ biến thiên từ $0.01$ đến $0.5$; cường độ nội tương tác $g_A = 3, g_B = 0$; cường độ tương tác liên thành phần $g_{AB}$ khảo sát từ $0$ đến $3$.

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

Khóa luận khảo sát hệ boson hai thành phần $A$ và $B$ ($\sigma \in {A, B}$) chuyển động trong không gian một chiều dưới tác dụng của trường thế năng ngoài $V_{ext}(x)$ và tương tác tiếp xúc dạng hàm delta Olshanii (1998). Hamiltonian của hệ được mô tả bởi biểu thức:

$$H = \sum_{\sigma \in {A, B}} \sum_{i=1}^{N_\sigma} \left[ -\frac{\hbar^2}{2m} \frac{\partial^2}{\partial (x_i^\sigma)^2} + V_{ext}(x_i^\sigma) \right] + \sum_{\sigma} \sum_{i<j}^{N_\sigma} g_\sigma \delta(x_i^\sigma - x_j^\sigma) + \sum_{i=1}^{N_A} \sum_{j=1}^{N_B} g_{AB} \delta(x_i^A - x_j^B)$$

Thế năng mạng tinh thể quang học được biểu diễn theo số lượng giếng thế $p$:

  • Trường hợp số giếng thế lẻ ($p = 1, 3, 5...$): $V_{ext}(x) = V_0 \sin^2\left(\frac{p\pi x}{L}\right)$ với $x \in \left[-\frac{L}{2}; \frac{L}{2}\right]$
  • Trường hợp số giếng thế chẵn ($p = 2, 4, 6...$): $V_{ext}(x) = V_0 \cos^2\left(\frac{p\pi x}{L}\right)$ với $x \in \left[-\frac{L}{2}; \frac{L}{2}\right]$

Các đơn vị trong tính toán được chuẩn hóa không thứ nguyên: năng lượng tỉ lệ với $\mathcal{E}_0 = \frac{\hbar^2}{m L^2}$, chiều dài tỉ lệ với $L$, và cường độ tương tác tỉ lệ với $g_0 = \frac{\hbar^2}{m L}$.

Chế độ $g_A$ $g_B$ $g_{AB}$ Hiện tượng vật lý tương ứng
1 $0$ $0$ $0$ Hai thành phần ở trạng thái BEC (BEC-BEC)
2 $\infty$ $0$ $0$ Thành phần A ở trạng thái TG, thành phần B ở trạng thái BEC
3 $\infty$ $\infty$ $0$ Cả hai thành phần ở trạng thái TG (TG-TG)
4 $\infty$ $0$ $g_{AB} \to \infty$ Fermion hóa tổng hợp (Composite Fermionization - CF)
5 $\infty$ $0$ $\infty$ Hiện tượng tách pha (Phase Separation - PS)
6 $g_A = g_B = g_{AB} \to \infty$ Fermion hóa hoàn toàn (Full Fermionization - FF)

Để khảo sát trạng thái dừng, khóa luận sử dụng các đại lượng đặc trưng lượng tử:

  • Ma trận mật độ đơn hạt rút gọn (RSPDM): $\rho_1^\sigma(x, x') = N_\sigma \int \Psi^*(x, x_2^\sigma, ..., x_{N_A}^A, x_1^B, ...) \Psi(x', x_2^\sigma, ..., x_{N_A}^A, x_1^B, ...) , dx_2^\sigma ... dx_{N_B}^B$. Đường chéo chính của RSPDM biểu diễn hàm phân bố mật độ $\rho^\sigma(x) = \rho_1^\sigma(x, x)$.
  • Quỹ đạo tự nhiên (NO) và số chiếm đóng tự nhiên ($\lambda_j^\sigma$): Thu được từ việc chéo hóa RSPDM $\rho_1^\sigma$, thỏa mãn $\sum_j \lambda_j^\sigma = N_\sigma$.
  • Hàm phân bố cặp (PDF): $\rho_2^\sigma(x_1, x_2)$ xác định xác suất tìm thấy hai hạt cùng loại $\sigma$ tại hai vị trí $x_1$ và $x_2$.
  • Hàm phân bố cặp chéo (CTBDF): $\rho_2^{AB}(x_1^A, x_2^B)$ xác định tương quan vị trí giữa một hạt thành phần $A$ và một hạt thành phần $B$.

Phương pháp tính toán: Phương pháp chéo hóa ma trận chính xác (ED) được triển khai qua biểu diễn lượng tử hóa lần hai. Không gian Fock $|F_w\rangle$ được xây dựng từ tập $M_\sigma$ trạng thái đơn hạt $|\phi_i\rangle$, vốn được giải chính xác từ phương trình Schrödinger không phụ thuộc thời gian (TISE) một hạt bằng phương pháp biến số rời rạc (DVR). Số chiều không gian Hilbert toàn phần được xác định theo tổ hợp:

$$D = \prod_{\sigma \in {A, B}} \frac{(N_\sigma + M_\sigma - 1)!}{N_\sigma!(M_\sigma - 1)!}$$

Nhằm xử lý hạn chế về dung lượng bộ nhớ khi số chiều $D$ tăng nhanh, khóa luận ứng dụng thuật toán ED cải tiến của Trần Dương Anh Tài và cộng sự (2023). Thuật toán khai thác tính đối xứng hoán vị của hàm sóng boson để chỉ giữ lại các trạng thái Fock chẵn, áp dụng phương pháp tương tác hiệu dụng (effective interaction approach) và giới hạn năng lượng cắt ngưỡng $E_{max}$, giúp rút giảm số chiều không gian ma trận nhưng vẫn đảm bảo độ hội tụ chính xác. Toàn bộ chương trình tính toán số được lập trình và thực thi trên môi trường MATLAB.

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

Chương 1. Tổng quan

Chương 1 hệ thống hóa cơ sở lý thuyết về khí boson siêu lạnh và các hệ một chiều. Tác giả trình bày lịch sử từ thí nghiệm BEC năm 1995 đến thí nghiệm tạo ngưng tụ hai thành phần $^{87}\text{Rb}$ ở hai trạng thái spin của Myatt (1996), cùng các phương pháp kiểm soát bẫy thế và cộng hưởng Feshbach. Chương này phân tích chi tiết 6 cơ chế tương tác giới hạn giữa các thành phần và thiết lập các công thức toán học định nghĩa RSPDM, phân bố mật độ, hàm phân bố cặp (PDF, CTBDF), quỹ đạo tự nhiên và số chiếm đóng tự nhiên.

Chương 2. Phương pháp chéo hóa ma trận chính xác

Chương 2 trình bày chi tiết quy trình giải số hệ nhiều hạt:

  • Thiết lập toán tử sinh/hủy trường boson $\hat{\psi}_\sigma(x)$ và khai triển Hamiltonian nhiều hạt sang dạng ma trận trong không gian Fock.
  • Ứng dụng phương pháp DVR để tìm hệ hàm riêng đơn hạt $|\phi_i\rangle$ từ phương trình TISE một hạt.
  • Trình bày thuật toán cải tiến của Trần Dương Anh Tài và cộng sự (2023): phân tích tính đối xứng hoán vị để loại bỏ các trạng thái lẻ, kết hợp phương pháp tương tác hiệu dụng và đặt ngưỡng năng lượng $E_{max}$ nhằm giảm số chiều ma trận $D$ mà không làm mất độ chính xác của các trạng thái năng lượng thấp.

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

Chương 3 trình bày các kết quả tính toán số cho hệ 1 thành phần và hệ 2 thành phần:

  • Kiểm chứng hệ 1 thành phần: Tính toán năng lượng trạng thái cơ bản khi tăng cường độ nội tương tác $g$ từ $0$ đến vô cùng. Kết quả giải số tiệm cận chính xác nghiệm giải tích TG (đạt tương tác rất mạnh tại $g=1$ với 2 giếng và $g=2.5$ với 3 giếng). Tác giả chọn $g=3$ làm giá trị đại diện cho giới hạn TG trong tính toán số. Đường phân bố mật độ giải số tại $g=3$ trùng khớp hoàn toàn với đường giải tích TG.
  • Hệ 2 thành phần trong mạng 3 giếng thế ($p=3$):
    • Trường hợp $g_A=3, g_B=0, g_{AB}=0$: Thành phần B ở trạng thái BEC tập trung tại giếng giữa; thành phần A ở trạng thái TG. Khi $V_0=0.01$, hiệu ứng xuyên ngầm lượng tử mạnh; khi tăng lên $V_0=0.5$, hiệu ứng xuyên ngầm bị triệt tiêu, các hạt tập trung tại đáy giếng. Do tỉ số lấp đầy $N_A/p = 2/3$ không nguyên, quá trình chuyển sang pha cách điện Mott bị cản trở, dẫn đến sự tồn tại của phần phi định xứ (các đỉnh ngoài đường chéo của RSPDM).
    • Trường hợp $g_A=3, g_B=0, g_{AB}=3$: Xảy ra hiện tượng tách pha (PS). Hai boson loại B định xứ ở giếng trung tâm (duy trì trạng thái BEC với $\lambda_1^B/N_B = 1$), trong khi hai boson loại A bị đẩy sang định xứ đối xứng ở hai giếng ngoài cùng (trạng thái TG cục bộ với $\lambda_1^A/N_A = \lambda_2^A/N_A = 0.5$). Quá trình chuyển pha diễn ra nhanh trong khoảng $g_{AB} \in [0; 0.1]$ và đạt trạng thái ổn định khi $g_{AB} \ge 0.5$.
  • Hệ 2 thành phần trong mạng 2 giếng thế ($p=2$):
    • Trường hợp $g_A=3, g_B=0, g_{AB}=0$: Vì $N_A/p = 2/2 = 1$ là số nguyên, thành phần A chuyển hoàn toàn sang pha cách điện Mott (mỗi giếng chứa đúng 1 hạt định xứ, triệt tiêu đường chéo chính của PDF và không có đỉnh ngoài đường chéo của RSPDM).
    • Trường hợp $g_A=3, g_B=0, g_{AB}=3$: Hiện tượng tách pha diễn ra trong từng giếng thế: tại mỗi giếng, boson loại B định xứ tại chính giữa đáy giếng, còn boson loại A định xứ đối xứng ở hai bên rìa của giếng đó. Do tính đối xứng chẵn của 2 giếng, số chiếm đóng tự nhiên của thành phần B chia đều cho hai quỹ đạo tự nhiên đầu tiên ($\lambda_1^B/N_B = \lambda_2^B/N_B = 0.5$). Khi quét $g_{AB}$ tại $V_0=0.1$, sự chuyển pha diễn ra mạnh quanh vùng $g_{AB} \approx 0.5$.

Chương 4. Kết luận và hướng phát triển

Chương 4 đúc kết các phát hiện chính về ảnh hưởng của số lượng giếng thế và cường độ tương tác đến trạng thái dừng của hệ, đồng thời vạch ra các định hướng nghiên cứu mở rộng trong tương lai.

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

Khóa luận đã định lượng và so sánh chi tiết các trạng thái dừng của hệ boson hai thành phần ($2+2$ hạt) thông qua bảng tổng hợp kết quả:

Thông số khảo sát Mạng 3 giếng thế ($p=3$) Mạng 2 giếng thế ($p=2$)
Cấu hình không tương tác liên thành phần ($g_A=3, g_B=0, g_{AB}=0$) - Thành phần B ở trạng thái BEC (đỉnh cực đại tại giếng giữa).
- Thành phần A ở trạng thái TG, xuất hiện đỉnh phi định xứ ngoài đường chéo RSPDM do $N_A/p = 2/3$ không nguyên.		 - Thành phần B ở trạng thái BEC.
- Thành phần A chuyển sang pha cách điện Mott hoàn chỉnh do $N_A/p = 2/2 = 1$ nguyên (mỗi giếng 1 hạt, PDF đường chéo triệt tiêu).
Cấu hình tách pha mạnh ($g_A=3, g_B=0, g_{AB}=3$) - Thành phần B định xứ hoàn toàn tại giếng giữa ($\lambda_1^B/N_B = 1.0$).
- Thành phần A định xứ ở 2 giếng ngoài ($\lambda_1^A/N_A = \lambda_2^A/N_A = 0.5$).		 - Tách pha nội giếng: Hạt B ở tâm mỗi giếng, hạt A ở hai rìa của từng giếng.
- Số chiếm đóng thành phần B chia đều: $\lambda_1^B/N_B = \lambda_2^B/N_B = 0.5$.
Động học chuyển pha theo $g_{AB}$ Chuyển pha rất nhanh trong khoảng $g_{AB} \in [0; 0.1]$, đạt bão hòa ổn định từ $g_{AB} \ge 0.5$. Chuyển pha diễn ra quanh vùng lân cận $g_{AB} \approx 0.5$ (khi $V_0=0.1$).
Ảnh hưởng của rào thế $V_0$ Xuyên ngầm lượng tử giảm dần từ $V_0=0.01$ và biến mất hoàn toàn khi $V_0 \ge 0.4 - 0.5$. Xuyên ngầm lượng tử biến mất khi $V_0 \ge 0.4 - 0.5$; các đỉnh mật độ tách biệt rõ rệt.

Đề xuất và kiến nghị của tác giả:

  • Áp dụng phương pháp chéo hóa ma trận chính xác cải tiến với cắt ngưỡng $E_{max}$ và tương tác hiệu dụng làm công cụ giải số tiêu chuẩn cho các bài toán hệ ít hạt có tương tác.
  • Thiết lập chiều cao giếng thế ở mức $V_0 \ge 0.4 - 0.5$ và cường độ tương tác $g_{AB} \ge 0.5$ trong các bố trí thực nghiệm mạng quang học nhằm đạt được trạng thái tách pha bão hòa ổn định.

Đóng góp khoa học tự nêu:

  • Là công trình đầu tiên khảo sát chi tiết tính chất dừng tại trạng thái cơ bản của hệ boson siêu lạnh ít hạt hai thành phần trong mạng tinh thể quang học một chiều (thay vì bẫy điều hòa truyền thống).
  • Chứng minh rõ vai trò của tỉ số lấp đầy $N_\sigma/p$ trong việc hình thành pha cách điện Mott và chỉ ra sự khác biệt về cơ chế tách pha giữa mạng giếng thế chẵn ($p=2$) và mạng giếng thế lẻ ($p=3$).

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

Hạn chế của nghiên cứu:

  • Nghiên cứu giới hạn ở bài toán hệ ít hạt với số lượng hạt cố định là $N_A = 2, N_B = 2$ (tổng 4 hạt).
  • Mô hình không gian chỉ xét trên mạng tinh thể quang học một chiều với số lượng giếng thế nhỏ ($p=2$ và $p=3$).
  • Chỉ tập trung khảo sát tính chất dừng tại trạng thái cơ bản đối với các giá trị tương tác cụ thể ($g_A=3, g_B=0, g_{AB} \in {0, 3}$), chưa khảo sát toàn diện các trạng thái fermion hóa tổng hợp (CF) và fermion hóa hoàn toàn (FF).

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

  • Khảo sát các hiện tượng tương tác mạnh khác như fermion hóa tổng hợp (CF) và fermion hóa hoàn toàn (FF) cho hệ boson hai thành phần $2+2$ hạt.
  • Tiếp tục cải tiến thuật toán ED nhằm nâng cao tốc độ và tối ưu hóa dung lượng tính toán, mở rộng khảo sát cho các hệ có số hạt lớn hơn ($N_\sigma = 3, 4...$ hạt).
  • Nghiên cứu các quá trình động lực học phi cân bằng phụ thuộc thời gian và mở rộng bài toán sang không gian hai chiều (2D) và ba chiều (3D).

Giá trị tham khảo

Khóa luận là tài liệu tham khảo cho sinh viên, học viên cao học và nghiên cứu sinh thuộc các chuyên ngành Vật lý lý thuyết, Vật lý tính toán, Quang lượng tử và Vật lý chất rắn:

  • Tham khảo về phương pháp tính toán: Quy trình tích hợp phương pháp biểu diễn biến số rời rạc (DVR) để giải TISE đơn hạt kết hợp với thuật toán chéo hóa ma trận chính xác (ED) cải tiến (sử dụng cơ sở Fock chẵn, tương tác hiệu dụng và ngưỡng năng lượng $E_{max}$) tại Chương 2 là tài liệu hữu ích cho việc xây dựng các chương trình mô phỏng lượng tử nhiều hạt.
  • Tham khảo về phân tích vật lý: Phương pháp sử dụng các đại lượng tương quan lượng tử (RSPDM, PDF, số chiếm đóng tự nhiên) tại Chương 3 cung cấp cơ sở đối sánh cho các nghiên cứu về chuyển pha lượng tử, hiện tượng tách pha và hiệu ứng giam hãm trong bẫy quang học.

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

1. Vì sao tác giả lựa chọn giá trị cường độ tương tác $g=3$ để mô tả giới hạn Tonks - Girardeau trong tính toán số?

Trong lý thuyết giải tích, trạng thái Tonks - Girardeau (TG) đạt được khi tương tác tiến tới vô cùng ($g \to \infty$). Tuy nhiên, kết quả giải số ED tại Chương 3 cho thấy năng lượng trạng thái cơ bản của hệ 1 thành phần tăng nhanh và tiệm cận mức năng lượng giải tích TG ngay từ $g=1$ (với mạng 2 giếng thế) và $g=2.5$ (với mạng 3 giếng thế). Đường phân bố mật độ tại $g=3$ trùng khớp hoàn toàn với nghiệm giải tích TG, do đó giá trị $g=3$ được chọn làm mốc tương tác rất mạnh xấp xỉ chính xác trạng thái TG trong tính toán số.

2. Sự khác biệt về trạng thái của thành phần A giữa mạng 2 giếng thế và 3 giếng thế khi không có tương tác liên thành phần ($g_A=3, g_B=0, g_{AB}=0$) là gì?

Sự khác biệt bắt nguồn từ tỉ số lấp đầy $N_A/p$:

  • Ở mạng 2 giếng thế ($p=2$), tỉ số $N_A/p = 2/2 = 1$ là số nguyên, dẫn đến sự hình thành pha cách điện Mott hoàn chỉnh (mỗi giếng thế chứa đúng 1 hạt loại A định xứ hoàn toàn, triệt tiêu đường chéo chính của hàm PDF).
  • Ở mạng 3 giếng thế ($p=3$), tỉ số $N_A/p = 2/3$ không phải số nguyên, làm cản trở quá trình chuyển pha từ siêu lỏng sang cách điện Mott, khiến mật độ hạt xuất hiện thành phần phi định xứ (thể hiện qua các đỉnh xác suất nằm ngoài đường chéo chính của RSPDM).

3. Hiện tượng tách pha trong mạng 3 giếng thế khác với mạng 2 giếng thế như thế nào khi $g_A=3, g_B=0, g_{AB}=3$?

  • Trong mạng 3 giếng thế ($p=3$): Các boson loại B định xứ hoàn toàn tại giếng thế trung tâm và duy trì trạng thái BEC ($\lambda_1^B/N_B = 1$). Các boson loại A bị đẩy ra định xứ đối xứng tại hai giếng ngoài cùng ở trạng thái TG cục bộ ($\lambda_1^A/N_A = \lambda_2^A/N_A = 0.5$).
  • Trong mạng 2 giếng thế ($p=2$): Hiện tượng tách pha diễn ra ngay trong từng giếng thế (tách pha nội giếng). Hạt loại B tập trung tại tâm của mỗi giếng, còn hạt loại A phân bố đối xứng ở hai rìa của giếng đó. Do tính đối xứng của 2 giếng, số chiếm đóng của thành phần B chia đều cho hai quỹ đạo tự nhiên đầu tiên ($\lambda_1^B/N_B = \lambda_2^B/N_B = 0.5$).

4. Chiều cao rào thế $V_0$ ảnh hưởng như thế nào đến hiệu ứng xuyên ngầm lượng tử trong mạng tinh thể quang học?

Khi chiều cao giếng thế nhỏ ($V_0 = 0.01$), hiệu ứng xuyên ngầm lượng tử diễn ra mạnh, mật độ xác suất tìm thấy hạt phân bố liên tục giữa các giếng thế và đạt cực đại ở vùng trung tâm. Khi tăng dần $V_0$, hiệu ứng xuyên ngầm suy yếu dần. Khi $V_0 \ge 0.4 - 0.5$, hiệu ứng xuyên ngầm bị triệt tiêu hoàn toàn, mật độ xác suất nối liền giữa các giếng thế biến mất và các hạt định xứ hoàn toàn tại các vị trí đáy giếng thế.

5. Thuật toán chéo hóa ma trận chính xác (ED) cải tiến giải quyết bài toán bùng nổ số chiều không gian Hilbert bằng cách nào?

Thuật toán cải tiến (Trần Dương Anh Tài và cộng sự, 2023) tận dụng tính đối xứng hoán vị của hàm sóng boson để chỉ giữ lại các trạng thái Fock chẵn trong không gian cơ sở. Đồng thời, thuật toán kết hợp phương pháp tương tác hiệu dụng và áp dụng ngưỡng năng lượng cắt $E_{max}$, cho phép tính toán chính xác các trạng thái năng lượng thấp mà không cần mở rộng toàn bộ kích thước không gian Hilbert $D$.

Kết luận

Khóa luận của tác giả Lê Minh Khang đã nghiên cứu thành công tính chất dừng tại trạng thái cơ bản của hệ boson siêu lạnh hai thành phần ($2+2$ hạt) trong mạng tinh thể quang học một chiều bằng phương pháp chéo hóa ma trận chính xác cải tiến. Kết quả nghiên cứu làm sáng tỏ ảnh hưởng của chiều cao rào thế $V_0$, các chế độ tương tác ($g_A, g_B, g_{AB}$) và số lượng giếng thế ($p=2, 3$) lên cấu trúc phân bố mật độ, ma trận RSPDM, hàm phân bố cặp PDF và số chiếm đóng tự nhiên. Công trình chứng minh rõ nét sự hình thành pha cách điện Mott phụ thuộc vào hệ số lấp đầy nguyên và cơ chế tách pha không gian đặc thù theo tính chẵn/lẻ của số giếng thế, tạo cơ sở cho các nghiên cứu tiếp theo về khí lượng tử ít hạt.