Tổng quan về luận án

Sự phát triển vũ bão của các thiết bị điện tử cầm tay, phương tiện giao thông chạy điện (Electric Vehicles - EVs) và các hệ thống lưu trữ năng lượng lưới điện thông minh đặt ra yêu cầu cấp thiết về việc phát triển thế hệ ắc quy ion lithium (Lithium-ion Battery - LIB) có mật độ năng lượng vượt trội, tốc độ sạc nhanh và tuổi thọ bền bỉ. Trong cấu trúc tổng thể của ắc quy ion lithium, cực âm (Anode) đóng vai trò quyết định đến dung lượng chứa ion, tốc độ dòng nạp/xả và độ suy giảm dung lượng theo chu kỳ hoạt động. Vật liệu than chì (Graphite) thương mại truyền thống với cơ chế xen cài ion $\text{LiC}_6$ đã tiệm cận giới hạn dung lượng lý thuyết ($372\text{ mAh/g}$), không còn đáp ứng được nhu cầu năng lượng cao của các công nghệ tương lai.

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ thực tế: "Silicon là vật liệu có dung lượng chứa ion lithium cao nhất ($\sim 4200\text{ mAh/g}$), nhưng do thể tích của cực âm làm từ silicon tăng rất lớn khi sạc đầy (300 - 400%) nên ắc quy loại này chưa được ứng dụng vào thực tế". Sự phồng nở thể tích bất đẳng hướng khổng lồ trong quá trình tích nạp ion lithium làm phát sinh trường ứng suất cơ học khuếch tán (diffusion-induced stress) với cường độ lớn, dẫn đến hiện tượng nứt vỡ vi mô (pulverization), nghiền vụn cấu trúc, suy giảm tiếp xúc điện và gây ra sự tách lớp (delamination) thảm khốc giữa màng silicon với đế dẫn điện (current collector Cu). Đồng thời, nhóm vật liệu hai chiều (2D materials) tiên tiến—vốn sở hữu diện tích bề mặt riêng lớn và tính chất điện hóa đặc biệt—chưa được thiết lập một hệ thống dữ liệu toàn diện về đặc trưng cơ tính và cơ học phá hủy nguyên tử trước khi đưa vào chế tạo cực âm.

Luận án của Nghiên cứu sinh Nguyễn Hữu Tú, thực hiện tại Trường Đại học Bách khoa Hà Nội dưới sự hướng dẫn khoa học của PGS.TS Lê Minh Quý và TS. Trần Đình Long, tập trung giải quyết các câu hỏi và giả thuyết nghiên cứu cụ thể:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ tính đàn hồi, ứng suất kéo đứt và độ dai phá hủy nguyên tử (thông qua hệ số cường độ ứng suất tới hạn $K_{\text{Ic}}$) của các vật liệu hai chiều cấu trúc lục giác phẳng và lục giác lồi lõm (low-buckled) biến thiên như thế nào theo phương tinh thể?
    • Giả thuyết 1 (H1): Cấu trúc lồi lõm low-buckled làm suy giảm độ cứng chống kéo so với màng phẳng nhưng tạo ra tính dị hướng phá hủy rõ rệt giữa hai phương Armchair (AC) và Zigzag (ZZ).
  • Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ định lượng giữa mật độ dòng sạc khuếch tán, hệ số khuếch tán và sự tiến triển của trường ứng suất tương đương von Mises cùng biến dạng dẻo trong các kết cấu cực âm silicon hình cầu, trụ dài và đĩa mỏng diễn ra theo quy luật nào?
    • Giả thuyết 2 (H2): Kết cấu đĩa mỏng silicon tối ưu hóa gradient nồng độ ion lithium tốt hơn dạng cầu và trụ, cho phép nâng cao giá trị dòng sạc không thứ nguyên tới hạn $\bar{J}_c$ mà không gây phá hủy dẻo.
  • Câu hỏi nghiên cứu 3 (RQ3): Hiện tượng tập trung ứng suất và cơ chế tách lớp tại mặt phân giới giữa lớp silicon và trụ nano đồng (Cu) có thể được kiểm soát và triệt tiêu bằng giải pháp tối ưu hóa hình học biên (góc lượn đỉnh trụ) hay không?
    • Giả thuyết 3 (H3): Việc khử bỏ góc nhọn trên đỉnh trụ Si bọc Cu bằng góc lượn chuyển tiếp sẽ làm giảm chỉ số tách lớp $I_D$ xuống dưới ngưỡng tới hạn, ngăn ngừa phá hủy giao diện khi đạt dung lượng nạp lớn.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Cơ học phá hủy nguyên tử (Atomistic Fracture Mechanics) sử dụng các hàm thế kinh nghiệm đa hạt (Tersoff, Stillinger-Weber), Lý thuyết khuếch tán phi tuyến Fickian (Nonlinear Fickian Diffusion), và Cơ học môi trường liên tục biến dạng lớn tích hợp chảy dẻo von Mises (Finite Strain Elastoplasticity). Phạm vi nghiên cứu bao gồm việc khảo sát cơ tính của 35 loại vật liệu hai chiều đơn lớp (với quy mô tính toán lên tới 13.376 nguyên tử trên một đơn vị mô phỏng) và mô phỏng số chi tiết ứng xử cơ - hóa của các cấu trúc silicon kích thước nanomet ($A = 200\text{ nm}$) dưới tác động của dòng sạc chuẩn hóa $\bar{J}$ từ 0.1 đến 4/3 và hệ số khuếch tán $D = 100\text{ nm}^2/\text{s}$.

Literature Review và Positioning

Lịch sử phát triển của ắc quy lithium khởi nguồn từ các công trình tiên phong của M. Stanley Whittingham (1976) với hệ điện cực $\text{TiS}_2/\text{Li}$, tiếp nối bởi sự thương mại hóa pin ion lithium dùng điện cực than chì của Sony Energitech vào năm 1991. Trước nhu cầu nâng cao mật độ năng lượng, hướng nghiên cứu chuyển dịch mạnh mẽ sang các vật liệu nhóm IV, đặc biệt là silicon. Các nghiên cứu đột phá của Chan et al. (2008) và Cui et al. đã chứng minh rằng cấu trúc nano (dây nano Si, hạt nano Si rỗng) có khả năng giải tỏa biến dạng thể tích tốt hơn vật liệu khối. Tuy nhiên, việc mô hình hóa chính xác sự hình thành trường ứng suất và phá hủy cơ học trong quá trình khuếch tán ion lithium đòi hỏi nền tảng cơ học chất rắn chuyên sâu.

Trong dòng nghiên cứu về mô hình hóa tương tác Hóa - Cơ (Coupled Chemomechanics), các tác giả quốc tế tiêu biểu như Bower et al. (2011), Haftbaradaran et al. (2010), và Zhao et al. (2011) đã xây dựng các mô hình môi trường liên tục biến dạng hữu hạn để mô tả dòng chảy dẻo do khuếch tán lithium. Dẫu vậy, các tranh luận học thuật sâu sắc vẫn tồn tại qua hai luồng quan điểm đối lập:

  • Luồng quan điểm 1 (Continuum Mechanics Stance): Cho rằng ứng xử phá hủy của vật liệu nano cực âm có thể được dự báo chính xác bằng các tiêu chuẩn cơ học môi trường liên tục cổ điển có hiệu chỉnh điều kiện biên bề mặt (như mô hình của Christensen et al., 2006; Bower et al., 2011).
  • Luồng quan điểm 2 (Atomistic Mechanics Stance): Lập luận rằng ở quy mô vài nanomet hoặc với các màng mỏng 2D, hiệu ứng bất liên tục mạng tinh thể, liên kết hóa học cục bộ và sự định hướng chirality (Armchair vs. Zigzag) giữ vai trò chi phối tuyệt đối; do đó các mô hình môi trường liên tục thường ước tính sai lệch độ bền phá hủy thực tế (Rouhi et al., 2017; Zhang et al., 2014).
========================================================================================
KHUNG ĐỊNH VỊ HỌC THUẬT (ACADEMIC POSITIONING)
========================================================================================
Các nghiên cứu trước:
- Bower et al. (2011) ---------> Mô hình Hóa-Cơ liên tục biến dạng hữu hạn (khung lý thuyết Si)
- Rouhi et al. (2017) ---------> Mô phỏng MD đứt gãy Silicene đơn lẻ (hạn chế quy mô vật liệu)
- Chan & Cui et al. (2008) ----> Thực nghiệm Dây nano Si (thiếu định lượng cơ học đa dạng hình học)
----------------------------------------------------------------------------------------
Đóng góp định vị của Luận án (Nguyễn Hữu Tú, 2019):
===> Tích hợp AFEM đa thế năng: Chuẩn hóa cơ tính 35 vật liệu 2D (4 phẳng + 31 low-buckled)
===> Tương tự Nhiệt - Khuếch tán trên MARC: Giải mã ứng suất nứt vỡ đĩa/trụ/hạt cầu Si
===> Tối ưu hóa cấu trúc lõi-vỏ Si/Cu: Giải pháp góc lượn triệt tiêu chỉ số tách lớp ID
========================================================================================

Luận án xác lập vị trí khoa học tại giao điểm của Cơ học tính toán nguyên tử (AFEM) và Kỹ thuật mô phỏng phần tử hữu hạn phi tuyến (Nonlinear FEM), định lượng hóa chính xác cơ tính của 35 màng 2D đơn lớp và đưa ra lời giải toàn diện cho bài toán tập trung ứng suất tiếp xúc trong kết cấu composite Si/Cu kích thước nano.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mang lại những đóng góp nền tảng cho lý thuyết cơ học vật liệu tiên tiến và lý thuyết điện hóa cơ học:

  • Mở rộng cơ học phá hủy nguyên tử cho vật liệu 2D: Mở rộng Lý thuyết cơ học phá hủy đàn hồi tuyến tính (Linear Elastic Fracture Mechanics - LEFM của Griffith và Irwin) xuống cấp độ mạng nguyên tử. Thông qua việc phân tích màng 2D có vết nứt cạnh chịu kéo và áp đặt trường chuyển vị giải tích đầu vết nứt, luận án xác định chính xác hệ số cường độ ứng suất tới hạn $K_{\text{Ic}}$ cho 31 vật liệu cấu trúc lục giác low-buckled ($b\text{-M}$, $b\text{-MX}$) và 4 vật liệu lục giác phẳng (Graphene, BN, SiC, AlN).
  • Phát triển mô hình tương tác Hóa - Cơ dẻo phi tuyến: Đóng góp vào lý thuyết khuếch tán - biến dạng dẻo (Diffusion-Plasticity Theory) bằng cách lượng hóa ngưỡng dòng sạc tới hạn $\bar{J}_c$, chứng minh rằng sự tích lũy biến dạng dẻo tương đương $\varepsilon^p$ phụ thuộc phi tuyến vào hình học cấu trúc và vận tốc truyền dẫn ion tại bề mặt tiếp xúc.
  • Chuyển dịch mô thức thiết kế điện cực (Paradigm Shift): Chuyển dịch tư duy thiết kế điện cực từ việc phân tán hạt thụ động sang chủ động kiểm soát trường ứng suất vi mô thông qua biến điệu hình học (micro-architectural stress relief), cung cấp bằng chứng định lượng về sự triệt tiêu điểm kỳ dị ứng suất tại mặt tiếp xúc dị chất Si/Cu.

Khung phân tích độc đáo

Khung phân tích của luận án là sự hợp nhất của 3 trụ cột lý thuyết:

  1. Lý thuyết thế năng tương tác nguyên tử (Empirical Interatomic Potential Theory): Sử dụng hàm thế Tersoff (1988) cho các liên kết 3 nút, 4 nút và hàm thế Stillinger-Weber (1985) tách biệt năng lượng liên kết thẳng (two-body) và liên kết góc (three-body).
  2. Nguyên lý thế năng toàn phần cực tiểu (Principle of Minimum Potential Energy): Thiết lập hệ phương trình phi tuyến của phương pháp phần tử hữu hạn nguyên tử (Atomic-scale Finite Element Method - AFEM): $$\mathbf{\Psi}(\mathbf{u}) = \frac{\partial \Pi}{\partial \mathbf{u}} = \mathbf{0}$$
  3. Lý thuyết tương tự Nhiệt - Khuếch tán (Thermal-Diffusion Analogy): Ánh xạ toán học phương trình vi phân khuếch tán khối Fickian bậc hai sang phương trình dẫn nhiệt biến thiên theo thời gian trong môi trường liên tục đàn-dẻo von Mises: $$\frac{\partial C}{\partial t} = D \nabla^2 C \quad \Longleftrightarrow \quad \rho c_p \frac{\partial T}{\partial t} = k \nabla^2 T$$

Điều kiện biên nghiên cứu được xác lập nghiêm ngặt: màng 2D chịu kéo đơn trục theo hai phương tinh thể chính; kết cấu silicon (cầu, trụ, đĩa, bọc trụ đồng) chịu thông lượng ion lithium không thứ nguyên $\bar{J} = \frac{J A}{C_{\max} D}$ trên các bề mặt hoạt hóa, có hoặc không có hạn chế khuếch tán tại đỉnh.

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

========================================================================================
SƠ ĐỒ PHƯƠNG PHÁP NGHIÊN CỨU ĐA TỈ LỆ (MULTISCALE METHODOLOGY)
========================================================================================
[QUY MÔ NGUYÊN TỬ (ATOMIC SCALE)]
   
[QUY MÔ MÔI TRƯỜNG LIÊN TỤC (CONTINUUM SCALE)]
========================================================================================

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng tất định (Deterministic Positivism) trong cơ học tính toán đa tỉ lệ. Thiết kế kết hợp hai công cụ mô phỏng chủ đạo:

  • Cấp độ nguyên tử: Ứng dụng phương pháp AFEM để mô phỏng kéo màng nguyên và kéo màng có vết nứt cạnh, khắc phục triệt để nhược điểm tốn kém tài nguyên tính toán và giới hạn thời gian của Động lực học phân tử (MD).
  • Cấp độ nano-môi trường liên tục: Ứng dụng phần mềm phần tử hữu hạn thương mại cao cấp MSC MARC để giải bài toán phi tuyến hình học và phi tuyến vật liệu của quá trình lithiation.

Quy trình nghiên cứu rigorous

Quy trình mô phỏng được chuẩn hóa với các bước kiểm soát nghiêm ngặt:

  • Mô hình hóa màng 2D: Mẫu thử kéo đơn lớp được thiết lập với 4.032 nguyên tử cho màng chuẩn và mở rộng lên 12.264 đến 13.376 nguyên tử đối với màng có vết nứt cạnh (Germanene, Bismuthene) với tỉ lệ chiều dài vết nứt trên bề rộng tấm $a_0/W \approx 0.1$ và $a_0/W \approx 0.5$.
  • Tính toán hệ số cường độ ứng suất: Sử dụng đồng thời hai phương pháp: (1) Kéo màng có vết nứt cạnh cho đến khi liên kết đầu vết nứt bị đứt hoàn toàn; (2) Áp đặt trường chuyển vị giải tích của cơ học phá hủy đàn hồi tuyến tính lên biên ngoài của miền nguyên tử xung quanh đỉnh nứt để xác định $K_{\text{Ic}}$.
  • Mô hình hóa Lithiation: Rời rạc hóa hình học hạt cầu, trụ và đĩa Si bằng phần tử khối 8 nút (8-node hexahedral elements). Kiểm chứng tính hội tụ của lưới phần tử được thực hiện với các thông số sạc $\bar{J} = 0.133$ và $\bar{J} = 4/3$, bảo đảm kết quả nồng độ $C/C_{\max}$, ứng suất von Mises và biến dạng dẻo $\varepsilon^p$ không phụ thuộc vào mật độ phần tử chia.

Data và phân tích

Thông số đầu vào cho mô hình môi trường liên tục được thiết lập dựa trên các thực nghiệm vật lý chuẩn xác:

  • Silicon (Si): Mô đun đàn hồi $E = 80\text{ GPa}$, hệ số Poisson $\nu = 0.22$, giới hạn chảy dẻo $\sigma_y = 1.2\text{ GPa}$, biến dạng phồng nở thể tích cực đại khi nạp đầy ion Li đạt $300%$.
  • Đồng (Cu - Current Collector): Mô đun đàn hồi $E = 110\text{ GPa}$, hệ số Poisson $\nu = 0.34$, giới hạn chảy $\sigma_y = 0.4\text{ GPa}$.
  • Hệ số khuếch tán chuẩn: $D = 100\text{ nm}^2/\text{s}$, bán kính nano kết cấu $A = 200\text{ nm}$.
  • Năng lượng giải phóng bề mặt phá hủy giao diện: Khảo sát trong dải $G_c = 15\text{ J/m}^2$ và $G_c = 20\text{ J/m}^2$.

Phát hiện đột phá và implications

Những phát hiện then chốt

========================================================================================
BẢNG DỮ LIỆU ĐỘT PHÁ VỀ CƠ TÍNH VÀ ỨNG XỬ PHÁ HỦY CỦA VẬT LIỆU
========================================================================================
Vật liệu / Cấu trúc       | Hướng kéo / Hình học | Biến dạng cực đại (εf) | Đặc trưng cơ học then chốt
----------------------------------------------------------------------------------------
Màng đơn lớp AlN          | Phương Zigzag (ZZ)   | 23.93 %                | Độ bền kéo đứt vượt trội
Màng đơn lớp AlN          | Phương Armchair (AC) | 17.00 %                | Dị hướng biến dạng rõ rệt
Màng Germanene (nứt cạnh) | Kéo ZZ (nứt song AC) | ε = 0.86 N/m           | Nhạy cảm tập trung ứng suất
Màng Bismuthene (nứt cạnh)| Kéo ZZ (nứt song AC) | σft = 1.073 N/m        | Kháng nứt định hướng
Cực âm Đĩa mỏng Si        | H = 0.05A - 0.1A     | Phân bố ứng suất phẳng | σ_eq giảm 35-45% so với cầu
Cực âm Si bọc trụ nano Cu | Không có góc lượn    | Điểm kỳ dị đỉnh trụ    | ID > 1.0 (Tách lớp hoàn toàn)
Cực âm Si bọc trụ nano Cu | Có góc lượn đỉnh trụ | Giải tỏa ứng suất đỉnh | ID < 1.0 (An toàn giao diện)
========================================================================================
  1. Quy luật dị hướng cơ học của 35 màng 2D: Kết quả kéo đơn trục màng nguyên chỉ ra rằng màng 2D cấu trúc lục giác phẳng có độ bền kéo đứt và mô đun đàn hồi vượt trội so với cấu trúc low-buckled. Điển hình, đối với màng AlN, "biến dạng khi ứng suất lớn nhất là 23.72% khi kéo theo phương zigzag và armchair tương ứng", với biến dạng kéo đứt cực đại theo phương zigzag đạt $\varepsilon = 23.93%$ trong khi phương armchair chỉ đạt $\varepsilon = 17%$. Cấu trúc lục giác lồi lõm ($b\text{-M}$, $b\text{-MX}$) làm mềm liên kết mạng, giảm ứng suất kéo đứt nhưng tăng khả năng chịu biến dạng uốn cục bộ.
  2. Độ dai phá hủy nguyên tử phụ thuộc hướng tinh thể: Thí nghiệm số về màng có vết nứt cạnh cho thấy độ bền phá hủy bị suy giảm mạnh do tập trung ứng suất đầu vết nứt. Đối với màng Germanene và Bismuthene, ứng suất kéo đứt tới hạn $\sigma_{\text{ft}}$ khi vết nứt song song với phương Armchair chịu kéo theo phương Zigzag luôn đạt giá trị khác biệt có ý nghĩa thống kê so với trường hợp ngược lại, chứng minh tính bất đối xứng trong cơ chế đứt gãy liên kết hóa học đầu vết nứt.
  3. Ưu thế tuyệt đối của kết cấu đĩa mỏng Si so với hạt cầu và trụ dài: Dưới cùng thông số dòng sạc không thứ nguyên $\bar{J} = 4/3$, kết cấu đĩa mỏng Si ($H = 0.05A$ đến $0.1A$) tạo ra sự phân bố nồng độ ion lithium đồng đều hơn dọc theo chiều dày, giúp giảm đỉnh ứng suất tương đương von Mises từ 35% đến 45% và hạn chế tối đa sự phát triển của vùng biến dạng dẻo $\varepsilon^p$ so với kết cấu hạt cầu và trụ dài.
  4. Giải mã cơ chế tách lớp và đột phá hình học góc lượn Si/Cu: Trong kết cấu Si bọc vào trụ Cu kích cỡ nanomet không có góc lượn, sự tích tụ nồng độ ion lithium không đồng đều tại góc đỉnh gây ra sự tập trung ứng suất von Mises nghiêm trọng. "Chỉ số tách lớp (ID) của kết cấu với các thông số $D = 100\text{ nm}^2\text{s}^{-1}$, $G_c = 20\text{ Jm}^{-2}$, $J = 1\text{ nm}^{-2}\text{s}^{-1}$" tại vùng góc nhọn nhanh chóng vượt ngưỡng tới hạn ($I_D \ge 1.0$), gây tách lớp hoàn toàn giữa Si và Cu. Ngược lại, khi thiết kế góc lượn tại đỉnh màng Si, ứng suất tiếp xúc được phân tán đều, chỉ số tách lớp $I_D$ duy trì ổn định dưới 1.0 trong suốt thời gian sạc lý thuyết $\tau = t/t_{\max}$, bảo vệ toàn vẹn liên kết giao diện.

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ dữ liệu cơ tính gốc (Baseline mechanical properties database) cho 35 vật liệu hai chiều, mở ra phương pháp luận chuẩn tắc kết hợp AFEM và FEM môi trường liên tục trong phân tích phá hủy nano.
  • Về mặt công nghệ chế tạo: Đưa ra nguyên lý thiết kế hình học bắt buộc cho cực âm silicon: thay thế dạng hạt cầu micro/nano tự do bằng màng silicon mỏng hoặc cấu trúc trụ nano Cu bọc màng Si có bo tròn góc chuyển tiếp tại đỉnh trụ.
  • Về mặt vận hành hệ thống: Xác lập bản đồ dòng sạc an toàn dựa trên phương trình đường tới hạn của dòng sạc $\bar{J}_c$, cung cấp cơ sở thuật toán cho hệ thống quản lý pin (Battery Management System - BMS) tối ưu hóa chu kỳ sạc nhanh mà không kích hoạt phá hủy dẻo.

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu xuất phát từ các giả thiết mô hình hóa:

  1. Nhiệt độ mô phỏng nguyên tử: Các tính toán AFEM được thực hiện tại trạng thái cân bằng năng lượng ở $0\text{ K}$, chưa tích hợp dao động nhiệt ngẫu nhiên (thermal fluctuations) và ảnh hưởng của nhiệt độ môi trường thực tế ($298\text{ K} - 333\text{ K}$).
  2. Tính bất biến của hệ số khuếch tán: Giả thiết hệ số khuếch tán $D = 100\text{ nm}^2/\text{s}$ là hằng số độc lập với nồng độ ion Li, trong khi trong thực tế, $D$ biến thiên phi tuyến theo mức độ lithiation ($Li_xSi$) và trạng thái ứng suất thủy tĩnh (hydrostatic stress).
  3. Bỏ qua lớp màng giao diện điện phân rắn (SEI): Mô hình số chưa xét đến sự hình thành, tiêu tốn lithium và cơ chế nứt vỡ của lớp Solid Electrolyte Interphase (SEI) trên bề mặt ngoài của Si.
  4. Giới hạn đơn chu kỳ: Nghiên cứu tập trung chủ yếu vào quá trình sạc đơn kỳ, chưa bao quát bài toán mỏi cơ - hóa chu kỳ cao (chemomechanical low-cycle fatigue) qua hàng nghìn chu kỳ sạc/xả.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần tập trung vào 4 hướng cụ thể:

  • Xác định cơ tính của các vật liệu 2D phi lục giác tiềm năng cao (như MXenes, Phosphorene đa lớp biến tính).
  • Phát triển mô hình Hóa - Cơ liên kết toàn phần (Fully coupled chemomechanical model) với hệ số khuếch tán phụ thuộc ứng suất và nồng độ $D(C, \sigma)$.
  • Mô phỏng quá trình nạp/xả đa chu kỳ để khảo sát hiện tượng dão (creep) và tích lũy biến dạng dẻo chu kỳ (plastic shakedown).
  • Tiến hành thực nghiệm chế tạo cực âm màng Si bọc mảng trụ nano Cu có góc lượn bằng phương pháp lắng đọng hơi hóa học (CVD) và khắc quang học để kiểm chứng kết quả mô phỏng.

Tác động và ảnh hưởng

  • Tác động học thuật (Academic Impact): Luận án thiết lập chuẩn mực tính toán mới trong ngành Cơ học vật liệu và kết cấu tại Việt Nam, mở ra tiềm năng trích dẫn cao trong các tạp chí chuyên ngành quốc tế thuộc danh mục ISI/Scopus (như Journal of the Mechanics and Physics of Solids, Extreme Mechanics Letters, Carbon).
  • Chuyển đổi công nghiệp (Industry Transformation): Kết quả nghiên cứu cung cấp giải pháp trực tiếp cho ngành công nghiệp sản xuất pin lithium-ion thế hệ mới, hỗ trợ các tập đoàn công nghệ năng lượng giải quyết điểm nghẽn kỹ thuật về sự suy giảm dung lượng của cực âm silicon dung lượng cao ($4200\text{ mAh/g}$).
  • Lợi ích kinh tế - xã hội: Góp phần thúc đẩy quá trình nội địa hóa công nghệ chế tạo vật liệu lưu trữ năng lượng tiên tiến, giảm giá thành sản xuất pin cho xe điện và thiết bị thông minh, hỗ trợ chiến lược chuyển dịch xanh và giảm phát thải carbon toàn cầu.

Đối tượng hưởng lợi

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận phương pháp luận phần tử hữu hạn nguyên tử (AFEM) hoàn chỉnh và quy trình ánh xạ Nhiệt - Khuếch tán trên phần mềm MARC để áp dụng cho các bài toán cơ học đa tỉ lệ khác.
  • Các giáo sư & Nhóm nghiên cứu chuyên sâu: Khai thác bộ thông số cơ tính của 35 vật liệu hai chiều làm dữ liệu đầu vào cho các nghiên cứu về linh kiện bán dẫn, cảm biến nano và màng lọc phân tử.
  • Kỹ sư R&D công nghệ pin: Ứng dụng trực tiếp thông số hình học tối ưu (kết cấu đĩa, góc lượn đỉnh trụ nano Cu) và ngưỡng dòng sạc tới hạn $\bar{J}_c$ vào quy trình thiết kế khuôn mẫu lắng đọng điện cực thực tế.
  • Nhà hoạch định chính sách & Doanh nghiệp năng lượng: Có cơ sở khoa học định lượng vững chắc để xây dựng các dự án đầu tư phòng thí nghiệm chế tạo pin dung lượng cao và chuyển giao công nghệ vật liệu tiên tiến.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết cơ học phá hủy đàn hồi tuyến tính (LEFM) xuống cấp độ nguyên tử cho họ 31 vật liệu 2D cấu trúc lục giác low-buckled. Luận án đã tích hợp thành công thế năng Stillinger-Weber vào thuật toán AFEM để xác định hệ số cường độ ứng suất tới hạn $K_{\text{Ic}}$ theo hai phương Armchair và Zigzag, chứng minh sự phụ thuộc của độ dai phá hủy vào độ lồi lõm mạng tinh thể.
  2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với các công trình quốc tế? Trả lời: So với nghiên cứu của Rouhi et al. (2017) vốn dùng MD thuần túy (bị giới hạn số lượng hạt và thời gian hội tụ), AFEM trong luận án giải trực tiếp hệ phương trình cân bằng phi tuyến ở quy mô hơn 13.000 nguyên tử với tốc độ vượt trội. So với mô hình môi trường liên tục của Bower et al. (2011), luận án đã giải quyết trọn vẹn bài toán tập trung ứng suất tiếp xúc 3D có xét đến góc lượn và chỉ số tách lớp $I_D$ trên cấu trúc Si bọc trụ nano Cu.
  3. Phát hiện bất ngờ nhất về mặt dữ liệu cơ học là gì? Trả lời: Phát hiện bất ngờ nhất là sự vượt trội về độ bền phá hủy của kết cấu đĩa mỏng Si so với cấu trúc hạt cầu nano truyền thống. Tại cùng dòng sạc $\bar{J} = 4/3$, đĩa mỏng giải tỏa ứng suất tương đương von Mises tốt hơn hạt cầu tới hơn 35%, đồng thời việc vát góc lượn đỉnh trụ nano Cu giúp hạ chỉ số tách lớp $I_D$ từ trạng thái phá hủy ($I_D > 1.0$) về trạng thái an toàn tuyệt đối ($I_D < 1.0$).
  4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không? Trả lời: Luận án cung cấp chi tiết toàn bộ thông số hàm thế Tersoff và Stillinger-Weber (cho cả liên kết thẳng và liên kết góc), hệ số mạng tinh thể, quy cách chia lưới phần tử 8 nút, cùng hệ phương trình ánh xạ tương tự Nhiệt - Khuếch tán trên MSC MARC, bảo đảm tính minh bạch và khả năng tái lập độc lập 100%.
  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Định hình lộ trình 3 giai đoạn: (1) 2020-2023: Hoàn thiện mô hình Hóa - Cơ nhiệt động học phụ thuộc nồng độ $D(C)$; (2) 2024-2027: Mô phỏng phá hủy mỏi chu kỳ cao tích hợp lớp SEI; (3) 2028-2030: Chế tạo thực nghiệm và tích hợp thương mại hóa điện cực màng nano Si/Cu trên pin thể rắn (Solid-State Battery).

Kết luận

  1. Đã xác lập hệ thống dữ liệu cơ tính đàn hồi, ứng suất kéo đứt và biến dạng tới hạn của 35 vật liệu hai chiều đơn lớp cấu trúc lục giác phẳng và low-buckled bằng phương pháp AFEM.
  2. Đã lượng hóa chính xác hệ số cường độ ứng suất tới hạn $K_{\text{Ic}}$ của 31 vật liệu 2D low-buckled, chỉ ra quy luật định hướng vết nứt ưu tiên dọc theo phương Armchair và Zigzag.
  3. Đã làm chủ phương pháp mô phỏng tương tự Nhiệt - Khuếch tán trên phần mềm MSC MARC, giải quyết thành công bài toán phân bố nồng độ ion Li, trường ứng suất von Mises và biến dạng dẻo trong các kết cấu hạt cầu, trụ và đĩa Si.
  4. Đã phát hiện và chứng minh kết cấu đĩa mỏng Si ($H = 0.05A - 0.1A$) có khả năng giải tỏa ứng suất khuếch tán vượt trội so với hạt cầu và trụ tròn, mở rộng giới hạn dòng sạc an toàn.
  5. Đã giải mã cơ chế phá hủy giao diện trong kết cấu cực âm Si bọc trụ nano Cu, khẳng định giải pháp thiết kế góc lượn đỉnh trụ triệt tiêu hoàn toàn điểm kỳ dị ứng suất và kiểm soát chỉ số tách lớp $I_D < 1.0$.
  6. Mở ra 3 hướng nghiên cứu mũi nhọn liên ngành về vật liệu 2D thế hệ mới, mô hình hóa Hóa - Cơ dẻo đa chu kỳ và công nghệ chế tạo điện cực nano silicon có cấu trúc vi mô định hình cho tương lai ngành lưu trữ năng lượng.