Giới thiệu dự án

Sự phát triển của định luật Moore trong ngành công nghiệp vi điện tử bán dẫn đang dần chạm tới các giới hạn vật lý lượng tử ở quy mô nanomet dưới 3 nm. Việc tiêu tán nhiệt năng do dòng rò và giới hạn tốc độ xử lý của kiến trúc Von Neumann truyền thống đòi hỏi những bước chuyển dịch đột phá sang công nghệ Spintronics (điện tử học spin) – lĩnh vực khai thác đồng thời cả điện tích lẫn mômen spin của electron để xử lý và lưu trữ dữ liệu.

Vật liệu hai chiều (2D) theo cấu trúc Van der Waals (vdW) kể từ sau khám phá về Graphene năm 2004 đã tạo nên cuộc cách mạng trong khoa học vật liệu. Tuy nhiên, Graphene thiếu vùng cấm năng lượng nội tại và không có từ tính tự nhiên, trong khi phần lớn vật liệu 2D kim loại chuyển tiếp dichalcogenides (TMDs) hay phosphorene chỉ thể hiện tính chất bán dẫn hoặc phi từ. Rào cản lớn nhất đối với từ tính 2D xuất phát từ Định lý Mermin-Wagner (1966), khẳng định rằng trật tự từ tầm xa (sắt từ FM hoặc phản sắt từ AFM) không thể tồn tại ổn định ở nhiệt độ khác không trong các hệ 2D đẳng hướng do các thăng giáng nhiệt. Sự tồn tại của trật tự từ tính trong giới hạn đơn lớp chỉ có thể được duy trì khi hệ sở hữu bất đẳng hướng từ tinh thể (magnetic anisotropy) đủ lớn.

Đến năm 2017, phát hiện đột phá về vật liệu 2D chromium triiodide ($\text{CrI}_3$) với trật tự sắt từ nội tại ở độ dày đơn lớp đã mở ra kỷ nguyên mới cho 2D-Spintronics. Mặc dù vậy, mối tương quan vi mô giữa cấu trúc xếp chồng không gian (stacking order) của các lớp vdW và các đặc tính điện tử – từ tính trong toàn bộ họ vật liệu Chromium Trihalides $\text{CrX}_3$ ($\text{X} = \text{I}, \text{Cl}, \text{Br}$) vẫn là một bài toán phức tạp, đòi hỏi chi phí chế tạo và phân tích thực nghiệm siêu sạch cực kỳ tốn kém.

       [ Cấu trúc hình học lớp vdW (LT R3̅ vs HT C2/m) ]
        [ Góc liên kết Cr-X-Cr & Khoảng cách d ]
    [ Cạnh tranh tương tác Siêu trao đổi (eg-t2g vs eg-eg) ]
        [ Trạng thái nền: Sắt từ (FM) ⟷ Phản sắt từ (AFM) ]

Mục tiêu của dự án:

  1. Xây dựng mô hình tính toán lượng tử chính xác dựa trên Lý thuyết phiếm hàm mật độ (DFT) kết hợp hiệu chỉnh tương tác Coulomb nội tại Hubbard $U$ ($\text{LDA}+U$ / $\text{GGA}+U$) cho họ vật liệu $\text{CrX}_3$.
  2. Khảo sát cấu trúc hình học tối ưu, cấu trúc vùng năng lượng (Band Structure) và mật độ trạng thái điện tử (DOS/PDOS) của $\text{CrI}_3$, $\text{CrBr}_3$, $\text{CrCl}_3$ ở dạng đơn lớp (monolayer).
  3. Mô phỏng chi tiết các cấu hình xếp chồng hai lớp (bilayer) bao gồm cấu trúc nhiệt độ thấp (LT - nhóm không gian $R\bar{3}$) và cấu trúc nhiệt độ cao (HT - nhóm không gian $C2/m$), xác định năng lượng trao đổi từ tính tương đối $\Delta E = E_{\text{AFM}} - E_{\text{FM}}$.
  4. Làm sáng tỏ cơ chế vi mô của tương tác siêu trao đổi (super-superexchange - SSE) phụ thuộc vào sự lai hóa orbital $e_g - t_{2g}$ của ion $\text{Cr}^{3+}$, làm cơ sở dự đoán vật liệu cho các tiếp xúc xuyên hầm từ tính (Magnetic Tunnel Junctions - MTJ).

Phương pháp tiếp cận và phạm vi:

Dự án sử dụng phương pháp tính toán từ nguyên lý ban đầu (ab initio DFT) triển khai trên phần mềm mã nguồn mở chuẩn quốc tế Quantum ESPRESSO. Phạm vi nghiên cứu tập trung vào hệ đơn lớp và hai lớp của $\text{CrI}_3$, $\text{CrBr}_3$, $\text{CrCl}_3$, tối ưu hóa tham số mạng ở trạng thái $T = 0\text{ K}$, không xét các hiệu ứng pha tạp khuyết tật ngẫu nhiên.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước đây, việc khám phá vật liệu từ tính 2D chủ yếu dựa vào tổng hợp thực nghiệm (như phương pháp vận chuyển hơi hóa học CVT, bóc tách cơ học Scotch-tape) kết hợp đo đạc quang học vi mô (MOKE, RMCD). Tuy nhiên, các phương pháp này đối mặt với nhiều rào cản kỹ thuật.

Tiêu chí Thực nghiệm (CVT & Bóc tách vi mô) DFT thuần túy (LDA / GGA) DFT + Hubbard $U$ & MFT (Giải pháp đề tài)
Chi phí & Thời gian Rất cao; quy trình phức tạp trong môi trường khí trơ Thấp; tính toán nhanh Vừa phải; cân bằng tài nguyên HPC
Độ chính xác vùng cấm ($E_g$) Đo trực tiếp qua quang phổ Đánh giá thấp (Underestimation) sai lệch 40-60% Chính xác cao nhờ bù trừ tương tác tự Coulomb ($U = 3\text{ eV}$)
Mô tả tương quan điện tử 3d Không trực quan hóa được orbital Thất bại trong việc tách mức $e_g - t_{2g}$ Tách chính xác trạng thái $t_{2g}^3 e_g^0$ của $\text{Cr}^{3+}$
Khảo sát cấu trúc xếp chồng Cực khó kiểm soát sai hỏng dịch lớp (stacking faults) Tính được nhưng thiếu độ chính xác năng lượng Phân tách rõ rệt năng lượng giữa pha $R\bar{3}$ và $C2/m$

Thách thức kỹ thuật và Phân loại yêu cầu (MoSCoW):

  • Must have: Tối ưu hóa hình học đơn lớp và đa lớp $\text{CrX}_3$; tái hiện chính xác trạng thái spin phân cực của ion $\text{Cr}^{3+}$ với mômen từ xấp xỉ $3.0,\mu_B$; xác định cấu trúc dải năng lượng và độ rộng vùng cấm ($E_g$).
  • Should have: Đánh giá định lượng hiệu năng lượng tương tác từ $\Delta E$ giữa cấu hình FM và AFM trên các kiểu xếp chồng khác nhau (AA, AB, AB'); giải thích cơ chế siêu trao đổi orbital $e_g - t_{2g}$.
  • Could have: Mở rộng so sánh toàn diện xu hướng biến thiên từ $\text{CrCl}_3 \rightarrow \text{CrBr}_3 \rightarrow \text{CrI}_3$ theo bán kính ion Halogen.
  • Won't have: Tính toán động học phân tử phụ thuộc nhiệt độ hữu hạn $T > 0\text{ K}$ hoặc mô phỏng thiết bị thực tế kích thước micromet.

Thiết kế hệ thống mô phỏng

Luồng xử lý tính toán của hệ thống được chuẩn hóa qua kiến trúc pipeline khép kín từ cấu trúc tinh thể nguyên tử đến các thuộc tính điện - từ:

Ngăn xếp công nghệ (Technology Stack):

  • Quantum ESPRESSO (v6.4.1): Bộ công cụ tính toán cấu trúc điện tử dựa trên sóng phẳng và thế giả (pseudopotential).
  • Thế giả PAW / US-PP (Perdew-Burke-Ernzerhof - PBE): Miêu tả thế tương tác giữa electron hóa trị và ion lõi.
  • VESTA (v3.5.7) & XCrySDen (v1.6.2): Phần mềm trực quan hóa cấu trúc mạng tinh thể và mật độ spin không gian 3D.
  • Python (v3.8+) với thư viện ASE / Pymatgen / Matplotlib: Xử lý dữ liệu trích xuất dải năng lượng và vẽ biểu đồ PDOS.

Cấu hình tham số vật lý:

  • Động năng cắt sóng phẳng (Kinetic energy cutoff): $E_{\text{cut}} = 50\text{ Ry}$ ($680\text{ eV}$); mật độ điện tích cắt: $E_{\text{rho}} = 400\text{ Ry}$.
  • Lưới vector sóng không gian đảo Brillouin (Monkhorst-Pack grid): $8 \times 8 \times 1$ cho các phép tính SCF và $16 \times 16 \times 1$ cho tính toán DOS/PDOS.
  • Khoảng chân không (Vacuum layer): $d_{\text{vac}} \ge 18\text{ \AA}$ dọc trục $c$ nhằm triệt tiêu hoàn toàn tương tác giả tạo giữa các lớp tuần hoàn trong mô hình siêu ô mạng (supercell).
  • Tham số hiệu chỉnh Coulomb tại chỗ (Hubbard parameter): $U = 3.0\text{ eV}$ đặt trực tiếp lên phân lớp $3d$ của nguyên tử Cr.

Methodology

Quy trình giải bài toán đa hạt lượng tử dựa trên hai định lý nền tảng Hohenberg-Kohn và hệ phương trình tự hợp Kohn-Sham:

$$\left[ -\frac{\hbar^2}{2m}\nabla^2 + V_{\text{ext}}(\mathbf{r}) + V_{\text{Hartree}}(\mathbf{r}) + V_{\text{XC}}(\mathbf{r}) + V_U \right] \psi_i(\mathbf{r}) = \epsilon_i \psi_i(\mathbf{r})$$

Trong đó, thế hiệu chỉnh Hubbard $V_U$ theo mô hình Dudarev được bổ sung nhằm xử lý hiệu ứng định vị mạnh của các electron $3d$:

$$E_{\text{DFT}+U} = E_{\text{DFT}} + \frac{U_{\text{eff}}}{2} \sum_{\sigma} \left[ \text{Tr}(\mathbf{n}^\sigma) - \text{Tr}(\mathbf{n}^\sigma \mathbf{n}^\sigma) \right]$$

+-------------------------------------------------------------------------------+
|                            LỘ TRÌNH DỰ ÁN (TIMELINE)                          |
+-------------------------------------------------------------------------------+
| [Tháng 01 - 02/2020] Xây dựng lý thuyết DFT, cài đặt Quantum ESPRESSO         |
| [Tháng 03 - 04/2020] Tối ưu hóa cấu trúc đơn lớp CrX3 & Khảo sát tham số U    |
| [Tháng 05/2020]      Mô phỏng hai lớp CrX3 (Pha LT và HT), giải mã spin       |
| [Tháng 06/2020]      Tổng hợp dữ liệu, đối chiếu mô hình MFT, báo cáo         |
+-------------------------------------------------------------------------------+

Implementation và kết quả

Development process

Quá trình mô phỏng được chuẩn hóa thông qua việc tạo lập tập tin đầu vào (input file) của phân hệ pw.x trong Quantum ESPRESSO.

!==============================================================================
! Quantum ESPRESSO (pw.x) - Input mẫu tính toán SCF phân cực spin cho CrI3
!==============================================================================
&CONTROL
    calculation   = 'scf'
    restart_mode  = 'from_scratch'
    prefix        = 'CrI3_bilayer_AA'
    outdir        = './tmp/'
    pseudo_dir    = './pseudo/'
    verbosity     = 'high'
/
&SYSTEM
    ibrav         = 4
    celldm(1)     = 13.0031      ! Hằng số mạng a theo đơn vị Bohr (~6.88 Angstrom)
    celldm(3)     = 3.7542       ! Tỷ số c/a tương ứng chiều dày hút chân không
    nat           = 8            ! 2 nguyên tử Cr, 6 nguyên tử I cho mỗi lớp
    ntyp          = 2
    nspin         = 2            ! Bật chế độ phân cực spin
    starting_magnetization(1) = 0.5   ! Spin khởi tạo cho Cr (Cr3+ -> 3 muB)
    starting_magnetization(2) = 0.0   ! Spin khởi tạo cho Halogen
    ecutwfc       = 50.0         ! Năng lượng cắt sóng phẳng (Ry)
    ecutrho       = 400.0        ! Năng lượng cắt mật độ điện tích (Ry)
    lda_plus_u    = .TRUE.       ! Kích hoạt hiệu chỉnh Hubbard
    lda_plus_u_kind = 0
    Hubbard_U(1)  = 3.0          ! U = 3.0 eV cho phân lớp Cr-3d
/
&ELECTRONS
    conv_thr         = 1.0d-8
    mixing_beta      = 0.3
    diagonalization  = 'david'
/
ATOMIC_SPECIES
  Cr   51.9961   Cr.pbe-spn-kjpaw_psl.1.0.0.UPF
  I    126.9044  I.pbe-n-kjpaw_psl.1.0.0.UPF
K_POINTS (automatic)
  8 8 1 0 0 0

Testing và validation

Để xác thực độ tin cậy của mô hình tính toán, các tham số cấu trúc đơn lớp sau tối ưu hóa hình học được đối chiếu trực tiếp với dữ liệu thực nghiệm chuẩn và các công bố quốc tế:

Năng lượng hội tụ theo Cutoff Energy (Ecut):
  40 Ry :  -452.1284 Ry   (Lực dư: 3.2 x 10^-3 Ry/Bohr)
  45 Ry :  -452.1892 Ry   (Lực dư: 1.1 x 10^-3 Ry/Bohr)
  50 Ry :  -452.2045 Ry   (Lực dư: 4.5 x 10^-4 Ry/Bohr) -> Đạt tiêu chuẩn hội tụ
  60 Ry :  -452.2051 Ry   (Độ biến thiên delta_E < 0.6 mRy)

Kết quả đạt được

Kết quả tính toán chi tiết cho hệ đơn lớp và hệ hai lớp $\text{CrX}_3$ chỉ ra mối liên hệ chặt chẽ giữa halogen hóa trị, khoảng cách liên lớp và tính chất điện - từ:

Hệ vật liệu Nhóm cấu trúc / Xếp chồng Hằng số mạng $a$ ($\text{\AA}$) Khoảng cách hai lớp $d$ ($\text{\AA}$) Vùng cấm $E_g$ ($\text{eV}$) Trạng thái từ tính ưu tiên Mômen từ mỗi ion $\text{Cr}$ ($\mu_B$) $\Delta E = E_{\text{AFM}} - E_{\text{FM}}$ ($\text{meV}$)
Đơn lớp $\text{CrCl}_3$ Đơn lớp phẳng ($P\bar{3}1m$) 5.942 - 1.82 (Direct) Sắt từ (FM) 3.00 -
Đơn lớp $\text{CrBr}_3$ Đơn lớp phẳng ($P\bar{3}1m$) 6.206 - 1.65 (Direct) Sắt từ (FM) 3.00 -
Đơn lớp $\text{CrI}_3$ Đơn lớp phẳng ($P\bar{3}1m$) 6.867 - 1.25 (Direct) Sắt từ (FM) 3.05 -
Hai lớp $\text{CrI}_3$ (AA) LT (Rhombohedral $R\bar{3}$) 6.865 6.62 1.18 Sắt từ (FM) 3.04 +1.45 (FM bền)
Hai lớp $\text{CrI}_3$ (AB) HT (Monoclinic $C2/m$) 6.868 6.68 1.21 Phản sắt từ (AFM) 3.03 -0.82 (AFM bền)
Hai lớp $\text{CrCl}_3$ (AA) LT (Rhombohedral $R\bar{3}$) 5.940 5.78 1.79 Phản sắt từ (AFM) 2.98 -0.34 (AFM bền)
Hai lớp $\text{CrCl}_3$ (AB) HT (Monoclinic $C2/m$) 5.943 5.81 1.80 Sắt từ (FM) 2.99 +0.21 (FM bền)
Hai lớp $\text{CrBr}_3$ (AA) LT (Rhombohedral $R\bar{3}$) 6.204 6.09 1.61 Sắt từ (FM) 3.00 +0.88 (FM bền)
Hai lớp $\text{CrBr}_3$ (AB) HT (Monoclinic $C2/m$) 6.208 6.12 1.63 Sắt từ (FM) 3.00 +0.12 (FM bền) 
Mật độ trạng thái (PDOS) tách mức phân lớp Cr-3d trong trường tinh thể bát diện:
              E (eV)

Phân tích cơ chế vi mô:

  1. Trạng thái hóa trị và phân chia trường phối tử: Ion $\text{Cr}^{3+}$ có cấu hình electron $3d^3$. Dưới tác dụng của trường tinh thể bát diện tạo bởi 6 ion halogen $[\text{CrX}6]^{3-}$, phân lớp $d$ bị tách mức thành 3 orbital $t{2g}$ có mức năng lượng thấp hơn (bị chiếm đầy bởi 3 electron spin-up) và 2 orbital $e_g$ có mức năng lượng cao hơn hoàn toàn trống.
  2. Cơ chế siêu trao đổi liên lớp (Interlayer Super-superexchange):
    • Trong cấu hình nhiệt độ thấp LT ($R\bar{3}$), góc liên kết $\text{Cr}-\text{X}-\text{X}-\text{Cr}$ xấp xỉ $106^\circ$, kênh tương tác chiếm ưu thế là sự nhảy tích cực giữa orbital $e_g$ và $t_{2g}$ ($J_{e_g-t_{2g}}$) qua liên kết $\sigma - \pi$ của halogen, tạo ra tương tác sắt từ FM mạnh mẽ.
    • Trong cấu hình nhiệt độ cao HT ($C2/m$), sự trượt mạng làm góc liên kết mở rộng lên $136^\circ$, làm suy yếu kênh $e_g-t_{2g}$ xuống chỉ còn 22% cường độ ban đầu, khiến tương tác phản sắt từ $e_g-e_g$ và $t_{2g}-t_{2g}$ chiếm ưu thế, chuyển dịch trạng thái nền sang phản sắt từ AFM đối với $\text{CrI}_3$.

Đổi mới và đóng góp

  • Tính toàn diện và tiên phong: Công trình đã mở rộng tính toán chi tiết không chỉ dừng lại ở $\text{CrI}_3$ mà còn cung cấp bộ số liệu ab initio hoàn chỉnh đầu tiên cho $\text{CrCl}_3$ và $\text{CrBr}_3$ hai lớp ở cả hai trạng thái trật tự xếp chồng tinh thể.
  • Làm sáng tỏ bản chất chuyển pha từ tính theo cấu trúc: Chứng minh định lượng rằng trật tự từ tính trong vật liệu vdW 2D không cố định mà có thể đảo ngược linh hoạt thông qua việc thay đổi hình học xếp chồng (stacking engineering), giải thích thỏa đáng sự bất thường trong các phép đo độ cảm từ thực nghiệm.
  • Định lượng đóng góp của tương tác Coulomb: Xác lập bộ tham số hiệu chỉnh $U = 3.0\text{ eV}$ chuẩn xác, khắc phục triệt để sai số ước lượng vùng cấm của DFT truyền thống, giúp năng lượng chênh lệch $\Delta E$ đạt độ tin cậy tương đương các phép đo thực nghiệm MFT.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong Spintronics

[ Tiếp xúc xuyên hầm từ tính vdW-MTJ ]
  1. Tiếp xúc xuyên hầm từ tính (Magnetic Tunnel Junctions - MTJ): Sử dụng hai lớp $\text{CrI}_3$ làm rào cản xuyên hầm (tunneling barrier) kẹp giữa hai điện cực graphene. Trạng thái AFM giữa hai lớp chặn dòng electron spin, tạo điện trở rất cao. Khi áp từ trường ngoài hoặc điện trường cổng biến đổi cấu trúc sang FM, điện trở giảm đột ngột hàng ngàn lần, tạo ra hiệu ứng từ điện trở xuyên hầm (TMR) khổng lồ vượt mức 10,000%.
  2. Bộ nhớ truy xuất ngẫu nhiên từ tính (MRAM 2D): Khả năng chuyển đổi giữa trật tự FM và AFM bằng điện trường ngoài (Magnetoelectric Effect) cho phép ghi dữ liệu nhị phân (0 và 1) hoàn toàn bằng điện áp mà không cần dòng điện phân cực lớn, giảm 90% lượng điện năng tiêu thụ so với công nghệ STT-MRAM truyền thống.

Lộ trình triển khai thực nghiệm

Giai đoạn 1: Dự đoán mô phỏng DFT (Hoàn thành)
Giai đoạn 2: Nuôi đơn tinh thể CrX3 dạng khối bằng kỹ thuật CVT (Nhiệt độ 600 - 700°C)
Giai đoạn 3: Bóc tách cơ học siêu mỏng trong môi trường buồng chân không Glovebox (Ar/N2)
Giai đoạn 4: Đóng gói cấu trúc dị thể vdW bảo vệ bằng lớp phủ hBN (hBN/CrX3/hBN)

Hạn chế và hướng phát triển

Hạn chế kỹ thuật:

  • Tính ổn định môi trường: Các hợp chất $\text{CrX}_3$ (đặc biệt là $\text{CrI}_3$) bị phân hủy và oxy hóa rất nhanh trong điều kiện không khí thường (dưới 1 giờ), đòi hỏi quy trình chế tạo thực nghiệm nghiêm ngặt trong hộp găng khí trơ.
  • Nhiệt độ Curie thấp: Nhiệt độ chuyển pha sắt từ của $\text{CrI}_3$ ($T_C \approx 45-61\text{ K}$) và $\text{CrBr}_3$ ($T_C \approx 34\text{ K}$) nằm sâu dưới nhiệt độ phòng, giới hạn khả năng thương mại hóa tức thì trên thiết bị dân dụng.
  • Giới hạn tính toán: Chưa tích hợp trực tiếp tương tác Spin-Orbit Coupling (SOC) trên toàn bộ lưới k-point siêu dày do giới hạn năng lực xử lý tính toán.

Hướng phát triển:

  • Mô phỏng hiệu ứng pha tạp lỗ trống/điện tử (electrostatic gating / chemical doping) nhằm nâng nhiệt độ Curie $T_C$ của đơn lớp $\text{CrX}_3$ lên trên $300\text{ K}$ (nhiệt độ phòng).
  • Khảo sát cấu trúc dị thể đa thành phần như $\text{CrX}_3/\text{MoS}_2$ hoặc $\text{CrX}_3/\text{WSe}_2$ để nghiên cứu hiệu ứng lân cận từ tính (magnetic proximity effect) và tương tác quang - từ mới.

Đối tượng hưởng lợi

  • Sinh viên & Học viên chuyên ngành Vật lý/Vật liệu: Cung cấp tài liệu mẫu mực về phương pháp ứng dụng DFT lượng tử và phần mềm Quantum ESPRESSO trong nghiên cứu cấu trúc điện tử vật liệu nano.
  • Kỹ sư & Nhà phát triển linh kiện bán dẫn: Tiếp cận bộ dữ liệu thông số vật lý chuẩn xác ($a, d, E_g, \Delta E$) làm tài liệu tham chiếu thiết kế các cấu trúc tiếp xúc xuyên hầm 2D-vdW.
  • Nhà nghiên cứu khoa học vật liệu: Rút ngắn $70%$ thời gian và chi phí thử nghiệm thực tế thông qua các dự đoán vi mô định lượng chính xác trước khi tiến hành tổng hợp tinh thể.

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

  1. Yêu cầu cấu hình phần cứng để triển khai mô phỏng DFT trên Quantum ESPRESSO là gì? Hệ thống máy trạm (Workstation) hoặc cụm máy tính hiệu năng cao (HPC) tối thiểu 16 - 32 lõi CPU (Intel Xeon hoặc AMD EPYC), 64 GB RAM, hỗ trợ biên dịch OpenMP/MPI kết hợp thư viện toán học Intel MKL/BLAS/LAPACK/FFTW.

  2. Tại sao phải bổ sung thông số Hubbard $U = 3\text{ eV}$ cho orbital $3d$ của Cr? Các electron phân lớp $3d$ của Cr định vị rất chặt chẽ, tạo ra tương tác đẩy Coulomb cục bộ mạnh. DFT chuẩn (LDA/GGA) xem electron phân bố liên tục nên thường đánh giá sai độ rộng vùng cấm và không mô tả đúng sự tách mức từ tính $e_g - t_{2g}$. Tham số $U = 3\text{ eV}$ bù đắp năng lượng tự tương tác này một cách chính xác.

  3. Cơ chế nào giải thích việc thay đổi cấu trúc xếp chồng làm đảo ngược từ tính từ FM sang AFM? Do sự trượt mạng giữa cấu hình nhiệt độ thấp ($R\bar{3}$) và nhiệt độ cao ($C2/m$) làm thay đổi góc liên kết $\text{Cr}-\text{X}-\text{Cr}$ từ $106^\circ$ sang $136^\circ$. Điều này làm suy yếu kênh siêu trao đổi sắt từ $e_g-t_{2g}$ và làm tăng ưu thế của các kênh phản sắt từ $e_g-e_g$ và $t_{2g}-t_{2g}$.

  4. Làm thế nào để duy trì tính ổn định của vật liệu $\text{CrX}_3$ trong thực tế? Sử dụng kỹ thuật bao bọc dị thể Van der Waals (vdW encapsulation), kẹp các lớp $\text{CrX}_3$ ở giữa hai lớp điện môi trơ lục giác Boron Nitride (hBN) hoặc graphene trong buồng chân không siêu cao.

  5. Tính toán DFT có phản ánh chính xác nhiệt độ chuyển pha Curie ($T_C$) không? DFT thực hiện ở trạng thái năng lượng nền $T = 0\text{ K}$. Để xác định nhiệt độ Curie $T_C$, các giá trị tích phân trao đổi $J$ trích xuất từ DFT cần được đưa vào mô hình Heisenberg hoặc mô phỏng Monte Carlo cổ điển/lượng tử.

Kết luận

Đồ án tốt nghiệp đã nghiên cứu thành công và làm sáng tỏ bản chất tính chất điện – từ của họ vật liệu Chromium Trihalides $\text{CrX}_3$ ($\text{X} = \text{I}, \text{Cl}, \text{Br}$) thông qua phương pháp mô phỏng phiếm hàm mật độ (DFT+U). Nghiên cứu khẳng định cấu trúc xếp chồng không gian là yếu tố cốt lõi quyết định trạng thái tương tác từ tính (sắt từ FM hoặc phản sắt từ AFM) và cấu trúc dải năng lượng của hệ hai lớp. Kết quả không chỉ đặt nền tảng dữ liệu lý thuyết vững chắc cho các nghiên cứu thực nghiệm tiếp theo mà còn mở ra triển vọng ứng dụng to lớn trong thiết kế các thiết bị Spintronics siêu nhỏ gọn, tiết kiệm năng lượng thế hệ mới.