Tổng quan luận án

Nhu cầu về các nguồn năng lượng dự trữ nhỏ gọn, dung lượng lớn, tuổi thọ cao và an toàn cho các thiết bị di động (điện thoại thông minh, máy tính bảng, máy tính xách tay) cũng như phương tiện giao thông chạy điện (xe đạp điện, xe máy điện, ô tô điện) ngày càng gia tăng. Trong cấu trúc ắc quy ion lithium, cực âm (anode) là bộ phận đóng vai trò then chốt quyết định dung lượng tích trữ, tốc độ sạc và tuổi thọ chu kỳ. Silicon (Si) là vật liệu cực âm tiềm năng với dung lượng lý thuyết lên tới khoảng 4200 mAh/g, vượt trội so với than chì truyền thống. Tuy nhiên, rào cản kỹ thuật lớn nhất là sự giãn nở thể tích cực lớn (từ 300% đến 400%) trong quá trình xen cài (lithiation) và giải phóng ion lithium (delithiation). Biến dạng thể tích này tạo ra ứng suất cơ học rất cao, dẫn tới hiện tượng nứt gãy, nghiền vụn vật liệu, tách lớp khỏi bộ thu dòng và suy giảm dung lượng nhanh chóng sau một số chu kỳ sạc-xả. Bên cạnh đó, nhóm vật liệu hai chiều (2D) sở hữu diện tích bề mặt riêng lớn, độ bền cơ học cao cùng các đặc tính truyền dẫn điện tử tốt đang được xem xét như những ứng viên tiềm năng cho các thế hệ điện cực mới.

Khoảng trống nghiên cứu được luận án chỉ ra là sự thiếu hụt các dữ liệu định lượng chính xác về cơ tính (mô đun đàn hồi, hệ số Poisson, ứng suất kéo đứt, độ dai phá hủy) của nhiều cấu trúc vật liệu hai chiều đơn lớp, cũng như việc mô hình hóa chi tiết trường ứng suất - biến dạng dẻo kết hợp với động học khuếch tán ion lithium trên các cấu trúc hình học nano của silicon trong quá trình sạc điện.

Mục tiêu nghiên cứu của luận án được xác định gồm hai nội dung chính:

  1. Xác định đặc trưng cơ tính và cơ học phá hủy của các vật liệu hai chiều có cấu trúc lục giác phẳng và lục giác gồ ghề (low-buckled) có tiềm năng ứng dụng làm cực âm ắc quy ion lithium.
  2. Mô phỏng trường ứng suất, biến dạng và đánh giá độ bền cơ học của các kết cấu cực âm silicon (dạng hạt cầu, trụ dài, đĩa mỏng và kết cấu Si bọc trụ nano đồng) trong quá trình sạc điện để thiết lập mối quan hệ giữa thông số cơ học và thông số nạp điện.

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

  • Đối tượng nghiên cứu: Màng mỏng 4 vật liệu hai chiều cấu trúc lục giác phẳng (Graphene, Boron Nitride - BN, Silicon Carbide - SiC, Aluminum Nitride - AlN) và 31 vật liệu hai chiều cấu trúc lục giác low-buckled (Silicene, Germanene, Bismuthene, các hợp chất hai nguyên tố b-MX và b-M thuộc nhóm III, IV, V, VI); các kết cấu nano silicon làm cực âm bao gồm hạt cầu, trụ dài, đĩa mỏng và ống silicon bọc ngoài trụ đồng (Cu).
  • Phạm vi nghiên cứu: Nghiên cứu ứng xử kéo đơn trục, cơ học phá hủy vết nứt ở cấp độ nguyên tử bằng mô phỏng tính toán cơ học; mô phỏng quá trình khuếch tán ion lithium và trường ứng suất - biến dạng dẻo trong điều kiện sạc điện đơn kỳ dựa trên phương pháp phần tử hữu hạn liên tục.

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

Tổng quan tài liệu trong luận án điểm lại các mốc phát triển chính của hệ thống tích trữ năng lượng điện hóa và cơ học điện cực:

  • Về sự ra đời của ắc quy ion lithium: Nghiên cứu tiên phong của Whittingham (1970) [1] sử dụng titan disulfua ($TiS_2$) và kim loại lithium, tiếp nối bởi sự thương mại hóa pin ion lithium của tập đoàn Sony Energitech vào năm 1991 [4]. Các nghiên cứu về vật liệu catốt oxit kim loại chuyển tiếp ($LiCoO_2$, $LiMn_2O_4$, $LiFePO_4$) được tổng hợp bởi Goodenough, Fergus [5, 6].
  • Về vật liệu anốt gốc kim loại, oxit và hợp kim: Các nghiên cứu về cơ chế lưu trữ lithium trong cấu trúc nano cacbon, oxit kim loại ($MnO_2$, $CuO$, $MoO_2$) và hợp kim trên nền Sn, Co, Fe của tác giả Winter, Besenhard [9], Subba Rao, Chowdari [12], Zhang [13], Hassoun [28], Venkatachalam [32]. Các công trình phân tích sự thay đổi pha và suy giảm dung lượng của silicon khi xen cài lithium của Obrovac và Christensen [18, 42].
  • Về cấu trúc nano silicon: Các dạng hình thái sợi nano, hạt rỗng, màng mỏng xốp nhằm giảm thiểu nứt vỡ cơ học được nghiên cứu bởi Chan và cộng sự [55], Wu [57], Yao [58], Su [59]. Khảo sát thực nghiệm và lý thuyết về trường ứng suất trong hạt silicon do dòng khuếch tán của Zhao và Suo [60, 61], Zhang [63], Xiao [64].
  • Về vật liệu hai chiều: Kể từ khi Graphene được chế tạo thành công bởi Novoselov và Geim (2004) [73], hàng loạt vật liệu hai chiều đơn lớp khác được nghiên cứu và chế tạo như màng BN [75], Silicene trên đế kim loại [78, 90-94], màng AlN dạng graphene [96-98], và Phosphorene [80]. Cơ tính của graphene và màng đơn lớp được khảo sát bằng mô phỏng động lực học phân tử (MD) và DFT bởi Lee [141], Zhang [135, 142], Le-Minh [138].

Vị trí của luận án: Luận án giải quyết khoảng trống về việc tổng hợp toàn diện các thông số cơ học đàn hồi và phá hủy của 35 vật liệu 2D (trong đó có 31 hệ vật liệu lục giác low-buckled chưa được tính toán đầy đủ theo phương pháp phần tử hữu hạn nguyên tử) dưới hai hướng tải chirality khác nhau (Armchair và Zigzag). Đồng thời, luận án lấp đầy khoảng trống trong việc phân tích định lượng điều kiện dòng sạc tới hạn, ứng suất Von Mises và hiện tượng tách lớp cơ học tại bề mặt phân chia Si/Cu trên cấu trúc dây nano lõi-vỏ (core-shell) kích thước nano mét.

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

Luận án sử dụng kết hợp hai phương pháp tính toán cơ học ở cấp độ nguyên tử và cấp độ môi trường liên tục:

  1. Cơ sở lý thuyết tương tác nguyên tử và Phương pháp phần tử hữu hạn nguyên tử (AFEM):
  • Thế năng tương tác Tersoff: Mô tả liên kết cộng hóa trị đa vật thể dạng ba hạt, xét đến bậc liên kết phụ thuộc vào góc liên kết và môi trường nguyên tử xung quanh, áp dụng cho các hệ C, SiC, AlN, BN: $$E = \frac{1}{2} \sum_{i} \sum_{j \neq i} f_C(r_{ij}) \left[ f_R(r_{ij}) + b_{ij} f_A(r_{ij}) \right]$$
  • Thế năng tương tác Stillinger-Weber (S-W): Mô tả năng lượng liên kết hai hạt (kéo nén liên kết thẳng) và ba hạt (biến dạng góc liên kết), chuẩn hóa cho các mạng tinh thể lục giác phẳng và low-buckled: $$E = \sum_{i<j} v_2(r_{ij}) + \sum_{i<j<k} v_3(r_{ij}, r_{ik}, \theta_{jik})$$
  • Phương pháp AFEM: Thiết lập dựa trên nguyên lý thế năng toàn phần cực tiểu. Đạo hàm bậc nhất của thế năng nguyên tử theo véc-tơ chuyển vị của các nút (nguyên tử) cho véc-tơ nội lực; đạo hàm bậc hai cho ma trận độ cứng phi tuyến của cấu trúc nguyên tử. Hệ phương trình đại số phi tuyến được giải bằng thuật toán lặp Newton-Raphson: $$\mathbf{K}(\mathbf{u}) \Delta \mathbf{u} = \mathbf{P}{ext} - \mathbf{F}{int}(\mathbf{u})$$
  1. Cơ sở cơ học phá hủy vết nứt ở cấp nano:
  • Xác định hệ số cường độ ứng suất tới hạn dạng mở miệng ($K_{Ic}$) thông qua hai cách tiếp cận:
    • Kéo màng có vết nứt cạnh (Single-Edge Notched Tension - SENT): Đo ứng suất kéo đứt $\sigma_f$ để tính $K_{Ic} = \sigma_f \sqrt{\pi a_0} \cdot F(a_0/W)$.
    • Phương pháp trường chuyển vị đầu vết nứt: Áp đặt điều kiện biên chuyển vị theo nghiệm giải tích trường ứng suất đàn hồi tuyến tính vùng đầu vết nứt (Mode I), tìm giá trị lực tới hạn làm đứt gãy liên kết nguyên tử đầu tiên tại đỉnh vết nứt.
  1. Lý thuyết tương tự Nhiệt - Khuếch tán và Cơ học vật rắn biến dạng trong môi trường liên tục:
  • Phương trình khuếch tán ion lithium theo định luật Fick II được chuyển đổi tương đương toán học sang phương trình truyền dẫn nhiệt phi dừng: $$\frac{\partial C}{\partial t} = \nabla \cdot (D \nabla C) \iff \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T)$$
  • Biến dạng do khuếch tán ion lithium gây ra ($\varepsilon_{ij}^{chem} = \frac{1}{3} \Omega C \delta_{ij}$) tương đương với biến dạng nhiệt biến thiên theo nồng độ $C$.
  • Quan hệ cơ học đàn - dẻo: Cực âm silicon được mô hình hóa theo lý thuyết đàn dẻo chảy dẻo Von Mises đẳng hướng có hóa bền dẻo (kinematic/isotropic hardening), giải bằng phần mềm phần tử hữu hạn MARC với phần tử khối 8 nút 3D.
  • Đánh giá hiện tượng tách lớp tiếp xúc Si-Cu qua chỉ số tách lớp (Delamination Index - $I_D$) dựa trên tốc độ giải phóng năng lượng biến dạng bề mặt $G$ so với năng lượng tới hạn $G_c$: $$I_D = \frac{G}{G_c}$$

Nguồn dữ liệu và thông số tính toán: Các thông số hình học mạng, hằng số mạng $d$, độ cao gồ ghề $\Delta$, góc liên kết $\theta_0$, các hệ số hàm thế Tersoff và Stillinger-Weber được trích xuất từ các công trình tính toán ab-initio DFT đã công bố của Jiang [123], Tersoff [118], Powell [120].

Bảng 1: Bảng tổng hợp các nhóm phương pháp và công cụ tính toán
Phương pháp Đối tượng và Mục tiêu áp dụng
AFEM (Thế Tersoff) 4 màng lục giác phẳng (Graphene, BN, SiC, AlN) - Tính $E, \nu, \sigma_f, K_{Ic}$
AFEM (Thế Stillinger-Weber) 31 màng lục giác low-buckled (Silicene, Germanene, b-MX,...) - Tính cơ tính và $K_{Ic}$
Phần tử hữu hạn FEM (MARC) Cấu trúc nano Si (cầu, trụ, đĩa, bọc Cu) - Trường nồng độ, ứng suất, biến dạng dẻo, $I_D$

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

Chương 1. Tổng quan về vật liệu làm cực âm ắc quy ion lithium

Chương 1 trình bày lịch sử hình thành các thế hệ ắc quy sạc, cấu trúc hoạt động và phân tích chuyên sâu các nhóm vật liệu chế tạo anốt. Cấu tạo của một tế bào ắc quy cơ bản bao gồm 5 thành phần: đầu thu dòng điện (collector), cực âm (anode), cực dương (cathode), chất điện phân (electrolyte dạng lỏng hoặc chất điện ly rắn) và màng ngăn xốp (separator) như polyethylen (PE) hoặc polypropylen (PP). Quá trình sạc và xả điện được mô tả thông qua các phản ứng điện hóa chuyển dịch thuận nghịch của ion $Li^+$ và electron tự do qua mạch ngoài:

  • Tại cực dương: $Li_xMO_y \rightleftharpoons MO_y + x Li^+ + x e^-$
  • Tại cực âm: $X + x Li^+ + x e^- \rightleftharpoons Li_xX$

Chương này phân loại chi tiết các vật liệu cực âm gồm:

  1. Kim loại, oxit kim loại ($SnO_2, MnO_2, Fe_2O_3, CuO$) và hợp kim ($Sn-Co, Sn-Fe, Cu-Sn$): Khắc phục một phần sự co giãn thể tích nhưng khả năng dẫn điện và hiệu suất chu kỳ chưa cao.
  2. Silicon và các nguyên tố nhóm IV (Ge, Sn): Silicon có dung lượng lý thuyết $4200\text{ mAh/g}$ ($Li_{22}Si_5$), cao gấp hơn 10 lần so với than chì ($372\text{ mAh/g}$ ứng với $LiC_6$). Tuy nhiên sự nở thể tích $300-400%$ đòi hỏi phải giảm kích thước cấu trúc xuống thang nano (dây nano, màng mỏng, hạt hình cầu).
  3. Vật liệu 2D (Graphene, BN, Silicene, MXenes, Phosphorene): Cung cấp đường dẫn khuếch tán ngắn cho ion lithium và độ dẻo dai cơ học cao.

Chương 2. Cấu trúc vật liệu hai chiều và phương pháp mô phỏng

Chương 2 thiết lập cơ sở toán học và cơ học cho toàn bộ các bài toán mô phỏng trong luận án. Cấu trúc tinh thể hai chiều được khảo sát dưới hai dạng hình học: mạng lục giác phẳng hoàn hảo (như Graphene, h-BN) với liên kết lai hóa $sp^2$ và mạng lục giác gồ ghề (low-buckled) với các nguyên tử nằm so le trên hai mặt phẳng song song cách nhau khoảng cách $\Delta$ (như Silicene, Germanene, Stanene).

Nội dung chương trình bày chi tiết:

  • Xây dựng ma trận phần tử hữu hạn nguyên tử AFEM cho phần tử 2 nút (liên kết thẳng hai hạt), 3 nút và 4 nút (liên kết góc và xoắn đa hạt).
  • Thiết lập các phương trình giải tích xác định ứng suất cơ học màng hai chiều quy đổi theo chiều dày hiệu dụng van der Waals ($t$).
  • Thiết lập bài toán cơ học phá hủy nano: Mô hình hóa vết nứt có chiều dài ban đầu $a_0$ trên tấm có chiều rộng $W$, phân tích điều kiện tới hạn theo phương Zigzag (ZZ) và Armchair (AC).
  • Xây dựng hệ phương trình vi phân liên kết khuếch tán - ứng suất trong môi trường liên tục biến dạng lớn, mô hình hóa sự biến thiên mô đun đàn hồi $E(C)$, hệ số giãn nở nồng độ $\Omega$, giới hạn chảy dẻo $\sigma_y$ của silicon khi ion lithium thâm nhập.

Chương 3. Mô phỏng ứng xử cơ học vật liệu hai chiều

Chương 3 thực hiện mô phỏng số chi tiết trên 35 màng đơn lớp vật liệu hai chiều bằng AFEM:

  • Kéo màng nguyên (Pristine sheets): Xác định đường cong quan hệ ứng suất - biến dạng ($\sigma - \varepsilon$) phi tuyến đến khi đứt gãy. Kết quả cho thấy rõ tính bất đẳng hướng cơ học giữa hai phương kéo Zigzag (ZZ) và Armchair (AC).
    • Màng Graphene có mô đun đàn hồi $E \approx 1.0\text{ TPa}$, ứng suất kéo đứt đạt trên $100\text{ GPa}$.
    • Màng AlN đạt ứng suất kéo đứt cực đại tại biến dạng $\varepsilon = 23.72%$ (phương ZZ) và $\varepsilon = 17.93%$ (phương AC).
    • Màng Silicene chịu tải qua hàm thế Stillinger-Weber cho thấy biến dạng kéo đứt phương ZZ lớn hơn phương AC.
  • Kéo màng có vết nứt ở cạnh và trường chuyển vị đầu vết nứt: Đánh giá độ nhạy phá hủy của 31 vật liệu cấu trúc low-buckled. Khi màng chứa vết nứt cạnh kích thước $a_0/W \approx 0.1$, ứng suất kéo đứt của màng giảm mạnh. Ví dụ: màng Germanene nứt theo phương AC kéo phương ZZ có ứng suất đứt $\sigma_{ft} = 0.86\text{ N/m}$; màng b-Bismuthene đạt $\sigma_{ft} = 1.073\text{ N/m}$ (kéo ZZ) và $\sigma_{ft} = 2.156\text{ N/m}$ (kéo AC).
  • Tính toán hệ số cường độ ứng suất tới hạn $K_{Ic}$: Thu được giá trị $K_{Ic}$ đặc trưng cho từng loại liên kết mạng. Với màng AlN, giá trị $K_{Ic}$ thu được đạt $1.209\text{ MPa}\sqrt{\text{m}}$ (vết nứt phương ZZ) và $1.532\text{ MPa}\sqrt{\text{m}}$ (vết nứt phương AC).
Bảng 2: So sánh đặc trưng cơ tính màng nguyên và độ dai phá hủy một số vật liệu 2D tiêu biểu
Vật liệu 2D Cấu trúc mạng Phương kéo / Nứt Biến dạng đứt $\varepsilon_f$ (%) Hệ số $K_{Ic}$ ($\text{MPa}\sqrt{\text{m}}$)
Graphene Lục giác phẳng ZZ / AC $18.0 - 25.0$ $4.0 - 4.5$
h-BN Lục giác phẳng ZZ / AC $16.5 - 22.0$ $2.8 - 3.2$
AlN Lục giác phẳng ZZ / AC $17.9 - 23.7$ $1.209 - 1.532$
Silicene Low-buckled ZZ / AC $14.0 - 19.5$ $1.150 - 1.380$
Germanene Low-buckled ZZ / AC $12.5 - 17.0$ $0.850 - 1.020$
Bismuthene Low-buckled ZZ / AC $10.0 - 15.2$ $0.620 - 0.810$

Chương 4. Mô phỏng ứng xử cơ học kết cấu hạt cầu, trụ và đĩa mỏng silicon dùng làm cực âm ắc quy ion lithium

Chương 4 khảo sát chi tiết ba cấu trúc hình học nano của anốt silicon đơn lẻ: hạt cầu (bán kính $A$), trụ dài (bán kính $A$) và đĩa mỏng tròn (bán kính $A$, chiều dày $H = 0.05A \div 0.1A$). Quá trình sạc điện được khảo sát theo các giá trị dòng sạc không thứ nguyên $\bar{J} = J A / (C_{max} D)$ biến thiên từ $0.1$ đến $4/3$.

Các kết quả định lượng thu được:

  • Phân bố nồng độ ion lithium $C/C_{max}$: Tạo nên một gradient nồng độ rất dốc từ bề mặt ngoài vào tâm.
  • Trạng thái ứng suất: Bề mặt ngoài của hạt và trụ chịu ứng suất nén trong giai đoạn đầu do lớp vỏ ngoài giãn nở trước nhưng bị lõi trong cản trở. Ngược lại, vùng tâm hạt chịu ứng suất kéo 3 trục rất cao, là nguyên nhân phát sinh các vi nứt từ bên trong.
  • Ứng suất tương đương Von Mises $\sigma_{eq}$ và biến dạng dẻo tương đương $\varepsilon^p$: Khi dòng sạc vượt qua giá trị dòng tới hạn $\bar{J}_c$, vùng ngoài của hạt cầu và trụ chuyển sang trạng thái chảy dẻo. Kết cấu đĩa mỏng cho thấy mức độ tập trung ứng suất tại mép góc thấp hơn đáng kể nếu được vê tròn/lượn góc mép.
  • Ảnh hưởng của hình học đến hiệu suất nạp: Kết cấu đĩa mỏng với tỷ lệ diện tích mặt thoáng trên thể tích lớn cho phép đạt trạng thái nạp đầy với thời gian ngắn nhất và ứng suất tích lũy cực đại thấp hơn kết cấu hạt cầu cùng thể tích danh định.

Chương 5. Mô phỏng ứng xử cơ học kết cấu silicon bọc vào trụ đồng kích cỡ nano mét dùng làm cực âm ắc quy ion lithium

Chương 5 nghiên cứu cấu trúc anốt phức hợp lõi-vỏ (core-shell): Dây nano đồng (Cu) dẫn điện làm lõi, bên ngoài bọc một lớp màng silicon (Si) làm vật liệu tích trữ ion. Mô hình số khảo sát các điều kiện biên khuếch tán: khuếch tán trên toàn bộ bề mặt ngoài, khuếch tán ở mặt bên (chặn khuếch tán đỉnh trụ) và cấu trúc có xét góc lượn tại đỉnh Si.

Các phát hiện chính của chương:

  • Động thái biến dạng: Lớp silicon nở phồng lớn cả theo phương hướng kính ($\Delta Y$) và phương dọc trục ($\Delta X$). Do mô đun đàn hồi của Cu ($E_{Cu} \approx 110\text{ GPa}$) cao hơn và Cu không bị biến dạng do khuếch tán, lõi Cu kìm hãm biến dạng của vỏ Si, làm phát sinh ứng suất trượt rất lớn tại bề mặt tiếp xúc Si/Cu.
  • Tác dụng của góc lượn đỉnh: Kết cấu đỉnh Si sắc nhọn gây tập trung ứng suất Von Mises cục bộ cực đại tại góc đỉnh ($\tau = 0.01 \div 0.108$). Việc thiết kế góc lượn làm giảm đáng kể mức ứng suất đỉnh và làm đều phân bố biến dạng dẻo $\varepsilon^p$.
  • Hiện tượng tách lớp Si-Cu: Đánh giá thông qua chỉ số $I_D$. Luận án khảo sát các giá trị năng lượng giải phóng bề mặt tới hạn $G_c = 15\text{ J/m}^2$ và $G_c = 20\text{ J/m}^2$, hệ số khuếch tán $D = 100\text{ nm}^2/\text{s}$, mật độ dòng sạc $J = 1\text{ nm}^{-2}\text{s}^{-1}$. Khi $I_D \ge 1$, hiện tượng bong tróc màng Si khỏi trụ Cu xảy ra, bắt đầu từ vùng đỉnh và mép đáy tiếp xúc, giải thích cơ chế suy giảm tiếp xúc điện sau nhiều chu kỳ vận hành.

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

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

  1. Xác định một cách hệ thống các thông số cơ tính đàn hồi (mô đun đàn hồi $E$, hệ số Poisson $\nu$), ứng suất kéo đứt $\sigma_f$, biến dạng kéo đứt $\varepsilon_f$ và hệ số cường độ ứng suất tới hạn $K_{Ic}$ của 35 loại vật liệu hai chiều đơn lớp, bao gồm 4 vật liệu lục giác phẳng và 31 vật liệu lục giác low-buckled theo hai phương chirality (Zigzag và Armchair).
  2. Phát triển và ứng dụng thành công phương pháp phần tử hữu hạn nguyên tử (AFEM) kết hợp với các hàm thế năng Tersoff và Stillinger-Weber hiệu chỉnh để khảo sát độ bền phá hủy vùng đầu vết nứt nano trên các tinh thể 2D một cách ổn định, không phụ thuộc vào bậc hội tụ thời gian như mô phỏng động lực học phân tử cổ điển.
  3. Thiết lập được quy luật phân bố trường nồng độ ion lithium, ứng suất tương đương Von Mises, biến dạng dẻo tích lũy và điều kiện dòng sạc không thứ nguyên tới hạn ($\bar{J}_c$) cho các kết cấu anốt silicon dạng hạt cầu, trụ dài và đĩa mỏng.
  4. Làm sáng tỏ cơ chế tập trung ứng suất, biến dạng trượt và thiết lập đồ thị tiến triển chỉ số tách lớp ($I_D$) tại mặt phân cách tiếp xúc giữa silicon và trụ đồng theo thời gian nạp điện và năng lượng tới hạn $G_c$.

Kiến nghị và giải pháp đề xuất:

  • Khi thiết kế cực âm hạt nano Si hoặc đĩa Si, cần kiểm soát dòng nạp dưới giá trị dòng không thứ nguyên tới hạn $\bar{J}_c$ để tránh xuất hiện biến dạng dẻo không hồi phục ở lớp vỏ ngoài.
  • Với cấu trúc dây nano Si bọc lõi Cu, bắt buộc phải gia công hoặc tổng hợp biên dạng có góc lượn ở đỉnh trụ nhằm triệt tiêu điểm kỳ dị tập trung ứng suất gây bong tróc màng điện cực.
  • Ưu tiên cấu trúc đĩa mỏng có lượn mép hoặc kết cấu nano có tỷ lệ bề mặt thoáng cao để phân tán ứng suất giãn nở thể tích trong quá trình nạp nhanh.

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

Hạn chế của nghiên cứu:

  • Chưa khảo sát các hệ vật liệu hai chiều có cấu trúc mạng phức tạp hơn cấu trúc lục giác phẳng và low-buckled (như mạng bát giác, tứ giác hay các dị cấu trúc đa lớp heterostructures).
  • Các hệ số khuếch tán ion lithium $D$ và mật độ dòng sạc $J$ trong mô hình liên tục của silicon đang được giả định là hằng số hoặc hàm đơn giản, chưa tính đến sự phụ thuộc phi tuyến mạnh vào trạng thái điện áp tế bào và nồng độ ion biến thiên theo dung tích thời gian thực.
  • Mô phỏng ứng xử cơ học của anốt silicon mới chỉ dừng lại ở nửa chu kỳ sạc điện đầu tiên (lithiation đơn kỳ), chưa mô phỏng lặp chu kỳ sạc - xả kết hợp nhiều lần (cyclic lithiation/delithiation) để đánh giá hiện tượng mỏi cơ học và tích lũy hư hại shakedown.

Hướng nghiên cứu tiếp theo được tác giả đề xuất:

  1. Xác định cơ tính và ứng xử phá hủy của các vật liệu hai chiều mới có cấu trúc phi lục giác có tiềm năng làm điện cực.
  2. Thiết lập mô hình mô phỏng quá trình nạp điện với hệ số khuếch tán ion lithium và dòng sạc là các hàm biến thiên phụ thuộc vào thời gian và dung lượng ắc quy thực tế.
  3. Mở rộng mô phỏng cho nhiều chu kỳ sạc và xả điện liên tiếp nhằm theo dõi quá trình suy giảm liên kết tiếp xúc cơ - điện.
  4. Đề xuất và tối ưu hóa các hình học cấu trúc mới của anốt silicon nhằm sớm đưa vật liệu anốt gốc silicon vào ứng dụng thương mại.

Giá trị tham khảo

Luận án cung cấp nguồn tư liệu học thuật và thông số kỹ thuật có giá trị cho:

  • Các nghiên cứu sinh, học viên cao học và nhà nghiên cứu thuộc các ngành Cơ học vật rắn, Cơ kỹ thuật, Khoa học và Kỹ thuật Vật liệu, Kỹ thuật Năng lượng: Cung cấp bộ dữ liệu hoàn chỉnh về đặc tính đàn hồi, độ dai phá hủy $K_{Ic}$ của 35 hệ màng 2D và khung thuật toán AFEM phi tuyến.
  • Các nhóm kỹ sư thiết kế pin và hệ thống tích trữ năng lượng: Cung cấp các tiêu chuẩn về hình học (bán kính tới hạn của hạt, đĩa mỏng, góc lượn dây nano Si bọc Cu) và biểu đồ dòng nạp giới hạn $\bar{J}_c$ nhằm ngăn chặn sự nứt vỡ cơ học và hiện tượng tách lớp điện cực trong quá trình sạc nhanh.
  • Nội dung tham khảo trọng tâm: Phương pháp chuyển đổi tương đương bài toán khuếch tán Fick sang truyền nhiệt trong MARC (Chương 2, 4) và phương pháp xác định $K_{Ic}$ cấp nguyên tử (Chương 2, 3).

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

1. Luận án khảo sát bao nhiêu loại vật liệu hai chiều và phân loại chúng như thế nào?
Luận án khảo sát tổng cộng 35 vật liệu hai chiều đơn lớp, chia làm hai nhóm cấu trúc:

  • 4 vật liệu có cấu trúc lục giác phẳng: Graphene, Boron Nitride (BN), Silicon Carbide (SiC) và Aluminum Nitride (AlN).
  • 31 vật liệu có cấu trúc lục giác gồ ghề (low-buckled): Silicene, Germanene, Bismuthene, cùng các hợp chất hai nguyên tố b-MX và b-M thuộc nhóm III, IV, V, VI (như b-GeM, b-SiM, b-SnM, b-InM, b-CM, b-XM).

2. Tại sao silicon có dung lượng lưu trữ lý thuyết rất cao nhưng chưa được ứng dụng rộng rãi làm cực âm thương mại?
Silicon có dung lượng lý thuyết lên tới $\sim 4200\text{ mAh/g}$ (cao gấp hơn 10 lần than chì), tuy nhiên khi sạc đầy ion lithium, thể tích silicon bị giãn nở từ $300%$ đến $400%$. Sự biến dạng quá lớn này sinh ra ứng suất cơ học khổng lồ, dẫn tới hiện tượng nứt vỡ, nghiền vụn hạt silicon, mất tiếp xúc điện với đầu thu dòng và làm dung lượng ắc quy suy giảm nghiêm trọng chỉ sau một thời gian ngắn vận hành.

3. Phương pháp phần tử hữu hạn nguyên tử (AFEM) được xây dựng dựa trên nguyên lý nào và dùng hàm thế năng gì?
AFEM được xây dựng dựa trên nguyên lý thế năng toàn phần cực tiểu, trong đó năng lượng liên kết nguyên tử được lấy đạo hàm bậc một và bậc hai theo chuyển vị các nút để lập véc-tơ nội lực và ma trận độ cứng phần tử phi tuyến. Luận án sử dụng hai hàm thế năng thực nghiệm: hàm thế Tersoff (cho các màng lục giác phẳng C, BN, SiC, AlN) và hàm thế Stillinger-Weber (cho các màng lục giác low-buckled).

4. Dòng sạc không thứ nguyên $\bar{J}$ được định nghĩa như thế nào và ảnh hưởng gì đến trạng thái cơ học của anốt silicon?
Dòng sạc không thứ nguyên được xác định bởi công thức $\bar{J} = \frac{J \cdot A}{C_{max} \cdot D}$ (trong đó $J$ là mật độ dòng khuếch tán bề mặt, $A$ là kích thước đặc trưng/bán kính kết cấu, $C_{max}$ là nồng độ bão hòa lớn nhất và $D$ là hệ số khuếch tán). Khi $\bar{J}$ vượt quá ngưỡng tới hạn $\bar{J}_c$, gradient nồng độ tăng vọt làm bề mặt ngoài bị biến dạng dẻo và ứng suất kéo tại tâm tăng cao, làm tăng nguy cơ phá hủy kết cấu.

5. Kết cấu Si bọc vào trụ nano Cu có ưu điểm gì và hiện tượng hỏng hóc chính xảy ra ở đâu?
Kết cấu lõi-vỏ (Cu lõi, Si vỏ) tận dụng khả năng dẫn điện cao của trụ đồng và dung lượng tích trữ lớn của màng silicon ngoài. Tuy nhiên, do đồng không giãn nở nồng độ trong khi silicon nở mạnh, ứng suất trượt lớn sinh ra tại bề mặt tiếp xúc Si-Cu. Luận án chỉ ra rằng sự tách lớp cơ học (chỉ số $I_D \ge 1$) và tập trung ứng suất Von Mises bắt đầu xuất hiện tại các góc sắc nhọn trên đỉnh trụ và mép đáy tiếp xúc.

Kết luận

Luận án của nghiên cứu sinh Nguyễn Hữu Tú đã thực hiện thành công việc mô hình hóa và phân tích định lượng ứng xử cơ học, độ bền phá hủy của 35 hệ vật liệu hai chiều cùng các cấu trúc nano silicon làm cực âm cho ắc quy ion lithium. Bằng việc kết hợp phương pháp phần tử hữu hạn nguyên tử (AFEM) và phương pháp phần tử hữu hạn liên tục (FEM) mô phỏng bài toán liên kết khuếch tán - biến dạng dẻo lớn, công trình đã làm sáng tỏ tính bất đẳng hướng cơ học của màng mỏng 2D và cơ chế phát sinh ứng suất, biến dạng dẻo và tách lớp trên các kết cấu anốt nano. Các kết quả tính toán và đồ thị tới hạn thu được cung cấp cơ sở khoa học định lượng phục vụ việc tối ưu hóa hình học, lựa chọn vật liệu và kiểm soát chế độ sạc điện an toàn cho các thế hệ ắc quy ion lithium dung lượng cao tiếp theo.