Từ điển học thuật Khoa học tự nhiên

Mô hình hubbard là gì? Các nghiên cứu khoa học liên quan

Tiếng AnhHubbard model

Mô hình Hubbard là một mô hình gần đúng cơ bản trong vật lý chất rắn và vật lý lý thuyết, mô tả sự chuyển tiếp giữa trạng thái dẫn điện và cách điện (chuyển pha Mott) của các electron tương tác mạnh trên một mạng tinh thể.

298 lượt xem Cập nhật 11/9/2026

Mô hình Hubbard (Hubbard model) là thuật ngữ khoa học then chốt.

Giới thiệu mô hình Hubbard

Mô hình Hubbard là một trong những mô hình tối giản nhưng nền tảng nhất trong vật lý chất đậm đặc, được đề xuất để miêu tả tương tác giữa các electron trên mạng tinh thể. Mục đích chính của mô hình là nghiên cứu sự cạnh tranh giữa động học khuếch tán electron và tác động tương tác Coulomb mạnh tại cùng một vị trí mạng.

Mô hình này đã trở thành cơ sở lý thuyết để giải thích các hiện tượng phức tạp như pha Mott, chuyển pha kim loại-điện môi và cơ chế siêu dẫn bất thường. Tính tổng quát và khả năng mở rộng của mô hình cho phép nghiên cứu cả hệ electron spin tự do và có tương tác, đồng thời hỗ trợ phát triển các phương pháp giải gần đúng và mô phỏng số học.

  • Khởi nguồn: đề xuất bởi John Hubbard năm 1963 để giải thích từ tính của kim loại chuyển tiếp.
  • Áp dụng rộng rãi trong: chất lượng cao phân tách, vật liệu siêu dẫn, mạng quang học với nguyên tử lạnh.
  • Phát triển: mở rộng thành mô hình t–J, Extended Hubbard, và tích hợp với lý thuyết trường trung bình động (DMFT).

Phương trình Hamiltonian cơ bản

Hamiltonian của mô hình Hubbard bao gồm hai thành phần chính: điều hòa electron giữa các vị trí lân cận và tương tác Coulomb khi hai electron có spin trái ngược cùng chiếm một vị trí mạng. Phương trình tổng quát:

H=ti,j,σci,σcj,σ  +  Uini,ni, H = -t \sum_{\langle i,j\rangle,\sigma} c_{i,\sigma}^\dagger c_{j,\sigma} \;+\; U \sum_{i} n_{i,\uparrow} n_{i,\downarrow}

Trong đó:

Ký hiệuÝ nghĩa
tHệ số hopping, xác định xác suất electron di chuyển giữa hai vị trí lân cận
UNăng lượng thế Coulomb khi hai electron ngược spin chiếm cùng một site
ci,σc_{i,\sigma}^\daggerToán tử tạo electron spin σ tại vị trí i
ni,σn_{i,\sigma}Toán tử số electron spin σ tại vị trí i (n = ccc^\dagger c)

Biểu thức Hamiltonian tóm gọn bản chất tương tác và khuếch tán của hệ electron, tạo nền tảng cho các phân tích pha và tính chất từ tính của vật liệu.

Ý nghĩa vật lý của tham số t và U

Tham số t kiểm soát độ rộng băng dẫn của hệ electron: giá trị t lớn làm tăng khả năng di chuyển của electron, dẫn đến trạng thái kim loại. Ngược lại, khi t nhỏ, electron bị cản trở chuyển động, dễ hình thành pha cách ly điện môi nếu tương tác U đủ mạnh.

Tham số U mô tả tương tác đẩy Coulomb tĩnh tại cùng một site, tạo ra hiệu ứng cản trở kép electron (coulomb blockade). Khi U vượt ngưỡng U_c so với băng dẫn W~zt (z: số lân cận), hệ có thể chuyển sang pha Mott:

  • Pha kim loại (U/W ≪ 1): electron phân bố rộng, dẫn điện tốt.
  • Pha Mott (U/W ≫ 1): electron bị khóa tại site, dẫn điện kém, tạo điện môi tương tác mạnh.

Sự cạnh tranh giữa t và U cũng quyết định mức độ từ tương ứng: từ tính ferromagnet khi U lớn tại bán đầy, đến antiferromagnet gỗ-số khi U trung bình và độ lấp đầy n=1.

Giải pháp lý thuyết: phương pháp gần đúng

Giải chính xác mô hình Hubbard trong chiều lớn hơn 1 là bài toán phi tuyến, vì vậy nhiều phương pháp gần đúng đã được phát triển:

  • Phương pháp trường trung bình (Mean–Field): Tách tương tác U theo giá trị trung bình ⟨n⟩, đơn giản hóa Hamiltonian thành dạng không tương tác hiệu dụng.
  • DMFT (Dynamical Mean-Field Theory): Áp dụng lý thuyết trường trung bình động, thay thế lattice bằng một impurity model gắn kết với bồn tắm điện tử tự điều chỉnh (link.aps.org).
  • Phương pháp biến phân (Variational): Sử dụng hàm sóng Gutzwiller hoặc Jastrow để ước lượng trạng thái nền, tối ưu hóa tham số biến phân.
  • Mô phỏng Monte Carlo: QMC (Quantum Monte Carlo) cho lattice hữu hạn, tính các đại lượng nhiệt động học và phổ excitations.

Bảng tóm tắt hiệu quả của từng phương pháp:

Phương phápƯu điểmHạn chế
Mean-Field Đơn giản, ít tốn tài nguyên Bỏ qua dao động lượng tử, không thích hợp pha Mott mạnh
DMFT Bao gồm dao động thời gian, mô tả tốt pha Mott Không tính được tương quan không gian, tốn kém tính toán
Variational Linh hoạt hàm sóng, hiệu quả cho trạng thái nền Phụ thuộc hàm sóng chọn, khó mở rộng kích thước lớn
QMC Kết quả gần chính xác cho lattice nhỏ Vấn đề dấu trừ (sign problem), hạn chế kích thước

Các phương pháp này thường được kết hợp hoặc so sánh để đạt được cái nhìn toàn diện về pha và tính chất của mô hình Hubbard.

Giải pháp chính xác trên mạng một chiều (Bethe ansatz)

Mô hình Hubbard một chiều có thể giải chính xác bằng phương pháp Bethe ansatz, như John Sutherland và Elliott Lieb – F. Y. Wu đã phát triển. Phương pháp này biến bài toán tương tác mạnh thành hệ phương trình tích phân cho các moment lượng tử (quasimomenta) của electron, cho phép xác định phổ năng lượng và trạng thái nền.

Phương trình Bethe ansatz cho hệ N electron trên L site với mật độ n = N/L thường được viết dưới dạng:

eikjL=ljsinkjsinkl+iU/2tsinkjsinkliU/2t,j=1,,N e^{ik_jL} = \prod_{l\neq j} \frac{\sin k_j - \sin k_l + iU/2t}{\sin k_j - \sin k_l - iU/2t}, \quad j=1,\dots,N

Giải hệ này cho phép tính chính xác năng lượng nền và các phổ kích thích (spinon, holon), làm rõ hiện tượng tách spin – charge đặc trưng cho một chiều .

Pha Mott và pha kim loại

Pha Mott xuất hiện khi tương tác Coulomb U lớn hơn băng dẫn W (~4t ở một chiều, ~~zt với z là số lân cận). Ở n=1 (bán đầy), khi U/W vượt ngưỡng U_c/W ∼1, hệ đóng mở khe năng lượng, trở thành điện môi Mott với khe Mott Δ ≈ U−W.

Tỷ số U/WTrạng tháiĐặc điểm
<0.5Kim loạiĐiện dẫn tốt, phổ liên tục
≈1Cận phaDao động mạnh, khởi phát Mott
>1.5Mott điện môiKhe năng lượng mở, tính cách ly

Chuyển pha kim loại–Mott quan sát được qua đo độ dẫn điện, quang phổ Raman hoặc quang phổ quang điện (ARPES) trên vật liệu như V2O3 và các oxit chuyển tiếp .

Mở rộng mô hình: t–J và Extended Hubbard

Từ giới hạn U≫t, mô hình Hubbard có thể quy về mô hình t–J qua khai triển perturbation bậc hai trong t/U, thu được Hamiltonian:

Ht-J=ti,j,σc~i,σc~j,σ+Ji,j(SiSj14ninj),J=4t2U H_{t\text{-}J} = -t \sum_{\langle i,j\rangle,\sigma} \tilde{c}_{i,\sigma}^\dagger \tilde{c}_{j,\sigma} + J\sum_{\langle i,j\rangle} \Bigl(\mathbf{S}_i\cdot\mathbf{S}_j - \frac{1}{4}n_i n_j\Bigr), \quad J=\frac{4t^2}{U}

Extended Hubbard bổ sung các tương tác ngoại-site V giữa electron lân cận, mang lại khả năng nghiên cứu pha điện mật độ (charge-density wave) và các pha trật tự dài hạn khác:

  • Thành phần V: Vi,jninjV\sum_{\langle i,j\rangle}n_i n_j
  • Pha CDW, SDW (spin-density wave) xuất hiện khi V cạnh tranh với U và t .

Ứng dụng trong vật lý chất đậm đặc và siêu dẫn

Mô hình Hubbard và t–J là nền tảng lý thuyết cho cơ chế siêu dẫn bất thường ở các chất động điện tử dày như cuprate (La2–xSrxCuO4) và sắt pnictide. Các tính toán DMFT, QMC và DCA (Dynamical Cluster Approximation) đã mô phỏng được pha siêu dẫn dx²–y² và gợi ý vai trò của tương tác điện tử trong cặp electron (pairing) .

Trong vật liệu nặng fermion (heavy fermion), phiên bản mở rộng của mô hình Hubbard kết hợp bồn tắm Kondo (Periodic Anderson Model) mô tả hòa trộn giữa electron địa phương và electron dẫn truyền, giải thích hiệu ứng lượng tử gần điểm tới hạn (quantum criticality).

  • Thí nghiệm quang phổ ARPES và neutron scattering kiểm chứng các dự đoán lý thuyết.
  • Tính chất nhiệt động: hệ số nhiệt, đối lưu Pauli biểu hiện tương tác mạnh.

Thực nghiệm và mô phỏng trên mạng quang học với nguyên tử lạnh

Với nguyên tử lạnh, hệ Rubidium hoặc Lithium trong bẫy quang học (optical lattice) tái tạo chính xác Hamiltonian Hubbard với khả năng điều chỉnh t, U độc lập. Thí nghiệm của I. Bloch và cộng sự đã quan sát chuyển pha kim loại–Mott qua đo mật độ lấp đầy và dao động Bragg .

Ưu điểm của hệ nguyên tử lạnh:

  1. Điều khiển tham số động học và tương tác với độ chính xác cao.
  2. Đọc tín hiệu không phá hủy qua ảnh học fluorescence và quang phổ gióng sóng (noise correlations).

Những tiến bộ gần đây bao gồm quan sát tách spin–charge và thực nghiệm Bethe ansatz bước sóng:các site 1D bằng kĩ thuật quét quang học phân giải cao.

Danh mục tài liệu tham khảo

  1. Lieb E. H., Wu F. Y. “Absence of Mott transition in an exact solution of the short-range, one-band model in one dimension.” Phys. Rev. Lett. 1968;20(25):1445–1448. link.aps.org
  2. Imada M., Fujimori A., Tokura Y. “Metal-insulator transitions.” Rev. Mod. Phys. 1998;70(4):1039–1263. link.aps.org
  3. Dagotto E. “Correlated electrons in high-temperature superconductors.” Rev. Mod. Phys. 1994;66(3):763–840. link.aps.org
  4. Gutzwiller M. C. “Effect of correlation on the ferromagnetism of transition metals.” Phys. Rev. Lett. 1963;10(5):159–162. link.aps.org
  5. Georges A., Kotliar G., Krauth W., Rozenberg M. J. “Dynamical mean‐field theory of strongly correlated fermion systems and the limit of infinite dimensions.” Rev. Mod. Phys. 1996;68(1):13–125. link.aps.org
  6. Bloch I., Dalibard J., Zwerger W. “Many-body physics with ultracold gases.” Rev. Mod. Phys. 2008;80(3):885–964. link.aps.org
  7. Jaksch D., Zoller P. “The cold atom Hubbard toolbox.” Ann. Phys. 2005;315(1):52–79. sciencedirect.com
  8. Essler F. H. L., Frahm H., Göhmann F., Klümper A., Korepin V. E. The One-Dimensional Hubbard Model. Cambridge University Press; 2005.
  9. Scalapino D. J. “A common thread: The pairing interaction for unconventional superconductors.” Rev. Mod. Phys. 2012;84(4):1383–1417. link.aps.org
  10. Schauss P., et al. “Observation of spatially ordered structures in a two-dimensional Rydberg gas.” Nature. 2015;491:87–91. nature.com

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

Hai tham số then chốt t (thông số nhảy) và U (năng lượng đẩy Coulomb tại chỗ) trong Hamilton Hubbard biểu thị ý nghĩa vật lý gì?

Tham số t mô tả động năng và xác suất electron nhảy giữa các nút mạng lân cận (thúc đẩy tính chất dẫn điện kim loại); còn tham số U mô tả năng lượng đẩy tĩnh điện khi hai electron trái spin cùng chiếm giữ một nút mạng (thúc đẩy tính chất định vị cách điện).

Cơ chế chuyển pha cách điện Mott (Mott insulator) phát sinh từ tương tác electron-electron được giải thích thế nào qua mô hình Hubbard?

Khi tỷ số U/t vượt quá một giá trị tới hạn tại mật độ lấp đầy một nửa (half-filling), lực đẩy Coulomb cản trở sự chuyển động của electron, dải dẫn bị tách thành dải Hubbard dưới và trên cách nhau một khe năng lượng, biến vật liệu thành chất cách điện dù dải năng lượng chưa lấp đầy.

Mô hình Hubbard đóng vai trò gì trong việc giải thích cơ chế vi mô của hiện tượng siêu dẫn nhiệt độ cao (cuprate)?

Mô hình Hubbard 2D mô tả chính xác mặt phẳng đồng-oxy (CuO2) trong chất siêu dẫn cuprate; khi pha tạp lỗ trống (hole doping) vào chất cách điện Mott phản sắt từ, sự cạnh tranh giữa động năng và tương tác từ dẫn đến sự bắt cặp electron siêu dẫn d-wave.

Tài liệu tham khảo

  1. Kawabata (1981). Electron Correlation in the Small Hubbard Model. Springer Series in Solid-State Sciences Electron Correlation and Magnetism in Narrow-Band Systems. DOI: 10.1007/978-3-642-81639-0_19
  2. Ovchinnikov, Shneider (2005). Effective Hamiltonian for HTSC cuprates taking into account electron-phonon interaction in the strong-correlation regime. Journal of Experimental and Theoretical Physics. DOI: 10.1134/1.2149064
  3. Chen, Wang (1998). Mapping the electron correlation onto a model Hamiltonian for Cs/GaAs(110): a Mott-Hubbard insulator at quarter filling. Journal of Physics: Condensed Matter. DOI: 10.1088/0953-8984/10/4/003