Tổng quan luận án

Luận án tiến sĩ "Mô phỏng Monte Carlo cho hệ vi cầu từ tính $\text{Fe}_3\text{O}_4\text{/Poly(Glycidyl Methacrylate)}$" được thực hiện bởi nghiên cứu sinh Nguyễn Thanh Hoàng, dưới sự hướng dẫn khoa học của TS. Nguyễn Mạnh Tuấn và TS. Trần Hoàng Hải, chuyên ngành Vật liệu quang học, quang điện tử và quang tử (mã số: 9440127) tại Học viện Khoa học và Công nghệ, Viện Hàn lâm Khoa học và Công nghệ Việt Nam (2024).

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

Trong những năm gần đây, vật liệu tổ hợp (nanocomposite) gồm các hạt nano từ tính phân tán trong vỏ bọc polymer tương thích sinh học được ứng dụng trong y sinh học (như hệ phân phối thuốc hướng đích, tăng thân nhiệt cục bộ, hệ vi lưu lab-on-chip) và xử lý môi trường (hấp phụ ion kim loại nặng). Để đưa vật liệu vào thực tế có sử dụng từ trường ngoài, việc phân tích và đánh giá các đặc trưng vi từ nội tại như dị hướng từ hiệu dụng, mô-men xoắn, nhiệt độ khóa ($T_B$) và tương tác lưỡng cực - lưỡng cực đóng vai trò chi phối.

Tuy nhiên, các phương pháp thực nghiệm đo đạc vĩ mô gặp nhiều rào cản: chi phí cao, tốn thời gian, độ lặp lại bị giới hạn bởi quá trình chế tạo, khó thực hiện ở các điều kiện giới hạn (nhiệt độ $0\text{ K}$, từ trường cực lớn hoặc xung laser siêu nhanh), đồng thời khó định lượng trực tiếp các đại lượng ở cấp độ vi mô của từng hạt. Ngược lại, các nghiên cứu lý thuyết hoặc mô phỏng thuần túy thường phải dựa trên các giả định xấp xỉ đơn giản hóa, bỏ qua cấu trúc thực hoặc thiếu sự đối chiếu với mẫu thực nghiệm cụ thể. Khoảng trống nghiên cứu được tác giả xác định là sự thiếu hụt một quy trình nghiên cứu tích hợp đa chiều, kết nối chặt chẽ giữa thực nghiệm chế tạo, tính toán lý thuyết số và mô phỏng máy tính đa quy mô nhằm định lượng các thông số vi từ nội tại của hệ vi cầu tổ hợp.

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

  1. Xây dựng một mô hình nghiên cứu đa chiều kết hợp giữa thực nghiệm, lý thuyết từ tính và mô phỏng máy tính nhằm khảo sát, định lượng và đánh giá chi tiết các tính chất vi từ nội tại của hệ vật liệu tổ hợp.
  2. Áp dụng mô hình trên đối tượng cụ thể là hệ vi cầu từ tính $\text{Fe}_3\text{O}_4\text{/Poly(glycidyl methacrylate)}$ ($\text{Fe}_3\text{O}_4\text{/PGMA}$) chứa các hạt nano từ $\text{Fe}_3\text{O}_4$.
  3. Khai thác các kết quả tính toán và mô phỏng số để cung cấp các thông tin tiên nghiệm cho thực nghiệm chế tạo và ứng dụng định hướng y sinh cũng như xử lý môi trường.

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

  • Đối tượng nghiên cứu: Hệ vật liệu vi cầu tổ hợp từ tính $\text{Fe}_3\text{O}_4\text{/Poly(glycidyl methacrylate)}$ gồm hạt nano oxit sắt từ $\text{Fe}_3\text{O}_4$ phân tán trong ma trận polymer nền poly(glycidyl methacrylate) (PGMA) đã biến tính amine.
  • Phạm vi không gian: Khảo sát cấu trúc và tính chất từ từ cấp độ vi mô nguyên tử (spin nguyên tử, ô mạng cơ sở $\text{Fe}_3\text{O}_4$), cấp độ hạt nano đơn lẻ (đường kính dưới $50\text{ nm}$), cụm cluster hạt phân bố trong vi cầu polymer (đường kính vi cầu cấp độ micromet), đến mẫu vĩ mô.
  • Phạm vi điều kiện vật lý: Nhiệt độ khảo sát từ $0\text{ K}$ đến nhiệt độ phòng ($300\text{ K}$) và các dải nhiệt độ lân cận; từ trường ngoài biến thiên phục vụ khảo sát từ hóa bão hòa, đường cong từ nhiệt ZFC/FC và mô-men xoắn từ.

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

Các hướng nghiên cứu trước được tác giả điểm lại

Tác giả hệ thống hóa các hướng nghiên cứu trong nước và quốc tế liên quan đến vật liệu tổ hợp polymer/hạt nano từ và các công cụ mô phỏng số:

  • Vật liệu nền polymer thông minh và PGMA: Các nghiên cứu về hydrogel thông minh (IHP) đáp ứng pH, nhiệt độ, từ trường trong phân phối thuốc (dẫn xuất các công trình về PLGA tải Doxorubicin). PGMA được xác định là polymer nhiệt dẻo họ poly(acrylic acid) có monomer là glycidyl methacrylate (GMA), chứa vòng epoxy dễ biến tính với các nhóm chức như $-\text{NH}_2$, $-\text{COOH}$, $-\text{SH}$. Tác giả điểm lại các công trình của Zohreh và cộng sự (2021) về vi cầu hai lõi - vỏ kép $\text{Fe}_3\text{O}_4\text{/PGMA@PEG}$ tải DOX đạt hiệu suất giải phóng trên $80%$ kháng tế bào MCF-7; công trình vi cầu nanocomposite lai hóa silica xốp/PGMA-CD đáp ứng đa kích thích.
  • Ứng dụng y sinh và môi trường của hạt nano oxit sắt (SPION/MNPs): Điểm lại các nghiên cứu ứng dụng hạt nano từ trong tăng thân nhiệt cục bộ tiêu diệt tế bào u não ác tính ở dải nhiệt độ $39 - 45^\circ\text{C}$; nghiên cứu về hệ vi lưu (microfluidics/lab-on-chip) kết hợp nhíp mô-men xoắn từ tính để định lượng Protein G, kháng thể IgG, troponin tim; nghiên cứu của các nhóm tác giả về vi cầu PMDA-PGMA hấp phụ ion $\text{Pd}^{2+}$ và vi cầu PPA-$\text{Fe}_3\text{O}_4\text{/PGMA}$ hấp phụ Uranium từ nước rỉ quặng.
  • Lý thuyết từ và mô hình khử từ: Osborn (1945) với hệ phương trình giải tích xác định hệ số khử từ $N_d$ cho các hạt từ hình ellipsoid.
  • Mô phỏng máy tính trong khoa học vật liệu: Điểm lại các phương pháp đa quy mô gồm lý thuyết phiếm hàm mật độ (DFT), động học phân tử (MD), mô phỏng mô hình spin nguyên tử (Atomistic Spin Model - ASM) dựa trên phương trình Landau-Lifshitz-Gilbert (LLG) và mô phỏng Monte Carlo (MC).

Khoảng trống luận án lựa chọn giải quyết

Các công trình trước đây thường phân tách riêng biệt: nhóm nghiên cứu thực nghiệm tập trung chế tạo và chỉ đo đạc các đường cong từ hóa vĩ mô (VSM); nhóm tính toán lý thuyết thường mô phỏng các hệ hạt đơn lẻ hoặc phân bố tuần hoàn lý tưởng, bỏ qua ảnh hưởng của ma trận polymer và tương tác lưỡng cực thực tế giữa các hạt trong vi cầu. Luận án lựa chọn giải quyết khoảng trống này bằng cách: dùng dữ liệu thực nghiệm (kết quả đo VSM, XRD, TEM) kết hợp với các lý thuyết từ số (định luật tiếp cận bão hòa, hàm Langevin) để xác định chính xác các thông số đầu vào thực (phân bố kích thước hạt logarit-chuẩn, hằng số dị hướng hiệu dụng $K_{eff}$). Từ đó, đưa các thông số này vào các thuật toán mô phỏng Monte Carlo và mô hình spin nguyên tử trên hệ tính toán hiệu năng cao để tái hiện sát thực tế hành vi vi từ của hệ vi cầu.

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

Cơ sở lý thuyết và các khung năng lượng vi từ

Hạt nano $\text{Fe}_3\text{O}_4$ có cấu trúc tinh thể ferrite spinel nghịch. Trong mỗi ô mạng cơ sở, các spin từ của 8 ion $\text{Fe}^{3+}$ ở vị trí tứ diện ngược chiều và triệt tiêu với spin của 8 ion $\text{Fe}^{3+}$ ở vị trí bát diện; mô-men từ thuần của phân tử do các ion $\text{Fe}^{2+}$ tại vị trí bát diện đóng góp (độ lớn $4\mu_B$).

Hệ vi từ được mô tả thông qua năng lượng tổng cộng $E_{total}$, là điều kiện biên để cực tiểu hóa trạng thái cân bằng spin: $$E_{total} = E_{ex} + E_k + E_d + E_Z + E_\sigma$$

Trong đó:

  1. Năng lượng trao đổi ($E_{ex}$): Xác định theo mô hình Heisenberg: $$E_{ex} = -2 \sum_{i<j} J_{ij} \vec{S}i \cdot \vec{S}j = -2 J S^2 \sum{i<j} \cos\psi{ij}$$
  2. Năng lượng dị hướng ($E_k$): Bao gồm dị hướng từ tinh thể và dị hướng hình dạng. Dị hướng đơn trục khai triển theo chuỗi chẵn: $$E_{uni} = K_0 + K_1\sin^2\theta + K_2\sin^4\theta + \dots$$ Với dị hướng lập phương của $\text{Fe}3\text{O}4$: $$E{cub} = K{c0}(c_x^2 + c_y^2 + c_z^2) + K_{c1}(c_x^2 c_y^2 + c_y^2 c_z^2 + c_z^2 c_x^2) + K_{c2}(c_x^2 c_y^2 c_z^2)$$ Tại $300\text{ K}$, giá trị hằng số dị hướng của $\text{Fe}_3\text{O}4$ khối được dẫn chiếu là $K_1 = -1{,}35 \times 10^5\text{ erg/cm}^3$ và $K_2 = -0{,}44 \times 10^5\text{ erg/cm}^3$. Năng lượng dị hướng hình dạng liên hệ qua thừa số khử từ: $$K{eff} = \frac{N_b - N_a}{2} M_s^2$$
  3. Năng lượng tĩnh từ / tương tác lưỡng cực - lưỡng cực ($E_d$): $$\vec{H}{dip}^i = \frac{\mu_0}{4\pi} \sum{j \neq i} \left[ \frac{3(\vec{m}j \cdot \vec{r}{ij})\vec{r}{ij}}{r{ij}^5} - \frac{\vec{m}j}{r{ij}^3} \right], \quad E_d = \sum_i \vec{m}i \cdot \vec{H}{dip}^i$$
  4. Năng lượng Zeeman ($E_Z$): Tương tác với từ trường ngoài $\vec{H}_{app}$: $$E_Z = -\mu_0 M_s V \vec{e}M \cdot \vec{H}{app}$$
  5. Năng lượng từ giảo ($E_\sigma$): Hình thành do tương tác spin - quỹ đạo dưới tác dụng của ứng suất cơ học $\sigma$: $$E_\sigma = \int_V \frac{3}{2}\lambda_s \sigma \sin^2\theta , dV$$

Phương pháp nghiên cứu thực nghiệm

  • Tổng hợp hóa học: Trùng hợp huyền phù PGMA từ monomer glycidyl methacrylate (GMA) có mặt chất khơi mào AIBN; biến tính bề mặt PGMA bằng ethylene diamine để tạo amino-PGMA mang nhóm chức hoạt tính $-\text{NH}_2$; kết tủa tại chỗ (in-situ co-precipitation) các muối sắt $\text{Fe}^{2+}/\text{Fe}^{3+}$ trong môi trường kiềm để tạo hạt nano $\text{Fe}_3\text{O}_4$ định hình trên và trong mạng lưới vi cầu polymer; tổng hợp đối chứng mẫu hạt nano $\text{Fe}_3\text{O}_4$ trần.
  • Kỹ thuật phân tích cấu trúc và hình thái: Nhiễu xạ tia X (XRD) xác định pha tinh thể spinel; phổ hồng ngoại biến đổi Fourier (FTIR) xác định liên kết hóa học và nhóm chức epoxy, amine; kính hiển vi điện tử quét (SEM) và truyền qua (TEM/HRTEM) khảo sát hình thái hạt, kích thước vi cầu và sự phân bố hạt nano từ.
  • Phép đo từ tính vĩ mô: Từ kế mẫu rung (VSM) đo đường cong từ hóa $M(H)$ tại nhiệt độ phòng ($300\text{ K}$) với từ trường quét lên đến vài Tesla, đo từ hóa phụ thuộc nhiệt độ theo quy trình làm lạnh không từ trường (ZFC) và làm lạnh có từ trường (FC).
  • Phép đo hấp phụ môi trường: Quang phổ hấp thụ nguyên tử (AAS) xác định dung lượng hấp phụ đẳng nhiệt ion kim loại nặng chì ($\text{Pb}^{2+}$) trong dung dịch nước.

Phương pháp tính toán lý thuyết số và mô phỏng máy tính

  • Định luật tiếp cận bão hòa (Law of Approach to Saturation - LAS): Dùng để xử lý dữ liệu vùng từ trường cao của đường cong $M(H)$ thực nghiệm, dựa trên mối quan hệ tuyến tính giữa $M/M_s$ và $1/H^2$ nhằm chiết xuất giá trị hằng số dị hướng từ hiệu dụng $K_{eff}$.
  • Làm khớp hàm Langevin (Langevin Fitting): Khớp đường cong từ hóa siêu thuận từ với hàm phân bố kích thước hạt logarit-chuẩn (log-normal distribution) được lập trình bằng Matlab, trích xuất đường kính trung bình và độ lệch chuẩn của hạt nano từ.
  • Mô phỏng Monte Carlo (MC): Sử dụng thuật toán Metropolis và mô phỏng Monte Carlo hằng số áp suất để tạo cấu hình phân bố không gian (clustering) của các hạt nano $\text{Fe}_3\text{O}_4$ bên trong vi cầu PGMA theo các tỷ số nén (mật độ không gian) khác nhau.
  • Mô phỏng mô hình spin nguyên tử (Atomistic Spin Model): Sử dụng gói phần mềm VAMPIRE, thiết lập cấu trúc tinh thể từ ô mạng cơ sở của $\text{Fe}_3\text{O}_4$, giải phương trình động học spin Landau-Lifshitz-Gilbert (LLG) sử dụng thuật toán tích phân Heun và phương pháp Monte Carlo cho spin nguyên tử để mô phỏng đường cong từ nhiệt $M(T)$, đường cong từ hóa $M(H)$ tại $0\text{ K}$ và $300\text{ K}$, đường cong ZFC/FC, cấu hình spin cục bộ và mô-men xoắn từ.
  • Hạ tầng tính toán: Thực hiện trên hệ máy tính hiệu năng cao (HPC) của Trung tâm Tin học và Tính toán (Viện Hàn lâm Khoa học và Công nghệ Việt Nam - VAST) và hệ điện toán lưới HCMIP.
Nhóm phương pháp Công cụ / Kỹ thuật cụ thể Mục tiêu dữ liệu trích xuất
Thực nghiệm chế tạo Trùng hợp dung dịch, amine hóa, kết tủa tại chỗ Mẫu vi cầu $\text{Fe}_3\text{O}_4\text{/PGMA}$, hạt $\text{Fe}_3\text{O}_4$ trần
Đặc trưng thực nghiệm XRD, FTIR, SEM, TEM, VSM, AAS Cấu trúc pha, liên kết hóa học, hình thái vi cầu, $M_s$, dung lượng hấp phụ $\text{Pb}^{2+}$
Tính toán lý thuyết số LAS ($M/M_s - 1/H^2$), Langevin fitting (Matlab) Hằng số dị hướng $K_{eff}$, hàm phân bố kích thước logarit-chuẩn
Mô phỏng máy tính Monte Carlo áp suất hằng số, VAMPIRE (LLG-Heun, MC) Cấu hình không gian cụm hạt, đường cong $M(T)$, $M(H)$, ZFC/FC, mô-men xoắn

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

Chương 1. Vi cầu từ tính Poly(glycidyl methacrylate)/Fe3O4

Chương 1 trình bày tổng quan cấu trúc, tính chất vật liệu thành phần, các cơ chế năng lượng từ, tiềm năng ứng dụng và cơ sở lý thuyết mô phỏng máy tính trong khoa học vật liệu.

  • Vật liệu nền PGMA: PGMA là polymer nhiệt dẻo họ poly(acrylic acid) tổng hợp từ monomer glycidyl methacrylate (GMA), công thức phân tử $(\text{C}7\text{H}{10}\text{O}_3)_n$, khối lượng phân tử trung bình khoảng $20.153\text{ g/mol}$, nhiệt độ nóng chảy $274 - 280^\circ\text{C}$, khối lượng riêng $0{,}805\text{ g/ml}$ tại $25^\circ\text{C}$, độ thấu quang đạt $92%$. PGMA có đặc tính tự phân hủy sinh học, nhạy pH (trương nở hoặc co lại khi thay đổi pH), nhóm chức epoxy dễ mở vòng để liên kết với các phối tử sinh học.
  • Hạt nano từ $\text{Fe}_3\text{O}_4$: Cấu trúc ferrite spinel nghịch với các cation $\text{Fe}^{3+}$ chiếm 8 vị trí tứ diện và 8 vị trí bát diện, ion $\text{Fe}^{2+}$ chiếm 8 vị trí bát diện tạo mô-men từ $4\mu_B$/phân tử. Đường kính tới hạn đơn ô-men $D_c$ được mô tả qua công thức: $$D_c = \frac{72\sqrt{A K}}{\mu_0 M_s^2}$$ Dưới kích thước tới hạn ($<50\text{ nm}$), hạt thể hiện tính chất siêu thuận từ. Phân tích chi tiết 5 dạng năng lượng vi từ: năng lượng trao đổi $E_{ex}$, năng lượng dị hướng tinh thể $E_{uni}$ và $E_{cub}$ (phương $[111]$ là trục dễ, $[100]$ là trục khó ở nhiệt độ $>130\text{ K}$), năng lượng dị hướng hình dạng dựa trên hệ số khử từ Osborn ($N_x + N_y + N_z = 4\pi$), năng lượng tĩnh từ lưỡng cực - lưỡng cực, năng lượng Zeeman và năng lượng từ giảo.
  • Ứng dụng tiềm năng:
    1. Hệ phân phối thuốc hướng đích (TDDS): Vi cầu PGMA mang thuốc Doxorubicin (DOX), điều hướng bằng gradient từ trường ngoài; tại khối u có pH axit ($6{,}2$ so với mô lành $7{,}4$), vỏ PGMA trương nở giải phóng thuốc.
    2. Tăng thân nhiệt cục bộ: Hạt nano từ sinh nhiệt dưới tác dụng của từ trường xoay chiều đạt dải nhiệt độ $39 - 45^\circ\text{C}$ tiêu diệt tế bào ác tính.
    3. Hệ vi lưu (microfluidics/lab-on-chip): Kênh dẫn PDMS gắn vi cầu từ tính nhận biết tế bào MCF7 và protein BSA thông qua nhíp mô-men xoắn từ tính.
    4. Hấp phụ ion kim loại nặng: Tách loại $\text{Pd}^{2+}$, Uranium và chì bằng cơ chế trao đổi ion và lực hút tĩnh điện với nhóm amine.
  • Mô phỏng máy tính: Phân loại hệ thống gián đoạn/liên tục; vai trò của mô hình vật lý và mô hình toán học; so sánh giải pháp giải tích với mô phỏng số trên máy tính.

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

Chương 2 mô tả chi tiết quy trình thực nghiệm tổng hợp mẫu, các phép đo phân tích, cùng thuật toán và công cụ tính toán số, mô phỏng máy tính.

  • Quy trình tổng hợp:
    • Trùng hợp PGMA: Monomer GMA được trùng hợp với chất khơi mào AIBN trong dung môi hữu cơ tạo polymer dạng hạt vi cầu.
    • Amine hóa bề mặt: Sử dụng ethylene diamine phản ứng với vòng epoxy của PGMA tạo amino-PGMA chứa các nhóm chức $-\text{NH}_2$.
    • Tổng hợp vi cầu $\text{Fe}_3\text{O}_4\text{/PGMA}$ và hạt $\text{Fe}_3\text{O}_4$ trần: Hòa tan muối $\text{Fe}^{2+}$ và $\text{Fe}^{3+}$ theo tỷ lệ thích hợp, kết tủa tại chỗ bằng dung dịch $\text{NH}_4\text{OH}$ tạo pha $\text{Fe}_3\text{O}_4$ cố định trong khung mạng amino-PGMA; mẫu hạt $\text{Fe}_3\text{O}_4$ trần được tổng hợp bằng phản ứng đồng kết tủa không có mặt polymer.
  • Phương pháp phân tích thực nghiệm: Đo XRD xác định cấu trúc mạng; FTIR kiểm tra nhóm chức; SEM/TEM quan sát kích thước vi cầu và sự phân bố hạt nano oxit sắt; VSM đo đường cong từ hóa $M(H)$ tại $300\text{ K}$ và đường cong ZFC/FC; AAS đo hàm lượng kim loại nặng sau hấp phụ.
  • Tính toán lý thuyết số:
    • Định luật tiếp cận bão hòa (LAS): Áp dụng phương trình cho vùng từ trường cao: $$M = M_s \left(1 - \frac{a}{H} - \frac{b}{H^2}\right) + \chi H$$ Ở vật liệu đẳng hướng ngẫu nhiên, hệ số $b$ liên hệ trực tiếp với hằng số dị hướng hiệu dụng: $b = \frac{8}{105}\frac{K_{eff}^2}{\mu_0^2 M_s^2}$.
    • Làm khớp hàm Langevin: Mô tả từ hóa siêu thuận từ kết hợp hàm phân bố kích thước logarit-chuẩn: $$P(D) = \frac{1}{\sqrt{2\pi}\sigma D} \exp\left(-\frac{\ln^2(D/D_0)}{2\sigma^2}\right)$$ Lập trình trên môi trường Matlab để tối ưu hóa việc xác định đường kính hạt trung bình $D_0$ và độ phân tán $\sigma$.
  • Mô phỏng máy tính:
    • Mô phỏng Monte Carlo áp suất hằng số: Định hình cụm hạt nano $\text{Fe}_3\text{O}_4$ bên trong không gian giới hạn của vi cầu polymer theo các mật độ không gian khác nhau.
    • Mô hình spin nguyên tử (ASM): Dựng tinh thể $\text{Fe}_3\text{O}_4$ từ ô mạng cơ sở lập phương, gắn các giá trị spin cục bộ và tích phân trao đổi giữa các vị trí bát diện - tứ diện; giải phương trình Landau-Lifshitz-Gilbert (LLG): $$\frac{\partial \vec{S}_i}{\partial t} = -\frac{\gamma}{1+\alpha^2} \left[ \vec{S}i \times \vec{H}{eff}^i + \alpha \vec{S}_i \times (\vec{S}i \times \vec{H}{eff}^i) \right]$$ Sử dụng giải thuật Heun và mô phỏng Monte Carlo spin trên gói phần mềm VAMPIRE. Thiết lập các kịch bản có và không có tương tác lưỡng cực - lưỡng cực, thay đổi khoảng cách liên hạt ($0{,}1 - 10\text{ nm}$) và kích thước hạt.

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

Chương 3 trình bày các phát hiện thực nghiệm, kết quả tính toán số và kết quả mô phỏng máy tính.

  • Kết quả thực nghiệm và tính toán số:
    • XRD: Giản đồ nhiễu xạ của $\text{Fe}_3\text{O}_4$ trần và vi cầu $\text{Fe}_3\text{O}_4\text{/PGMA}$ đều xuất hiện đầy đủ các đỉnh đặc trưng của pha spinel nghịch của $\text{Fe}_3\text{O}_4$ chuẩn, không có đỉnh phụ của tạp pha oxit sắt khác.
    • FTIR: Phổ của vi cầu $\text{Fe}_3\text{O}_4\text{/PGMA}$ xác nhận sự mở vòng epoxy của PGMA và sự xuất hiện của liên kết peptide/amine, đồng thời có đỉnh dao động liên kết $\text{Fe}-\text{O}$ tại vùng tần số thấp ($\sim 570\text{ cm}^{-1}$), chứng minh hạt nano từ đã liên kết chặt chẽ với khung polymer.
    • SEM và TEM: Ảnh SEM cho thấy vi cầu có dạng hình cầu tương đối đồng đều; ảnh TEM chỉ ra các hạt nano $\text{Fe}_3\text{O}_4$ có kích thước trung bình cỡ nanomet phân bố đều bên trong nền PGMA, có lớp nhiễu loạn từ tính trên bề mặt hạt.
    • Tính chất từ vĩ mô và LAS: Đường cong $M(H)$ tại $300\text{ K}$ cho thấy mẫu biểu hiện tính chất siêu thuận từ với độ từ dư và lực kháng từ xấp xỉ bằng không. Độ từ hóa bão hòa $M_s$ của vi cầu thấp hơn so với hạt $\text{Fe}_3\text{O}4$ trần do sự có mặt của lớp polymer không mang từ tính. Khớp tuyến tính $M/M_s$ theo $1/H^2$ theo định luật LAS xác định được hằng số dị hướng từ hiệu dụng $K{eff}$.
    • Khớp hàm Langevin: Kết quả khớp đường cong $M(H)$ thực nghiệm với hàm Langevin trên Matlab xác định được phân bố kích thước hạt nano thực tế trong vi cầu tuân theo quy luật logarit-chuẩn.
    • Ứng dụng xử lý nước: Đường hấp phụ đẳng nhiệt của vi cầu $\text{Fe}_3\text{O}_4\text{/PGMA}$ đối với ion chì ($\text{Pb}^{2+}$) khẳng định khả năng tách loại ion kim loại nặng nhờ các nhóm amine tự do trên bề mặt vi cầu, thu hồi vật liệu nhanh chóng bằng nam châm ngoài.
  • Kết quả mô phỏng máy tính:
    • Mô phỏng Monte Carlo cụm hạt: Tái tạo thành công phân bố không gian cụm hạt nano $\text{Fe}_3\text{O}_4$ định hình trong vi cầu PGMA theo các tỷ số nén khác nhau.
    • Đường cong từ nhiệt $M(T)$: So sánh hai phương pháp LLG-Heun và Monte Carlo spin cho thấy sự suy giảm độ từ hóa theo nhiệt độ phù hợp tốt giữa hai cách tiếp cận.
    • Ảnh hưởng của tương tác lưỡng cực - lưỡng cực: Tại $0\text{ K}$, khi có tương tác lưỡng cực, đường cong từ hóa có độ từ dư và lực kháng từ lớn hơn rõ rệt so với khi bỏ qua tương tác. Ở $300\text{ K}$, tương tác lưỡng cực làm chậm quá trình đạt tới bão hòa từ của hệ hạt.
    • Ảnh hưởng của khoảng cách liên hạt: Khảo sát các khoảng cách từ $0{,}1\text{ nm}$ đến $10\text{ nm}$. Khi khoảng cách rất nhỏ ($0{,}1\text{ nm}$), tương tác lưỡng cực rất mạnh chi phối định hướng spin; khi khoảng cách tăng lên $10\text{ nm}$, tương tác lưỡng cực suy giảm nhanh theo quy luật $1/r^3$, tính chất vi từ tiến dần về trạng thái hạt độc lập.
    • Ảnh hưởng của kích thước hạt: Tăng đường kính hạt làm tăng độ từ hóa và năng lượng dị hướng toàn phần của từng hạt, dẫn đến thay đổi trường bão hòa.
    • Nhiệt độ khóa ($T_B$) và đường cong ZFC: Mô phỏng đường cong ZFC chứng minh sự hiện diện của tương tác lưỡng cực làm dịch chuyển đỉnh $T_B$ về phía nhiệt độ cao hơn và mở rộng độ rộng phân bố năng lượng chắn. Snapshot cấu hình spin tại $0\text{ K}$ và tại $T_B$ thể hiện trực quan trạng thái đóng băng spin và sự mất định hướng nhiệt khi vượt qua rào cản năng lượng. Khảo sát tại các góc phương vị $\phi$ khác nhau giải thích rõ sự bất đối xứng của rào cản năng lượng trên giản đồ năng lượng 2D và 3D.
    • Mô-men xoắn từ tính: Mô phỏng mô-men xoắn của hệ hạt $\text{Fe}_3\text{O}_4$ ở các dải nhiệt độ: vùng nhiệt độ rất thấp, vùng dưới lân cận $300\text{ K}$ và vùng trên lân cận $300\text{ K}$. Kết quả chỉ ra sự phụ thuộc phi tuyến của biên độ mô-men xoắn vào nhiệt độ và góc từ hóa, cung cấp tham số cho việc điều khiển hạt trong hệ vi lưu.
Thông số / Đặc trưng khảo sát Phương pháp xác định Kết quả / Xu hướng ghi nhận trong văn bản
Pha tinh thể Phép đo XRD Pha spinel nghịch của $\text{Fe}_3\text{O}_4$, không lẫn tạp oxit sắt khác
Nhóm chức liên kết Đo phổ FTIR Xuất hiện liên kết amine/peptide, liên kết $\text{Fe}-\text{O}$ ($\sim 570\text{ cm}^{-1}$)
Hằng số dị hướng $K_{eff}$ Định luật LAS ($1/H^2$) Trích xuất từ độ dốc đường thẳng khớp vùng trường cao $M(H)$
Phân bố kích thước hạt Khớp hàm Langevin (Matlab) Tuân theo hàm phân bố logarit-chuẩn
Tương tác lưỡng cực Mô phỏng VAMPIRE (ASM) Tăng lực kháng từ tại $0\text{ K}$, làm dịch chuyển đỉnh $T_B$ của ZFC lên nhiệt độ cao
Khoảng cách liên hạt Mô phỏng MC / ASM Tương tác lưỡng cực suy giảm mạnh khi khoảng cách tăng từ $0{,}1\text{ nm}$ đến $10\text{ nm}$
Hấp phụ kim loại nặng Đo quang phổ AAS Tách loại hiệu quả ion $\text{Pb}^{2+}$ theo đường hấp phụ đẳng nhiệt

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

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

  1. Xây dựng quy trình tích hợp ba nhánh nghiên cứu: Thiết lập thành công mô hình nghiên cứu kết hợp đa chiều giữa thực nghiệm chế tạo, tính toán lý thuyết từ số và mô phỏng máy tính đa quy mô. Mô hình này khắc phục tính rời rạc của các phép đo vĩ mô truyền thống và hạn chế của các mô phỏng lý tưởng hóa.
  2. Sử dụng tham số thực nghiệm làm đầu vào mô phỏng: Ứng dụng thành công định luật tiếp cận bão hòa (LAS) và hàm Langevin để trích xuất trực tiếp hằng số dị hướng hiệu dụng $K_{eff}$ cùng hàm phân bố kích thước hạt thực tế từ mẫu chế tạo, đưa trực tiếp vào làm thông số đầu vào cho mô phỏng Monte Carlo và mô hình spin nguyên tử.
  3. Làm rõ bản chất vi từ của tương tác lưỡng cực - lưỡng cực: Định lượng chi tiết vai trò của tương tác lưỡng cực giữa các hạt nano $\text{Fe}_3\text{O}_4$ định hình trong nền polymer PGMA, chứng minh tương tác này làm biến đổi đường cong từ hóa, thay đổi lực kháng từ ở nhiệt độ thấp, nâng cao nhiệt độ khóa $T_B$ và làm biến dạng giản đồ năng lượng vi từ.
  4. Định lượng động học spin và mô-men xoắn: Tính toán chính xác hành vi mô-men xoắn và chuyển động hồi phục spin theo thời gian thực (thông qua phương trình LLG-Heun và Monte Carlo) ở các vùng nhiệt độ khác nhau xung quanh nhiệt độ phòng ($300\text{ K}$).

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

  1. Quy trình tổng hợp ổn định: Làm chủ quy trình tổng hợp vi cầu $\text{Fe}_3\text{O}_4\text{/PGMA}$ bằng phương pháp kết tủa tại chỗ trên nền PGMA amine hóa, giúp hạt nano từ phân tán đồng đều, hạn chế kết tụ không kiểm soát.
  2. Tiềm năng ứng dụng môi trường: Chứng minh thực nghiệm khả năng xử lý nước ô nhiễm thông qua phép đo hấp phụ đẳng nhiệt ion $\text{Pb}^{2+}$, cho phép tách và thu hồi chất hấp phụ nhanh chóng bằng từ trường ngoài.
  3. Cung cấp dữ liệu tiên nghiệm cho y sinh: Cung cấp các giá trị tính toán định lượng về mô-men từ, độ từ hóa bão hòa và mô-men xoắn định hướng, đóng vai trò số liệu tiền trạm cho việc thiết kế thiết bị vi lưu (lab-on-chip) và tối ưu hóa hệ phân phối thuốc hướng đích.

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

Hạn chế được ghi nhận

  • Đề tài tập trung chuyên sâu vào một hệ vật liệu đại diện cụ thể là $\text{Fe}_3\text{O}_4\text{/PGMA}$ nhằm đơn giản hóa và cụ thể hóa mô hình tính toán, do đó chưa khảo sát thực nghiệm mở rộng trên toàn bộ các cấu trúc mang khác (như liposome, micelles, dendrimer hay graphene).
  • Quá trình mô phỏng mô hình spin nguyên tử đòi hỏi tài nguyên tính toán rất lớn trên hệ HPC, do đó số lượng nguyên tử và kích thước cụm hạt mô phỏng vẫn bị giới hạn trong một thể tích vi mô nhất định, chưa thể mô phỏng đồng thời toàn bộ thể tích thực của một vi cầu polymer kích thước micromet chứa hàng triệu hạt nano từ.
  • Phép đo hấp phụ ion kim loại nặng mới dừng lại ở việc đánh giá dung lượng hấp phụ đối với ion chì ($\text{Pb}^{2+}$) trong điều kiện phòng thí nghiệm.

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

  • Mở rộng mô hình tính toán và mô phỏng số đã xây dựng cho các hệ vật liệu nanocomposite từ tính khác như vi nang, micelles, liposome, dendrimer và vật liệu nền carbon (CNT, graphene khử oxit - rGO).
  • Phát triển các thuật toán song song hóa nâng cao để tăng quy mô tính toán của mô hình spin nguyên tử lên các cụm hạt lớn hơn.
  • Triển khai các thử nghiệm lâm sàng hoặc in-vitro sâu hơn đối với khả năng tải/nhả thuốc Doxorubicin và hiệu ứng tăng thân nhiệt cục bộ dựa trên các thông số vi từ đã mô phỏng.

Giá trị tham khảo

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

  • Nghiên cứu sinh và học viên cao học ngành Khoa học vật liệu, Vật lý chất rắn, Vật lý tính toán: Tham khảo phương pháp thiết lập mô hình spin nguyên tử trên phần mềm VAMPIRE, thuật toán mô phỏng Monte Carlo áp suất hằng số, và cách giải phương trình Landau-Lifshitz-Gilbert.
  • Các nhà nghiên cứu thực nghiệm vật liệu từ nano và y sinh: Tham khảo quy trình chiết xuất tham số vi từ nội tại ($K_{eff}$, phân bố kích thước) từ dữ liệu VSM bằng định luật tiếp cận bão hòa (LAS) và khớp hàm Langevin bằng code Matlab, tránh việc chỉ mô tả dữ liệu từ hóa một cách định tính.
  • Kỹ sư và nhà công nghệ trong lĩnh vực vi lưu và xử lý môi trường: Tham khảo dữ liệu mô phỏng mô-men xoắn từ tính phục vụ việc thiết kế nhíp từ trong chip vi lưu (lab-on-chip) và quy trình chế tạo vi cầu polymer hấp phụ ion kim loại nặng.

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

1. Tại sao luận án lựa chọn hệ vật liệu vi cầu $\text{Fe}_3\text{O}_4\text{/PGMA}$ làm đối tượng nghiên cứu?

Hệ vật liệu này kết hợp được ưu điểm của hai thành phần: hạt nano từ $\text{Fe}_3\text{O}_4$ có độ từ hóa bão hòa cao, tương thích sinh học tốt, dễ tổng hợp và có tính siêu thuận từ ở kích thước dưới $50\text{ nm}$; polymer nền PGMA là polymer nhiệt dẻo thông minh, nhạy pH, có khả năng tự phân hủy sinh học và chứa các vòng epoxy dễ biến tính amine hóa ($-\text{NH}_2$) để phân tán giữ hạt và liên kết phối tử sinh học. Hệ này đại diện cho cấu trúc nanocomposite điển hình, thuận lợi để lượng hóa mô hình mà không làm mất tính tổng quát.

2. Định luật tiếp cận bão hòa (LAS) được tác giả sử dụng nhằm mục đích gì?

Định luật LAS được sử dụng để phân tích dữ liệu đường cong từ hóa $M(H)$ vùng từ trường cao thu được từ phép đo VSM thực nghiệm. Bằng cách vẽ đồ thị phụ thuộc tuyến tính giữa $M/M_s$ và $1/H^2$, tác giả chiết xuất được giá trị số của hằng số dị hướng từ hiệu dụng $K_{eff}$, cung cấp tham số đầu vào cho mô phỏng máy tính.

3. Tương tác lưỡng cực - lưỡng cực ảnh hưởng như thế nào đến tính chất từ của hệ vi cầu trong kết quả mô phỏng?

Mô phỏng cho thấy tương tác lưỡng cực - lưỡng cực đóng vai trò quan trọng: tại $0\text{ K}$, nó làm tăng độ từ dư và lực kháng từ; tại $300\text{ K}$, nó cản trở sự định hướng đồng thời của các spin, làm chậm quá trình đạt bão hòa từ. Trên đường cong ZFC, tương tác lưỡng cực làm dịch chuyển nhiệt độ khóa $T_B$ về vùng nhiệt độ cao hơn và mở rộng rào cản năng lượng dị hướng. Khi khoảng cách giữa các hạt tăng từ $0{,}1\text{ nm}$ lên $10\text{ nm}$, tương tác này suy giảm nhanh theo hàm bậc ba nghịch đảo của khoảng cách ($1/r^3$).

4. Hệ thống tính toán nào được sử dụng để thực thi các mô phỏng trong luận án?

Tác giả sử dụng hệ tính toán hiệu năng cao (HPC) của Trung tâm Tin học và Tính toán thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam (VAST) thông qua các tập lệnh quản lý tải và điều phối tác vụ (Job Scheduler), kết hợp với hệ điện toán lưới HCMIP để xử lý các mô phỏng song song quy mô lớn của gói phần mềm VAMPIRE và các chương trình tính toán số trên nền tảng Matlab.

Kết luận

Luận án của nghiên cứu sinh Nguyễn Thanh Hoàng đã thiết lập một phương pháp luận nghiên cứu đa chiều hoàn chỉnh kết hợp giữa thực nghiệm hóa học, tính toán lý thuyết từ số và mô phỏng số đa quy mô cho hệ vi cầu từ tính $\text{Fe}_3\text{O}_4\text{/PGMA}$. Công trình đã làm rõ cấu trúc vi từ, lượng hóa hằng số dị hướng hiệu dụng, phân bố kích thước hạt và xác định vai trò của tương tác lưỡng cực - lưỡng cực chi phối các hành vi từ nhiệt, đường cong từ hóa và mô-men xoắn. Các kết quả đạt được cung cấp cơ sở dữ liệu khoa học tin cậy và các tham số tiên nghiệm có giá trị cho việc định hướng ứng dụng vật liệu tổ hợp từ tính trong y sinh và kỹ thuật xử lý môi trường.