Tổng quan luận án

Vật liệu hai chiều (2D) là trọng tâm nghiên cứu của ngành khoa học vật liệu và vật lý chất rắn kể từ khi graphene được cô lập thành công. Graphene sở hữu các đặc tính vật lý vượt trội về sức bền cơ học, độ dẫn nhiệt và dẫn điện; tuy nhiên, rào cản lớn nhất của vật liệu này trong việc chế tạo linh kiện bán dẫn và chuyển mạch điện tử là độ rộng vùng cấm bằng không ($0\text{ eV}$). Nhằm khắc phục hạn chế này, các vật liệu đơn lớp thuộc nhóm IVA có cấu trúc tương tự graphene nhưng sở hữu độ rộng vùng cấm nội tại hoặc vùng cấm có thể điều chỉnh được đã được tập trung nghiên cứu, tiêu biểu là germanene (màng đơn lớp nguyên tử germanium). Germanene có cấu trúc mạng lục giác dạng tổ ong nhưng không phẳng hoàn toàn như graphene mà sở hữu độ vênh giữa hai phân lớp nguyên tử là $0{,}737\text{ \AA}$ cùng khoảng cách liên kết Ge–Ge xấp xỉ $2{,}443\text{ \AA}$. Nhờ độ vênh đặc thù này, germanene có khả năng mở vùng cấm khi đặt trong điện trường ngoài vuông góc và đạt độ linh động chất mang rất cao.

Dù đã được dự đoán lý thuyết từ năm 2009 và chế tạo thực nghiệm thành công từ năm 2014 trên nhiều loại chất nền khác nhau, các công bố khoa học liên quan đến quá trình chuyển pha của germanene vẫn còn rất hạn chế. Các yếu tố tác động đến quá trình hình thành vật liệu chưa được khảo sát một cách có hệ thống; nhiệt độ chuyển pha (nhiệt độ tinh thể hóa, nhiệt độ nóng chảy, nhiệt độ chuyển pha vô định hình) chưa được xác định; loại chuyển pha nhiệt động học chưa được làm sáng tỏ. Thêm vào đó, cấu trúc vi mô, các dạng sai hỏng mạng (khuyết tật điểm, khuyết tật Stone-Wales, biên hạt) và khả năng tồn tại trạng thái vô định hình của germanene ở cấp độ nguyên tử chưa được phân tích đầy đủ. Các nghiên cứu thực nghiệm gặp trở ngại do không thể quan sát chi tiết diễn biến vi mô trong quá trình chuyển pha ở nhiệt độ cao, trong khi các tính toán bằng phương pháp phiếm hàm mật độ (DFT) chỉ thực hiện được trên các mô hình kích thước nhỏ (vài chục nguyên tử) nên không phản ánh được đầy đủ tính chất của hệ thống kích thước lớn.

Luận án "Khảo sát cấu trúc và tính chất nhiệt động học của vật liệu hai chiều germanene bằng phương pháp mô phỏng" được thực hiện nhằm giải quyết các khoảng trống học thuật nêu trên thông qua bốn mục tiêu nghiên cứu cụ thể:

  1. Khảo sát quá trình hình thành màng tinh thể germanene từ trạng thái lỏng cân bằng thông qua quá trình làm lạnh; xác định nhiệt độ tinh thể hóa ($T_c$), loại chuyển pha và phân tích vi cấu trúc mô hình đa tinh thể thu được tại $300\text{ K}$.
  2. Khảo sát cơ chế và đặc trưng nhiệt động của quá trình nóng chảy đối với cả hai cấu hình màng germanene đơn tinh thể và đa tinh thể; xác định nhiệt độ nóng chảy ($T_m$), loại chuyển pha và vai trò của các dạng khuyết tật trong quá trình hóa lỏng.
  3. Đánh giá sự ảnh hưởng của tốc độ làm lạnh ($\gamma$) lên quá trình hình thành pha của màng germanene; xác định ngưỡng hình thành pha tinh thể và pha vô định hình; tìm nhiệt độ chuyển pha vô định hình ($T_g$) và phân tích cấu trúc của màng germanene vô định hình ở $300\text{ K}$.
  4. Khảo sát hiệu ứng kích thước hệ (số lượng nguyên tử từ $3600$ đến $10000$ nguyên tử) tác động đến vùng nhiệt độ chuyển pha, sự phân bố số phối vị và phân bố kích thước vòng cấu trúc trong quá trình tinh thể hóa.
  • Đối tượng nghiên cứu: Màng hai chiều tự do (free-standing) đơn lớp nguyên tử germanium (germanene) ở các trạng thái lỏng, đơn tinh thể, đa tinh thể và vô định hình.
  • Phạm vi nghiên cứu:
    • Không gian mô phỏng: Không gian hai chiều mở rộng có áp dụng biên tuần hoàn theo phương $x, y$ và giới hạn bằng hai biên phản xạ đàn hồi theo phương $z$ nhằm loại bỏ tương tác nền, với quy mô hệ từ $3600$, $6400$ đến $10000$ nguyên tử.
    • Nhiệt độ: Dải nhiệt độ khảo sát biến thiên liên tục từ $300\text{ K}$ (nhiệt độ phòng) đến $3000\text{ K}$ (trạng thái lỏng hoàn toàn) trong các quá trình gia nhiệt và làm lạnh đẳng áp - đẳng nhiệt (hệ NPT).

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

Luận án đã tổng hợp và hệ thống hóa các công trình nghiên cứu trong và ngoài nước về vật liệu hai chiều nói chung và germanene nói riêng theo ba nhánh trọng tâm:

Các nghiên cứu về tổng hợp thực nghiệm germanene

  • Bianco và cộng sự (2013): Khởi xướng việc tổng hợp màng hydro hóa germanane ($\text{GeH}$) thông qua chất trung gian $\text{CaGe}_2$. Hỗn hợp $\text{Ca}$ và $\text{Ge}$ tinh khiết $99{,}999%$ được ủ kín trong ống thạch anh ở nhiệt độ $1223\text{ K} - 1320\text{ K}$ trong $16 - 20\text{ giờ}$, làm lạnh về $300\text{ K}$ ($1 - 5\text{ ngày}$), sau đó khuấy trong dung dịch $\text{HCl}$ ở $223\text{ K} - 253\text{ K}$ trong $5 - 10\text{ ngày}$.
  • Farajian (2013): Sử dụng tính toán ab initio dự đoán $\text{GeH}$ có độ rộng vùng cấm xấp xỉ $3\text{ eV}$ và độ linh động điện tử bề mặt vượt trội so với germanium khối.
  • Tổng hợp trên đế Ag(111): Các thử nghiệm ban đầu cho thấy màng germanene bị kéo giãn mạng khoảng $10%$, phản ánh sự không tương thích mạng giữa màng và đế.
  • Li và cộng sự (tháng 7/2014): Lắng đọng chùm tia điện tử (electron-beam evaporator) nguyên tử $\text{Ge}$ lên đế $\text{Pt}(111)$ trong môi trường chân không siêu cao, ủ ở $873\text{ K} - 1023\text{ K}$ trong $30\text{ phút}$, xác nhận cấu trúc germanene bằng kính hiển vi quét xuyên hầm (STM) và nhiễu xạ electron.
  • Davila và cộng sự (2014): Công bố tổng hợp thành công màng germanene trên đế $\text{Au}(111)$ bằng kỹ thuật bốc bay lắng đọng, tương tự quy trình chế tạo silicene trên $\text{Ag}(111)$.
  • Bampoulis và cộng sự (2014): Quan sát sự hình thành màng germanene trên bề mặt hợp kim $\text{Ge}_2\text{Pt}$.
  • Derivaz và cộng sự (2015): Kết hợp thực nghiệm và tính toán lý thuyết phiếm hàm mật độ (DFT) chứng minh sự tồn tại của germanene trên đế nhôm $\text{Al}(111)$.
  • Zhang và cộng sự (2016): Lắng đọng thành công germanene trên đế vật liệu 2D $\text{MoS}_2$.

Các nghiên cứu về tính chất vật lý và tiềm năng ứng dụng

  • Q. Pang và cộng sự (2011): Tính toán DFT trên dải nano germanene biên ghế bành (armchair) và dải zigzag pha tạp, chỉ ra tiềm năng ứng dụng trong tranzito hiệu ứng trường (FET), điện trở vi sai âm (NDR) và bộ lọc spin (SF).
  • Houssa và cộng sự (2011): Dự đoán màng hydro hóa germanane có vùng cấm trực tiếp trung bình khoảng $3{,}2\text{ eV}$, phù hợp cho thiết bị quang điện tử vùng phổ xanh lam/tím.
  • Chen Si và cộng sự (2014): Nghiên cứu hệ vật liệu $\text{GeI}$ với độ rộng vùng cấm $0{,}3\text{ eV}$ và độ linh động điện tử cao.
  • Xue-Sheng Ye và cộng sự (2014): Kết hợp phương trình Boltzmann và phép gần đúng thời gian hồi phục từ tính toán ban đầu, tính được độ linh động chất mang nội tại (ICM) của tấm germanene đạt tới $6 \times 10^5\text{ cm}^2\text{V}^{-1}\text{s}^{-1}$, cao hơn đáng kể so với graphene.
  • Liang và cộng sự: Xác định tính chất bán dẫn phụ thuộc vào nguyên tố chức hóa bề mặt: $\text{Ge}_2\text{H}$ mang từ tính cao, trong khi $\text{Ge}_2\text{F}$ và $\text{Ge}_2\text{Cl}$ có tính chất chống sắt từ.
  • Minglei Sun và cộng sự: Khảo sát pha tạp kim loại chuyển tiếp dãy $3d$ ($\text{Sc}-\text{Ni}$), phát hiện từ tính xuất hiện khi pha tạp $\text{V}, \text{Cr}, \text{Mn}, \text{Fe}, \text{Co}$.
  • Nghiên cứu hấp phụ khí độc: Khảo sát màng germanene biên armchair hấp phụ $\text{NH}_3, \text{SO}_2, \text{NO}_2$; trong đó phân tử $\text{NO}_2$ có năng lượng hấp phụ thấp nhất ($273{,}72\text{ meV}$, cấu hình loại B) và mức độ truyền điện tích lớn nhất.
  • X. Chen và cộng sự: Cấu trúc dị thể hai chiều germanene/antimonene liên kết thông qua lai hóa quỹ đạo, độ rộng vùng cấm trực tiếp có thể điều chỉnh linh hoạt bằng hướng và cường độ điện trường ngoài.
  • P. Kaloni (2014): Nghiên cứu sự hấp phụ của kim loại $3d$ ($\text{Sc}-\text{Zn}$) trên germanene; tổng mômen từ biến thiên từ $0{,}97\text{ }\mu_B$ đến $4{,}95\text{ }\mu_B$ (từ $\text{Sc}$ đến $\text{Mn}$); $\text{Sc}$ mở vùng cấm $74\text{ meV}$; $\text{Ti}, \text{Cr}, \text{Mn}, \text{Fe}, \text{Co}$ cho trạng thái kim loại; trường hợp hấp phụ $\text{V}$ gây ra hiệu ứng Hall dị thường lượng tử.
  • Q. Pang và cộng sự: Khảo sát hấp phụ kim loại kiềm ($\text{Li}, \text{Na}, \text{K}$); tương tác mạnh do cấu trúc vênh, vị trí Dirac dịch chuyển xuống dưới mức Fermi thể hiện đặc tính pha tạp loại n.
  • Srimathi và cộng sự: Khảo sát ống nano germanene hấp phụ các hợp chất hữu cơ dễ bay hơi (VOC) thoát ra từ các giai đoạn chín của quả chuối để ứng dụng làm cảm biến kiểm tra độ tươi thực phẩm.
  • F. Bechstedt và cộng sự: Tính toán hàm điện môi quang học, chứng minh germanene có độ hấp thụ hồng ngoại tần số thấp tương tự graphene.
  • Tính chất cơ học và lưu trữ năng lượng:
    • Kaloni và cộng sự: Mạng germanene ổn định cơ học và duy trì tần số phonon dương dưới biến dạng kéo hai trục lên đến $16%$.
    • B. Mortazavi và cộng sự: Đánh giá đáp ứng ứng suất - biến dạng đơn trục; khảo sát khả năng lưu trữ ion kiềm, tính toán dung lượng sạc lý thuyết đạt $954\text{ mAh}\cdot\text{g}^{-1}$ cho silicene, $369\text{ mAh}\cdot\text{g}^{-1}$ cho germanene và $226\text{ mAh}\cdot\text{g}^{-1}$ cho stanene.
    • J. Zhu và cộng sự: Dị thể silicene/germanene kết hợp $\text{MgX}_2$ ($\text{X} = \text{Cl}, \text{Br}, \text{I}$) làm điện cực pin Li-ion; trên $\text{MgCl}_2$, dung lượng lưu trữ $\text{Li}$ đạt lần lượt $68%$ và $94%$ so với màng tự do.
    • P. Panigrahi và cộng sự: Ứng dụng germanene làm vật liệu neo giữ các phân tử trung gian lithium polysulfides ($\text{LPS}$) trong pin $\text{Li}-\text{S}$ với năng lượng liên kết $1{,}8\text{ eV} - 2{,}6\text{ eV}$, khắc phục hiệu ứng con thoi.
  • Ứng dụng y sinh:
    • Chan Feng và cộng sự: Chế tạo hạt nano germanene kết hợp hydrogel ($\text{Ge@hydrogel}$) hỗ trợ điều trị phẫu thuật ung thư và diệt khuẩn nhờ kích hoạt nhiệt cục bộ bởi tia cận hồng ngoại (NIR).
    • Ge và cộng sự: Tấm nano germanene phân hủy sinh học điều chế từ $\text{CaGe}_2$ cho hiệu suất chuyển đổi quang nhiệt đạt $38{,}4%$ ở bước sóng $1064\text{ nm}$ trong cửa sổ sinh học NIR-II.

Các nghiên cứu tính chất nhiệt động của vật liệu 2D tương tự

  • GS. Võ Văn Hoàng và cộng sự: Mô phỏng Động lực học phân tử quá trình nóng chảy của graphene, xác định nhiệt độ nóng chảy $T_m = 7750\text{ K}$ thuộc chuyển pha loại I, quá trình bắt đầu từ khuyết tật Stone-Wales tiếp nối bởi khuyết tật mất nguyên tử và phát triển vùng hóa lỏng cục bộ.
  • Khảo sát silicene: Xác định nhiệt độ tinh thể hóa cấu trúc lục giác tổ ong là $1775\text{ K}$, cấu trúc vòng 4 (tetra-silicene) là $1950\text{ K}$, và nhiệt độ chuyển pha vô định hình của silicene là $1350\text{ K}$.

Khoảng trống nghiên cứu luận án giải quyết: Mặc dù các đặc tính điện tử, cơ học và quang học của germanene đã được nghiên cứu rộng rãi, nhưng cho đến trước công trình này, chưa có nghiên cứu nào làm sáng tỏ bức tranh nhiệt động học tổng thể của màng germanene tự do: chưa xác định được nhiệt độ tinh thể hóa, nhiệt độ nóng chảy của pha đơn tinh thể và đa tinh thể, loại chuyển pha nhiệt động học, điều kiện động học hình thành trạng thái vô định hình và ảnh hưởng của quy mô kích thước hệ nguyên tử. Luận án định vị là công trình mô phỏng quy mô lớn giải quyết trọn vẹn các thông số nhiệt động học cơ bản này.

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

Khung lý thuyết Động lực học phân tử (MD)

Phương pháp mô phỏng Động lực học phân tử cổ điển tính toán các tính chất cân bằng nhiệt động và sự dịch chuyển cấu trúc của hệ $N$ hạt thông qua việc tích phân số phương trình chuyển động Newton:

$$\frac{d^2 \vec{r}i}{dt^2} = \frac{\vec{F}i}{m_i} = -\frac{1}{m_i} \sum{j \neq i}^N \nabla U{ij}(\vec{r})$$

Hệ phương trình vi phân chuyển động được giải bằng thuật toán tích phân Verlet với bước thời gian mô phỏng tích phân $\Delta t \approx 10^{-15}\text{ s}$ ($1\text{ fs}$). Vận tốc các nguyên tử tuân theo phân bố Maxwell-Boltzmann, và nhiệt độ tức thời $T$ của hệ được xác định theo biểu thức:

$$T = \frac{1}{3 N k_B} \sum_{i=1}^N m_i v_i^2$$

Thế tương tác Stillinger-Weber (SW)

Trong mô phỏng vật liệu 2D nhóm IV, hai dạng hàm thế thường được áp dụng là thế Tersoff và thế Stillinger-Weber (SW). Luận án lựa chọn thế tương tác Stillinger-Weber được tham số hóa bởi Jin-Wu Jiang và cộng sự (2017) cho màng hai chiều germanene, vì thế SW mô tả chính xác hơn sự phân tán phonon và các đặc tính cấu trúc dạng vênh ($0{,}737\text{ \AA}$) so với thế Tersoff.

Biểu thức thế Stillinger-Weber gồm năng lượng tương tác hai hạt ($U_2$) thể hiện sự giãn dài liên kết và tương tác ba hạt ($U_3$) thể hiện sự uốn góc liên kết:

$$U = \sum_{i} \sum_{j > i} U_2(\vec{r}{ij}) + \sum{i} \sum_{j \neq i} \sum_{k > j} U_3(\vec{r}{ij}, \vec{r}{ik}, \theta_{ijk})$$

Trong đó:

  • Số hạng tương tác hai hạt: $$U_2(r_{ij}) = A \epsilon \left[ B \left(\frac{\sigma}{r_{ij}}\right)^p - \left(\frac{\sigma}{r_{ij}}\right)^q \right] \exp\left( \frac{\sigma}{r_{ij} - a\sigma} \right)$$
  • Số hạng tương tác ba hạt: $$U_3(r_{ij}, r_{ik}, \theta_{ijk}) = \lambda \epsilon \exp\left[ \frac{\gamma\sigma}{r_{ij} - a\sigma} + \frac{\gamma\sigma}{r_{ik} - a\sigma} \right] (\cos\theta_{ijk} - \cos\theta_0)^2$$
Tham số $\epsilon\text{ (eV)}$ $\sigma\text{ (\AA)}$ $a$ $\lambda\text{ (eV)}$ $\gamma$ $\cos\theta_0$ $A$ $B$ $p$ $q$ $tol$
Giá trị $1{,}930$ $2{,}181$ $1{,}190$ $19{,}570$ $3{,}205$ $-0{,}380$ $7{,}049$ $19{,}500$ $4$ $0$ $0{,}0$

Điều kiện biên và công cụ phân tích

  • Điều kiện biên: Áp dụng kết hợp điều kiện biên tuần hoàn (Periodic Boundary Conditions - PBC) theo hai trục $x$ và $y$ để mô hình co giãn thể tích tự nhiên, cùng điều kiện biên phản xạ đàn hồi (Reflective Wall) cố định theo phương $z$ nhằm giữ màng vật liệu trong một mặt phẳng giới hạn mà không chịu tương tác bám dính từ chất nền.
  • Phần mềm thực hiện:
    • LAMMPS: Thực hiện toàn bộ tính toán mô phỏng Động lực học phân tử ở hệ điều kiện đẳng áp - đẳng nhiệt (NPT ensemble).
    • ISAACS: Phân tích thống kê phân bố kích thước vòng mạng theo tiêu chí con đường ngắn nhất của Guttman.
    • VMD (Visual Molecular Dynamics): Trực quan hóa cấu hình nguyên tử và biểu diễn phân bố vector độ dịch chuyển nguyên tử (Atomic Displacements - AD).
    • Origin & Fortran: Vẽ đồ thị nhiệt động học; lập trình tính toán hàm phân bố xuyên tâm $g(r)$, số phối vị $Z$, phân bố khoảng cách liên kết, phân bố góc liên kết và tạo lập cấu trúc tinh thể ban đầu.

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

Chương 1: Tổng quan

Chương 1 hệ thống hóa cơ sở dữ liệu về lịch sử phát triển vật liệu hai chiều từ graphene đến các vật liệu đơn lớp nhóm IVA như silicene, germanene, stanene. Phân tích chi tiết các quy trình tổng hợp thực nghiệm màng germanene trên các bề mặt nền kim loại ($\text{Pt}, \text{Au}, \text{Al}, \text{Ag}, \text{Ge}_2\text{Pt}$) và nền bán dẫn ($\text{MoS}_2$). Tổng hợp toàn diện các nghiên cứu tính toán về cấu trúc điện tử, tính chất quang, độ linh động chất mang, khả năng lưu trữ ion kim loại trong pin năng lượng và tiềm năng ứng dụng trong cảm biến hóa học cũng như y sinh học quang nhiệt. Đồng thời, chương này tổng kết các phát hiện về tính chất nhiệt động của graphene và silicene, chỉ rõ khoảng trống dữ liệu về chuyển pha của germanene.

Chương 2: Đối tượng và phương pháp tính toán

Chương 2 trình bày chi tiết cơ sở toán học của mô phỏng MD, phương trình thế Stillinger-Weber, các điều kiện biên hỗn hợp tuần hoàn - phản xạ đàn hồi và quy trình 4 giai đoạn khảo sát mô phỏng:

  1. Mô hình khảo sát tinh thể hóa: Cấu hình $10^4$ nguyên tử $\text{Ge}$, hồi phục ở $3000\text{ K}$ ($10^6$ bước MD), làm lạnh về $300\text{ K}$ với tốc độ $\gamma = 10^{10}\text{ K}\cdot\text{s}^{-1}$, hồi phục $10^6$ bước ở $300\text{ K}$.
  2. Mô hình khảo sát nóng chảy: Khảo sát đồng thời màng đơn tinh thể hoàn hảo ($10^4$ nguyên tử, khoảng cách $2{,}443\text{ \AA}$, độ vênh $0{,}737\text{ \AA}$) và màng đa tinh thể chứa khuyết tật; hồi phục $10^6$ bước tại $300\text{ K}$, đun nóng lên $3000\text{ K}$ với tốc độ $10^{11}\text{ K}\cdot\text{s}^{-1}$.
  3. Mô hình ảnh hưởng tốc độ làm lạnh: Khởi tạo mô hình lỏng $10^4$ nguyên tử ở $3000\text{ K}$ (kích thước $L_x = 338{,}2097\text{ \AA}$, $L_y = 196{,}0189\text{ \AA}$, $L_z = 0{,}7370\text{ \AA}$, mật độ $4{,}18 \times 10^{-6}\text{ atom}\cdot\text{\AA}^{-3}$); làm lạnh về $300\text{ K}$ với ba tốc độ: $10^{11}, 10^{12}$ và $10^{13}\text{ K}\cdot\text{s}^{-1}$; áp dụng bán kính cắt $R_c$ tương ứng là $2{,}80\text{ \AA}$, $2{,}90\text{ \AA}$ và $3{,}00\text{ \AA}$.
  4. Mô hình ảnh hưởng kích thước: So sánh các cấu hình kích thước khác nhau ($3600$, $6400$ và $10000$ nguyên tử) trong quá trình làm lạnh $\gamma = 10^{10}\text{ K}\cdot\text{s}^{-1}$.

Chương 3: Kết quả

Chương 3 trình bày các kết quả mô phỏng thu được, được chia thành 4 nội dung chuyên sâu:

1. Quá trình hình thành màng tinh thể vật liệu germanene

  • Đặc trưng nhiệt động học: Đường phụ thuộc năng lượng toàn phần theo nhiệt độ trong quá trình làm lạnh từ $3000\text{ K}$ xuống $300\text{ K}$ giảm tuyến tính đến $1380\text{ K}$, sau đó xuất hiện bước nhảy năng lượng đột ngột trong khoảng $1380\text{ K} - 1300\text{ K}$ và tiếp tục giảm tuyến tính về $300\text{ K}$. Điểm uốn năng lượng này xác định nhiệt độ tinh thể hóa của germanene là $T_c = 1340\text{ K}$. Sự thay đổi gián đoạn của năng lượng và thể tích khẳng định quá trình tinh thể hóa là chuyển pha loại I.
  • Đặc trưng cấu trúc tại $300\text{ K}$: Hàm phân bố xuyên tâm $g(r)$ thể hiện các đỉnh nhọn đặc trưng của trật tự xa. Màng thu được có cấu trúc mạng lục giác tổ ong chủ đạo với đa số nguyên tử có số phối vị $Z = 3$, góc liên kết tập trung quanh $120^\circ$ và kích thước vòng lục giác $R = 6$ chiếm ưu thế. Tuy nhiên, hệ tồn tại các khuyết tật Stone-Wales (SW - cặp vòng 5-7 hình thành do sự xoay liên kết $90^\circ$), khuyết tật mất nguyên tử (đơn khoảng trống SVs, đôi khoảng trống DVs) và chuỗi vòng 5-7 liên kết tạo thành biên hạt (Grain Boundary - GB), cấu thành trạng thái đa tinh thể germanene.

2. Quá trình nóng chảy của germanene đơn tinh thể và đa tinh thể

  • Nhiệt độ nóng chảy: Đường cong năng lượng và nhiệt dung riêng phụ thuộc nhiệt độ cho thấy nhiệt độ nóng chảy của hai mô hình có sự khác biệt rõ rệt:
    • Màng đơn tinh thể germanene nóng chảy tại $T_m = 1670\text{ K}$.
    • Màng đa tinh thể germanene nóng chảy tại $T_m = 1540\text{ K}$.
  • Cơ chế chuyển pha: Cả hai mô hình đều thể hiện bước nhảy năng lượng đặc trưng của chuyển pha loại I. Sự chênh lệch nhiệt độ nóng chảy ($130\text{ K}$) xuất phát từ năng lượng tự do dư thừa tại các vị trí khuyết tật trong màng đa tinh thể.
  • Phân tích độ dịch chuyển nguyên tử (AD): Phân tích vector độ dịch chuyển chỉ ra rằng quá trình nóng chảy không diễn ra đồng nhất trên toàn bộ bề mặt. Ở màng đơn tinh thể, sự nóng chảy khởi đầu từ các khuyết tật Stone-Wales ngẫu nhiên được kích hoạt nhiệt, tạo thành các tâm mất trật tự cục bộ rồi lan rộng. Ở màng đa tinh thể, quá trình nóng chảy bắt đầu tại các biên hạt (nơi tập trung khuyết tật vòng 5-7) rồi phá hủy dần các hạt tinh thể lân cận.
Mô hình Nhiệt độ chuyển pha ($T_t$) Loại chuyển pha Khuyết tật khởi đầu chủ đạo Cấu trúc thu được tại $300\text{ K}$
Tinh thể hóa ($\gamma = 10^{10}\text{ K/s}$) $T_c = 1340\text{ K}$ Loại I Vòng 5-7 (Stone-Wales), SVs, DVs Đa tinh thể dạng tổ ong
Nóng chảy đơn tinh thể $T_m = 1670\text{ K}$ Loại I Khuyết tật Stone-Wales do nhiệt Trạng thái lỏng 2D
Nóng chảy đa tinh thể $T_m = 1540\text{ K}$ Loại I Biên hạt, khuyết tật Stone-Wales Trạng thái lỏng 2D

3. Ảnh hưởng của tốc độ làm lạnh lên sự hình thành màng - Germanene vô định hình

  • Phân hóa pha theo tốc độ:
    • Với tốc độ làm lạnh chậm ($\gamma = 10^{11}\text{ K}\cdot\text{s}^{-1}$ và $10^{12}\text{ K}\cdot\text{s}^{-1}$), hệ có đủ thời gian hồi phục để sắp xếp vào mạng tổ ong, tạo thành màng đa tinh thể với tỷ lệ phối vị $Z = 3$ cao và cấu trúc vòng $R = 6$ chiếm ưu thế. Tốc độ làm lạnh càng chậm, mật độ khuyết tật càng thấp.
    • Với tốc độ làm lạnh nhanh ($\gamma = 10^{13}\text{ K}\cdot\text{s}^{-1}$), sự giảm nhiệt độ đột ngột làm đóng băng các nguyên tử ở trạng thái mất trật tự, ngăn cản quá trình kết tinh và dẫn đến sự hình thành pha vô định hình của germanene tại $300\text{ K}$.
  • Nhiệt độ chuyển pha vô định hình: Đồ thị năng lượng và số phối vị trung bình theo nhiệt độ xác định nhiệt độ chuyển hóa thủy tinh / vô định hình là $T_g = 1210\text{ K}$ (ứng với $\gamma = 10^{13}\text{ K}\cdot\text{s}^{-1}$).
  • Cấu trúc vô định hình: Hàm $g(r)$ của pha vô định hình bị tù ở các đỉnh thứ hai và thứ ba (mất trật tự tầm xa, chỉ còn trật tự tầm gần). Phân bố kích thước vòng trải rộng từ $R = 3$ đến $R = 10$ thay vì tập trung ở $R = 6$; phổ nhiễu xạ không còn các điểm nhiễu xạ sắc nét mà chuyển thành các vành khuếch tán đồng tâm.

4. Ảnh hưởng của kích thước mô hình lên sự hình thành màng germanene

  • Vùng nhiệt độ chuyển pha: So sánh các mô hình chứa $3600, 6400$ và $10000$ nguyên tử cho thấy kích thước mô hình ảnh hưởng đáng kể đến độ rộng vùng chuyển pha. Đối với các mô hình kích thước lớn ($\ge 6400$ nguyên tử), dải nhiệt độ diễn ra quá trình chuyển pha tinh thể hóa mở rộng hơn so với mô hình nhỏ ($3600$ nguyên tử).
  • Thống kê vi mô: Sự phân bố số phối vị $Z$, phân bố kích thước vòng $R$, khoảng cách liên kết và góc liên kết chịu sự chi phối của hiệu ứng kích thước và mật độ khuyết tật biên. Khi số lượng nguyên tử đạt từ $6400$ đến $10000$ nguyên tử, các đại lượng nhiệt động và cấu trúc vi mô đạt trạng thái ổn định, triệt tiêu sai số kích thước hữu hạn của điều kiện biên tuần hoàn.

Chương 4: Kết luận và hướng phát triển đề tài

Chương 4 tóm tắt các kết quả đạt được về các điểm chuyển pha nhiệt động học, cơ chế hình thành khuyết tật và điều kiện tạo pha vô định hình của vật liệu hai chiều germanene; đồng thời đưa ra các đề xuất nghiên cứu mở rộng trong tương lai.

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

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

  1. Xác định các thông số chuyển pha nhiệt động học cơ bản: Lần đầu tiên công bố đầy đủ các giá trị nhiệt độ chuyển pha chính xác của màng tự do germanene: nhiệt độ tinh thể hóa ($T_c = 1340\text{ K}$), nhiệt độ nóng chảy của màng đơn tinh thể ($T_m = 1670\text{ K}$), nhiệt độ nóng chảy của màng đa tinh thể ($T_m = 1540\text{ K}$) và nhiệt độ chuyển pha vô định hình ($T_g = 1210\text{ K}$ ứng với $\gamma = 10^{13}\text{ K}\cdot\text{s}^{-1}$).
  2. Làm sáng tỏ bản chất loại chuyển pha: Chứng minh bằng định lượng rằng cả hai quá trình tinh thể hóa và nóng chảy của màng 2D germanene đều tuân theo quy luật của chuyển pha loại I thông qua sự thay đổi nhảy bậc của năng lượng toàn phần và thể tích riêng.
  3. Mô tả cơ chế khuyết tật và động học vi mô: Xác định khuyết tật Stone-Wales (cặp vòng 5-7) là khuyết tật cấu trúc then chốt xuất hiện đầu tiên trong quá trình kết tinh và là tâm khởi đầu cho sự mất ổn định mạng trong quá trình nóng chảy; làm sáng tỏ vai trò của chuỗi vòng 5-7 trong việc hình thành biên hạt đa tinh thể.
  4. Phát hiện và đặc trưng hóa pha germanene vô định hình: Chứng minh bằng mô phỏng sự tồn tại bền vững của trạng thái germanene vô định hình ở $300\text{ K}$ khi tốc độ làm nguội đạt ngưỡng $\ge 10^{13}\text{ K}\cdot\text{s}^{-1}$, thiết lập bộ dữ liệu cấu trúc chi tiết về phân bố vòng ($R = 3 - 10$) và hàm phân bố xuyên tâm cho pha vô định hình này.
  5. Đánh giá hiệu ứng kích thước hệ: Chỉ ra rằng quy mô hệ từ $6400$ nguyên tử trở lên là điều kiện cần thiết để mô phỏng chính xác dải nhiệt độ chuyển pha và triệt tiêu ảnh hưởng nhân tạo của điều kiện biên.

Ý nghĩa thực tiễn và kiến nghị giải pháp

  • Dữ liệu nhiệt độ tinh thể hóa ($1340\text{ K}$) và nóng chảy ($1540\text{ K} - 1670\text{ K}$) cung cấp cơ sở kỹ thuật định lượng để thiết lập khoảng nhiệt độ ủ nhiệt, lắng đọng pha hơi hóa học (CVD) hoặc bốc bay chùm tia điện tử trong thực nghiệm chế tạo germanene chất lượng cao.
  • Việc kiểm soát tốc độ làm lạnh cho phép định hướng chế tạo chủ động màng đa tinh thể ít khuyết tật (với tốc độ làm lạnh chậm $\le 10^{10}\text{ K}\cdot\text{s}^{-1}$) hoặc màng vô định hình (với tốc độ làm lạnh siêu nhanh $\ge 10^{13}\text{ K}\cdot\text{s}^{-1}$) phục vụ các ứng dụng linh kiện điện tử và chuyển hóa quang nhiệt chuyên biệt.

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

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

  • Mô hình mô phỏng được thực hiện trên hệ màng mỏng đứng tự do (free-standing) đặt giữa hai biên phản xạ đàn hồi, chưa xét đến các hiệu ứng tương tác bám dính, ứng suất trượt và truyền nhiệt tại mặt tiếp xúc với các loại chất nền thực nghiệm phổ biến (như $\text{Pt}, \text{Au}, \text{Al}, \text{MoS}_2$).
  • Phương pháp Động lực học phân tử cổ điển sử dụng thế Stillinger-Weber không trực tiếp mô tả sự thay đổi của cấu trúc vùng năng lượng điện tử và sự tái phân bố mật độ điện tích trong suốt quá trình chuyển pha.

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

  • Khảo sát quá trình chuyển pha và động lực học màng germanene có đặt trên các lớp đế đơn tinh thể khác nhau nhằm đánh giá định lượng ảnh hưởng của tương tác màng - đế lên nhiệt độ tinh thể hóa và cấu trúc khuyết tật.
  • Khảo sát các tính chất cơ nhiệt dưới tác động của biến dạng kéo/nén hai trục và đánh giá sự thay đổi tính chất nhiệt động khi pha tạp nguyên tử kim loại chuyển tiếp hoặc kim loại kiềm.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và học viên sau đại học ngành Vật lý kỹ thuật, Vật lý tính toán: Luận án là tài liệu tham khảo chi tiết về quy trình xây dựng mô hình mô phỏng MD màng hai chiều trên phần mềm LAMMPS, kỹ thuật thiết lập điều kiện biên hỗn hợp (PBC kết hợp phản xạ đàn hồi) và phương pháp phân tích vi cấu trúc mạng bằng thuật toán vòng ISAACS.
  • Đối với các nhà khoa học thực nghiệm trong lĩnh vực Khoa học vật liệu và Công nghệ bán dẫn: Cung cấp dải thông số nhiệt độ giới hạn cho các quá trình ủ nhiệt và xử lý nhiệt màng mỏng germanene, hỗ trợ tối ưu hóa quy trình tổng hợp vật liệu đơn lớp và màng nano bán dẫn.

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

1. Nhiệt độ tinh thể hóa của vật liệu germanene là bao nhiêu và thuộc loại chuyển pha nào?

Theo kết quả mô phỏng làm lạnh từ trạng thái lỏng ở tốc độ $\gamma = 10^{10}\text{ K}\cdot\text{s}^{-1}$, nhiệt độ tinh thể hóa của germanene được xác định là $T_c = 1340\text{ K}$ (vùng chuyển pha xảy ra trong khoảng $1380\text{ K} - 1300\text{ K}$). Quá trình này là chuyển pha loại I do có sự thay đổi gián đoạn nhảy bậc của năng lượng toàn phần và thể tích.

2. Nhiệt độ nóng chảy của màng germanene đơn tinh thể và đa tinh thể khác nhau như thế nào?

Nhiệt độ nóng chảy của màng đơn tinh thể là $T_m = 1670\text{ K}$, trong khi màng đa tinh thể nóng chảy ở nhiệt độ thấp hơn là $T_m = 1540\text{ K}$. Sự chênh lệch $130\text{ K}$ này do màng đa tinh thể chứa sẵn các khuyết tật mạng và biên hạt, làm giảm độ bền liên kết và kích hoạt sự nóng chảy sớm hơn. Cả hai quá trình đều thuộc chuyển pha loại I.

3. Điều kiện động học nào dẫn đến sự hình thành pha germanene vô định hình?

Khi tốc độ làm lạnh đạt mức rất nhanh $\gamma = 10^{13}\text{ K}\cdot\text{s}^{-1}$, các nguyên tử $\text{Ge}$ bị đông cứng ở trạng thái mất trật tự và không kịp kết tinh thành mạng tổ ong, tạo ra màng germanene vô định hình tại $300\text{ K}$. Nhiệt độ chuyển pha vô định hình tương ứng được xác định là $T_g = 1210\text{ K}$.

4. Dạng khuyết tật nào đóng vai trò chủ đạo trong quá trình tinh thể hóa và nóng chảy của germanene?

Khuyết tật Stone-Wales (cặp vòng 5-7 tạo bởi sự xoay $90^\circ$ của liên kết $\text{Ge}-\text{Ge}$) là dạng khuyết tật chủ đạo. Trong quá trình kết tinh, sự liên kết của các chuỗi vòng 5-7 tạo nên biên hạt của màng đa tinh thể; trong quá trình nóng chảy, khuyết tật Stone-Wales xuất hiện đầu tiên do kích hoạt nhiệt và đóng vai trò tâm khởi phát phá hủy mạng tinh thể.

5. Hàm thế tương tác và các phần mềm chính nào được tác giả sử dụng trong mô phỏng?

Tác giả sử dụng hàm thế tương tác ba hạt Stillinger-Weber (tham số hóa bởi Jin-Wu Jiang và cộng sự, 2017). Các phần mềm chính bao gồm: LAMMPS (chạy mô phỏng MD hệ NPT), ISAACS (phân tích phân bố vòng cấu trúc theo tiêu chí Guttman), VMD (trực quan hóa cấu hình và độ dịch chuyển nguyên tử), Origin (vẽ đồ thị nhiệt động) và các chương trình viết bằng ngôn ngữ Fortran (tạo cấu hình, tính số phối vị và hàm phân bố xuyên tâm).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Hoàng Giang đã ứng dụng thành công phương pháp mô phỏng Động lực học phân tử cổ điển với thế Stillinger-Weber để khảo sát cấu trúc vi mô và tính chất nhiệt động học của vật liệu hai chiều germanene. Nghiên cứu đã xác định chính xác các thông số nhiệt động cơ bản bao gồm nhiệt độ tinh thể hóa ($1340\text{ K}$), nhiệt độ nóng chảy của màng đơn tinh thể ($1670\text{ K}$) và đa tinh thể ($1540\text{ K}$), cùng nhiệt độ chuyển pha vô định hình ($1210\text{ K}$), đồng thời khẳng định bản chất chuyển pha loại I của hệ. Bên cạnh đó, luận án đã làm sáng tỏ cơ chế hình thành và vai trò của các khuyết tật Stone-Wales, biên hạt cũng như tác động của tốc độ làm nguội và kích thước hệ lên cấu trúc màng. Những kết quả thu được bổ sung các dữ liệu học thuật quan trọng cho ngành vật lý màng mỏng 2D và định hướng trực tiếp cho các quy trình chế tạo thực nghiệm vật liệu germanene trong công nghệ bán dẫn.