Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Trong lĩnh vực khoa học vật liệu và công nghệ vi điện tử hiện đại, xu hướng thu nhỏ kích thước linh kiện, tăng tốc độ truyền dẫn và xử lý thông tin, giảm thiểu mức tiêu thụ năng lượng đã thúc đẩy mạnh mẽ các nghiên cứu về vật liệu thấp chiều. Khi kích thước không gian của vật liệu giảm xuống dưới 100 nm, năng lượng giam hãm lượng tử và các hiệu ứng bề mặt bắt đầu đóng vai trò chi phối, làm biến đổi căn bản các tính chất vật lý so với vật liệu dạng khối. Trong nhóm vật liệu này, các oxit phức hợp có cấu trúc perovskite $ABO_3$ và các biến thể cấu trúc của chúng giữ vị trí trung tâm nhờ sự đa dạng về các tính chất điện môi, sắt điện, áp điện, hỏa điện, trật tự từ và khả năng xuất hiện các trạng thái lượng tử mới.

Đặc trưng cốt lõi của oxit phức hợp là tương tác lai hóa giữa các trạng thái $2p$ của anion oxy và $3d$ của cation kim loại chuyển tiếp tại vị trí $B$ trong mạng phối trí $BO_x$. Mọi biến dạng nhỏ của mạng phối trí (méo mạng Jahn-Teller, méo mạng sắt điện, sự nghiêng khối bát diện) hoặc sự suy giảm số chiều không gian đều dẫn đến sự thay đổi mạnh mẽ của tương tác $p-d$, từ đó chi phối sự cạnh tranh giữa các cơ chế tương tác trao đổi spin, trao đổi điện tích và liên kết spin - quỹ đạo. Sự can thiệp vào mạng phối trí $BO_x$ mở ra khả năng điều khiển các hiệu ứng từ điện trở khổng lồ (CMR), tính đa phân cực (multiferroic), tính dẫn điện hai chiều tại bề mặt phân cách và trật tự từ một chiều. Tuy nhiên, cơ chế vi mô của tương tác $p-d$ trong các hệ cấu trúc thấp chiều chưa được làm sáng tỏ đầy đủ. Do đó, việc mô phỏng tính toán cấu trúc điện tử từ các nguyên lý ban đầu kết hợp với khảo sát thực nghiệm là nhiệm vụ cần thiết nhằm xây dựng cơ sở khoa học cho việc thiết kế và chế tạo các linh kiện spintronics thế hệ mới.

Mục tiêu nghiên cứu

Mục tiêu chính của luận án được xác định bao gồm:

  1. Tìm hiểu mối liên hệ định lượng giữa cấu trúc tinh thể, cấu trúc điện tử và các tính chất vật lý đặc trưng của các hệ oxit phức hợp thấp chiều.
  2. Tiên đoán, giải thích các cơ chế vi mô chi phối tính chất mới xuất hiện tại bề mặt phân cách và trong cấu trúc thấp chiều.
  3. Chế tạo thực nghiệm và khảo sát các đặc tính cấu trúc, dao động mạng trên hệ vật liệu một chiều tiêu biểu để kiểm chứng và bổ trợ cho các mô hình lý thuyết.

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

  • Đối tượng nghiên cứu:
    • Hệ màng mỏng tiếp xúc dị thể hai chiều $La_{1-x}Sr_xMnO_3/BaTiO_3$ (LSMO/BTO).
    • Hệ màng mỏng tiếp xúc dị thể hai chiều $LaAlO_3/SrTiO_3$ (LAO/STO).
    • Hệ chuỗi spin phản sắt từ một chiều $Ca_2CuO_3$ và hệ pha tạp uranium $Ca_2CuO_3:U_x$.
  • Phạm vi nghiên cứu: Luận án tập trung vào tính toán mô phỏng cấu trúc điện tử dựa trên lý thuyết phiếm hàm mật độ (DFT) và phương pháp trường tự hợp Hartree-Fock; đồng thời thực hiện quy trình chế tạo mẫu đa tinh thể $Ca_2CuO_3:U_x$ bằng phương pháp sol-gel và phân tích cấu trúc, phổ dao động tại Bộ môn Vật lý Chất rắn và Phòng thí nghiệm Tính toán trong Khoa học Vật liệu, Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (hoàn thành năm 2017).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước

Quá trình phát triển lý thuyết và thực nghiệm về oxit phức hợp được tác giả hệ thống hóa qua các công trình nền tảng:

  • Cấu trúc điện tử và phân loại trạng thái cơ bản: Zaanen, Sawatzky và Allen (1985) đã thiết lập giản đồ pha ZSA dựa trên mô hình đơn tạp Anderson, khảo sát tương quan giữa năng lượng đẩy Coulomb tại chỗ $U$, năng lượng chuyển điện tích $\Delta_{pd}$, độ rộng vùng O $2p$ ($W$) và tích phân truyền dịch $t_{pd}$ để phân loại các trạng thái cách điện Mott-Hubbard, cách điện chuyển điện tích và kim loại. Các phát triển tiếp theo mở rộng cho chất cách điện cộng hóa trị với khe năng lượng phụ thuộc nhiệt độ.
  • Tương tác trao đổi từ tính: Cơ chế tương tác trao đổi gián tiếp RKKY qua điện tử dẫn (Ruderman, Kittel, Kasuya, Yosida); cơ chế trao đổi kép (Double Exchange - DE) giữa các ion $Mn^{3+}$ và $Mn^{4+}$ giải thích tính sắt từ trong manganite; lý thuyết siêu trao đổi (Super Exchange - SE) của Anderson (1950, 1959) và hệ thống quy tắc bán thực nghiệm Goodenough – Kanamori – Anderson (GKA) (Goodenough, 1955, 1958; Kanamori, 1959) định hình dấu và độ lớn tương tác spin qua góc liên kết cation-anion-cation ($180^\circ$ và $90^\circ$).
  • Vật liệu đa phân cực và tiếp xúc dị thể LSMO/BTO: Newnham và cộng sự phân loại các cấu hình composite đa pha ($0-3, 2-2, 1-3$). Yau và cộng sự (2007) chế tạo tiếp xúc chui ngầm đa phân cực (MFTJ) $LSMO\text{-}0,3/BTO/NiFe$ hoạt động 4 trạng thái điện trở. Y. Yin và cộng sự khảo sát hệ $LCMO\text{-}0,3/BSTO\text{-}0,5/LCMO\text{-}0,3/Au$ và $LSMO\text{-}0,3/BTO/LCMO\text{-}0,5/LSMO\text{-}0,5$ ghi nhận hệ số điện trở chui ngầm sắt điện (TER) đạt mức $5000%$ ở $40\text{ K}$. J. Tsymbal đề xuất cơ chế liên kết từ - điện thông qua hiện tượng chắn tĩnh điện tại mặt phân cách manganite/sắt điện.
  • Hệ khí điện tử hai chiều tại mặt phân cách LAO/STO: Ohtomo và Hwang (2004) phát hiện tính dẫn hai chiều giả định tại tiếp xúc loại $n$ ($LaO/TiO_2$) giữa hai chất cách điện dải rộng. Brinkman và cộng sự ghi nhận hiệu ứng từ điện trở âm và trật tự từ ở nhiệt độ thấp ($0,3\text{ K}$). M. Stengel kết hợp phương trình Maxwell với tính toán ab initio chứng minh vai trò của méo mạng phân cực và tính phi tuyến điện môi trong STO đối với sự hình thành hố thế giam hãm điện tử. Z. Popovic và cộng sự phân loại ba nhóm điện tử Ti $3d$ ($d_{xy}, d_{xz}, d_{yz}$). N. Pavlenko và cộng sự nghiên cứu vai trò của nút khuyết oxy đối với từ tính và tính siêu dẫn. Lee và cộng sự áp dụng hiệu chỉnh LDA+U để tái hiện nồng độ hạt tải thực nghiệm.
  • Hệ spin một chiều $A_2CuO_3$ ($A = Ca, Sr$): Ami và cộng sự nghiên cứu tán xạ neutron trên $Sr_2CuO_3$, xác định cấu trúc trực giao bền vững từ nhiệt độ phòng xuống $11\text{ K}$. Zhang và cộng sự kiểm chứng độ bền cấu trúc của $Ca_2CuO_3$ dưới áp suất thủy tĩnh tới $34\text{ GPa}$ ($B_0 \approx 165\text{ GPa}$). K. Maiti và cộng sự sử dụng phổ XPS và BIS để xác định khe năng lượng $E_g$ của $Ca_2CuO_3$ ($1,7\text{ eV}$) và $Sr_2CuO_3$ ($1,5\text{ eV}$), đồng thời tham số hóa tương tác cụm $CuO_4$. Các phép đo cộng hưởng spin muon ($\mu$SR) xác định nhiệt độ chuyển pha phản sắt từ 3D $T_N = 5,41\text{ K}$ ($Sr_2CuO_3$) và $T_N = 8\text{ K}$ ($Ca_2CuO_3$) với mô men từ cục bộ rất nhỏ ($0,06\text{ }\mu_B$ và $0,09\text{ }\mu_B$).

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

  1. Thiếu các nghiên cứu định lượng chi tiết về ảnh hưởng của định hướng tinh thể ([001] và [111]) cũng như ứng suất biến dạng đế (đế BTO so với đế STO) lên liên kết từ - điện và sự chuyển đổi trật tự orbital tại mặt phân cách siêu mạng $La_{1-x}Sr_xMnO_3/BaTiO_3$.
  2. Cơ chế vi mô dẫn đến sự giam hãm mạnh của khí điện tử tự do trong tiếp xúc $LaAlO_3/SrTiO_3$ loại $n$ và sự vắng mặt của tính dẫn thực nghiệm trong tiếp xúc loại $p$ chưa được phân tích thấu đáo trên cơ sở tương tác trường tinh thể và méo mạng sắt điện nội tại.
  3. Sự mâu thuẫn giữa kết luận thực nghiệm quang phổ (xem $Ca_2CuO_3$ là chất cách điện cộng hóa trị) và các tính toán cấu trúc dải năng lượng thông thường (xem là cách điện chuyển điện tích), cùng với việc thiếu các dữ liệu thực nghiệm về ảnh hưởng của pha tạp ion hóa trị cao (uranium) lên cấu trúc vi mô và dao động phonon mạng của hệ chuỗi spin $Ca_2CuO_3$.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và khái niệm sử dụng

  • Cấu trúc tinh thể mạng Perovskite: Khảo sát mạng $ABO_3$ lý tưởng đối xứng lập phương $Pm\bar{3}m$, các dạng méo mạng hạ thấp đối xứng (méo Jahn-Teller do suy biến orbital, méo sắt điện do dịch chuyển cation khỏi tâm bát diện, sự nghiêng/quay bát diện $BO_6$), và các biến thể perovskite lớp (pha Ruddleson-Popper $A_{n+1}B_nO_{3n+1}$, pha Aurivillius $Bi_2O_2/A_{n-1}B_nO_{3n-1}$, pha Dion-Jacobson $M^{+1}A_{n-1}B_nO_{3n+1}$).
  • Lý thuyết trường tinh thể và cấu trúc dải năng lượng: Sự tách mức năng lượng của các orbital kim loại chuyển tiếp $3d$ trong các trường phối trí bát diện ($e_g$ và $t_{2g}$), kim tự tháp và vuông ($d_{x^2-y^2}, d_{z^2}, d_{xy}, d_{xz}, d_{yz}$). Mô hình Hubbard và mô hình ZSA mô tả sự cạnh tranh giữa $U$ và $\Delta_{pd}$.
  • Cơ chế tương tác từ và từ - điện: Tương tác trao đổi trực tiếp, siêu trao đổi qua cầu oxy theo quy tắc GKA, trao đổi kép phụ thuộc định hướng spin của điện tử $e_g$ chuyển động; hiện tượng phân cực điện cảm ứng, hiệu ứng chắn tĩnh điện tại mặt phân cách kim loại/sắt điện và tương tác qua liên kết hóa học.
  • Mô hình "Đột biến phân cực" (Polar Catastrophe): Sự tích tụ điện thế do sự xen kẽ của các lớp tích điện $(LaO)^+$ và $(AlO_2)^-$ trên đế trung hòa $(SrO)^0/(TiO_2)^0$, dẫn đến sự dịch chuyển điện tích $0,5|e|$ sang lớp tiếp xúc.

Phương pháp nghiên cứu

Luận án kết hợp giữa phương pháp mô phỏng lý thuyết và phương pháp thực nghiệm:

Phương pháp Công cụ / Kỹ thuật Đối tượng / Mục đích
Tính toán Ab Initio Lý thuyết phiếm hàm mật độ (DFT), phép gần đúng LDA, GGA, phiếm hàm lai hóa B3LYP, B-F:LYP, gần đúng Hartree-Fock (HF), thuật toán lặp tự hợp SCF. Mô phỏng cấu trúc dải năng lượng, mật độ trạng thái (DOS), mô men từ, phân bố điện tích và trật tự orbital của siêu mạng LSMO/BTO, dị thể LAO/STO, và hệ chuỗi spin $Ca_2CuO_3$.
Tổng hợp mẫu Phương pháp hóa học Sol-Gel. Chế tạo hệ vật liệu đa tinh thể $Ca_2CuO_3$ pha tạp uranium ($Ca_2CuO_3:U_x$ với $x = 0 \div 0,05$).
Phân tích nhiệt Phân tích nhiệt vi sai (DTA) và nhiệt trọng lượng (TGA). Xác định các giai đoạn phân hủy nhiệt của gel, xác lập chế độ nhiệt thiêu kết tạo pha tinh thể.
Phân tích cấu trúc Nhiễu xạ tia X (XRD) trên bột tinh thể, phân tích cấu trúc tinh thể theo phương pháp Rietveld sử dụng phần mềm WinMProf. Xác định độ tinh khiết pha, nhóm không gian, hằng số mạng ($a, b, c$), vị trí tọa độ nguyên tử Wyckoff, độ dài liên kết và hệ số dao động nhiệt.
Hình thái & Thành phần Kính hiển vi điện tử quét (SEM) và phổ tán sắc năng lượng tia X (EDX / EDS). Khảo sát vi cấu trúc bề mặt, kích thước hạt tinh thể và định lượng thành phần nguyên tố, kiểm tra sự phân bố của tạp chất U.
Quang phổ dao động Phổ tán xạ Raman và phổ hồng ngoại truyền qua (FTIR). Khảo sát các mode dao động phonon quang hoạt động Raman ($A_g$) và hồng ngoại ($B_{2u}$), đối chiếu với tính toán lý thuyết nhóm.

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

Chương 1 – Tổng quan về oxit phức hợp

Chương này cung cấp bức tranh toàn diện về các đặc tính cấu trúc, điện tử, từ tính và liên kết từ - điện của họ vật liệu oxit phức hợp:

  • Đặc điểm cấu trúc: Trình bày cấu trúc perovskite lập phương đối xứng cao ($Pm\bar{3}m$) và nguyên nhân gây bất ổn định cấu trúc dẫn đến các chuyển pha hình học (tứ giác, trực thoi, hộp thoi, một nghiêng). Phân tích chi tiết ba họ perovskite lớp: Ruddleson-Popper, Aurivillius và Dion-Jacobson.
  • Cấu trúc vùng năng lượng: Hệ thống hóa giản đồ pha ZSA và mở rộng mô hình Hubbard đa vùng, chỉ rõ sự chuyển đổi giữa các pha Mott-Hubbard, pha cách điện chuyển điện tích, pha kim loại và pha cách điện cộng hóa trị. Phân tích chi tiết sự tách mức trường tinh thể của orbital $3d$ trong các trường bát diện, kim tự tháp và mặt phẳng vuông.
  • Tương tác trao đổi và tương tác từ - điện: Trình bày bản chất cơ học lượng tử của tương tác trao đổi, cơ chế trao đổi trực tiếp, trao đổi gián tiếp qua điện tử dẫn RKKY, trao đổi kép trong manganite và siêu trao đổi theo quy tắc GKA. Phân tích các cơ chế đa phân cực đơn pha (tính sắt điện hình học trong $AMnO_3$, phân cực do spin xoắn qua tương tác Dzyaloshinskii-Moriya) và đa phân cực hai pha dạng composite màng mỏng.
  • Các hệ vật liệu cụ thể của luận án:
    • Hệ $La_{1-x}Sr_xMnO_3/BaTiO_3$: Phân tích giản đồ pha của LSMO-x, đặc tính sắt điện của BTO (chuyển pha từ hộp thoi sang một nghiêng ở $-90^\circ\text{C}$, sang tứ giác ở $5^\circ\text{C}$, sang lập phương ở $135^\circ\text{C}$ với độ phân cực $P = 26\text{ }\mu\text{C/cm}^2$), cấu trúc tiếp xúc chui ngầm MFTJ và cơ chế chắn tĩnh điện điều khiển trật tự từ.
    • Hệ $LaAlO_3/SrTiO_3$: So sánh đặc tính của hai chất cách điện thành phần (LAO có $E_g = 5,6\text{ eV}, \varepsilon \approx 24$; STO có $E_g = 3,2\text{ eV}, \varepsilon \approx 300$ ở nhiệt độ phòng). Phân tích mô hình phân cực lớp tiếp xúc loại $n$ ($LaO/TiO_2$) và loại $p$ ($AlO_2/SrO$), các thực nghiệm về khí điện tử hai chiều (độ linh động đạt $10^4\text{ cm}^2/\text{V.s}$ ở $4\text{ K}$), sự giam hãm trong phạm vi $4\text{ nm}$ ở $300\text{ K}$ và $12\text{ nm}$ ở $8\text{ K}$, cùng các giả thuyết về tán xạ Kondo và siêu dẫn.
    • Hệ spin một chiều $A_2CuO_3$: Phân tích cấu trúc trực giao $Immm$ gồm các chuỗi vuông phẳng $CuO_4$ chia sẻ góc, đối chiếu với cấu trúc tứ giác $I4/mmm$ của $La_2CuO_4$. Tổng hợp các kết quả thực nghiệm XPS/BIS và tham số hóa tương tác cụm.

Bảng tổng hợp tham số hóa phổ vùng dẫn và vùng hóa trị trong cụm $CuO_4$ (theo dữ liệu trích xuất từ Chương 1):

Vật liệu Số chiều mạng $Cu-O-Cu$ $U\text{ (eV)}$ $\Delta_{pd}\text{ (eV)}$ $t_{pd}\text{ (eV)}$
$Ca_2CuO_3$ 1D $7,2$ $-1,5$ $0,5^{(*)} ;/; 1,4^{(**)}$
$Sr_2CuO_3$ 1D $7,2$ $-1,5$ $0,0^{(*)} ;/; 0,8^{(**)}$
$La_2CuO_4$ 2D $7,2$ $1,0$ $-1,5^{(*)}$
$CuO$ 3D $7,8$ $2,2$ $-1,15$

Ghi chú: $^{()}$ Chưa xét độ rộng vùng O $2p$; $^{(**)}$ Đã xét độ rộng vùng O $2p$.*

Chương 2 – Phương pháp nghiên cứu

Chương này trình bày cơ sở toán học và quy trình kỹ thuật của các phương pháp được sử dụng:

  • Phương pháp tính toán từ các nguyên lý ban đầu: Thiết lập phương trình Hartree-Fock và các định lý Hohenberg-Kohn trong lý thuyết phiếm hàm mật độ (DFT), phương trình Kohn-Sham và quy trình giải lặp tự hợp SCF. Mô tả các phép gần đúng năng lượng trao đổi - tương quan địa phương (LDA), gradient suy rộng (GGA) và phiếm hàm lai hóa B3LYP.
  • Phương pháp thực nghiệm: Trình bày nguyên lý đo DTA/TGA; cơ chế phát xạ tia X đặc trưng trong phân tích EDX; định luật phản xạ Bragg ($2d\sin\theta = n\lambda$) và phương pháp Rietveld tinh chỉnh cấu trúc trên phần mềm WinMProf; nguyên lý quét chùm điện tử trong kính hiển vi SEM; giản đồ chuyển mức năng lượng phonon trong tán xạ đàn hồi Rayleigh và tán xạ không đàn hồi Raman.

Chương 3 – Kết quả và thảo luận

1. Kết quả tính toán trên hệ tiếp xúc dị thể $La_{1-x}Sr_xMnO_3/BaTiO_3$

  • Mô hình tính toán: Xây dựng ô siêu mạng LSMO/BTO theo hai định hướng tinh thể [001] và [111], khảo sát các trạng thái định hướng phân cực điện $\vec{P}$ của lớp BTO hướng về bề mặt phân cách thứ nhất ($IF_1$) hoặc bề mặt phân cách thứ hai ($IF_2$).
  • Cơ chế chắn tĩnh điện và mô men từ: Tính toán chỉ ra sự bất đối xứng của mô men từ ion Mn giữa hai bề mặt phân cách do hiện tượng tích tụ điện tích âm tại $IF_1$ và điện tích dương tại $IF_2$ nhằm chắn điện trường phân cực tự phát của BTO. Độ biến thiên mô men từ tổng cộng $\Delta M$ và hệ số liên kết từ - điện $S$ phụ thuộc rõ rệt vào định hướng màng.
  • Trật tự orbital và trạng thái từ:
    • Với định hướng [001] dưới ứng suất của đế BTO (tối ưu hóa $c = 39,35\text{ \AA}$ và kéo dãn $c = 40\text{ \AA}$), mật độ trạng thái Mn $3d$ cho thấy sự chiếm đóng ưu tiên của orbital $d_{x^2-y^2}$ tại bề mặt tích tụ điện tích, thúc đẩy tương tác siêu trao đổi phản sắt từ trong mặt phẳng (pha A-AFM). Ngược lại, tại bề mặt suy giảm điện tử, orbital $d_{z^2}$ chiếm ưu thế thúc đẩy pha C-AFM.
    • So sánh năng lượng giữa trạng thái sắt từ (FM) và các cấu hình phản sắt từ (AFM) khẳng định sự lật độ phân cực điện $\vec{P}$ của BTO có thể điều khiển trực tiếp quá trình chuyển pha từ tính tại lớp ranh giới LSMO.

2. Kết quả tính toán trên hệ tiếp xúc dị thể $LaAlO_3/SrTiO_3$

  • Mô hình tính toán: Sử dụng siêu ô mạng gồm 5 đơn lớp $LaAlO_3$ đặt trên 3 đơn lớp $SrTiO_3$ theo phương [001], bổ sung lớp chân không dày $15\text{ \AA}$ để khử tương tác tuần hoàn dọc trục $c$. Các lớp nguyên tử được đánh số từ 1 (SrO bề mặt STO) đến 6 ($TiO_2$ tại mặt phân cách) và từ 7 (LaO tại mặt phân cách) đến 16 ($AlO_2$ bề mặt LAO).
  • Phân bố mật độ trạng thái và giam hãm điện tử: Tính toán dải năng lượng và DOS từng lớp khẳng định sự chuyển điện tích từ bề mặt LAO sang các trạng thái Ti $3d$ tại mặt phân cách $TiO_2$.
  • Cơ chế giam hãm hố thế: Phân tích chỉ rõ méo mạng sắt điện nội tại trong các lớp bát diện $TiO_6$ lân cận mặt phân cách kết hợp với đặc tính phi tuyến điện môi của STO tạo ra một giếng thế hẹp và sâu. Điều này bẻ cong dải dẫn, giam hãm chặt các điện tử Ti $d_{xy}$ trong phạm vi vài đơn lớp tinh thể sát mặt tiếp xúc, giải thích nguồn gốc vi mô của khí điện tử hai chiều (2DEG).

3. Kết quả nghiên cứu trên hệ spin phản sắt từ một chiều $Ca_2CuO_3$

  • Tổng hợp và phân tích nhiệt: Đường cong DTA/TGA của gel $Ca_2CuO_3:U_{0,04}$ ghi nhận các hiệu ứng tỏa nhiệt và mất khối lượng liên quan đến quá trình cháy hợp chất hữu cơ và phân hủy muối nitrate, xác định nhiệt độ thiêu kết tối ưu để thu được pha đơn tinh khiết.
  • Cấu trúc tinh thể và hình thái vi mô:
    • Phổ XRD của chuỗi mẫu $Ca_2CuO_3:U_x$ ($x = 0; 0,01; 0,025; 0,05$) được làm khớp chính xác theo phương pháp Rietveld trên cấu trúc trực giao tâm khối, nhóm không gian $Immm$. Hằng số mạng thu được đạt độ chính xác cao: $a \approx 3,26\text{ \AA}, b \approx 3,78\text{ \AA}, c \approx 12,26\text{ \AA}$.
    • Xác định chi tiết các tọa độ vị trí Wyckoff và chỉ số dao động nhiệt của các nguyên tử Ca, Cu, O1, O2 trong ô mạng.
    • Ảnh SEM và phổ EDX xác nhận các hạt tinh thể phát triển đồng đều, không xuất hiện phân pha cục bộ của nguyên tố uranium; kích thước hạt có xu hướng biến đổi có quy luật theo nồng độ pha tạp $x$.
  • Cấu trúc điện tử: So sánh tính toán DFT và HF giữa mô hình không xét trật tự spin (ô mạng $Immm$) và có xét trật tự phản sắt từ (ô mạng $P1$). Kết quả khẳng định việc đưa trật tự spin AFM vào tính toán là điều kiện bắt buộc để mở ra khe năng lượng vùng cấm phù hợp với thực nghiệm, phản ánh bản chất tương quan điện tử mạnh dọc theo chuỗi $-Cu-O-Cu-$.
  • Quang phổ dao động phonon:
    • Phổ hồng ngoại truyền qua (FTIR) ghi nhận 3 mode dao động quang hoạt động hồng ngoại loại $B_{2u}$, hoàn toàn phù hợp với tính toán Hartree-Fock sử dụng hệ cơ sở Gaussian 3-21G.
    • Phổ tán xạ Raman xác định 2 mode hoạt động $A_g$. Khi tăng nồng độ pha tạp uranium, xuất hiện các đỉnh tán xạ mới có cường độ tăng dần theo $x$. Phân tích mô hình cụm chuỗi $(Cu-O)n$ với độ dài khác nhau chỉ ra rằng các đỉnh mới này bắt nguồn từ sự phá vỡ tính đối xứng tịnh tiến một chiều và dao động đối xứng của nguyên tử oxy dọc theo phương chuỗi ($O{sym,par}$).

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và mô phỏng

  1. Làm sáng tỏ cơ chế liên kết từ - điện tại bề mặt phân cách $La_{1-x}Sr_xMnO_3/BaTiO_3$: Chứng minh định lượng rằng tương tác từ - điện bắt nguồn từ sự kết hợp đồng thời giữa hiệu ứng chắn tĩnh điện bề mặt và sự tái phân bố trật tự orbital $3d$ ($d_{x^2-y^2}$ so với $d_{z^2}$) của ion Mn. Kết quả chỉ ra khả năng điều khiển trật tự từ bề mặt (chuyển đổi giữa FM và AFM) thông qua việc đảo chiều véc tơ phân cực điện tự phát của lớp sắt điện BTO.
  2. Xác lập cơ chế giam hãm hạt tải trong dị thể $LaAlO_3/SrTiO_3$: Chứng minh rằng sự kết hợp giữa méo mạng sắt điện phân cực tại lớp biên $TiO_2$ và tính phi tuyến của hằng số điện môi trong vật liệu STO tạo ra một hố thế năng sâu, đóng vai trò quyết định trong việc giam hãm khí điện tử tự do trong phạm vi hẹp ở mặt tiếp xúc loại $n$.
  3. Giải thích bản chất cách điện và dao động mạng của hệ một chiều $Ca_2CuO_3$: Làm rõ vai trò của trật tự spin phản sắt từ dọc theo chuỗi đối với sự mở rộng khe năng lượng dải dẫn - hóa trị; cung cấp lời giải thích nhất quán cho cơ chế mở rộng dải phổ Raman khi pha tạp dựa trên hiệu ứng suy giảm độ dài chuỗi tương quan $(Cu-O)_n$.

Đóng góp mới về mặt thực nghiệm và công nghệ

  1. Hoàn thiện quy trình sol-gel chế tạo $Ca_2CuO_3:U_x$: Xác lập thành công thông số công nghệ nhiệt phân và thiêu kết để tổng hợp hệ vật liệu một chiều $Ca_2CuO_3$ pha tạp uranium với độ đơn pha cao và phân bố tạp chất đồng nhất trong dải nồng độ $x = 0 \div 0,05$.
  2. Cung cấp bộ dữ liệu thực nghiệm chuẩn xác: Xác định hệ thống thông số mạng tinh thể Rietveld, tọa độ Wyckoff, hệ số dao động nhiệt và phổ tán xạ Raman/FTIR chuẩn của hệ $Ca_2CuO_3$ pha tạp uranium, đóng góp dữ liệu cấu trúc và phổ phonon cho cơ sở dữ liệu oxit phức hợp thấp chiều.

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

Hạn chế của luận án

  • Các mô phỏng cấu trúc điện tử dựa trên lý thuyết phiếm hàm mật độ (DFT) và Hartree-Fock chủ yếu thực hiện ở trạng thái cơ bản ($T = 0\text{ K}$) trên các cấu hình tuần hoàn lý tưởng, chưa mô phỏng đầy đủ ảnh hưởng của dao động nhiệt lượng tử ở nhiệt độ phòng và các dạng sai hỏng phức tạp (như sai hỏng phối trí ngẫu nhiên, nút khuyết oxy đa dạng).
  • Các khảo sát thực nghiệm trên hệ $Ca_2CuO_3:U_x$ mới dừng lại ở mẫu đa tinh thể dạng bột, chưa thực hiện nuôi đơn tinh thể kích thước lớn hoặc chế tạo màng mỏng định hướng để đo trực tiếp các tính chất dị hướng truyền tải spin và từ trở dọc theo trục chuỗi.

Hướng nghiên cứu tiếp theo

  • Mở rộng tính toán cấu trúc dải năng lượng và động học spin cho các siêu ô mạng kích thước lớn chứa nút khuyết oxy và sai hỏng bề mặt bằng các phương pháp hiệu chỉnh tương quan nâng cao (như DFT+DMFT).
  • Nghiên cứu công nghệ chế tạo màng mỏng dị thể tiếp xúc đa lớp chất lượng cao (LSMO/BTO và LAO/STO) bằng phương pháp lắng đọng xung laser (PLD), khảo sát thực nghiệm các đặc trưng chuyển mạch từ - điện của cấu trúc MFTJ ở nhiệt độ phòng.
  • Khảo sát sâu hơn các tính chất từ nhiệt độ thấp và từ điện trở của hệ $Ca_2CuO_3$ pha tạp trên các mẫu đơn tinh thể nuôi bằng phương pháp vùng nổi quang học (TSFZ).

Giá trị tham khảo

  • Đối với nghiên cứu sinh và học viên cao học ngành Vật lý chất rắn, Khoa học vật liệu: Luận án là tài liệu tham khảo chi tiết về phương pháp luận kết hợp giữa tính toán hóa lượng tử/vật lý chất rắn từ nguyên lý ban đầu (DFT, HF) với các phép phân tích thực nghiệm cấu trúc tinh thể chuyên sâu (tinh chỉnh Rietveld qua WinMProf, phân tích phổ tán xạ Raman và phổ hồng ngoại).
  • Đối với giảng viên và cán bộ nghiên cứu vật lý vật liệu: Cung cấp hệ thống lý thuyết chuẩn mực về oxit phức hợp, các cơ chế vi mô của tương tác $p-d$, chuyển pha kim loại - cách điện theo mô hình ZSA, quy tắc siêu trao đổi GKA và các mô hình tiếp xúc dị thể oxit.
  • Đối với các kỹ sư và chuyên gia phát triển công nghệ Spintronics: Các dữ liệu định lượng về liên kết từ - điện tại mặt phân cách manganite/sắt điện và cấu trúc vùng năng lượng tiếp xúc $LaAlO_3/SrTiO_3$ là cơ sở định hướng hữu ích cho việc thiết kế, tối ưu hóa các linh kiện nhớ không xóa đa trạng thái (MFTJ), van spin và transistor có độ linh động điện tử cao (HEMT).

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

1. Cơ chế chính dẫn đến tính dẫn điện hai chiều tại mặt tiếp xúc dị thể $LaAlO_3/SrTiO_3$ loại $n$ là gì?

Tính dẫn điện hai chiều tại mặt tiếp xúc loại $n$ ($LaO/TiO_2$) xuất phát từ hiện tượng "đột biến phân cực" (Polar Catastrophe). Do sự không liên tục của trường cảm ứng điện giữa các lớp tích điện $(LaO)^+/ (AlO_2)^-$ của LAO và các lớp trung hòa $(SrO)^0/ (TiO_2)^0$ của STO, điện thế tĩnh điện bị tích tụ tăng dần theo độ dày của LAO. Khi độ dày đạt tới giá trị tới hạn (khoảng 4 đơn lớp), có sự dịch chuyển điện tích khoảng $0,5|e|$ từ bề mặt LAO sang các trạng thái Ti $3d$ tại mặt phân cách STO để triệt tiêu sự phân cực. Các điện tử này bị giam hãm trong một hố thế hẹp tạo bởi méo mạng sắt điện của STO, hình thành nên lớp khí điện tử tự do giả hai chiều (2DEG).

2. Sự đảo chiều phân cực điện trong $BaTiO_3$ ảnh hưởng như thế nào đến tính chất từ của $La_{1-x}Sr_xMnO_3$ trong siêu mạng LSMO/BTO?

Khi véc tơ phân cực điện tự phát $\vec{P}$ của lớp sắt điện BTO đổi hướng, hiệu ứng chắn tĩnh điện tại mặt phân cách kim loại/sắt điện làm thay đổi cục bộ mật độ điện tử dẫn trong lớp LSMO:

  • Bề mặt mà $\vec{P}$ hướng tới ($IF_1$) sẽ tích tụ điện tích âm, dẫn đến sự chiếm đóng ưu tiên của orbital Mn $d_{x^2-y^2}$, thúc đẩy tương tác siêu trao đổi phản sắt từ trong mặt phẳng (pha A-AFM).
  • Bề mặt mà $\vec{P}$ hướng ra xa ($IF_2$) bị suy giảm điện tích, làm orbital Mn $d_{z^2}$ chiếm ưu thế, ưu tiên pha C-AFM. Sự thay đổi này làm biến đổi mô men từ ranh giới và cho phép sử dụng điện trường ngoài để chuyển đổi trạng thái từ tính của vật liệu.

3. Cấu trúc tinh thể của hệ spin một chiều $Ca_2CuO_3$ có đặc điểm gì khác biệt so với pha Ruddleson-Popper $La_2CuO_4$?

Trong khi $La_2CuO_4$ có cấu trúc tứ giác $I4/mmm$ với các mạng phối trí bát diện $CuO_6$ tạo thành mặt phẳng hai chiều, cấu trúc trực thoi $Immm$ của $Ca_2CuO_3$ bị mất đi các nguyên tử oxy liên kết giữa các chuỗi trong mặt phẳng. Cation Cu chỉ có số phối trí 4, tạo thành các hình vuông phẳng $CuO_4$ chia sẻ góc dọc theo trục tinh thể, ngăn cách nhau bởi các lớp oxit $CaO$. Cấu trúc này triệt tiêu tương tác trao đổi giữa các chuỗi, biến $Ca_2CuO_3$ thành một hệ chuỗi spin Heisenberg một chiều ($S = 1/2$) điển hình với tương tác siêu trao đổi mạnh dọc trục chuỗi.

4. Pha tạp uranium ($Ca_2CuO_3:U_x$) gây ra những biến đổi gì đối với phổ dao động mạng Raman của vật liệu?

Phép đo tán xạ Raman cho thấy ở mẫu $Ca_2CuO_3$ thuần khiết có 2 mode hoạt động đối xứng $A_g$. Khi pha tạp uranium ($x = 0,01 \div 0,05$), các ion U thay thế làm đứt đoạn tuần hoàn của các chuỗi $-Cu-O-Cu-$, làm giảm độ dài tương quan hiệu dụng của cụm $(Cu-O)n$. Sự phá vỡ đối xứng tịnh tiến một chiều này kích hoạt các mode dao động mới, đặc biệt là mode dao động đối xứng của nguyên tử oxy dọc theo phương chuỗi ($O{sym,par}$), làm xuất hiện các đỉnh tán xạ mới có cường độ tăng tuyến tính theo nồng độ pha tạp uranium.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thùy Trang đã thực hiện thành công việc kết hợp giữa phương pháp mô phỏng cấu trúc điện tử từ nguyên lý ban đầu và các kỹ thuật thực nghiệm nhiệt, cấu trúc, quang phổ để nghiên cứu có hệ thống các hệ oxit phức hợp thấp chiều ($La_{1-x}Sr_xMnO_3/BaTiO_3$, $LaAlO_3/SrTiO_3$ và $Ca_2CuO_3:U_x$). Công trình đã làm sáng tỏ cơ chế vi mô của liên kết từ - điện qua hiệu ứng chắn tĩnh điện và trật tự orbital tại mặt phân cách dị thể, giải thích nguồn gốc giam hãm của khí điện tử hai chiều tại tiếp xúc oxit, đồng thời làm rõ mối liên hệ giữa trật tự phản sắt từ và các mode dao động phonon trong chuỗi spin một chiều. Các kết quả đạt được cung cấp những luận cứ khoa học và bộ thông số thực nghiệm chuẩn xác, đóng góp thiết thực cho định hướng thiết kế và phát triển các vật liệu oxit chức năng trong công nghệ điện tử và spintronics.