Tổng quan luận án

Luận án tiến sĩ "Tính chất cơ lý của vật liệu sắt điện ở kích thước nano mét" do nghiên cứu sinh Trần Thế Quang thực hiện dưới sự hướng dẫn khoa học của PGS. Đỗ Văn Trường tại Trường Đại học Bách khoa Hà Nội (chuyên ngành Cơ học, mã số 9440109, bảo vệ năm 2022). Luận án thuộc hướng tiếp cận tính toán mô phỏng số nhiều quy mô, kết hợp giữa tính toán lượng tử nguyên lý đầu và mô phỏng nguyên tử sử dụng mô hình thế năng vỏ - lõi.

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

Vật liệu gốm sắt điện có cấu trúc perovskite $ABO_3$ (như $\text{BaTiO}_3$, $\text{PbTiO}_3$, $\text{Pb}(\text{Zr}x\text{Ti}{1-x})\text{O}3$ - PZT, $\text{PbLaZrTiO}3$ - PLZT, $\text{Pb}(\text{Mg}{1/3}\text{Nb}{2/3})\text{O}_3$ - PMN) sở hữu các đặc tính phân cực tự phát, hiệu ứng áp điện, nhiệt điện và quang điện. Trong đó, Chì Titanate ($\text{PbTiO}_3$ - PTO) có phân cực tự phát lớn, nhiệt độ chuyển pha Curie ($T_c$) cao và năng lượng chuyển đổi phân cực thấp, là ứng viên hàng đầu trong công nghệ chế tạo bộ nhớ truy cập ngẫu nhiên sắt điện không mất dữ liệu khi ngắt nguồn (NVFRAM/FRAM) và các vi hệ thống cơ điện tử (MEMS).

Xu hướng phát triển của linh kiện điện tử đòi hỏi việc thu nhỏ kích thước hình học xuống dưới 100 nm (kích thước nano mét). Tuy nhiên, khi độ dày lớp vật liệu sắt điện giảm xuống dưới 100 nm, do xuất hiện hiệu ứng bề mặt và hiệu ứng kích thước, các xoáy phân cực (polarization vortices) xuất hiện tự phát. Các cấu trúc xoáy này làm triệt tiêu các miền phân cực đồng hướng $180^\circ$ ($180^\circ$ domain), dẫn đến việc mất khả năng lưu trữ dữ liệu theo phương pháp ghi nhị phân truyền thống dựa trên sự đảo hướng phân cực đồng nhất. Đây là rào cản vật lý lớn trong việc gia tăng mật độ linh kiện của bộ nhớ FRAM.

Về mặt phương pháp luận, các nghiên cứu tính toán lượng tử nguyên lý đầu (First-principles/DFT) dù có độ chính xác cao nhưng bị giới hạn về kích thước mô hình (thường chỉ mô phỏng được hệ vài chục đến vài trăm nguyên tử) và tính toán tại 0 K. Ngược lại, phương pháp mô phỏng nguyên tử trên cơ sở mô hình thế năng vỏ - lõi (Core-shell model) cho phép mở rộng quy mô tính toán lên hàng nghìn đến hàng vạn nguyên tử ở các điều kiện nhiệt độ và biến dạng khác nhau, nhưng đòi hỏi một bộ thông số hàm thế năng chính xác khớp với dữ liệu lượng tử và thực nghiệm của $\text{PbTiO}_3$.

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

Luận án xác định 2 mục tiêu cụ thể:

  1. Xác định quy luật ảnh hưởng của biến dạng cơ học, nhiệt độ, kích thước, và ảnh hưởng đồng thời của biến dạng và nhiệt độ đến tính phân cực của vật liệu sắt điện $\text{PbTiO}_3$.
  2. Nghiên cứu cách tạo và điều khiển xoáy phân cực đơn trong sợi nano và hạt nano sắt điện $\text{PbTiO}_3$ nhằm ứng dụng vào lưu trữ dữ liệu, giảm kích thước vật lý và tăng dung lượng bộ nhớ.

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

  • Đối tượng nghiên cứu: Vật liệu gốm sắt điện Chì Titanate ($\text{PbTiO}_3$).
  • Phạm vi nghiên cứu:
    • Về không gian cấu trúc: Nghiên cứu từ cấu trúc ô mạng đơn vị (unit cell), cấu trúc khối ($3a \times 3a \times 3c$), đến cấu trúc màng nano ($14a \times 10c$ với các độ dày từ $5c$ đến $20c$), sợi nano ($20a \times 5b \times 20c$) và hạt nano.
    • Về yếu tố tác động: Khảo sát trường nhiệt độ ($T$ từ 0 K đến nhiệt độ chuyển pha), trường biến dạng đơn trục ($\varepsilon_{zz}$), biến dạng đồng thời hai trục ($\varepsilon_{xx} = \varepsilon_{yy}$), biến dạng cắt ($\varepsilon_{xy}$), tác động của điện trường ngoài ($E$) và các khuyết tật hình học (vết nứt tại tâm, vết nứt cạnh biên đối xứng/không đối xứng, sai lệch tường miền phân cực $180^\circ$ DW).

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

Lịch sử nghiên cứu và các công trình tiền đề

Luận án điểm lại tiến trình phát triển của vật lý và cơ học vật liệu sắt điện:

  • Hiện tượng sắt điện ban đầu được Valasek (1921, 1922) phát hiện trên muối Rochelle ($\text{KNaC}_4\text{H}_4\text{O}_6 \cdot 4\text{H}_2\text{O}$), qua đó đo được đường cong điện trễ phân cực - điện trường ($P-E$).
  • Về cấu trúc vi mô và tường miền phân cực: Känzig (1957), Megaw (1957), Jona (1962), Merz (1954), Lines & Glass (2001), Arlt (1990) đã thiết lập nền tảng về cấu trúc tinh thể perovskite, miền phân cực (domain) và các tường miền phân cực $180^\circ$ (180° DW), 90° (90° DW).
  • Về ứng dụng bộ nhớ sắt điện: Scott (1989, 2005), Paz De Araujo (1989, 1996), Ishiwara (2004), Arimoto (2004) làm rõ nguyên lý lưu trữ thông tin dựa trên hai trạng thái phân cực dư ($+P_r$ và $-P_r$) trong ô nhớ FRAM, cũng như các giới hạn khi thu nhỏ kích thước.
  • Về cấu trúc xoáy phân cực và hiệu ứng kích thước: Fu & Bellaiche (2003, 2004), Bellaiche (2005, 2006, 2007), Naumov (2004), Rodriguez (1999), Stachiotti (2001, 2004), Du (2009, 2010), Kitamura (2008, 2010, 2011, 2012, 2016, 2017), Yadav (2016), Tang (2015), Li (2017) đã phát hiện và phân tích các dạng xoáy phân cực (polarization vortices, flux-closure domains) trong các chấm lượng tử, dây nano, màng siêu mỏng của $\text{BaTiO}_3$, $\text{BiFeO}_3$ và $\text{PbTiO}_3$.
  • Về phương pháp mô phỏng: Cohen (1992, 2000, 2002, 2004, 2011), Krakauer (1990, 1992, 1998), Rabe (2005), Vanderbilt (1990, 2002, 2007), Migoni (2000, 2001, 2005, 2006), Sepliarsky (2011) đã phát triển các phương pháp phiếm hàm mật độ (DFT) và xây dựng mô hình vỏ - lõi cho vật liệu perovskite.
  • Tại Việt Nam: Các nghiên cứu của Vũ Ngọc Hùng (2017), Ngô Đức Quân (2016) và các nhóm tại Viện Khoa học Vật liệu, Đại học Bách khoa Hà Nội, Đại học Công nghệ (ĐHQGHN) đã khảo sát tính chất điện môi, áp điện trên các hệ gốm $\text{BaTiO}_3$, $\text{BiFeO}_3$, $\text{CoFe}_2\text{O}_4-\text{PZT}$, $\text{Pb}(\text{Zr,Ti})\text{O}_3$ nhưng chủ yếu tiếp cận ở dạng khối hoặc màng dày và pha tạp.

Khoảng trống nghiên cứu được luận án giải quyết

  1. Thiếu một bộ thông số hoàn chỉnh và chính xác cho hàm thế năng vỏ - lõi của $\text{PbTiO}_3$ có khả năng mô tả đồng thời cả cấu trúc mạng tinh thể, hằng số đàn hồi và sự biến đổi của đường cong điện trễ dưới tác động đa trường (cơ - nhiệt - điện).
  2. Chưa có quy trình cụ thể để khởi tạo và điều khiển trạng thái định hướng của một xoáy phân cực đơn (đơn xoáy) trong cấu trúc nano $\text{PbTiO}_3$ bằng điện trường ngoài phục vụ chuyển đổi trạng thái bit dữ liệu.
  3. Chưa có các phân tích định lượng về độ ổn định của xoáy phân cực đơn dưới ảnh hưởng của biến dạng cơ học lớn và các dạng khuyết tật hình học cụ thể như vết nứt và sự lệch vị trí tường miền.

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

Khung lý thuyết và mô hình tính toán

Luận án kết hợp hai phương pháp mô phỏng số trên thang vi mô và nguyên tử:

+-------------------------------------------------------------+
|               Tính toán Nguyên lý đầu (DFT)                 |
|   - Phần mềm Quantum ESPRESSO (QE)                          |
|   - Xác định: Hằng số mạng (a, c), tọa độ nguyên tử,        |
|     các hằng số đàn hồi C_ij (C11, C12, C13, C33, C44, C66) |
+------------------------------+------------------------------+
                               |
                               v (Dữ liệu chuẩn làm khớp)
+-------------------------------------------------------------+
|            Mô hình Thế năng Vỏ - Lõi (Core-Shell)           |
|   - Tối ưu hóa hàm thế bằng thuật toán BFGS trong GULP      |
|   - Tương tác: Tầm gần (Buckingham), Tầm xa (Coulomb),      |
|     Tương tác Vỏ - Lõi phi điều hòa (k2, k4)                |
+------------------------------+------------------------------+
                               |
                               v
+-------------------------------------------------------------+
|                Khảo sát Tính chất Cơ - Lý                   |
|   - Dựng đường cong điện trễ P-E (kích thước khối)          |
|   - Tạo và điều khiển xoáy phân cực đơn (màng, sợi, hạt)    |
|   - Khảo sát biến dạng kéo/nén và khuyết tật hình học       |
+-------------------------------------------------------------+

  1. Phương pháp phiếm hàm mật độ (Density Functional Theory - DFT):

    • Dựa trên định lý Hohenberg-Kohn (1964) và hệ phương trình Kohn-Sham (1965).
    • Sử dụng phép xấp xỉ mật độ địa phương (LDA) và xấp xỉ độ dốc tổng quát (GGA) của Perdew (1985), kết hợp chuẩn thế mềm Vanderbilt (Ultrasoft Pseudopotential - USPP).
    • Lưới lấy mẫu điểm $k$ trong vùng Brillouin thứ nhất theo chuẩn Monkhorst-Pack (1976).
    • Công cụ thực thi: Phần mềm mã nguồn mở Quantum ESPRESSO (QE).

  2. Mô hình nguyên tử Vỏ - Lõi (Core-Shell Model):

    • Mô hình do Dick và Overhauser (1958) đề xuất, xem mỗi ion gồm một "lõi" (core) mang điện tích $q_c$ tập trung khối lượng nguyên tử và một "vỏ" (shell) không khối lượng mang điện tích $q_e$ biểu diễn đám mây điện tử bao quanh. Tổng điện tích ion là $q = q_c + q_e$.
    • Tổng thế năng của hệ thống gồm 3 thành phần chính: $$V_{\text{total}} = V_{SR} + V_{LR} + V_{CS}$$
      • Thế tương tác tầm ngắn (Short-Range - $V_{SR}$): Sử dụng hàm thế Buckingham: $$V_{SR}(r_{ij}) = A \exp\left(-\frac{r_{ij}}{\rho}\right) - \frac{C}{r_{ij}^6}$$ trong đó $r_{ij}$ là khoảng cách giữa các vỏ của ion $i$ và $j$; $A, \rho, C$ là các thông số tương tác.
      • Thế tương tác tầm xa (Long-Range - $V_{LR}$): Tương tác tĩnh điện Coulomb giữa tất cả các điện tích (lõi - lõi, lõi - vỏ, vỏ - vỏ): $$V_{LR}(r_{ij}) = \sum_{i \neq j} \frac{q_i q_j}{4\pi \varepsilon_0 r_{ij}}$$ sử dụng kỹ thuật tổng Ewald có xét bán kính cắt $R_c$.
      • Thế tương tác giữa vỏ và lõi trong cùng một ion ($V_{CS}$): Lò xo phi điều hòa liên kết lõi và vỏ với chuyển vị tương đối $d$: $$V_{CS}(d) = \frac{1}{2}k_2 d^2 + \frac{1}{24}k_4 d^4$$ trong đó $k_2$ là độ cứng lò xo điều hòa, $k_4$ là hệ số phi điều hòa bậc bốn.

  3. Thuật toán tối ưu hóa:

    • Thuật toán Broyden-Fletcher-Goldfarb-Shanno (BFGS) được áp dụng để làm khớp đa mục tiêu, cực tiểu hóa hàm thế và tìm bộ thông số tương tác trong phần mềm GULP (General Utility Lattice Program).

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

Chương 1. Tổng quan về vật liệu sắt điện và phương pháp mô phỏng

Chương này hệ thống hóa các khái niệm cơ sở về vật liệu sắt điện perovskite $ABO_3$, cấu trúc tinh thể của $\text{PbTiO}_3$, cơ chế chuyển pha và các công cụ mô phỏng.

  • Cấu trúc tinh thể và chuyển pha: Ở nhiệt độ cao ($T > T_c = 763\text{ K}$), $\text{PbTiO}_3$ ở pha lập phương (Cubic) có tâm đối xứng, không có tính sắt điện (pha thuận điện). Khi hạ nhiệt độ xuống dưới $T_c$, ô mạng chuyển sang cấu trúc lăng trụ tứ giác đều (Tetragonal), ion $\text{Ti}^{4+}$ và $\text{Pb}^{2+}$ dịch chuyển tương đối so với mạng ion $\text{O}^{2-}$, phá vỡ tâm đối xứng và tạo ra phân cực tự phát ($P_s$) dọc theo trục [001].
  • Miền phân cực và tường miền: Cấu trúc sắt điện hình thành các miền phân cực (domain) đồng hướng để giảm thiểu năng lượng tĩnh điện và cơ học. Giữa các miền phân cực là các tường miền $180^\circ$ ($180^\circ\text{ DW}$) hoặc $90^\circ$ ($90^\circ\text{ DW}$).
  • Đường cong điện trễ ($P-E$): Dưới điện trường ngoài biến thiên tuần hoàn, vật liệu biểu diễn vòng lặp trễ đặc trưng bởi các đại lượng: độ phân cực bão hòa ($P_s$), độ phân cực dư ($P_r$) khi ngắt điện trường, và trường điện kháng ($E_c$) cần thiết để đưa độ phân cực về 0.
  • Ứng dụng trong bộ nhớ: Chương 1 phân tích cấu trúc hoạt động của bộ nhớ NVFRAM, bóng bán dẫn hiệu ứng trường sắt điện (FeFET), so sánh hiệu năng với DRAM, SRAM, Flash, làm rõ ưu thế về tốc độ ghi xóa, số chu kỳ đọc ghi cao và tính không mất dữ liệu.
Loại bộ nhớ Tốc độ đọc / ghi Điện áp hoạt động Độ bền chu kỳ đọc/ghi Tính không mất dữ liệu
SRAM Rất nhanh (<10 ns) Thấp (1 - 3.3 V) Vô hạn ($>10^{16}$) Mất dữ liệu khi ngắt điện
DRAM Nhanh (~20 - 50 ns) Thấp (1.2 - 1.8 V) Vô hạn ($>10^{16}$) Mất dữ liệu (cần làm tươi)
Flash (NAND/NOR) Đọc nhanh, ghi chậm ($\mu\text{s}-\text{ms}$) Cao (10 - 12 V khi ghi) Giới hạn ($10^4 - 10^5$) Không mất dữ liệu
FRAM (NVFRAM) Nhanh (~50 - 70 ns) Thấp (1.5 - 3.3 V) Rất cao ($10^{12} - 10^{14}$) Không mất dữ liệu

Chương 2. Xác định hàm thế năng của mô hình vỏ - lõi cho vật liệu $\text{PbTiO}_3$ và ứng dụng trong khảo sát sự phân cực điện

Chương 2 trình bày quy trình tính toán các đặc trưng cơ học của $\text{PbTiO}_3$ bằng DFT, xây dựng bộ thông số thế năng vỏ - lõi và mô phỏng đường cong điện trễ của vật liệu khối.

1. Tính toán nguyên lý đầu bằng Quantum ESPRESSO

Tác giả thực hiện khảo sát hội tụ năng lượng theo động năng cắt sóng phẳng ($E_{\text{cutwfc}}$), động năng cắt mật độ điện tích ($E_{\text{cutrho}}$) và lưới điểm $k$ ($nk_1 \times nk_2 \times nk_3$). Kết quả xác định các thông số mạng và hằng số đàn hồi của $\text{PbTiO}_3$ ở pha tứ giác (Tetragonal) làm dữ liệu chuẩn:

  • Thông số mạng: $a = b = 3.904\text{ \AA}$, $c = 4.152\text{ \AA}$, tỷ lệ co dãn $c/a = 1.063$.
  • Tọa độ nguyên tử trong ô đơn vị: $\text{Pb}\ (0, 0, 0)$; $\text{Ti}\ (0.5, 0.5, 0.538)$; $\text{O}_1\ (0.5, 0.5, 0.112)$; $\text{O}_2\ (0.5, 0, 0.631)$; $\text{O}_3\ (0, 0.5, 0.631)$.
  • Hằng số đàn hồi $C_{ij}$ (GPa): $C_{11} = 236.4$, $C_{12} = 101.5$, $C_{13} = 78.4$, $C_{33} = 104.2$, $C_{44} = 65.8$, $C_{66} = 103.1$.

2. Kết quả tối ưu hóa bộ thông số mô hình vỏ - lõi

Sử dụng GULP với giải thuật BFGS, tác giả thu được bộ thông số hàm thế năng vỏ - lõi cho $\text{PbTiO}_3$:

Ion Điện tích lõi $q_c$ (e) Điện tích vỏ $q_e$ (e) Cặp tương tác vỏ - vỏ $A$ (eV) $\rho$ (\AA) $C$ ($\text{eV}\cdot\text{\AA}^6$) Lò xo $k_2$ ($\text{eV}/\text{\AA}^2$) Lò xo $k_4$ ($\text{eV}/\text{\AA}^4$)
Pb $+1.84$ $+0.16$ $\text{Pb}{\text{shell}} - \text{O}{\text{shell}}$ $854.21$ $0.362$ $0.00$ $190.4$ $0.00$
Ti $+4.00$ $+0.00$ $\text{Ti}{\text{shell}} - \text{O}{\text{shell}}$ $1245.30$ $0.345$ $0.00$ $450.0$ $6500.0$
O $+0.86$ $-2.86$ $\text{O}{\text{shell}} - \text{O}{\text{shell}}$ $2150.50$ $0.285$ $24.50$ $85.6$ $42000.0$

3. Khảo sát ảnh hưởng của cơ - nhiệt đến đường cong điện trễ dạng khối

Mô hình cấu trúc khối kích thước $3a \times 3a \times 3c$ được đưa vào khảo sát:

  • Biến dạng đơn trục ($\varepsilon_{zz}$): Biến dạng kéo theo phương $z$ làm tăng độ phân cực dư $P_r$ và làm mở rộng vòng trễ (tăng trường điện kháng $E_c$). Ngược lại, biến dạng nén trục $z$ làm giảm $P_r$ và $E_c$.
  • Biến dạng đồng thời hai trục ($\varepsilon_{xx} = \varepsilon_{yy}$): Biến dạng nén hai trục trong mặt phẳng $(x, y)$ đóng vai trò tương đương như kéo trục $z$ do hiệu ứng Poisson, làm tăng mạnh độ phân cực dọc theo trục $z$.
  • Biến dạng cắt ($\varepsilon_{xy}$): Gây nghiêng ô mạng, làm suy giảm tính đối xứng phân cực, dẫn đến giảm phân cực dư theo phương $z$ và làm biến dạng hình học của vòng trễ $P-E$.
  • Ảnh hưởng nhiệt độ ($T$): Khi nhiệt độ tăng từ 0 K đến gần nhiệt độ Curie $T_c$, độ phân cực tự phát $P_s$, phân cực dư $P_r$ và trường điện kháng $E_c$ đều giảm dần theo hàm phi tuyến; khi $T \ge T_c$, vòng trễ khép lại thành một đường thẳng đi qua gốc tọa độ (trạng thái thuận điện).

Chương 3. Nghiên cứu cách tạo xoáy phân cực đơn cho vật liệu $\text{PbTiO}_3$ ở kích thước nano mét

Chương 3 tập trung giải quyết bài toán cốt lõi: tạo lập và điều khiển một xoáy phân cực đơn trong các cấu trúc có kích thước nano.

  • Hiệu ứng chiều dày màng nano: Khảo sát mô hình mặt cắt $(x, z)$ kích thước $14a \times H$ (với chiều dày $H$ thay đổi: $5c, 10c, 15c, 20c$). Kết quả cho thấy ở chiều dày lớn ($20c$), các miền phân cực $180^\circ$ phẳng chiếm ưu thế. Khi chiều dày giảm xuống dưới $10c$, tương tác khử phân cực bề mặt chiếm ưu thế làm các véc tơ phân cực tại biên uốn cong, tạo thành chuỗi các dòng xoáy phân cực xen kẽ nhau.
  • Xây dựng xoáy phân cực đơn (Single Vortex):
    • Trong cấu trúc sợi nano (nanowire) và hạt nano (nanoparticle) $\text{PbTiO}_3$, bằng cách thiết lập cấu hình ban đầu gồm một cặp miền phân cực $180^\circ$ ngược chiều nhau cạnh nhau, hệ thống tự động hồi phục về trạng thái cân bằng năng lượng tạo thành một xoáy phân cực duy nhất khép kín trong mặt phẳng $(x, z)$.
    • Tùy thuộc vào chiều phân cực ban đầu, hệ tạo ra xoáy phân cực theo chiều kim đồng hồ (Clockwise - CW) hoặc ngược chiều kim đồng hồ (Counter-Clockwise - CCW).
+-------------------------------------------------------------------------+
|                  Cơ chế Tạo và Điều khiển Xoáy Đơn                      |
|                                                                         |
|   Trạng thái Bit "0" (Xoáy CW)             Trạng thái Bit "1" (Xoáy CCW)|
|        +------>------+                          +------<------+         |
|        |             |                          |             |         |
|        ^     CW      v                          v     CCW     ^         |
|        |             |                          |             |         |
|        +------<------+                          +------>------+         |
|                                                                         |
|        [Điện trường E tác dụng]                 [Điện trường -E tác dụng|
|        Nửa trái: +z | Nửa phải: -z              Nửa trái: -z | Nửa phải:|
+-------------------------------------------------------------------------+
  • Cơ chế điều khiển bằng điện trường ngoài:
    • Để đảo chiều xoáy một cách chủ động, tác giả đề xuất áp dụng một cặp điện trường ngoài phản đối xứng (hai điện trường có hướng ngược nhau đặt vào hai nửa của hạt/sợi nano: một bên hướng lên $+z$, một bên hướng xuống $-z$).
    • Mối quan hệ giữa mô men hình xuyến phân cực (Toroidal Moment $G = \frac{1}{2V}\sum (\mathbf{r}_i \times \mathbf{p}_i)$) và điện trường ngoài $E$ biểu diễn một đường cong trễ dạng vòng xuyến ($G-E$ hysteresis loop). Qua đó, chiều xoáy CW và CCW có thể đảo nghịch tuần hoàn và duy trì ổn định sau khi ngắt điện trường.

Chương 4. Khảo sát ảnh hưởng của biến dạng và khuyết tật hình học đến xoáy phân cực đơn

Chương 4 kiểm chứng tính bền vững của xoáy phân cực đơn trong sợi nano $\text{PbTiO}_3$ (mô hình $20a \times 5b \times 20c$) dưới tác động của các điều kiện cơ học bất lợi.

1. Ảnh hưởng của biến dạng đơn trục kéo và nén dọc trục $z$ ($\varepsilon_z$)

  • Khảo sát biến dạng kéo từ $+2%$ đến $+10%$ và biến dạng nén từ $-2%$ đến $-10%$.
  • Biến dạng nén: Làm giảm tỉ lệ $c/a$, suy giảm liên kết phân cực dọc, cấu trúc xoáy trở nên dẹt và tâm xoáy dịch chuyển nhẹ nhưng cấu trúc vòng xoáy khép kín vẫn được duy trì đến mức nén $-6%$. Khi nén quá mức ($-8%$ đến $-10%$), xoáy bị phá vỡ chuyển sang trạng thái phân cực trong mặt phẳng.
  • Biến dạng kéo: Tăng cường độ lớn các véc tơ phân cực theo phương $z$, làm cấu trúc xoáy phân cực bền vững hơn dọc theo trục thẳng đứng.

2. Ảnh hưởng của khuyết tật vết nứt hình học

Tác giả khảo sát chi tiết ảnh hưởng của vị trí và kích thước vết nứt:

  • Vết nứt tại tâm (Center crack):
    • Phát triển theo phương $x$ (kích thước từ $2a \times 2c$ đến $10a \times 2c$): Khi kích thước vết nứt nhỏ ($2a \times 2c$), xoáy phân cực đi vòng qua mép vết nứt, cấu trúc xoáy tổng thể được bảo toàn. Khi chiều dài vết nứt theo phương $x$ tăng lên $\ge 6a$, dòng phân cực bị chặn lại, xoáy đơn bị phân tách thành hai xoáy nhỏ cục bộ ở hai phía vết nứt.
    • Phát triển theo phương $z$ (kích thước $2a \times 4c$ đến $2a \times 16c$): Tác động phá vỡ xoáy diễn ra chậm hơn; cấu trúc xoáy đơn duy trì được tính liên tục ở kích thước vết nứt nhỏ hơn.
  • Vết nứt tại cạnh biên (Edge crack):
    • Khảo sát vết nứt không đối xứng (một bên biên) và đối xứng (hai bên biên) theo cả hai phương $x$ và $z$.
    • Vết nứt cạnh biên phát triển theo phương $x$ làm suy giảm tính đối xứng cơ học, gây lệch tâm xoáy về phía không bị nứt. Vết nứt đối xứng làm co hẹp hành lang phân cực, làm tăng mật độ dòng xoáy tại vùng vật liệu nguyên vẹn còn lại.
Loại khuyết tật hình học Kích thước khảo sát Trạng thái cấu trúc xoáy phân cực
Vết nứt tâm theo phương $x$ $2a \times 2c$ đến $4a \times 2c$ Duy trì xoáy đơn, dòng phân cực uốn quanh vết nứt
Vết nứt tâm theo phương $x$ $6a \times 2c$ đến $10a \times 2c$ Tách thành hai xoáy nhỏ cục bộ ở hai nửa
Vết nứt tâm theo phương $z$ $2a \times 4c$ đến $2a \times 16c$ Tâm xoáy dãn dài, xoáy đơn duy trì ở mức nứt vừa phải
Vết nứt biên không đối xứng $2a \times 2c$ đến $6a \times 2c$ Tâm xoáy bị đẩy lệch sang vùng không khuyết tật
Vết nứt biên đối xứng $2a \times 2c$ đến $8a \times 2c$ Co thắt dòng phân cực; xoáy duy trì nếu hành lang giữa đủ lớn
Lệch tường miền ($180^\circ\text{ DW}$) Lệch 1, 2, 3, 4 ô mạng Gây bất đối xứng kích thước hai nửa xoáy, tâm xoáy dịch chuyển

3. Ảnh hưởng của khuyết tật sai lệch tường miền $180^\circ$ DW

Khi vị trí phân chia ban đầu của tường miền $180^\circ$ bị lệch khỏi mặt phẳng đối xứng trung tâm (khảo sát các vị trí lệch $1, 2, 3, 4$ ô đơn vị): Xoáy phân cực vẫn hình thành nhưng tâm xoáy bị dịch chuyển tỉ lệ thuận với độ lệch của tường miền, làm cho phân bố mô men lưỡng cực giữa hai nửa trở nên không đồng đều.

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

Đóng góp mới về mặt khoa học (Lý luận)

  1. Xây dựng thành công bộ thông số hàm thế năng tương tác mô hình vỏ - lõi cho vật liệu $\text{PbTiO}3$ có độ chính xác cao, khớp với các đặc tính cấu trúc mạng tinh thể ($a, c$) và tensor hằng số đàn hồi ($C{ij}$) từ tính toán lượng tử nguyên lý đầu (DFT/QE).
  2. Làm sáng tỏ quy luật ứng xử cơ - nhiệt - điện của $\text{PbTiO}_3$ dạng khối: Định lượng hóa mối quan hệ giữa các trạng thái biến dạng đơn trục, biến dạng hai trục, biến dạng cắt và nhiệt độ với các tham số đường cong điện trễ ($P_s, P_r, E_c$).
  3. Thiết lập cơ sở lý thuyết về sự chuyển tiếp pha và hình thành xoáy phân cực đơn trong cấu trúc màng siêu mỏng, sợi nano và hạt nano $\text{PbTiO}_3$.

Đóng góp mới về mặt thực tiễn (Ứng dụng)

  1. Đề xuất phương pháp lưu trữ dữ liệu mới: Đề xuất sử dụng hai chiều quay ngược nhau của một xoáy phân cực đơn trong sợi nano/hạt nano $\text{PbTiO}_3$ để mã hóa bit dữ liệu nhị phân:
    • Chiều kim đồng hồ (CW) đại diện cho bit dữ liệu "0".
    • Chiều ngược kim đồng hồ (CCW) đại diện cho bit dữ liệu "1".
  2. Giải pháp tăng dung lượng bộ nhớ FRAM: Thay thế màng sắt điện phẳng truyền thống bằng các mảng sợi nano hoặc hạt nano chứa xoáy đơn giúp khắc phục hoàn toàn vấn đề triệt tiêu phân cực do hiệu ứng bề mặt, cho phép giảm kích thước vật lý ô nhớ xuống thang nano mét và có tiềm năng tăng mật độ lưu trữ lên hàng trăm lần.
  3. Quy trình điều khiển đảo chiều bit: Chứng minh khả năng ghi/xóa dữ liệu bằng cặp điện trường ngoài phản đối xứng và xác định các ngưỡng bền vững cơ học của ô nhớ trước biến dạng kéo/nén và các dạng khuyết tật hình học (vết nứt, sai lệch biên miền).

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

Hạn chế của đề tài

  • Các kết quả trong luận án hoàn toàn thu được từ phương pháp mô phỏng số nhiều quy mô (tính toán nguyên lý đầu kết hợp mô hình vỏ - lõi trên máy tính), chưa được kiểm chứng trực tiếp bằng các phép đo thực nghiệm chế tạo linh kiện nano thực tế tại phòng thí nghiệm do giới hạn về trang thiết bị công nghệ nano.
  • Mô hình chưa xét đến ảnh hưởng phức hợp của môi trường điện cực thực tế, sự khuếch tán tạp chất hóa học hoặc sự phân bố ngẫu nhiên của các nút khuyết oxy/chì ở nhiệt độ cao.

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

Luận án chỉ ra hai định hướng phát triển cụ thể:

  1. Tiến hành thực nghiệm chế tạo mẫu vật liệu cấu trúc nano (sợi nano, hạt nano $\text{PbTiO}_3$), thực hiện tạo và đo kiểm việc điều khiển chiều xoáy phân cực đơn bằng kính hiển vi lực áp điện (PFM) hoặc kính hiển vi điện tử truyền qua độ phân giải cao (HRTEM).
  2. Khảo sát thực nghiệm các ảnh hưởng của biến dạng cơ học, nhiệt độ, các dạng khuyết tật cấu trúc tinh thể thực tế và điện trường ngoài đến tính ổn định của xoáy phân cực đơn.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Cơ học vật rắn, Cơ kỹ thuật, Vật lý chất rắn, Khoa học & Kỹ thuật vật liệu, Cơ điện tử. Có thể tham khảo phương pháp luận kết hợp tính toán DFT và mô hình thế năng vỏ - lõi (Core-shell), cũng như quy trình tối ưu hóa đa mục tiêu bằng GULP.
  • Kỹ sư và các nhà thiết kế vi mạch, linh kiện bán dẫn: Tham khảo các số liệu định lượng về ảnh hưởng của biến dạng cơ học, nhiệt độ và khuyết tật hình học đến độ bền hoạt động của màng mỏng sắt điện, làm cơ sở tối ưu hóa thiết kế ô nhớ FRAM, cảm biến áp điện MEMS và các bộ biến đổi năng lượng nano.
  • Phần đáng tham khảo nhất:
    • Chương 2: Bảng thông số hàm thế năng vỏ - lõi và quy trình xác định thông số từ DFT.
    • Chương 3 & 4: Cơ chế thiết lập xoáy phân cực đơn và đồ thị trễ $G-E$ biểu diễn quá trình chuyển đổi trạng thái xoáy bằng điện trường ngoài.

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

1. Tại sao các ô nhớ FRAM truyền thống bị mất khả năng lưu trữ khi giảm kích thước màng sắt điện xuống dưới 100 nm?

Khi chiều dày màng sắt điện $\text{PbTiO}_3$ giảm xuống dưới 100 nm, tỷ lệ diện tích bề mặt trên thể tích tăng mạnh, dẫn đến sự gia tăng của trường khử phân cực tĩnh điện tại bề mặt. Để cực tiểu hóa năng lượng toàn phần, hệ thống tự động bẻ cong các véc tơ phân cực để tạo thành các dòng xoáy phân cực khép kín (polarization vortices). Các xoáy này triệt tiêu miền phân cực đồng hướng $180^\circ$, làm mất hai trạng thái phân cực dư đồng nhất $+P_r$ và $-P_r$ vốn là cơ sở để đọc/ghi bit nhị phân truyền thống.

2. Luận án đề xuất giải pháp gì để biến nhược điểm của xoáy phân cực thành ưu điểm lưu trữ dữ liệu?

Tác giả đề xuất phương pháp lưu trữ dữ liệu trực tiếp dựa trên trạng thái định hướng của một xoáy phân cực đơn (Single Vortex) trong sợi nano hoặc hạt nano:

  • Xoáy phân cực quay theo chiều kim đồng hồ (CW) được quy ước là trạng thái bit 0.
  • Xoáy phân cực quay ngược chiều kim đồng hồ (CCW) được quy ước là trạng thái bit 1. Nhờ đó, không cần triệt tiêu xoáy phân cực mà sử dụng chính chiều quay của xoáy làm đơn vị mang thông tin, cho phép thu nhỏ ô nhớ xuống kích thước nano mét và tăng mật độ lưu trữ lên nhiều lần.

3. Quy trình điều khiển đảo chiều xoáy phân cực đơn được thực hiện như thế nào?

Việc đảo chiều xoáy phân cực giữa CW và CCW được thực hiện bằng cách áp dụng một cặp điện trường ngoài phản đối xứng tác dụng đồng thời vào hai nửa của cấu trúc nano (một nửa chịu điện trường $+E$ hướng lên, nửa còn lại chịu điện trường $-E$ hướng xuống). Mối quan hệ giữa mô men hình xuyến $G$ và điện trường $E$ tạo thành một vòng trễ $G-E$, đảm bảo tính đảo chiều có thể kiểm soát và duy trì trạng thái phân cực dư sau khi ngắt điện trường.

4. Vết nứt tại tâm sợi nano ảnh hưởng như thế nào đến sự ổn định của xoáy phân cực đơn?

  • Khi vết nứt tại tâm phát triển theo phương $x$ với kích thước nhỏ ($2a \times 2c$), xoáy phân cực uốn cong đi vòng quanh mép vết nứt, cấu trúc xoáy đơn khép kín vẫn được bảo toàn.
  • Khi chiều dài vết nứt theo phương $x$ đạt từ $6a \times 2c$ trở lên, dòng phân cực bị chặn hoàn toàn, cấu trúc xoáy đơn bị phá vỡ và phân tách thành hai xoáy phân cực nhỏ cục bộ nằm ở hai bên vết nứt.

5. Sự khác biệt giữa phương pháp DFT và mô hình thế năng vỏ - lõi trong luận án là gì?

  • Phương pháp DFT (Quantum ESPRESSO): Giải phương trình cơ học lượng tử cho hệ điện tử, có độ chính xác cao nhưng tiêu tốn tài nguyên tính toán lớn, chỉ áp dụng cho ô mạng đơn vị chứa 5 nguyên tử tại 0 K để xác định thông số mạng ($a, c$) và hằng số đàn hồi ($C_{ij}$).
  • Mô hình vỏ - lõi (GULP): Xem ion gồm lõi mang điện tích dương và vỏ đám mây điện tử liên kết qua lò xo phi điều hòa. Phương pháp này nhận các thông số từ DFT làm dữ liệu chuẩn để tối ưu hóa, cho phép mô phỏng hệ thống nano chứa hàng nghìn nguyên tử ở các mức nhiệt độ và trạng thái biến dạng thực tế với tốc độ tính toán nhanh.

Kết luận

Luận án của nghiên cứu sinh Trần Thế Quang đã giải quyết trọn vẹn bài toán xây dựng bộ thông số hàm thế năng vỏ - lõi cho vật liệu sắt điện $\text{PbTiO}_3$ từ tính toán nguyên lý đầu, đồng thời phân tích chi tiết ảnh hưởng của các trường cơ học, nhiệt độ và điện trường đến đường cong điện trễ ở dạng khối. Đóng góp nổi bật của công trình là đề xuất và kiểm chứng thành công quy trình tạo lập, điều khiển xoáy phân cực đơn trong sợi và hạt nano $\text{PbTiO}_3$ bằng cặp điện trường ngoài phản đối xứng, mở ra giải pháp cấu trúc mới cho bộ nhớ NVFRAM thế hệ tiếp theo. Các phân tích về biến dạng kéo/nén và khuyết tật hình học đã cung cấp luận cứ khoa học về độ bền và giới hạn làm việc của cấu trúc xoáy phân cực ở kích thước nano mét.