Tổng quan luận án

Hiệu ứng đốt nóng cảm ứng từ (Magnetic Inductive Heating – MIH) của chất lỏng hạt nano từ đang thu hút sự quan tâm lớn trong các ứng dụng y sinh như nhiệt từ trị ung thư, giải phóng thuốc có kiểm soát bằng kích nhiệt từ và rã đông sinh học. Trong các ứng dụng này, yêu cầu bắt buộc là lượng hạt nano từ đưa vào cơ thể phải ở mức tối thiểu nhằm đảm bảo tính tương thích sinh học, nhưng đồng thời nhiệt lượng sinh ra dưới tác dụng của từ trường xoay chiều (Alternative Magnetic Field – AMF) phải đủ lớn để đạt hiệu quả điều trị. Để tối ưu hóa công suất tỏa nhiệt, nghiên cứu thực nghiệm gặp nhiều trở ngại do công suất này phụ thuộc đồng thời vào rất nhiều tham số vật lý: kích thước hạt ($D$), độ rộng phân bố kích thước ($\sigma$), từ độ bão hòa ($M_S$), dị hướng từ ($K$), độ nhớt môi trường ($\eta$), cũng như tần số ($f$) và cường độ ($H$) của từ trường ngoài. Hơn nữa, sự đóng góp và cạnh tranh giữa hai cơ chế hồi phục Néel và Brown chưa được phân tích một cách chi tiết và hệ thống trên các hệ vật liệu khác nhau. Luận án lựa chọn giải quyết khoảng trống này bằng cách kết hợp tính toán mô phỏng số và thực nghiệm kiểm chứng.

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

  1. Nghiên cứu các đặc trưng đốt nóng của chất lỏng hạt nano từ và các yếu tố ảnh hưởng dựa trên kỹ thuật tính toán bằng số theo mô hình lý thuyết đáp ứng tuyến tính (LRT).
  2. Thực hiện nghiên cứu thực nghiệm trên hai hệ mẫu hạt nano ferit spinel $\text{CoFe}_2\text{O}_4$ và $\text{MnFe}_2\text{O}_4$ để đánh giá ảnh hưởng của các yếu tố vật lý lên đặc trưng hiệu ứng MIH, đối chứng và kiểm chứng các kết quả tính toán lý thuyết.

Hệ thống câu hỏi nghiên cứu mà luận án tập trung giải quyết bao gồm:

  1. Kích thước hạt tới hạn nào tối ưu hóa hiệu ứng MIH đối với từng loại vật liệu cụ thể?
  2. Từ độ bão hòa, đường kính động học và đặc biệt là dị hướng từ ảnh hưởng như thế nào đến các tham số đặc trưng của hiệu ứng MIH?
  3. Có thể phân loại vật liệu dựa trên tham số dị hướng từ để xác định cơ chế tổn hao sinh nhiệt chiếm ưu thế hay không?
  4. Sự thay đổi của độ phân tán kích thước hạt và độ nhớt môi trường tác động ra sao đến công suất tổn hao riêng?

Đối tượng nghiên cứu của luận án là hệ chất lỏng hạt nano từ, bao gồm các hạt nano từ dạng kim loại, hợp kim và oxit kim loại ($\text{FeCo}$, $\text{La}{0.7}\text{Sr}{0.3}\text{MnO}_3$ – $\text{LSMO}$, $\text{MnFe}_2\text{O}_4$, $\text{Fe}_3\text{O}_4$, $\text{CoFe}_2\text{O}_4$, $\text{FePt}$) được phân tán trong dung môi lỏng. Phạm vi không gian và thực nghiệm được tiến hành tại Phòng Vật lý vật liệu từ - siêu dẫn và Phòng Vật liệu Nano Y sinh thuộc Viện Khoa học Vật liệu – Viện Hàn lâm Khoa học và Công nghệ Việt Nam.

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

Hiệu ứng đốt nóng cảm ứng từ trong y sinh được mở đầu từ công trình của Gilchrist và cộng sự (1957), sử dụng hạt nano $\gamma\text{-Fe}_2\text{O}_3$ kích thước 20–100 nm dưới từ trường xoay chiều 1,2 MHz để nâng nhiệt độ lên $42\text{ }^\circ\text{C} \div 46\text{ }^\circ\text{C}$ nhằm tiêu diệt tế bào ung thư. Kể từ đó, nhiều hướng nghiên cứu đã mở rộng ứng dụng của hiệu ứng MIH:

  • Nghiên cứu giải hấp và tái tạo vật liệu hấp phụ thu hồi hợp chất hữu cơ dễ bay hơi (VOC) thông qua khả năng tự khống chế nhiệt độ Curie ($T_C$) của hạt nano ferit do N. Kikukawa và cộng sự đề xuất.
  • Kỹ thuật gia công tế bào điều chỉnh nồng độ insulin thông qua cổng nhiệt tế bào TRPV1 gắn hạt nano oxit sắt dưới từ trường tần số thấp do S. Stanley và cộng sự công bố.
  • Phương pháp đánh giá nồng độ hạt nano từ tích tụ tại các cơ quan nội tạng in vitro dựa trên nhiệt độ gia tăng của hiệu ứng MIH do P. Nam và cộng sự thực hiện.
  • Kỹ thuật rã đông sinh học cho các mô mỏng và tế bào đông lạnh bằng hạt nano từ của Etheridge và cộng sự.
  • Các cơ chế phân phối và nhả thuốc kích nhiệt từ qua phá vỡ liên kết polymer nhạy nhiệt (DBB) hoặc qua độ thẩm từ tăng cường (DEP).

Về mặt lý thuyết cơ bản, Weiss đề xuất giả thuyết đô men từ, tiếp theo là mô hình hồi phục siêu thuận từ của Néel và Brown. Về lý thuyết tổn hao sinh nhiệt, các mô hình chính bao gồm: mô hình tổn hao từ trễ Rayleigh cho từ trường thấp, mô hình Stoner – Wohlfarth (SWMBTs) cho hạt đơn đô men dị hướng đơn trục trong từ trường lớn ($H > 2H_C$), và mô hình lý thuyết đáp ứng tuyến tính (LRT) do Rosensweig và Carrey phát triển cho vùng siêu thuận từ hoặc từ trường nhỏ. Các nghiên cứu thực nghiệm của Hergt và cộng sự trên $\text{Fe}_3\text{O}_4$ đã xác định sự chuyển đổi cơ chế tổn hao từ trễ sang tổn hao hồi phục Néel khi kích thước hạt giảm xuống dưới 18 nm. Đồng thời, Brezovich và cộng sự đã xác lập ngưỡng giới hạn an toàn đối với cơ thể người khi chịu tác dụng của từ trường xoay chiều ($H \cdot f \le 4,85 \times 10^8 \text{ A}/(\text{m}\cdot\text{s})$ đến $5 \times 10^9 \text{ A}/(\text{m}\cdot\text{s})$).

Tại Việt Nam, nhóm nghiên cứu của GS. Nguyễn Xuân Phúc tại Viện Khoa học Vật liệu đã chế tạo các hệ hạt nano kim loại $\text{Fe}$, oxit $\text{Fe}3\text{O}4$, manganite pha tạp $\text{La}{0.7}\text{Sr}{0.3}\text{MnO}_3$ và các hệ cấu trúc lõi - vỏ polymer $\text{Fe}_3\text{O}_4\text{@P(St-co-AA)}$, $\text{Fe}_3\text{O}_4\text{@P(NIPAM-co-AA)}$.

Khoảng trống nghiên cứu mà luận án lựa chọn giải quyết là: Các công trình trước đây chủ yếu khảo sát thực nghiệm đơn lẻ trên một vài hệ mẫu hoặc tính toán lý thuyết đơn giản, thiếu sự phân tích định lượng toàn diện về vai trò của từng tham số vật lý ($D, \sigma, M_S, K, \eta, H, f$). Luận án xác lập ranh giới phân tách cơ chế tổn hao hồi phục Néel và Brown dựa trên hằng số dị hướng từ phân tách ($K_C$), đồng thời chế tạo thực nghiệm hai hệ mẫu có dị hướng từ tương phản ($\text{MnFe}_2\text{O}_4$ và $\text{CoFe}_2\text{O}_4$) để kiểm chứng trực tiếp các quy luật lý thuyết.

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

Cơ sở lý thuyết và khung phân tích

Luận án sử dụng mô hình Lý thuyết đáp ứng tuyến tính (Linear Response Theory – LRT). Điều kiện áp dụng mô hình là khi tham số không thứ nguyên $\xi = \frac{\mu_0 M_S V H}{k_B T} < 1$ hoặc khi từ trường ngoài thỏa mãn $H \ll H_C$, trong đó $\mu_0$ là độ từ thẩm chân không, $k_B$ là hằng số Boltzmann, $V$ là thể tích hạt và $T$ là nhiệt độ tuyệt đối.

Các biểu thức toán học nền tảng trong mô hình LRT bao gồm:

  • Thời gian hồi phục Néel: $$\tau_N = \frac{\sqrt{\pi}}{2} \tau_0 \frac{\exp\left(\frac{K V}{k_B T}\right)}{\sqrt{\frac{K V}{k_B T}}}$$ (với $\tau_0 \approx 10^{-9} \div 10^{-13}\text{ s}$).
  • Thời gian hồi phục Brown: $$\tau_B = \frac{3 \eta V_H}{k_B T} = \frac{\pi \eta D_H^3}{2 k_B T}$$ (với $D_H = D + 2\delta$ là đường kính động học, $\delta$ là bề dày lớp bọc bề mặt, $\eta$ là độ nhớt dung môi).
  • Thời gian hồi phục hiệu dụng $\tau$: $$\frac{1}{\tau} = \frac{1}{\tau_N} + \frac{1}{\tau_B} \implies \tau = \frac{\tau_N \tau_B}{\tau_N + \tau_B}$$
  • Công suất đốt nóng cảm ứng từ trung bình theo mô hình LRT: $$P_{\text{LRT}} = \mu_0 \pi \chi''(f) H^2 f$$ với phần ảo của độ cảm từ xoay chiều $\chi''(f)$: $$\chi''(f) = \chi_0 \frac{2\pi f \tau}{1 + (2\pi f \tau)^2}$$
  • Công suất tổn hao riêng ($\text{SLP}$) và Tốc độ hấp thụ riêng ($\text{SAR}$): $$\text{SLP} = \frac{P_{\text{LRT}}}{\rho}$$ trong đó $\rho$ là khối lượng riêng của vật liệu từ. Hiệu ứng đạt cực đại ($\text{SLP}_{\max}$) khi hệ thỏa mãn điều kiện cộng hưởng: $\omega\tau = 2\pi f \tau = 1$.
  • Năng lượng dị hướng từ hiệu dụng có tính đến đóng góp của bề mặt: $$K_{\text{eff}} = K_V + \frac{6 K_S}{D}$$

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

Luận án kết hợp giữa phương pháp mô phỏng số và thực nghiệm đối chứng:

  • Phương pháp tính toán số: Xây dựng các thuật toán mô phỏng trên phần mềm MATLAB để tính toán $P_{\text{LRT}}$, $\text{SLP}$, khảo sát sự phụ thuộc của $\text{SLP}$ vào các biến số vật lý và tìm các giá trị tối ưu ($D_{cp}, \Delta D_{cp}, \text{SLP}_{\max}, K_C$).
  • Phương pháp thực nghiệm:
    • Tổng hợp hạt nano từ: Chế tạo hai hệ hạt nano ferit spinel $\text{MnFe}_2\text{O}_4$ và $\text{CoFe}_2\text{O}_4$ bằng phương pháp thủy nhiệt ở các nhiệt độ phản ứng khác nhau ($100\text{ }^\circ\text{C}, 120\text{ }^\circ\text{C}, 140\text{ }^\circ\text{C}, 160\text{ }^\circ\text{C}, 180\text{ }^\circ\text{C}$).
    • Phân tích cấu trúc và hình thái: Nhiễu xạ tia X (XRD) trên hệ máy Siemens D5000 để xác định cấu trúc pha, kích thước tinh thể ($D_{\text{XRD}}$) và hằng số mạng ($a_{\exp}$); kính hiển vi điện tử quét phát xạ trường (FESEM) để quan sát hình thái và kích thước hạt.
    • Phân tích kích thước động học và độ nhớt: Đo đường kính động học ($D_H$) bằng kỹ thuật tán xạ ánh sáng động (DLS); đo độ nhớt chất lỏng từ ($\eta$) bằng máy đo độ nhớt rung sóng sin Sine wave Vibro Viscometer SV 10.
    • Phân tích tính chất từ: Sử dụng từ kế mẫu rung (VSM) tại nhiệt độ phòng để đo đường cong từ trễ $M(H)$, xác định từ độ bão hòa ($M_S$), từ dư ($M_r$), lực kháng từ ($H_C$) và dị hướng từ hiệu dụng ($K_{\text{eff}}$).
    • Khảo sát đặc trưng đốt nóng cảm ứng từ: Sử dụng hệ máy phát xoay chiều Model RDO-HFI công suất 5 kW để đo đường cong nhiệt độ theo thời gian $T(t)$, từ đó xác định tốc độ gia nhiệt ban đầu $\left(\frac{\Delta T}{\Delta t}\right)$ và tính toán giá trị $\text{SAR}$ thực nghiệm: $$\text{SAR} = \frac{C \cdot m_s}{m_i} \left(\frac{\Delta T}{\Delta t}\right)_{t \to 0}$$ trong đó $C$ là nhiệt dung riêng của chất lỏng từ, $m_s$ là khối lượng dung dịch và $m_i$ là khối lượng hạt nano từ trong mẫu.

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

Chương 1. Hiệu ứng đốt nóng cảm ứng từ của chất lỏng hạt nano từ

Chương 1 hệ thống hóa cơ sở lý thuyết về từ học hạt nano và tổng quan các kết quả nghiên cứu về hiệu ứng MIH:

  • Trình bày sự hình thành cấu trúc đơn đô men khi kích thước hạt nhỏ hơn bán kính tới hạn $r_C$, trạng thái siêu thuận từ và sự thăng giáng nhiệt của mô men từ theo định luật Néel – Brown. Phân tích sự phụ thuộc của nhiệt độ khóa $T_B$ và kích thước siêu thuận từ tới hạn $D_{\text{SPM}}$ vào thời gian đặc trưng của phép đo ($\tau_m$).
  • Phân tích nguồn gốc dị hướng từ trong vật liệu nano, chỉ rõ dị hướng từ hiệu dụng $K_{\text{eff}}$ nhận đóng góp quan trọng từ dị hướng bề mặt ($K_S$) khi kích thước hạt giảm, do sự phá vỡ tính đối xứng tinh thể và bất trật tự spin ở lớp biên.
  • Trình bày cấu trúc ba thành phần của chất lỏng từ (lõi từ vô cơ, lớp vỏ hoạt hóa bề mặt và dung môi phân tán) cùng các phương pháp tổng hợp phổ biến (đồng kết tủa, sol-gel, thủy nhiệt, phân hủy nhiệt).
  • Tổng thuật các cơ chế sinh nhiệt trong từ trường xoay chiều: tổn hao dòng điện xoáy Joule (không đáng kể ở thang nano), tổn hao từ trễ (mô hình Rayleigh và SWMBTs), và tổn hao hồi phục (Néel và Brown theo mô hình LRT).
  • Nêu rõ các thách thức kỹ thuật trong thực nghiệm tối ưu hóa hiệu ứng MIH: khống chế độ hẹp phân bố kích thước hạt, hiện tượng suy giảm từ độ do lớp chết từ bề mặt, kiểm soát độ dị hướng từ và ảnh hưởng của độ nhớt môi trường trong điều kiện giới hạn an toàn Brezovich.

Chương 2. Kết quả tính toán công suất đốt nóng cảm ứng từ theo mô hình đáp ứng tuyến tính

Chương 2 trình bày các kết quả mô phỏng số chi tiết bằng MATLAB dựa trên mô hình LRT cho sáu hệ hạt nano: $\text{FeCo}$, $\text{LSMO}$, $\text{MnFe}_2\text{O}_4$, $\text{Fe}_3\text{O}_4$, $\text{CoFe}_2\text{O}_4$ và $\text{FePt}$. Các phát hiện chính gồm:

  • Cạnh tranh giữa cơ chế hồi phục Néel và Brown: Với hạt kích thước nhỏ, thời gian hồi phục Néel $\tau_N$ ngắn hơn $\tau_B$ nên cơ chế Néel chiếm ưu thế. Khi kích thước hạt tăng, $\tau_N$ tăng theo hàm mũ trong khi $\tau_B$ tăng theo hàm bậc ba thể tích; do đó có sự chuyển dịch sang cơ chế Brown.
  • Đặc trưng đường $\text{SLP}(D)$: Đồ thị $\text{SLP}$ theo kích thước hạt $D$ xuất hiện đỉnh cực đại tại kích thước tới hạn $D_{cp}$ với độ rộng bán vạch $\Delta D_{cp}$. Vị trí đỉnh $D_{cp}$ tương ứng với điều kiện cộng hưởng $\omega\tau = 1$.
  • Ảnh hưởng của các tham số vật lý:
    • Từ trường ($H$): $\text{SLP}$ tỉ lệ thuận với $H^2$; vị trí $D_{cp}$ và độ rộng $\Delta D_{cp}$ không đổi theo $H$.
    • Tần số ($f$): Khi tần số tăng, $\text{SLP}{\max}$ tăng phi tuyến, đồng thời vị trí $D{cp}$ dịch chuyển về phía kích thước nhỏ hơn do thời gian hồi phục cần thiết để cộng hưởng giảm xuống.
    • Từ độ bão hòa ($M_S$): $\text{SLP}{\max}$ tăng tuyến tính với $M_S$. Độ dốc $\frac{\Delta\text{SLP}{\max}}{\Delta M_S}$ phản ánh mức độ nhạy nhiệt của từng loại vật liệu.
    • Độ nhớt dung môi ($\eta$): Ảnh hưởng phụ thuộc vào cơ chế tổn hao. Ở các vật liệu có cơ chế Néel chiếm ưu thế, sự thay đổi độ nhớt ít tác động lên $\text{SLP}$; ngược lại, ở các vật liệu có cơ chế Brown chiếm ưu thế, tăng độ nhớt làm giảm mạnh $\text{SLP}{\max}$ và dịch chuyển $D{cp}$.
    • Phân bố kích thước ($\sigma$): Khi độ lệch chuẩn $\sigma$ của hàm phân bố log-normal tăng lên, giá trị $\text{SLP}_{\max}$ giảm nhanh chóng, đỉnh phổ bị mở rộng và dịch chuyển về phía kích thước hạt nhỏ hơn.
  • Khái niệm dị hướng từ phân tách $K_C$ và phân nhóm vật liệu: Luận án phát hiện sự tồn tại của giá trị dị hướng từ giới hạn $K_C$. Dựa vào $K_C$, các vật liệu nano từ được phân làm hai nhóm rõ rệt:
    • Nhóm A ($K < K_C$, dị hướng từ thấp: $\text{MnFe}_2\text{O}_4$, $\text{Fe}_3\text{O}_4$, $\text{FeCo}$, $\text{LSMO}$): Cơ chế hồi phục Néel chiếm ưu thế tại đỉnh cộng hưởng; $D_{cp}$ phụ thuộc rất mạnh vào hằng số dị hướng từ $K$ nhưng hầu như không phụ thuộc vào độ nhớt $\eta$.
    • Nhóm B ($K > K_C$, dị hướng từ cao: $\text{CoFe}_2\text{O}_4$, $\text{FePt}$): Cơ chế hồi phục Brown chiếm ưu thế; $D_{cp}$ không thay đổi theo dị hướng từ $K$ nhưng phụ thuộc chặt chẽ vào độ nhớt $\eta$ và bề dày lớp bọc $\delta$.
Hệ vật liệu Hằng số dị hướng $K$ ($\text{J/m}^3$) Từ độ bão hòa $M_S$ ($\text{emu/g}$) Phân nhóm theo $K_C$ Cơ chế tổn hao chiếm ưu thế tại $D_{cp}$
$\text{FeCo}$ $1,5 \times 10^4$ $200$ Nhóm A Hồi phục Néel
$\text{LSMO}$ $0,8 \times 10^4$ $45$ Nhóm A Hồi phục Néel
$\text{MnFe}_2\text{O}_4$ $0,5 \times 10^4$ $60$ Nhóm A Hồi phục Néel
$\text{Fe}_3\text{O}_4$ $2,0 \times 10^4$ $80$ Nhóm A Hồi phục Néel
$\text{CoFe}_2\text{O}_4$ $20,0 \times 10^4$ $70$ Nhóm B Hồi phục Brown
$\text{FePt}$ $20,0 \times 10^4$ $50$ Nhóm B Hồi phục Brown

Chương 3. Thực nghiệm kiểm chứng kết quả tính toán lý thuyết

Chương 3 trình bày các kết quả thực nghiệm chế tạo và đo đạc trên hai hệ mẫu đại diện cho hai nhóm dị hướng từ: $\text{MnFe}_2\text{O}_4$ (Nhóm A) và $\text{CoFe}_2\text{O}_4$ (Nhóm B).

  • Chế tạo mẫu: Tổng hợp các dãy mẫu $\text{MnFe}_2\text{O}_4$ (ký hiệu MFT100 đến MFT180) và $\text{CoFe}_2\text{O}_4$ (ký hiệu CFT100 đến CFT180) bằng phương pháp thủy nhiệt ở các nhiệt độ từ $100\text{ }^\circ\text{C}$ đến $180\text{ }^\circ\text{C}$.
  • Đặc trưng cấu trúc và từ tính:
    • Giản đồ XRD chứng minh các mẫu đều có cấu trúc spinel đơn pha. Kích thước hạt tinh thể $D_{\text{XRD}}$ tăng dần theo nhiệt độ thủy nhiệt (MFT: từ 7,5 nm đến 14,2 nm; CFT: từ 6,8 nm đến 12,5 nm).
    • Ảnh FESEM xác nhận các hạt có dạng hình cầu, độ đồng đều tương đối cao.
    • Phân tích DLS cho thấy đường kính động học $D_H$ lớn hơn kích thước tinh thể $D_{\text{XRD}}$ do sự xuất hiện của lớp bọc hoạt hóa bề mặt.
    • Đường cong từ hóa VSM tại nhiệt độ phòng xác định: Hệ $\text{MnFe}_2\text{O}4$ thể hiện tính chất siêu thuận từ rõ nét với lực kháng từ $H_C \approx 0$ và dị hướng từ hiệu dụng thấp ($K{\text{eff}} \approx 0,45 \times 10^4 \div 0,78 \times 10^4\text{ J/m}^3$). Hệ $\text{CoFe}_2\text{O}4$ có lực kháng từ $H_C$ lớn ($150 \div 850\text{ Oe}$) và dị hướng từ hiệu dụng cao ($K{\text{eff}} \approx 12 \times 10^4 \div 22 \times 10^4\text{ J/m}^3$).
  • Kết quả đo đốt nóng cảm ứng từ và kiểm chứng lý thuyết:
    • Khảo sát theo từ trường và tần số: Giá trị $\text{SAR}$ của cả hai hệ mẫu MFT100 và CFT100 đều tỉ lệ thuận với bình phương cường độ từ trường ($H^2$) và tăng theo tần số $f$, hoàn toàn phù hợp với biểu thức lý thuyết của mô hình LRT.
    • Kiểm chứng vai trò kích thước hạt: Đối với hệ $\text{MnFe}_2\text{O}_4$ (Nhóm A), $\text{SAR}$ phụ thuộc rất mạnh vào kích thước hạt lõi tinh thể và đạt cực đại tại vùng kích thước tối ưu quanh 10–12 nm. Với hệ $\text{CoFe}_2\text{O}_4$ (Nhóm B), $\text{SAR}$ phụ thuộc chủ yếu vào đường kính động học $D_H$.
    • Kiểm chứng ảnh hưởng của độ nhớt: Khi thay đổi độ nhớt của dung môi phân tán, giá trị $\text{SAR}$ của hệ $\text{MnFe}_2\text{O}_4$ hầu như không bị suy giảm (chứng minh cơ chế Néel đóng vai trò quyết định, không phụ thuộc vào chuyển động quay cơ học của hạt). Ngược lại, giá trị $\text{SAR}$ của hệ $\text{CoFe}_2\text{O}_4$ giảm mạnh khi độ nhớt môi trường tăng (chứng minh cơ chế Brown chiếm ưu thế, do chuyển động quay của khối hạt bị cản trở bởi ma sát nhớt).
Dãy mẫu Kích thước $D_{\text{XRD}}$ (nm) Hằng số mạng $a_{\exp}$ (Å) Từ độ $M_S$ (emu/g) Lực kháng từ $H_C$ (Oe) Dị hướng $K_{\text{eff}}$ ($10^4\text{ J/m}^3$) Cơ chế sinh nhiệt thực nghiệm
MFT100–180 $7,5 \div 14,2$ $8,38 \div 8,42$ $42 \div 65$ $\approx 0$ $0,45 \div 0,78$ Tổn hao hồi phục Néel
CFT100–180 $6,8 \div 12,5$ $8,36 \div 8,39$ $38 \div 58$ $150 \div 850$ $12,0 \div 22,0$ Tổn hao hồi phục Brown

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

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

  1. Luận án đã làm sáng tỏ bản chất cơ chế cạnh tranh giữa hai quá trình tổn hao hồi phục Néel và Brown trong mô hình lý thuyết đáp ứng tuyến tính (LRT) đối với chất lỏng hạt nano từ.
  2. Đã xác lập và đưa ra khái niệm hằng số dị hướng từ phân tách $K_C$. Dựa vào $K_C$, lần đầu tiên các hệ vật liệu từ được phân loại một cách định lượng thành hai nhóm đặc trưng: Nhóm A ($K < K_C$, dị hướng từ thấp, cơ chế Néel chiếm ưu thế) và Nhóm B ($K > K_C$, dị hướng từ cao, cơ chế Brown chiếm ưu thế).
  3. Đã mô hình hóa và làm rõ sự phụ thuộc định lượng của kích thước hạt tới hạn ($D_{cp}$), độ rộng bán vạch ($\Delta D_{cp}$) và công suất tổn hao cực đại ($\text{SLP}_{\max}$) vào các tham số $H, f, M_S, K, \eta$ và độ rộng phân bố kích thước $\sigma$.

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

  1. Cung cấp quy trình công nghệ chế tạo hạt nano ferit spinel $\text{MnFe}_2\text{O}_4$ và $\text{CoFe}_2\text{O}_4$ bằng phương pháp thủy nhiệt có khả năng kiểm soát kích thước hạt và tính chất từ thông qua việc điều chỉnh nhiệt độ phản ứng.
  2. Cung cấp bộ số liệu thực nghiệm chuẩn xác về các tham số đốt nóng cảm ứng từ ($\text{SAR}, \text{SLP}$) của hai hệ vật liệu đại diện có dị hướng từ tương phản, qua đó kiểm chứng thành công các dự đoán của mô hình tính toán số.
  3. Đề xuất các giải pháp định hướng chế tạo vật liệu trong ứng dụng y sinh:
    • Đối với vật liệu nhóm A (như $\text{MnFe}_2\text{O}_4, \text{Fe}_3\text{O}4$): Cần tập trung khống chế chính xác kích thước hạt vô cơ trong vùng lân cận $D{cp}$ và thu hẹp phân bố kích thước hạt $\sigma$; không cần quá chú trọng đến độ nhớt của môi trường xung quanh.
    • Đối với vật liệu nhóm B (như $\text{CoFe}_2\text{O}_4$): Cần kiểm soát kích thước động học $D_H$, tối ưu hóa độ dày lớp vỏ bọc $\delta$ và lưu ý sự suy giảm công suất tỏa nhiệt khi độ nhớt môi trường sinh học tăng cao.

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

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

  • Mô hình tính toán lý thuyết LRT chỉ áp dụng chính xác trong vùng từ trường xoay chiều biên độ nhỏ thỏa mãn điều kiện $\xi < 1$ hoặc $H \ll H_C$, chưa tính đến đóng góp của tổn hao từ trễ khi từ trường lớn.
  • Các mô phỏng số chưa tính đến ảnh hưởng của tương tác lưỡng cực từ giữa các hạt nano khi mật độ nồng độ hạt trong dung dịch tăng cao.
  • Các phép đo kiểm chứng thực nghiệm hiệu ứng MIH mới được thực hiện trong môi trường dung dịch in vitro, chưa tiến hành đánh giá trên mô sống hoặc điều kiện in vivo phức tạp.

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

  • Mở rộng mô hình tính toán lý thuyết cho vùng từ trường lớn kết hợp mô hình Stoner – Wohlfarth và lý thuyết đáp ứng tuyến tính có xét đến năng lượng tương tác lưỡng cực giữa các hạt.
  • Thực hiện các thử nghiệm chức năng hóa bề mặt hạt nano ferit bằng các polymer tương thích sinh học gắn kháng thể hướng đích và khảo sát hiệu ứng MIH trong môi trường mô tế bào ung thư thực tế.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo hữu ích cho:

  • Nghiên cứu sinh và học viên cao học chuyên ngành Vật liệu điện tử, Khoa học Vật liệu, Vật lý chất rắn và Công nghệ Nano: Cung cấp phương pháp luận tính toán "thực nghiệm số" trên MATLAB và quy trình thực nghiệm thủy nhiệt chế tạo hạt nano ferit spinel.
  • Giảng viên và các nhà nghiên cứu trong lĩnh vực Vật lý vật liệu từ và Vật liệu Y sinh: Tham khảo Chương 2 về các bảng số liệu, đường cong phụ thuộc và tiêu chí phân nhóm vật liệu $K_C$ để định hướng thiết kế các cấu trúc hạt nano sinh nhiệt tối ưu.
  • Các nhóm nghiên cứu ứng dụng y dược học: Tham khảo số liệu về ngưỡng an toàn Brezovich ($H \cdot f$), mức độ ảnh hưởng của độ nhớt môi trường và kích thước động học nhằm tối ưu hóa liều lượng hạt nano đưa vào cơ thể trong nhiệt từ trị ung thư và dẫn thuốc hướng đích.

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

1. Tham số dị hướng từ phân tách $K_C$ trong luận án có ý nghĩa gì?

$K_C$ là giá trị ngưỡng dị hướng từ được xác lập từ tính toán mô hình LRT để phân định cơ chế tổn hao nào chiếm ưu thế tại đỉnh công suất cực đại $\text{SLP}{\max}$. Khi $K < K_C$ (Nhóm A), tổn hao hồi phục Néel chiếm ưu thế, kích thước tới hạn $D{cp}$ phụ thuộc chặt chẽ vào $K$ và không phụ thuộc vào độ nhớt $\eta$. Khi $K > K_C$ (Nhóm B), tổn hao hồi phục Brown chiếm ưu thế, $D_{cp}$ không đổi theo $K$ nhưng phụ thuộc mạnh vào $\eta$.

2. Hai hệ vật liệu $\text{MnFe}_2\text{O}_4$ và $\text{CoFe}_2\text{O}_4$ được chọn trong thực nghiệm đại diện cho cơ chế sinh nhiệt nào?

  • Hệ $\text{MnFe}_2\text{O}4$ có dị hướng từ hiệu dụng thấp ($K{\text{eff}} \approx 0,45 \times 10^4 \div 0,78 \times 10^4\text{ J/m}^3 < K_C$), thuộc Nhóm A, sinh nhiệt chủ yếu qua cơ chế tổn hao hồi phục Néel.
  • Hệ $\text{CoFe}_2\text{O}4$ có dị hướng từ hiệu dụng cao ($K{\text{eff}} \approx 12 \times 10^4 \div 22 \times 10^4\text{ J/m}^3 > K_C$), thuộc Nhóm B, sinh nhiệt chủ yếu qua cơ chế tổn hao hồi phục Brown.

3. Điều kiện để công suất đốt nóng cảm ứng từ đạt cực đại theo mô hình LRT là gì?

Hiệu ứng MIH đạt cực đại ($\text{SLP}{\max}$ hoặc $\text{SAR}{\max}$) khi hệ chất lỏng hạt nano từ thỏa mãn điều kiện cộng hưởng tần số: $\omega\tau = 2\pi f \tau = 1$, trong đó $\omega$ là tần số góc của từ trường xoay chiều và $\tau$ là thời gian hồi phục hiệu dụng của hạt nano từ.

4. Độ rộng phân bố kích thước hạt $\sigma$ ảnh hưởng như thế nào đến công suất đốt nóng?

Khi độ rộng phân bố kích thước $\sigma$ tăng (hạt kém đồng đều), giá trị công suất tổn hao riêng cực đại $\text{SLP}{\max}$ suy giảm rất nhanh, đồng thời đỉnh cực đại bị mở rộng và vị trí kích thước tới hạn $D{cp}$ dịch chuyển dần về phía các giá trị kích thước nhỏ hơn so với hệ đơn phân tán ($\sigma = 0$).

5. Giới hạn Brezovich đối với từ trường xoay chiều trong y sinh được xác định ra sao?

Giới hạn an toàn Brezovich quy định tích số giữa cường độ từ trường và tần số từ trường xoay chiều tác dụng lên cơ thể người không được vượt quá ngưỡng an toàn: $H \cdot f \le 4,85 \times 10^8\text{ A}/(\text{m}\cdot\text{s})$ (hoặc mức mở rộng $5 \times 10^9\text{ A}/(\text{m}\cdot\text{s})$), nhằm tránh gây ra hiện tượng sinh nhiệt không mong muốn do dòng điện xoáy cảm ứng trong các mô lành của cơ thể.

Kết luận

Luận án đã nghiên cứu một cách có hệ thống các đặc trưng đốt nóng cảm ứng từ của chất lỏng hạt nano từ thông qua việc kết hợp mô phỏng số theo mô hình đáp ứng tuyến tính và thực nghiệm kiểm chứng. Kết quả nghiên cứu đã làm sáng tỏ sự cạnh tranh giữa hai cơ chế tổn hao hồi phục Néel và Brown, đồng thời xác lập tiêu chí hằng số dị hướng từ phân tách $K_C$ để phân loại vật liệu thành hai nhóm sinh nhiệt đặc trưng. Các kết quả thực nghiệm trên hai hệ hạt nano ferit $\text{MnFe}_2\text{O}_4$ và $\text{CoFe}_2\text{O}_4$ chế tạo bằng phương pháp thủy nhiệt đã kiểm chứng đầy đủ các quy luật lý thuyết, cung cấp cơ sở khoa học và định hướng công nghệ quan trọng cho việc thiết kế, tối ưu hóa vật liệu nano từ ứng dụng trong nhiệt từ trị y sinh.