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

Trong lĩnh vực khoa học và kỹ thuật vật liệu hiện đại, vật lý tính toán đóng vai trò đặc biệt quan trọng. Phương pháp tính toán mô phỏng lượng tử dựa trên lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) không chỉ mang tính tiên phong trong việc dự đoán các cấu trúc tinh thể và tính chất vật lý mới, mà còn giúp giải quyết bài toán hạn chế về kinh phí thực nghiệm đối với các hệ vật liệu thấp chiều, đặc biệt là vật liệu đơn lớp (2D monolayer), nơi trang thiết bị chế tạo và đo đạc đòi hỏi chi phí đầu tư rất lớn.

Sau sự phát hiện về các đặc tính vận chuyển của graphene, một loạt vật liệu hai chiều khác đã được nghiên cứu như phosphorene, silicene hay các dichalcogenide kim loại chuyển tiếp (TMDs). Tuy nhiên, phần lớn các tinh thể hai chiều này đều thiếu trật tự từ tính nội tại. Sự khám phá về từ tính trong các vật liệu hai chiều như đơn lớp $\text{CrI}_3$ (nhiệt độ Curie $T_C = 45\text{ K}$), $\text{Cr}_2\text{Ge}_2\text{Te}_6$, $\text{Fe}_3\text{GeTe}_2$ hay $\text{MnSe}_2$ đã kích thích mạnh mẽ hướng nghiên cứu vật liệu 2D từ tính. Trong đó, họ vật liệu chalcogenide phốt pho kim loại ($\text{MPX}_n$), đặc biệt là trichalcogenide phốt pho kim loại chuyển tiếp ($\text{MPX}_3$ với $\text{M}$ là kim loại chuyển tiếp và $\text{X} = \text{S}, \text{Se}$), sở hữu cấu trúc mạng tinh thể dạng lớp liên kết van der Waals yếu, cho phép bóc tách thành các lớp đơn. Trong số các hợp chất này, đơn lớp $\text{CrPSe}_3$ nhận được sự chú ý lớn nhờ đặc tính từ tính nội tại cùng tiềm năng ứng dụng trong các thiết bị quang điện tử, spintronics và xúc tác chuyển đổi năng lượng. Khóa luận tốt nghiệp của tác giả Võ Thị Sanh tập trung giải quyết bài toán khảo sát cấu trúc, tính chất điện và tính chất từ của đơn lớp $\text{CrPSe}_3$, đồng thời tìm hiểu sự thay đổi các tính chất này khi pha tạp nguyên tố sắt (Fe).

Mục tiêu và nhiệm vụ nghiên cứu của đề tài được xác định cụ thể như sau:

  1. Nghiên cứu tổng quan cơ sở lý thuyết phiếm hàm mật độ (DFT), các xấp xỉ tương quan trao đổi, kỹ thuật giả thế và phương pháp hiệu chỉnh Hubbard ($U$).
  2. Tổng hợp các dữ liệu thực nghiệm và lý thuyết đã công bố về cấu trúc hình học, tính chất từ tính và cấu trúc điện tử của vật liệu hai chiều $\text{CrPSe}_3$.
  3. Thiết lập mô hình mô phỏng, tối ưu hóa cấu trúc và tính toán năng lượng toàn phần của các cấu hình từ tính (FM, nAFM, zAFM, sAFM, NM) đối với đơn lớp $\text{CrPSe}_3$ nguyên bản bằng phần mềm Quantum ESPRESSO.
  4. Xây dựng mô hình siêu ô mạng của đơn lớp $\text{CrPSe}_3$ pha tạp $25%\text{ Fe}$ thay thế vị trí Cr; khảo sát sự biến đổi về độ dài liên kết, góc liên kết, cấu trúc vùng năng lượng, mật độ trạng thái (DOS/PDOS), điện tích Bader và hàm định xứ electron (ELF).

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

  • Đối tượng nghiên cứu: Cấu trúc hình học, cấu trúc vùng năng lượng điện tử và trật tự từ tính của vật liệu hai chiều đơn lớp $\text{CrPSe}_3$ thuần túy và hệ $\text{CrPSe}_3$ pha tạp $25%\text{ Fe}$.
  • Phạm vi nghiên cứu: Nghiên cứu được thực hiện hoàn toàn bằng phương pháp mô phỏng ab-initio dựa trên DFT với điều kiện biên tuần hoàn 2D (mạng tinh thể monoclinic, bổ sung lớp chân không $22\text{ \AA}$ theo phương trục $z$ để triệt tiêu tương tác giữa các lớp), không tiến hành chế tạo thực nghiệm trực tiếp mà sử dụng số liệu thực nghiệm từ các công bố quốc tế để đối chiếu.

Cơ sở lý thuyết và phương pháp tính toán mô phỏng

Khóa luận được xây dựng trên nền tảng cơ học lượng tử và phương pháp mô phỏng tính toán vật liệu từ nguyên tắc đầu (ab-initio):

  • Cơ sở cơ học lượng tử đa hạt: Bài toán trạng thái dừng của hệ đa hạt electron và hạt nhân được quy về phương trình Schrödinger. Tác giả áp dụng phép gần đúng Born - Oppenheimer dựa trên sự chênh lệch khối lượng giữa hạt nhân và electron (khoảng 1800 lần) để cố định vị trí các ion hạt nhân và chỉ giải chuyển động của hệ electron trong trường thế ngoài.
  • Lý thuyết phiếm hàm mật độ (DFT): Dựa trên hai định lý nền tảng Hohenberg - Kohn (năng lượng trạng thái cơ bản là phiếm hàm duy nhất của mật độ electron $n(\mathbf{r})$; mật độ cực tiểu hóa năng lượng toàn phần chính là mật độ trạng thái nền) và hệ phương trình tự hợp Kohn - Sham. Cách tiếp cận này chuyển bài toán hàm sóng $3N$ biến phức tạp thành bài toán biến phân của mật độ electron chỉ phụ thuộc 3 tọa độ không gian.
  • Phiếm hàm tương quan trao đổi: Luận văn xem xét xấp xỉ mật độ địa phương (LDA) dựa trên khí electron đồng nhất và xấp xỉ gradient tổng quát (GGA) kết hợp thông tin về gradient cục bộ của mật độ electron. Cụ thể, phiếm hàm Perdew - Burke - Ernzerhof (PBE) thuộc họ GGA được lựa chọn cho tính toán.
  • Hiệu chỉnh tương quan Coulomb cục bộ Hubbard ($DFT+U$ / $GGA+U$): Do các phiếm hàm gần đúng LDA/GGA thường đánh giá thấp tương quan electron mạnh ở các phân lớp $d$ của kim loại chuyển tiếp, lý thuyết Hubbard được bổ sung với thông số $U = 4\text{ eV}$ đặt tại vị trí phân lớp $d$ của nguyên tử Cr (và Fe).
  • Phương pháp giả thế sóng tăng cường xạ ảnh (PAW - Projector Augmented Wave): Kỹ thuật do Peter Blöchl phát triển giúp biến đổi hàm sóng vật lý gần hạt nhân thành hàm sóng phụ trơn tru, cho phép giảm năng lượng cắt sóng phẳng mà vẫn bảo toàn độ chính xác của các electron hóa trị.

Phương pháp nghiên cứu và điều kiện tính toán cụ thể:

  • Công cụ phần mềm: Phần mềm Quantum ESPRESSO (QE) phiên bản mã nguồn mở dùng để thực hiện các phép giải tự hợp (SCF), tối ưu hóa hình học, tính toán cấu trúc dải năng lượng và mật độ trạng thái; phần mềm XcrySDen và VESTA dùng để mô phỏng hình ảnh không gian tinh thể, đo đạc khoảng cách và góc liên kết, vẽ phân bố mật độ điện tích và hàm ELF; phần mềm Origin dùng để xử lý số liệu và kết xuất đồ thị.
  • Các thông số kỹ thuật mô phỏng:
    • Tiêu chuẩn hội tụ năng lượng SCF: $10^{-6}\text{ Ry}$.
    • Năng lượng cắt sóng phẳng (energy cutoff): $48\text{ Ry}$.
    • Lưới điểm $k$ không gian mạng đảo theo phương pháp Monkhorst - Park: $6 \times 6 \times 1$.
    • Đường dẫn vùng Brillouin vẽ dải năng lượng: Đi qua các điểm đối xứng đặc biệt $\Gamma - X - K$ với tổng số 21 điểm $k$.
    • Lớp đệm chân không (vacuum): Chiều dày $22\text{ \AA}$ dọc theo trục $z$.
    • Phân tích điện tích bằng phương pháp tích phân điện tích Bader và hàm định xứ electron (Electron Localization Function - ELF).

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

Chương 1: Cơ sở lý thuyết phiếm hàm mật độ (DFT)

Chương 1 hệ thống hóa cơ sở toán - lý lượng tử của phương pháp phiếm hàm mật độ. Tác giả trình bày phương trình Schrödinger tổng quát cho hệ $N$ electron chuyển động trong trường hạt nhân, qua đó phân tích rào cản tính toán khi số chiều không gian lên tới $3N$ biến. Để giải quyết, tác giả trình bày nguyên lý loại trừ Pauli, thế tương tác tĩnh điện Hartree $V_H(\mathbf{r})$ cùng vấn đề tự tương tác phi vật lý cần hiệu chỉnh, và khái niệm lỗ trống tương quan trao đổi (exchange - correlation hole).

Tiếp đó, chương luận giải chi tiết hai định lý Hohenberg - Kohn và phương trình Kohn - Sham, dẫn dắt thuật toán giải trường tự hợp (Self-Consistent Field - SCF) qua 4 bước: khởi tạo mật độ $n(\mathbf{r})$, giải phương trình tìm hàm sóng đơn electron $\Psi_i(\mathbf{r})$, tính lại mật độ $n_{KS}(\mathbf{r})$ và kiểm tra điều kiện hội tụ. Phần cuối chương so sánh các phiếm hàm tương quan trao đổi từ LDA đến GGA (nhấn mạnh hai dạng PW91 và PBE), đồng thời trình bày nguyên lý phương pháp giả thế, đặc biệt là ưu thế của giả thế sóng tăng cường xạ ảnh PAW trong việc xử lý lớp electron hóa trị và giảm năng lượng cắt sóng phẳng.

Chương 2: Các tính chất cơ bản của vật liệu đơn lớp CrPSe3

Chương 2 tổng quan toàn diện các kết quả thực nghiệm và mô phỏng lý thuyết đã được công bố quốc tế về hợp chất $\text{CrPSe}_3$ nói riêng và họ $\text{MPX}_3$ nói chung. Về mặt thực nghiệm, các dữ liệu tổng hợp tinh thể khối $\text{MPSe}_3$ ($\text{M} = \text{Cd}, \text{Cr}, \text{Fe}, \text{Mn}, \text{Sn}, \text{Zn}$) được trình bày thông qua ảnh kính hiển vi điện tử quét SEM (thể hiện cấu trúc phân phiến nhiều tấm xếp chồng), kết quả phân tích quang phổ tán sắc năng lượng tia X (EDS) và giản đồ nhiễu xạ tia X bột (XRD). Tỷ lệ thành phần nguyên tố xác định từ phổ EDS của $\text{CrPSe}_3$ cho thấy tỷ lệ $\text{Se}/\text{M} = 2.7$ và $\text{P}/\text{M} = 1.2$. Phân tích XRD khẳng định cấu trúc tinh thể thuộc nhóm không gian $C2/m$ đối với $\text{CrPSe}_3$.

Về mặt lý thuyết mô phỏng, chương trình bày cấu tạo mạng lưới tổ ong tương tự graphene của kim loại chuyển tiếp trong đơn lớp $\text{MPX}_3$, được liên kết bởi các đơn vị lưỡng tháp $(\text{P}_2\text{X}_6)^{4-}$. Các thông số mạng tối ưu hóa từ tài liệu tham khảo được tổng hợp: hằng số mạng $a = b = 6.364\text{ \AA}$ (hoặc $6.352\text{ \AA}$), độ dày lớp $h = 3.394\text{ \AA}$ (hoặc $3.534\text{ \AA}$), khoảng cách $\text{Cr}-\text{Cr} = 3.674\text{ \AA}$, $\text{Cr}-\text{Se} = 2.711\text{ \AA}$, $\text{P}-\text{Se} = 2.211\text{ \AA}$, khoảng cách $\text{P}-\text{P} = 2.15 - 2.7\text{ \AA}$, góc $\text{Cr}-\text{Se}-\text{Cr} = 85.3^\circ$ và $\text{Se}-\text{P}-\text{Se} = 113.5^\circ$.

Đối với tính chất từ, tác giả phân tích bốn cấu hình sắp xếp spin: sắt từ (FM), phản sắt từ Néel (nAFM), phản sắt từ Zigzag (zAFM) và phản sắt từ Stripy (sAFM). Mômen từ định vị chủ yếu tại các ion kim loại chuyển tiếp thông qua cơ chế siêu trao đổi gián tiếp $\text{Cr}-\text{Se}\cdots\text{Se}-\text{Cr}$. Đơn lớp $\text{CrPSe}3$ thể hiện trạng thái cơ bản là phản sắt từ Néel (nAFM) với nhiệt độ Néel $T_N \approx 105 - 136\text{ K}$, mômen từ hiệu dụng $\mu{eff} = 3.11\mu_B$, hằng số Weiss $\theta = -323\text{ K}$, năng lượng dị hướng từ trường $\text{MAE} = 650.3\text{ meV}$ và vận tốc Fermi $v_F \approx 16.11 \times 10^5\text{ m/s}$. Về cấu trúc điện tử, độ rộng vùng cấm của đơn lớp biến thiên từ $0.25\text{ eV}$ đến $1.99\text{ eV}$ phụ thuộc vào giá trị tham số Coulomb $U$ và phép xấp xỉ thế trao đổi (PBE, HSE hay DFT+U). Cuối chương, tác giả điểm qua các tiềm năng ứng dụng của vật liệu trong xúc tác phản ứng tạo hydro (HER), phản ứng tạo oxy (OER), phản ứng khử oxy (ORR), pin lưu trữ năng lượng lithium/natri và linh kiện spintronics.

Chương 3: Sử dụng lý thuyết phiếm hàm mật độ mô phỏng cấu trúc, tính chất điện, tính chất từ của vật liệu 2D CrPSe3 pha tạp Fe

Chương 3 trình bày các kết quả tính toán mô phỏng gốc do chính tác giả thực hiện bằng phần mềm Quantum ESPRESSO. Phần đầu chương mô tả chi tiết quy trình xây dựng siêu ô mạng, thiết lập các điều kiện biên và hội tụ thông số tính toán.

Trong phần khảo sát hệ đơn lớp $\text{CrPSe}_3$ nguyên bản, tác giả tiến hành tối ưu hóa cấu trúc tinh thể mạng monoclinic. Khi so sánh năng lượng toàn phần giữa các trạng thái từ tính khác nhau (NM, FM, nAFM, zAFM, sAFM), kết quả khẳng định cấu hình phản sắt từ Néel (nAFM) là trạng thái có mức năng lượng thấp nhất. Tính toán cấu trúc vùng năng lượng và mật độ trạng thái riêng (PDOS) khi không có và có bổ sung thông số Hubbard ($U = 4\text{ eV}$) chỉ ra rõ rệt vai trò của hiệu chỉnh Coulomb: việc bổ sung $U$ giúp tách biệt các dải obitan $3d$ của Cr, hiệu chỉnh mức Fermi và phản ánh chính xác bản chất điện tử của hệ. Tác giả cũng phân tích sự phân bố mật độ điện tử của các trạng thái spin-up và spin-down.

Trong phần khảo sát hệ đơn lớp $\text{CrPSe}_3$ pha tạp $25%\text{ Fe}$, tác giả thay thế một nguyên tử Cr bằng nguyên tử Fe trong siêu ô mạng hình chữ nhật. Kết quả tối ưu hóa hình học cho thấy các thông số mạng tinh thể và độ dài liên kết có sự biến dạng nhẹ xung quanh tâm pha tạp Fe do bán kính ion và cấu hình vỏ electron của Fe khác biệt. Khi tính toán năng lượng toàn phần của các trạng thái từ tính cho hệ pha tạp:

  • Khác biệt hoàn toàn với hệ nguyên bản, cấu hình sắt từ (FM) trở thành trạng thái cơ bản có năng lượng thấp nhất, vượt qua trạng thái nAFM.
  • Phân tích cấu trúc dải năng lượng và PDOS cho thấy vật liệu pha tạp $25%\text{ Fe}$ chuyển sang tính chất bán kim loại (half-metal): kênh spin-up có các dải năng lượng cắt qua mức Fermi (thể hiện tính dẫn điện kim loại), trong khi kênh spin-down tồn tại vùng cấm năng lượng rõ rệt (thể hiện tính cách điện).
  • Phân tích điện tích Bader cho thấy sự dịch chuyển điện tích cục bộ giữa Fe, Cr và các nguyên tử Se lân cận. Hàm định xứ electron (ELF) ở cả định dạng 2D và 3D xác nhận mức độ xen phủ và bản chất liên kết cộng hóa trị - ion phân cực trong cấu trúc sau khi pha tạp.

Kết quả và đóng góp khoa học

Kết quả chính kèm số liệu trong văn bản

Khóa luận đã thu thập, tính toán và lượng hóa các đại lượng vật lý đặc trưng của đơn lớp $\text{CrPSe}_3$ trước và sau pha tạp Fe:

Thông số / Đại lượng Hệ đơn lớp $\text{CrPSe}_3$ nguyên bản Hệ đơn lớp $\text{CrPSe}_3$ pha tạp $25%\text{ Fe}$
Cấu trúc mạng tinh thể Monoclinic, $a = b = 6.364\text{ \AA}$, độ dày $h = 3.394\text{ \AA}$ Biến dạng cục bộ tại vị trí thay thế Fe trong siêu ô mạng
Độ dài liên kết đặc trưng $\text{Cr}-\text{Cr} = 3.674\text{ \AA}$, $\text{Cr}-\text{Se} = 2.711\text{ \AA}$, $\text{P}-\text{Se} = 2.211\text{ \AA}$ Xuất hiện các liên kết mới $\text{Fe}-\text{Se}$, $\text{Fe}-\text{Cr}$
Góc liên kết $\text{Cr}-\text{Se}-\text{Cr} = 85.3^\circ$, $\text{Se}-\text{P}-\text{Se} = 113.5^\circ$ Biến thiên góc liên kết cục bộ quanh nguyên tử Fe
Trạng thái từ tính cơ bản Phản sắt từ Néel (nAFM) Sắt từ (FM)
Tính chất điện tử Bán dẫn hoặc kim loại tùy cấu hình spin và mức hiệu chỉnh $U$ Bán kim loại (Half-metal) với phân cực spin $100%$ tại mức Fermi
Năng lượng hiệu chỉnh Hubbard $U = 4\text{ eV}$ cho phân lớp $3d$ của Cr $U = 4\text{ eV}$ cho phân lớp $3d$ của Cr và Fe
Đặc trưng điện tích Điện tích tích tụ quanh các nguyên tử Se Dịch chuyển điện tích Bader quanh vị trí Fe; hàm ELF 2D/3D xác nhận liên kết

Đối với hệ thuần túy, các nghiên cứu trước được tổng hợp cho thấy năng lượng dị hướng từ trường $\text{MAE} = 650.3\text{ meV}$, vận tốc Fermi đạt $16.11 \times 10^5\text{ m/s}$ (hoặc $7.98 \times 10^5\text{ m/s}$), nhiệt độ chuyển pha Néel $T_N = 105 - 136\text{ K}$, độ rộng vùng cấm dao động từ $0.25\text{ eV}$ ($U=3\text{ eV}$) đến $1.99\text{ eV}$ ($U=4\text{ eV}$).

Giải pháp và kiến nghị kỹ thuật

  • Tác giả khẳng định việc đưa thông số tương tác Coulomb tại chỗ Hubbard $U = 4\text{ eV}$ vào mô hình DFT là bắt buộc để mô tả chính xác tương quan điện tử mạnh của các electron $3d$ trên ion kim loại chuyển tiếp trong hệ vật liệu $\text{MPX}_3$.
  • Đề xuất phương pháp pha tạp kim loại chuyển tiếp Fe với tỷ lệ $25%$ như một giải pháp điều khiển trật tự từ tính và cấu trúc điện tử, biến đổi vật liệu từ trạng thái phản sắt từ sang sắt từ và tạo ra trạng thái bán kim loại phục vụ linh kiện spin.

Đóng góp so với các nghiên cứu trước

  • Tác giả tự nêu rõ rằng trong các năm gần đây, các nghiên cứu về $\text{CrPSe}_3$ đơn lớp hầu như chỉ tập trung vào hệ thuần túy hoặc các hợp chất $\text{MPX}_3$ khác, hoàn toàn chưa có công trình nào khảo sát hệ $\text{CrPSe}_3$ pha tạp Fe.
  • Khóa luận là công trình đầu tiên tính toán khảo sát một cách có hệ thống ảnh hưởng của việc pha tạp nguyên tố sắt (Fe) $25%$ vào vị trí Cr trong đơn lớp $\text{CrPSe}_3$ bằng phương pháp phiếm hàm mật độ, cung cấp bộ dữ liệu lượng tử chi tiết về sự chuyển pha từ tính từ nAFM sang FM và sự hình thành trạng thái bán kim loại (half-metal).

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

Dựa trên các nhận xét phản biện và nội dung trình bày trong văn bản:

  • Hạn chế của nghiên cứu:
    • Nghiên cứu hoàn toàn dựa trên phương pháp tính toán mô phỏng lượng tử ở điều kiện $0\text{ K}$, chưa tiến hành tổng hợp thực nghiệm để kiểm chứng trực tiếp tính chất của đơn lớp $\text{CrPSe}_3$ pha tạp $25%\text{ Fe}$.
    • Luận văn mới chỉ dừng lại ở một tỷ lệ pha tạp cố định duy nhất là $25%\text{ Fe}$ thay thế nguyên tử Cr, chưa khảo sát dải nồng độ pha tạp rộng hơn (chẳng hạn $12.5%$, $50%$) để thiết lập giản đồ pha từ tính và điện tử hoàn chỉnh.
    • Theo nhận xét của giáo viên phản biện (TS. Đỗ Huy Bình), một số phát biểu trong văn bản ban đầu còn thiếu hình ảnh mô phỏng đi kèm để củng cố tính chắc chắn của luận điểm, đồng thời cần giải thích sâu hơn về nguyên nhân dẫn đến tính chất bán dẫn khi không tính đến tương tác từ cũng như sự khác biệt về độ dẫn điện giữa nghiên cứu này và tài liệu tham khảo [27].
  • Hướng phát triển tiếp theo:
    • Mở rộng tính toán mô phỏng cho các tỷ lệ pha tạp Fe khác nhau và khảo sát thêm các nguyên tố kim loại chuyển tiếp khác nhằm tìm kiếm các trạng thái từ tính và điện tử mới.
    • Nghiên cứu độ ổn định nhiệt động học và động học phân tử (ab-initio molecular dynamics - AIMD) của hệ pha tạp ở nhiệt độ phòng.
    • Khảo sát ảnh hưởng của các yếu tố ngoại vi như điện trường ngoài hoặc biến dạng cơ học (strain engineering) lên tính chất từ và điện của hệ vật liệu $\text{CrPSe}_3$ sau pha tạp.

Giá trị tham khảo

Khóa luận là nguồn tài liệu tham khảo hữu ích cho các nhóm đối tượng sau:

  • Sinh viên, học viên cao học chuyên ngành Công nghệ Vật liệu, Vật lý chất rắn, Hóa học lý thuyết và tính toán, Khoa học vật liệu: Tài liệu cung cấp bức tranh chi tiết từ lý thuyết nền tảng (DFT, phương trình Kohn - Sham, thế Coulomb Hubbard, giả thế PAW) đến kỹ thuật tính toán thực hành.
  • Người nghiên cứu về mô phỏng lượng tử sử dụng phần mềm Quantum ESPRESSO: Phần 3.2 cung cấp các thông số kỹ thuật chuẩn đã được tối ưu hóa cho hệ vật liệu hai chiều $\text{MPX}_3$, bao gồm năng lượng cắt sóng phẳng ($48\text{ Ry}$), lưới điểm $k$ Monkhorst - Park ($6 \times 6 \times 1$), điều kiện biên chân không ($22\text{ \AA}$) và thông số $U = 4\text{ eV}$.
  • Các nhà nghiên cứu về vật liệu hai chiều và spintronics: Chương 2 tổng hợp một khối lượng lớn dữ liệu thực nghiệm (SEM, EDS, XRD) và số liệu mô phỏng từ nhiều công trình quốc tế, giúp tiết kiệm thời gian tra cứu và đối chiếu các thông số mạng, hằng số tương tác trao đổi $j_1, j_2, j_3$ và nhiệt độ chuyển pha Néel.

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

1. Tại sao bắt buộc phải bổ sung tham số tương tác Coulomb tại chỗ Hubbard ($U$) trong các phép tính toán DFT cho hệ vật liệu $\text{MPX}_3$?

Trong các phép tính DFT sử dụng phiếm hàm xấp xỉ tiêu chuẩn như LDA hay GGA, năng lượng tự tương tác không được triệt tiêu hoàn toàn và tương quan điện tử mạnh giữa các electron định xứ trên phân lớp $3d$ của kim loại chuyển tiếp (như Cr, Fe) bị đánh giá thấp. Việc này thường dẫn đến sai lệch lớn trong việc xác định độ rộng vùng cấm hoặc mô tả sai trạng thái dẫn điện. Việc bổ sung tham số Hubbard ($U = 4\text{ eV}$) giúp mô tả chính xác lực đẩy Coulomb nội tại giữa các electron $3d$, tách biệt các dải năng lượng và phản ánh đúng trạng thái nền điện tử cũng như cấu hình từ tính của vật liệu.

2. Trạng thái từ tính cơ bản của đơn lớp $\text{CrPSe}_3$ nguyên bản là gì và nó thay đổi ra sao khi pha tạp $25%\text{ Fe}$?

Ở trạng thái nguyên bản, các tính toán so sánh năng lượng toàn phần giữa các cấu hình từ tính (FM, nAFM, zAFM, sAFM) chỉ ra rằng đơn lớp $\text{CrPSe}_3$ có trạng thái cơ bản là phản sắt từ Néel (nAFM). Tuy nhiên, khi pha tạp $25%\text{ Fe}$ thay thế vị trí của Cr, trật tự từ tính trong hệ bị thay đổi căn bản: trạng thái sắt từ (FM) có mức năng lượng thấp nhất và trở thành trạng thái cơ bản bền vững của vật liệu.

3. Sự khác biệt về tính chất điện tử giữa đơn lớp $\text{CrPSe}_3$ thuần túy và hệ pha tạp $25%\text{ Fe}$ được thể hiện như thế nào?

Đối với đơn lớp $\text{CrPSe}_3$ thuần túy, vật liệu thể hiện tính chất bán dẫn ở trạng thái phản sắt từ nAFM (với độ rộng vùng cấm phụ thuộc vào thông số $U$) hoặc tính chất kim loại ở trạng thái sắt từ. Trong khi đó, khi pha tạp $25%\text{ Fe}$, vật liệu chuyển hóa thành chất bán kim loại (half-metal): kênh phân cực spin-up thể hiện tính dẫn điện kim loại do các dải năng lượng cắt qua mức Fermi, trong khi kênh spin-down có khe năng lượng rõ rệt thể hiện tính cách điện, tạo ra độ phân cực spin $100%$ tại mức Fermi.

4. Các điều kiện và thông số kỹ thuật chính được thiết lập trong phần mềm Quantum ESPRESSO cho nghiên cứu này là gì?

Các điều kiện tính toán bao gồm: sử dụng phiếm hàm GGA kết hợp hàm giả thế sóng tăng cường xạ ảnh (PAW); bổ sung hiệu chỉnh Hubbard $U = 4\text{ eV}$ cho lớp $3d$ của Cr; tiêu chuẩn hội tụ năng lượng đạt $10^{-6}\text{ Ry}$; năng lượng cắt sóng phẳng (cutoff) là $48\text{ Ry}$; lưới lấy mẫu điểm $k$ trong vùng Brillouin là $6 \times 6 \times 1$ theo chuẩn Monkhorst - Park; đường dẫn tính vùng năng lượng gồm 21 điểm qua $\Gamma - X - K$; và lớp chân không dày $22\text{ \AA}$ dọc trục $z$ để ngăn chặn tương tác tuần hoàn giữa các lớp giả định.

5. Kết quả phân tích thành phần nguyên tố qua phổ EDS của tinh thể $\text{CrPSe}_3$ trong tài liệu thực nghiệm tham chiếu cho các giá trị nào?

Theo số liệu phân tích quang phổ tán sắc năng lượng tia X (EDS) được trích dẫn trong Bảng 2.1 của luận văn, thành phần phần trăm nguyên tố của mẫu $\text{CrPSe}_3$ được xác định với các tỷ lệ: tỷ lệ nguyên tử $\text{Se}/\text{M}$ (ở đây $\text{M} = \text{Cr}$) đạt giá trị $2.7$ và tỷ lệ $\text{P}/\text{M}$ đạt giá trị $1.2$ (với các tỷ lệ phần trăm được liệt kê gồm Cr: $7.1%$, Fe: $8.5%$, P: $19.5%$, Se: $34.7%$, Se/M: $3.3$ - $2.7$).

Kết luận

Khóa luận tốt nghiệp của tác giả Võ Thị Sanh đã ứng dụng thành công phương pháp tính toán từ nguyên tắc đầu dựa trên lý thuyết phiếm hàm mật độ kết hợp hiệu chỉnh Hubbard ($DFT+U$, $U=4\text{ eV}$) thông qua bộ phần mềm Quantum ESPRESSO để khảo sát chi tiết cấu trúc, tính chất điện và tính chất từ của vật liệu hai chiều $\text{CrPSe}_3$ đơn lớp. Nghiên cứu đã chứng minh rằng việc pha tạp $25%\text{ Fe}$ thay thế vị trí Cr tạo ra sự chuyển pha trật tự từ tính từ phản sắt từ Néel (nAFM) sang sắt từ (FM), đồng thời biến đổi vật liệu thành trạng thái bán kim loại (half-metal) có độ phân cực spin tuyệt đối tại mức Fermi. Kết quả tính toán cung cấp cơ sở dữ liệu định lượng có giá trị khoa học về họ trichalcogenide kim loại chuyển tiếp, đóng góp định hướng hữu ích cho việc thiết kế các vật liệu 2D thế hệ mới ứng dụng trong lĩnh vực spintronics và công nghệ điện tử nano.