Tổng quan luận án

Luận án tiến sĩ khoa học vật liệu với đề tài "Các đặc trưng đốt nóng cảm ứng của chất lỏng hạt nano từ và các yếu tố ảnh hưởng" do nghiên cứu sinh Lưu Hữu Nguyên thực hiện tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, dưới sự hướng dẫn khoa học của GS. Nguyễn Xuân Phúc và PGS. Phạm Thanh Phong (Hà Nội, 2019). Luận án thuộc chuyên ngành Vật liệu điện tử, mã số đào tạo 9.23.

+-------------------------------------------------------------------------+
|                          THÔNG TIN TỔNG QUAN                             |
+---------------------+---------------------------------------------------+
| Đề tài luận án      | Các đặc trưng đốt nóng cảm ứng của chất lỏng hạt  |
|                     | nano từ và các yếu tố ảnh hưởng                   |
| Tác giả             | Lưu Hữu Nguyên                                    |
| Người hướng dẫn     | GS. Nguyễn Xuân Phúc, PGS. Phạm Thanh Phong       |
| Chuyên ngành        | Vật liệu điện tử (Mã số: 9.23)                    |
| Cơ sở đào tạo       | Học viện Khoa học và Công nghệ                    |
| Cơ sở thực hiện     | Viện Khoa học Vật liệu (VKHVL - VHLKHCNVN)        |
| Năm hoàn thành      | 2019                                              |
+---------------------+---------------------------------------------------+

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

Trong những thập kỷ gần đây, hiệu ứng đốt nóng cảm ứng từ (Magnetic Inductive Heating - MIH) của các chất lỏng hạt nano từ đang thu hút sự quan tâm lớn nhờ các tiềm năng ứng dụng trong y sinh học (như nhiệt từ trị ung thư, nhả thuốc có kiểm soát thông qua kích nhiệt, rã đông tế bào và mô đông lạnh, điều chỉnh thụ thể ion tế bào) và trong xử lý môi trường (giải hấp tái tạo chất hấp phụ). Đối với các ứng dụng y sinh, yêu cầu thực tiễn đặt ra là liều lượng hạt nano từ đưa vào cơ thể phải ở mức tối thiểu nhưng vẫn đảm bảo khả năng sinh nhiệt đủ lớn tại vùng đích, đồng thời hệ chất lỏng từ phải giữ được độ ổn định phân tán lâu dài và có tính tương thích sinh học.

Để đáp ứng các yêu cầu trên, việc tối ưu hóa công suất tỏa nhiệt (đại diện bởi các đại lượng như công suất tổn hao riêng SLP hoặc tốc độ hấp thụ riêng SAR) là vấn đề cốt lõi. Tuy nhiên, công suất đốt nóng cảm ứng từ là một hàm phức tạp, phụ thuộc đồng thời vào rất nhiều tham số:

  • Tham số nội tại của hạt: kích thước hạt ($D$), độ rộng phân bố kích thước ($\sigma$), từ độ bão hòa ($M_S$), hằng số dị hướng từ hiệu dụng ($K, K_{eff}$), cấu trúc pha tinh thể và hiệu ứng bề mặt.
  • Tham số của môi trường phân tán: độ nhớt chất lỏng ($\eta$), bề dày lớp vỏ bọc hoạt hóa ($\delta$), đường kính thủy động học ($D_H$).
  • Tham số của từ trường ngoài xoay chiều (AMF): cường độ từ trường ($H$) và tần số ($f$), vốn bị giới hạn bởi ngưỡng an toàn sinh học Brezovich ($H \cdot f \le C$).

Do có quá nhiều biến số đan xen, việc khảo sát thực nghiệm đơn thuần gặp rất nhiều khó khăn, tốn kém và thiếu tính hệ thống. Khoảng trống nghiên cứu được tác giả xác định là: còn thiếu các nghiên cứu lý thuyết chi tiết, toàn diện đóng vai trò "thực nghiệm số" nhằm phân tích định lượng ảnh hưởng độc lập cũng như tương hỗ của từng tham số vật lý lên công suất đốt nóng cảm ứng từ; đồng thời cần có các thực nghiệm đối chứng có kiểm soát chặt chẽ trên các hệ hạt có dị hướng từ khác biệt rõ rệt để kiểm chứng tính xác thực của mô hình lý thuyết.

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

Luận án xác định hai mục tiêu chính:

  1. Mục tiêu 1 (Lý thuyết và mô phỏng số): Nghiên cứu các đặc trưng đốt nóng cảm ứng của chất lỏng hạt nano từ và khảo sát có hệ thống ảnh hưởng của các tham số vật lý ($D, \sigma, M_S, K, \eta, H, f$) 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 (Linear Response Theory - LRT).
  2. Mục tiêu 2 (Thực nghiệm kiểm chứng): Từ các định hướng và dự đoán của kết quả tính toán số, tiến hành chế tạo thực nghiệm hai hệ hạt nano từ ferit đại diện có dị hướng từ khác biệt là $CoFe_2O_4$ (dị hướng từ cao) và $MnFe_2O_4$ (dị hướng từ thấp), phân tán thành chất lỏng từ và đo đạc các đặc trưng nhiệt từ để so sánh, kiểm chứng mô hình lý thuyết.

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

  • Đối tượng nghiên cứu:
    • Về mặt lý thuyết: Sáu hệ chất lỏng hạt nano từ đại diện bao gồm $FeCo$, $LSMO$ ($La_{0.7}Sr_{0.3}MnO_3$), $MnFe_2O_4$, $Fe_3O_4$, $CoFe_2O_4$ và $FePt$.
    • Về mặt thực nghiệm: Hệ hạt nano spinel ferit $MnFe_2O_4$ và $CoFe_2O_4$ chế tạo bằng phương pháp thủy nhiệt ở các nhiệt độ phản ứng khác nhau (100, 120, 140, 160, 180 °C), cùng các hệ chất lỏng từ phân tán trong dung môi với độ nhớt thay đổi.
  • Phạm vi nghiên cứu: Khảo sát trong vùng từ trường xoay chiều biên độ thấp và trung bình thỏa mãn điều kiện áp dụng của mô hình LRT ($H \ll H_C$, hoặc tham số $\xi = \frac{\mu_0 M_S V H}{k_B T} < 1$) và phù hợp với giới hạn an toàn sinh học Brezovich ($f$ trong khoảng 100 kHz – 300 kHz, $H$ trong khoảng vài chục đến vài trăm Oe). Quá trình thực nghiệm được triển khai tại 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

Luận án đã hệ thống hóa các mốc phát triển và các hướng nghiên cứu lớn về hiệu ứng đốt nóng cảm ứng từ trong và ngoài nước:

Các công trình và tác giả tiền tiêu

  • Khởi nguồn ứng dụng nhiệt từ trị: Gilchrist và cộng sự (1957) là những người đầu tiên sử dụng các hạt nano $\gamma-Fe_2O_3$ kích thước 20–100 nm dưới dòng điện xoay chiều tần số 1,2 MHz để nâng nhiệt độ lên 42–46 °C nhằm tiêu diệt tế bào ung thư.
  • Các ứng dụng y sinh học hiện đại:
    • Nhả thuốc kích nhiệt: Nghiên cứu của Derfus, Hu, Ho (2006–2010) về hai cơ chế nhả thuốc: phá vỡ liên kết polymer (DBB) và tăng cường độ thấm qua khe hở ma trận polymer (DEP).
    • Gia công tế bào: S. Stanley và cộng sự (2012) gắn hạt nano oxit sắt vào cổng thụ thể nhạy nhiệt TRPV1 trên màng tế bào, kích hoạt dòng ion $Ca^{2+}$ dưới từ trường xoay chiều để điều khiển biểu hiện gen và giải phóng insulin.
    • Đánh giá nồng độ hạt tích tụ in vivo: P. Nam và cộng sự (2014) sử dụng tín hiệu gia nhiệt của hiệu ứng MIH để định lượng phân bố hạt nano trong các cơ quan nội tạng.
    • Rã đông sinh học: Etheridge và cộng sự (2014) sử dụng MIH để làm ấm nhanh, đồng đều các mô tế bào đông lạnh sâu, hạn chế tổn thương kết tinh lại của nước đá.
  • Ứng dụng công nghiệp và môi trường: N. Kikukawa và cộng sự (2004, 2009) nghiên cứu cơ chế đốt nóng tự khống chế nhiệt độ dựa trên nhiệt độ Curie ($T_C$) của hạt ferit để giải hấp thu hồi các hợp chất hữu cơ dễ bay hơi (VOC) trên chất hấp phụ than hoạt tính, zeolit, SBA.
  • Cơ sở lý thuyết và mô hình tính toán:
    • Lý thuyết siêu thuận từ và thời gian hồi phục: Néel (1949), Brown (1963).
    • Mô hình Stoner – Wohlfarth (SWMBTs) cho hạt đơn đô men sắt từ có tổn hao trễ (Stoner & Wohlfarth, 1948).
    • Lý thuyết đáp ứng tuyến tính (LRT): Rosensweig (2002), Hergt (1998, 2004, 2006), Carrey (2011). Hergt chỉ ra rằng với hạt nano $Fe_3O_4$ có đường kính $D < 18\text{ nm}$, tổn hao trễ biến mất và tổn hao hồi phục trở thành nguồn sinh nhiệt chủ đạo.
    • Giới hạn sinh học: Brezovich (1988) thiết lập ngưỡng chịu đựng an toàn của cơ thể người đối với tích số $H \cdot f \le 4,85 \times 10^8\text{ A/(m}\cdot\text{s)}$ (tương đương $6 \times 10^7\text{ Oe}\cdot\text{Hz}$).
+-----------------------------------------------------------------------------------+
|                 CÁC MÔ HÌNH LÝ THUYẾT TỔN HAO NHIỆT TỪ                           |
+-------------------+---------------------------+-----------------------------------+
| Mô hình           | Điều kiện áp dụng         | Cơ chế sinh nhiệt chính           |
+-------------------+---------------------------+-----------------------------------+
| Rayleigh          | Hạt đa đô men, từ trường  | Chuyển động vách đô men,          |
|                   | nhỏ ($H \ll H_C$)         | tổn hao tỉ lệ $H^3$               |
+-------------------+---------------------------+-----------------------------------+
| Stoner-Wohlfarth  | Hạt đơn đô men sắt từ,    | Đảo mô men từ bất thuận nghịch    |
| (SWMBTs)          | từ trường lớn ($H > 2H_C$)| qua chu trình trễ, tỉ lệ $2K$     |
+-------------------+---------------------------+-----------------------------------+
| Linear Response   | Hạt siêu thuận từ hoặc    | Tổn hao hồi phục Néel ($\tau_N$)  |
| Theory (LRT)      | $H \ll H_C$ ($\xi < 1$)   | và tổn hao hồi phục Brown ($\tau_B$)|
+-------------------+---------------------------+-----------------------------------+

Hướng nghiên cứu trong nước

Tại Việt Nam, các nghiên cứu về vật liệu nano từ được triển khai tại Viện Khoa học Vật liệu, Viện Kỹ thuật Nhiệt đới, Trường Đại học Bách khoa Hà Nội, Trường Đại học Khoa học Tự nhiên (ĐHQGHN). 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 đã tiên phong chế tạo hạt nano kim loại Fe, oxit $Fe_3O_4$, ferit pha tạp $Mn_{0.3}Zn_{0.7}Fe_2O_4$, manganite $La_{0.7}Sr_{0.3}MnO_3$, và các hệ cấu trúc lõi-vỏ polymer nhằm định hướng nhiệt từ trị.

Vị trí và đóng góp của luận án

Luận án của Lưu Hữu Nguyên giải quyết một mắt xích then chốt: kết nối giữa lý thuyết vi mô LRT và thực nghiệm chế tạo thông qua một chương trình "thực nghiệm số" quy mô trên MATLAB. Luận án phân định tường minh vai trò của dị hướng từ $K$ trong việc định hướng cơ chế hồi phục (Néel hay Brown), xác lập khái niệm dị hướng tới hạn $K_C$ để phân nhóm vật liệu, đồng thời chứng minh bằng thực nghiệm trên hai hệ mẫu đại diện ($MnFe_2O_4$ đại diện cho nhóm Néel và $CoFe_2O_4$ đại diện cho nhóm Brown).

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

Cơ sở lý thuyết

  1. Trạng thái siêu thuận từ và các thời gian hồi phục:

    • Khi kích thước hạt nhỏ hơn đường kính tới hạn $D_{SPM}$, năng lượng dị hướng $KV$ trở nên tương đương năng lượng nhiệt $k_B T$. Mô men từ dao động nhiệt giữa các hướng dễ từ hóa với thời gian hồi phục Néel ($\tau_N$): $$\tau_N = \frac{\sqrt{\pi}}{2} \tau_0 \frac{\exp\left(\frac{KV}{k_B T}\right)}{\sqrt{\frac{KV}{k_B T}}}$$ (trong đó $\tau_0 \approx 10^{-9} - 10^{-13}\text{ s}$, $V = \frac{\pi D^3}{6}$)
    • Khi hạt nano lơ lửng trong chất lỏng, toàn bộ hạt có thể quay cơ học theo sự định hướng của trường xoay chiều, đặc trưng bởi thời gian hồi phục Brown ($\tau_B$): $$\tau_B = \frac{3\eta V_H}{k_B T} = \frac{\pi \eta (D + 2\delta)^3}{2 k_B T}$$ (trong đó $\eta$ là độ nhớt động học, $V_H$ là thể tích thủy động học, $\delta$ là bề dày lớp chất hoạt hóa bề mặt)
    • Thời gian hồi phục hiệu dụng ($\tau$) của toàn hệ là sự kết hợp song song của hai quá trình: $$\frac{1}{\tau} = \frac{1}{\tau_N} + \frac{1}{\tau_B} \implies \tau = \frac{\tau_N \tau_B}{\tau_N + \tau_B}$$

  2. Công suất tổn hao theo Lý thuyết Đáp ứng Tuyến tính (LRT):

    • Công suất sinh nhiệt theo mô hình LRT ($P_{LRT}$) được tính thông qua phần ảo của độ cảm từ xoay chiều $\chi''(f)$: $$P_{LRT} = \mu_0 \pi \chi''(f) H^2 f = \mu_0 \pi \left[ \chi_0 \frac{2\pi f \tau}{1 + (2\pi f \tau)^2} \right] H^2 f$$
    • Trong đó độ cảm từ tĩnh $\chi_0$ phụ thuộc vào từ độ bão hòa đô men $M_D$ qua hàm Langevin: $$\chi_0 = \frac{M_S}{H} \left[ \coth\left(\frac{\mu_0 H M_D V}{k_B T}\right) - \frac{k_B T}{\mu_0 H M_D V} \right]$$
    • Điều kiện cộng hưởng để công suất đạt cực đại ($SLP_{max}$) xảy ra khi tần số góc của từ trường ngoài phù hợp với thời gian hồi phục: $$\omega \tau = 2\pi f \tau = 1$$
    • Khi tính đến sự phân bố kích thước không đồng đều của hạt thực nghiệm, hàm phân bố Log-normal $g(D, D_0, \sigma)$ được tích hợp để tính công suất tổn hao trung bình $\overline{P_{LRT}}$.

  3. Dị hướng từ hiệu dụng ($K_{eff}$):

    • Kết hợp đóng góp của dị hướng thể tích ($K_V$) và dị hướng bề mặt ($K_S$): $$K_{eff} = K_V + \frac{6}{D} K_S$$
+------------------------------------------------------------------------------------+
|               BẢNG TỔNG HỢP CÁC PHƯƠNG PHÁP THỰC NGHIỆM                            |
+--------------------+----------------------------------+----------------------------+
| Phương pháp / Đo   | Thiết bị sử dụng                 | Thông số thu nhận          |
+--------------------+----------------------------------+----------------------------+
| Tổng hợp mẫu       | Bình thủy nhiệt áp suất          | Hệ hạt $MnFe_2O_4$,        |
|                    | nhiệt độ 100 - 180 °C            | $CoFe_2O_4$                |
+--------------------+----------------------------------+----------------------------+
| Nhiễu xạ tia X     | Siemens D5000                    | Cấu trúc pha, hằng số mạng |
| (XRD)              |                                  | $a_{exp}$, kích thước tinh |
|                    |                                  | thể $D_{XRD}$              |
+--------------------+----------------------------------+----------------------------+
| Hiển vi điện tử    | Field Emission Scanning Electron | Hình thái học, kích thước, |
| phát xạ trường     | Microscope (FESEM)               | phân bố kích thước hạt     |
+--------------------+----------------------------------+----------------------------+
| Tán xạ ánh sáng    | Dynamic Light Scattering         | Đường kính thủy động học   |
| động (DLS)         | (DLS)                            | $D_H$                      |
+--------------------+----------------------------------+----------------------------+
| Đo độ nhớt         | Sine wave Vibro Viscometer SV-10 | Độ nhớt chất lỏng từ $\eta$|
+--------------------+----------------------------------+----------------------------+
| Từ kế mẫu rung     | Vibrating Sample Magnetometer    | Đường từ hóa $M(H), M_S,$  |
| (VSM)              | (VSM)                            | $H_C, M_r, K_{eff}, D_{SPM}$|
+--------------------+----------------------------------+----------------------------+
| Đốt nóng cảm ứng   | Máy phát RDO-HFI                 | Đường cong $T(t)$, tốc độ  |
| từ (MIH)           | (Công suất 5 kW)                 | gia nhiệt $\Delta T/\Delta t$,|
|                    |                                  | SAR, SLP                   |
+--------------------+----------------------------------+----------------------------+

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

Mở đầu

Tác giả trình bày tổng quan về tầm quan trọng của vật liệu nano và chất lỏng từ trong bối cảnh cuộc cách mạng khoa học công nghệ thế kỷ 21. Nêu rõ vấn đề tối ưu hóa công suất tỏa nhiệt cho ứng dụng nhiệt từ trị y sinh, các thách thức thực nghiệm do sự phụ thuộc đa tham số, từ đó xác lập mục tiêu, đối tượng, phương pháp và bố cục 3 chương của luận án.

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 là chương tổng quan chuyên sâu về các hiện tượng từ học kích thước nanomet và các cơ chế vật lý của hiệu ứng MIH:

  • Cấu trúc đô men và trạng thái siêu thuận từ: Giải thích sự chuyển đổi từ hạt đa đô men sang đơn đô men dưới bán kính tới hạn $r_C$ (theo phương trình Brown-Morrish). Phân tích hàng rào năng lượng dị hướng $E_a = KV$, phương trình Néel-Brown và nhiệt độ khóa $T_B$, kích thước giới hạn siêu thuận từ $D_{SPM}$.
  • Bản chất của dị hướng từ: Trình bày chi tiết về dị hướng từ tinh thể đối với cấu trúc lập phương (Fe, Ni) và đơn trục (Co); phân tích sự bất trật tự spin ở bề mặt dẫn tới dị hướng bề mặt $K_S$ và sự phụ thuộc kích thước của dị hướng hiệu dụng $K_{eff} = K_V + \frac{6}{D}K_S$.
  • Cấu trúc và quy trình chế tạo chất lỏng từ: Mô hình 3 thành phần gồm lõi nano từ (kim loại, hợp kim, oxit ferit), lớp vỏ bọc chức năng hóa/chống kết tụ (dextran, PEG, polyme, silica, citrate) với bề dày $\delta$, và môi trường phân tán (nước, dung môi hữu cơ).
  • Ứng dụng của hiệu ứng MIH: Khảo sát chi tiết 6 nhóm ứng dụng:
    1. Giải hấp và tái tạo chất hấp phụ VOC tự khống chế nhiệt độ qua điểm Curie $T_C$.
    2. Kích nhiệt cổng TRPV1 điều khiển dòng $Ca^{2+}$ giải phóng insulin.
    3. Đánh giá nồng độ tích tụ hạt nano trong cơ quan nội tạng.
    4. Rã đông mô và tế bào bảo quản lạnh sâu.
    5. Nhả thuốc có kiểm soát qua cơ chế phá vỡ liên kết (DBB) hoặc tăng tính thấm ma trận (DEP).
    6. Nhiệt từ trị ung thư (nâng nhiệt độ khối u lên 42–46 °C gây hoại tử chọn lọc).
  • Các cơ chế sinh nhiệt và mô hình lý thuyết:
    • Phân tích 4 thành phần tổn hao: tổn hao dòng điện Foucault bề mặt (không đáng kể ở hạt nano oxit), tổn hao từ trễ $P_{hys}$, tổn hao hồi phục Néel $P_N$ và tổn hao hồi phục Brown $P_B$.
    • Phân tích ranh giới của các mô hình lý thuyết: mô hình Rayleigh ($H \ll H_C$, hạt đa đô men), mô hình Stoner-Wohlfarth SWMBTs ($H > 2H_C$, hạt đơn đô men sắt từ), và mô hình LRT (hạt siêu thuận từ hoặc $H \ll H_C$).
    • Khái quát các thách thức trong thực nghiệm khống chế phân bố kích thước hạt, hiện tượng lớp chết từ tính ở bề mặt làm suy giảm $M_S$, và các yêu cầu khắt khe về độ nhớt môi trường sinh học.

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 là phần "thực nghiệm số" cốt lõi của luận án, sử dụng phần mềm MATLAB để mô phỏng định lượng phương trình LRT trên 6 hệ vật liệu: $FeCo, LSMO, MnFe_2O_4, Fe_3O_4, CoFe_2O_4, FePt$.

+------------------------------------------------------------------------------------------------+
|            ĐẶC TRƯNG HỒI PHỤC VÀ KÍCH THƯỚC TỚI HẠN Dcp TẠI TẦN SỐ f = 236 kHz                 |
+-------------------+-----------------+-----------------------+------------------+---------------+
| Hệ vật liệu       | Dị hướng K      | Cơ chế hồi phục       | Kích thước tới   | Bề rộng bán   |
|                   | ($J/m^3$)       | đóng góp chính        | hạn $D_{cp}$ (nm)| vạch $\Delta D_{cp}$ (nm)|
+-------------------+-----------------+-----------------------+------------------+---------------+
| $FeCo$            | $1,0 \times 10^4$ | Néel ($\tau_N < \tau_B$) | 17,2             | 3,8           |
| $LSMO$            | $1,5 \times 10^3$ | Néel ($\tau_N < \tau_B$) | 25,6             | 5,6           |
| $MnFe_2O_4$       | $3,0 \times 10^3$ | Néel ($\tau_N < \tau_B$) | 22,0             | 4,8           |
| $Fe_3O_4$         | $2,3 \times 10^4$ | Néel ($\tau_N < \tau_B$) | 14,8             | 3,2           |
| $CoFe_2O_4$       | $2,0 \times 10^5$ | Brown ($\tau_B < \tau_N$)| 7,5              | 12,5          |
| $FePt$            | $2,0 \times 10^5$ | Brown ($\tau_B < \tau_N$)| 7,8              | 13,0          |
+-------------------+-----------------+-----------------------+------------------+---------------+

Các phát hiện quy luật quan trọng thu được từ tính toán:

  1. Sự cạnh tranh giữa hai cơ chế Néel và Brown:
    • Khi kích thước hạt nhỏ hơn một giá trị phân tách $D_{NB}$, thời gian hồi phục Néel $\tau_N$ nhỏ hơn Brown $\tau_B$ nhiều bậc, cơ chế Néel hoàn toàn chi phối.
    • Khi $D > D_{NB}$, $\tau_N$ tăng vọt theo hàm mũ với thể tích $V$, trong khi $\tau_B$ chỉ tăng bậc nhất với $V_H$; do đó thời gian hồi phục Brown trở nên ngắn hơn và đóng vai trò quyết định.
  2. Đặc trưng của đỉnh cực đại $SLP(D)$:
    • Đường cong $SLP$ theo kích thước $D$ luôn xuất hiện đỉnh cực đại tại kích thước tới hạn $D_{cp}$ (tương ứng điều kiện cộng hưởng $2\pi f \tau = 1$).
    • Nhóm vật liệu có $K$ thấp ($FeCo, LSMO, MnFe_2O_4, Fe_3O_4$ - Nhóm A): Đỉnh $SLP$ rất hẹp ($\Delta D_{cp}$ nhỏ, từ 3,2 đến 5,6 nm), cơ chế tổn hao Néel chi phối hoàn toàn đỉnh phát nhiệt.
    • Nhóm vật liệu có $K$ cao ($CoFe_2O_4, FePt$ - Nhóm B): Đỉnh $SLP$ trải rộng trên dải kích thước lớn ($\Delta D_{cp}$ rộng, từ 12,5 đến 13,0 nm), cơ chế tổn hao Brown chi phối hoàn toàn.
  3. Ảnh hưởng của các tham số vật lý:
    • Cường độ từ trường ($H$): $SLP$ tỉ lệ thuận với $H^2$; vị trí $D_{cp}$ không đổi khi thay đổi $H$.
    • Tần số từ trường ($f$): $SLP_{max}$ tăng tuyến tính theo $f$; khi tần số tăng thì $D_{cp}$ dịch chuyển về phía kích thước nhỏ hơn.
    • Từ độ bão hòa ($M_S$): $SLP_{max}$ phụ thuộc tuyến tính với $M_S$; độ dốc $\frac{\Delta SLP_{max}}{\Delta M_S}$ phụ thuộc vào bản chất vật liệu.
    • Độ nhớt dung môi ($\eta$): Đối với nhóm Néel (Nhóm A), độ nhớt hầu như không làm thay đổi $SLP$; ngược lại, đối với nhóm Brown (Nhóm B), khi độ nhớt tăng, $D_{cp}$ dịch chuyển mạnh về phía hạt nhỏ hơn và giá trị $SLP_{max}$ suy giảm rõ rệt.
    • Độ rộng phân bố kích thước ($\sigma$): Khi hạt có độ phân tán kích thước càng rộng ($\sigma$ tăng), đỉnh $SLP_{max}$ bị suy giảm nhanh chóng theo hàm suy giảm hàm mũ/phi tuyến. Mức độ suy giảm ở nhóm A nghiêm trọng hơn nhiều so với nhóm B.
  4. Xác lập tham số dị hướng từ phân tách $K_C$:
    • Luận án đã phát hiện sự tồn tại của một giá trị dị hướng tới hạn $K_C$. Khi $K < K_C$, hệ thuộc vùng cơ chế Néel; khi $K > K_C$, hệ chuyển hoàn toàn sang vùng cơ chế Brown. Giá trị $K_C$ phụ thuộc trực tiếp vào tần số $f$ và độ nhớt $\eta$.

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

Để kiểm chứng các dự đoán lý thuyết ở Chương 2, tác giả đã tiến hành thực nghiệm tổng hợp và khảo sát hai hệ mẫu điển hình: $MnFe_2O_4$ (đại diện cho nhóm $K < K_C$, cơ chế Néel) và $CoFe_2O_4$ (đại diện cho nhóm $K > K_C$, cơ chế Brown).

+----------------------------------------------------------------------------------------------------+
|               BẢNG THÔNG SỐ CẤU TRÚC VÀ TỪ TÍNH CỦA CÁC HỆ MẪU THỰC NGHIỆM                         |
+----------+--------------------+----------------+----------------+----------------+-----------------+
| Ký hiệu  | Nhiệt độ thủy      | Kích thước     | Từ độ bão hòa  | Lực kháng từ   | Hằng số dị      |
| mẫu      | nhiệt (°C)         | XRD $D_{XRD}$  | $M_S$ (emu/g)  | $H_C$ (Oe)     | hướng $K_{eff}$ |
|          |                    | (nm)           |                |                | ($erg/cm^3$)    |
+----------+--------------------+----------------+----------------+----------------+-----------------+
| MFT100   | 100                | 7,5            | 38,2           | 0              | $3,5 \times 10^4$|
| MFT120   | 120                | 8,6            | 44,1           | 0              | $3,8 \times 10^4$|
| MFT140   | 140                | 10,2           | 49,8           | 0              | $4,1 \times 10^4$|
| MFT160   | 160                | 12,4           | 55,4           | 0              | $4,5 \times 10^4$|
| MFT180   | 180                | 15,1           | 62,0           | 0              | $5,0 \times 10^4$|
+----------+--------------------+----------------+----------------+----------------+-----------------+
| CFT100   | 100                | 5,2            | 35,6           | 45             | $1,8 \times 10^5$|
| CFT120   | 120                | 6,8            | 42,3           | 120            | $2,1 \times 10^5$|
| CFT140   | 140                | 8,4            | 48,7           | 260            | $2,5 \times 10^5$|
| CFT160   | 160                | 10,5           | 54,2           | 480            | $2,9 \times 10^5$|
| CFT180   | 180                | 13,8           | 59,5           | 750            | $3,4 \times 10^5$|
+----------+--------------------+----------------+----------------+----------------+-----------------+

Các kết quả thực nghiệm đạt được:

  1. Đặc trưng cấu trúc và hình thái: Giản đồ XRD xác nhận cấu trúc pha spinel đơn pha cho cả hai hệ mẫu; kích thước tinh thể $D_{XRD}$ tăng dần từ ~5 nm đến ~15 nm khi tăng nhiệt độ thủy nhiệt từ 100 °C lên 180 °C. Ảnh FESEM cho thấy các hạt có dạng hình cầu đồng đều. Đo tán xạ ánh sáng động DLS xác định đường kính thủy động học $D_H$ của các hạt sau khi phân tán.
  2. Đặc trưng từ tính:
    • Hệ $MnFe_2O_4$ thể hiện tính chất siêu thuận từ điển hình tại nhiệt độ phòng (lực kháng từ $H_C = 0$, từ dư $M_r = 0$). Đường từ hóa $M(H)$ khớp tốt với hàm Langevin.
    • Hệ $CoFe_2O_4$ có dị hướng từ lớn, thể hiện tính chất sắt từ yếu với lực kháng từ $H_C$ tăng dần theo kích thước hạt (từ 45 Oe đến 750 Oe). Dị hướng từ hiệu dụng $K_{eff}$ của $CoFe_2O_4$ lớn hơn $MnFe_2O_4$ xấp xỉ một bậc độ lớn.
  3. Đo đạc đặc trưng đốt nóng cảm ứng từ (MIH):
    • Khảo sát sự phụ thuộc vào từ trường: SAR tỉ lệ bậc hai với biên độ từ trường $H^2$ tại vùng trường thấp, hoàn toàn phù hợp với biểu thức lý thuyết LRT.
    • Khảo sát sự phụ thuộc vào tần số: SAR tăng tuyến tính với tần số kích thích $f$ (đo tại 3 tần số khác nhau dưới cùng biên độ 80 Oe).
    • Khảo sát sự phụ thuộc vào kích thước hạt: Khi kích thước hạt tăng, SAR chuẩn hóa theo từ độ ($SAR/M_S$) của cả hai hệ mẫu đều có xu hướng tăng và đạt giá trị tối ưu tương ứng với vùng kích thước dự đoán của lý thuyết.
    • Khảo sát ảnh hưởng của độ nhớt môi trường: Khi thay đổi độ nhớt $\eta$ của chất lỏng từ bằng cách pha thêm glycerol, kết quả thực nghiệm cho thấy SAR của hệ chất lỏng $MnFe_2O_4$ hầu như không bị ảnh hưởng (khẳng định cơ chế Néel là chủ đạo), trong khi SAR của hệ chất lỏng $CoFe_2O_4$ bị suy giảm mạnh khi độ nhớt tăng (khẳng định sự đóng góp chiếm ưu thế của cơ chế Brown).

Kết luận chung

Tóm lược toàn bộ các kết quả chính đã đạt được: xây dựng thành công bộ công cụ tính toán số khảo sát toàn diện các tham số của mô hình LRT; làm sáng tỏ quy luật cạnh tranh giữa hai cơ chế hồi phục Néel và Brown; xác lập tham số dị hướng tới hạn $K_C$; và kiểm chứng thực nghiệm thành công các quy luật này trên hai hệ hạt $MnFe_2O_4$ và $CoFe_2O_4$.

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

Luận án đã đạt được các đóng góp mới về mặt khoa học lý thuyết và thực tiễn ứng dụng như sau:

  1. Về mặt lý thuyết và mô phỏng số:

    • Đã khảo sát và làm sáng tỏ một cách toàn diện, có hệ thống ảnh hưởng độc lập cũng như tương hỗ của toàn bộ các tham số vật lý ($D, \sigma, M_S, K, \eta, H, f$) đến công suất đốt nóng cảm ứng từ theo mô hình lý thuyết LRT.
    • Phát hiện quy luật phân nhóm vật liệu dựa trên hằng số dị hướng từ: đưa ra khái niệm và phương pháp xác định giá trị dị hướng từ tới hạn $K_C$. Phân định rõ hai nhóm vật liệu: nhóm dị hướng thấp ($K < K_C$, tổn hao Néel chi phối, đỉnh $SLP$ nhọn, nhạy cảm với phân bố kích thước $\sigma$ nhưng không phụ thuộc độ nhớt $\eta$) và nhóm dị hướng cao ($K > K_C$, tổn hao Brown chi phối, đỉnh $SLP$ rộng, phụ thuộc mạnh vào độ nhớt $\eta$).
    • Thiết lập các tiêu chí tối ưu hóa tham số cho từng hệ vật liệu cụ thể, cung cấp bộ dữ liệu tham chiếu ("thực nghiệm số") chuẩn xác cho việc thiết kế chất lỏng từ.

  2. Về mặt thực nghiệm:

    • Chế tạo thành công bằng phương pháp thủy nhiệt hai hệ hạt nano từ ferit spinel $MnFe_2O_4$ và $CoFe_2O_4$ có kích thước kiểm soát được trong dải 5–15 nm, cấu trúc tinh thể hoàn hảo và từ tính xác định.
    • Chứng minh bằng thực nghiệm sự đúng đắn của mô hình lý thuyết: khẳng định quy luật phụ thuộc của SAR vào $H^2$, vào tần số $f$, và đặc biệt đã kiểm chứng trực tiếp vai trò của độ nhớt $\eta$ lên hai cơ chế hồi phục Néel (ở mẫu $MnFe_2O_4$) và Brown (ở mẫu $CoFe_2O_4$).

  3. Ý nghĩa thực tiễn:

    • Kết quả nghiên cứu đưa ra định hướng rõ ràng cho việc lựa chọn vật liệu trong nhiệt từ trị y sinh: đối với các mô có độ nhớt thay đổi hoặc môi trường nội bào phức tạp, nên ưu tiên chọn các hạt nano thuộc nhóm Néel (như $MnFe_2O_4$ hoặc $Fe_3O_4$ với kích thước tối ưu ~15–22 nm) để công suất nhiệt sinh ra không bị suy giảm do độ nhớt môi trường xung quanh.

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

Dựa trên phạm vi và nội dung trình bày trong luận án:

  • Hạn chế:
    • Mô hình tính toán tập trung chủ yếu vào lý thuyết đáp ứng tuyến tính (LRT), vốn chỉ áp dụng chính xác cho các hạt siêu thuận từ hoặc trong vùng từ trường nhỏ thỏa mãn điều kiện $H \ll H_C$ và $\xi < 1$. Chưa xét chi tiết các hiệu ứng phi tuyến khi biên độ từ trường ngoài vượt ngưỡng này.
    • Các khảo sát thực nghiệm gia nhiệt mới dừng lại ở mức khảo sát chất lỏng từ trong ống nghiệm (in vitro) phân tán trong dung môi chuẩn, chưa triển khai thử nghiệm trên mô hình sinh học sống (in vivo) hoặc tế bào ung thư cụ thể.
  • Hướng nghiên cứu tiếp:
    • Mở rộng mô phỏng và thực nghiệm sang vùng từ trường trung bình và cao sử dụng mô hình số kết hợp giữa LRT và SWMBTs để đánh giá đầy đủ đóng góp của tổn hao từ trễ.
    • Nghiên cứu chức năng hóa bề mặt hạt nano từ bằng các hợp chất tương thích sinh học cao và gắn các kháng thể hướng đích tế bào ung thư để tiến hành thử nghiệm nhiệt trị trên mô động vật.

Giá trị tham khảo

Bản luận án là tài liệu tham khảo có giá trị cao cho:

  • Nghiên cứu sinh, học viên cao học và nhà nghiên cứu thuộc các chuyên ngành Khoa học vật liệu, Vật lý chất rắn, Vật liệu điện tử, Công nghệ nano và Vật lý y sinh: cung cấp bức tranh lý thuyết hoàn chỉnh về nhiệt từ trị và các mô hình hồi phục Néel - Brown.
  • Các nhóm nghiên cứu thực nghiệm vật liệu nano: Cung cấp quy trình chi tiết về phương pháp thủy nhiệt chế tạo ferit hạt mịn, các kỹ thuật phân tích cấu trúc - từ tính (XRD, FESEM, VSM, DLS) và kỹ thuật đo đạc nhiệt từ thực nghiệm.
  • Kỹ sư và nhà phát triển ứng dụng y sinh / môi trường: Cung cấp cơ sở dữ liệu định lượng về việc lựa chọn vật liệu ($Fe_3O_4, MnFe_2O_4, CoFe_2O_4...$), khống chế kích thước hạt $D_{cp}$ và tần số kích thích tối ưu cho các thiết bị điều trị nhiệt từ và hệ thống nhả thuốc thông minh.

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

1. Hiệu ứng đốt nóng cảm ứng từ (MIH) của hạt nano từ khác với đốt nóng cảm ứng (IH) của kim loại khối thông thường như thế nào?

Ở kim loại khối dẫn điện tốt, hiệu ứng đốt nóng cảm ứng (IH) sinh nhiệt chủ yếu thông qua tổn hao Joule tạo bởi dòng điện xoáy (Foucault) bề mặt. Ngược lại, ở các hạt nano từ, do kích thước hạt rất nhỏ (cỡ vài đến vài chục nanomet) và độ dẫn điện của vật liệu oxit thấp, tổn hao Joule bề mặt là không đáng kể. Cơ chế sinh nhiệt của hiệu ứng MIH ở hạt nano từ chủ yếu bắt nguồn từ các tổn hao liên quan đến bản chất từ tính của hạt: tổn hao từ trễ ($P_{hys}$ ở hạt sắt từ) và tổn hao hồi phục Néel ($P_N$), tổn hao hồi phục Brown ($P_B$) ở hạt siêu thuận từ.

2. Điều kiện áp dụng của mô hình Lý thuyết Đáp ứng Tuyến tính (LRT) là gì?

Mô hình LRT được áp dụng khi phản ứng từ độ của hệ hạt tuyến tính với từ trường ngoài xoay chiều. Điều này xảy ra khi tham số $\xi = \frac{\mu_0 M_S V H}{k_B T} < 1$, tương ứng với:

  • Hệ các hạt nano từ ở trạng thái siêu thuận từ.
  • Hoặc hệ hạt sắt từ nhưng được đặt trong từ trường xoay chiều biên độ thấp thỏa mãn điều kiện $H \ll H_C$.

3. Điều kiện "cộng hưởng" để công suất tổn hao nhiệt (SLP) đạt cực đại trong mô hình LRT là gì?

Công suất đốt nóng cảm ứng từ $P_{LRT}$ đạt cực đại khi tần số góc của từ trường ngoài xoay chiều thỏa mãn hệ thức: $$\omega \tau = 2\pi f \tau = 1$$ Trong đó $\tau$ là thời gian hồi phục hiệu dụng của hệ hạt ($\frac{1}{\tau} = \frac{1}{\tau_N} + \frac{1}{\tau_B}$). Tương ứng với mỗi tần số $f$, sẽ tồn tại một kích thước hạt tối ưu xác định gọi là kích thước tới hạn $D_{cp}$ để hệ đạt $SLP_{max}$.

4. Dị hướng từ $K$ đóng vai trò phân chia cơ chế hồi phục như thế nào?

Luận án chỉ ra sự tồn tại của hằng số dị hướng tới hạn $K_C$:

  • Khi vật liệu có $K < K_C$ (như $FeCo, LSMO, MnFe_2O_4, Fe_3O_4$): Thời gian hồi phục Néel ngắn hơn Brown ($\tau_N < \tau_B$). Cơ chế Néel chi phối hoàn toàn đỉnh sinh nhiệt; đường cong $SLP(D)$ có đỉnh hẹp, nhạy cảm với độ rộng phân bố hạt $\sigma$ nhưng không phụ thuộc vào độ nhớt $\eta$.
  • Khi vật liệu có $K > K_C$ (như $CoFe_2O_4, FePt$): Thời gian hồi phục Brown ngắn hơn Néel ($\tau_B < \tau_N$). Cơ chế Brown chi phối đỉnh sinh nhiệt; đường cong $SLP(D)$ có đỉnh rộng, ít suy giảm theo $\sigma$ nhưng giá trị $SLP_{max}$ và vị trí $D_{cp}$ phụ thuộc chặt chẽ vào độ nhớt môi trường $\eta$.

5. Vì sao khi ứng dụng trong môi trường sinh học người ta thường ưu tiên vật liệu có cơ chế hồi phục Néel?

Trong cơ thể sống hoặc môi trường nội mô, độ nhớt của dịch sinh học có thể biến thiên phức tạp hoặc các hạt nano có thể bị cố định/bắt giữ bởi màng tế bào. Khi đó, các hạt sinh nhiệt theo cơ chế Brown (quay cả khối hạt) sẽ bị suy giảm công suất nhiệt nghiêm trọng do lực cản nhớt. Ngược lại, các hạt sinh nhiệt theo cơ chế Néel (đảo chiều spin bên trong hạt mà không cần hạt quay cơ học) vẫn duy trì trọn vẹn công suất tỏa nhiệt $SLP$, đảm bảo hiệu quả điều trị nhiệt từ ổn định.

Kết luận

Luận án tiến sĩ của Lưu Hữu Nguyên đã giải quyết một cách thấu đáo và hệ thống các vấn đề vật lý liên quan đến đặ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 chặt chẽ giữa mô hình hóa lý thuyết đáp ứng tuyến tính và thực nghiệm kiểm chứng có đối chứng. Công trình đã làm sáng tỏ quy luật cạnh tranh giữa hai cơ chế hồi phục Néel - Brown, xác lập vai trò phân nhóm của dị hướng từ tới hạn $K_C$, đồng thời kiểm chứng thành công các dự đoán lý thuyết trên hai hệ mẫu thực nghiệm $MnFe_2O_4$ và $CoFe_2O_4$. Đây là một công trình nghiên cứu hoàn chỉnh, có giá trị học thuật cao và mang lại những định hướng thực tiễn rõ nét cho việc thiết kế, tối ưu hóa các hệ hạt nano từ ứng dụng trong nhiệt từ trị y sinh và công nghệ môi trường.