Dampak Interaksi Magneto-Dipol terhadap Laju Serapan Spesifik Assembly Nanopartikel Besera Oksida
Abstrak
Kami menghitung laju serapan spesifik (SAR) dari assembly nanopartikel besera oksida yang diatur dalam kluster acak yang diterapkan pada medan magnet bergantian. Perhitungan menggunakan persamaan stokastik Landau–Lifshitz, yang menyertakan fluktuasi termal momen magnet dan interaksi magneto‑dipol di antara partikel. Untuk kluster 3‑D konvensional, kekuatan interaksi diatur terutama oleh kepadatan pengemasan η = NpV/Vcl, di mana Np adalah jumlah rata‑rata partikel per kluster, V adalah volume partikel, dan Vcl volume kluster. Kami mengamati penurunan yang signifikan dalam area histeresis frekuensi rendah dan dalam SAR seiring meningkatnya η dari 0,005 hingga <0,4. Ketergantungan SAR pada diameter rata‑rata terus ada tetapi melemah dengan kepadatan pengemasan yang lebih tinggi. Dalam kluster fraktal—umum dalam media biologis—SAR menurun lebih jauh, dengan puncak bergeser ke diameter yang lebih kecil. Peningkatan ketebalan cangkang non‑magnetik di permukaan secara parsial memulihkan nilai SAR.
Latar Belakang
Hipertemia magnetik adalah strategi terkemuka untuk pengobatan kanker, yang mengandalkan panas yang dihasilkan oleh nanopartikel magnetik (MNP) dalam medan magnet bergantian. Efisiensi pemanasan, yang dinyatakan sebagai SAR, bergantung pada ukuran partikel, bahan, amplitudo medan, frekuensi, dan properti medium di sekitarnya. Nanopartikel besera oksida sangat menarik karena biokompatibilitas, biodegradabilitas, dan visibilitas dalam MRI.
In vivo, MNP sering kali terikat erat dengan jaringan, yang membatasi rotasi fisik mereka. Akibatnya, relaksasi Brownian dapat diabaikan, dan hanya relaksasi Néel—yaitu, reorientasi momen magnet dalam partikel—yang berkontribusi pada pemanasan. Selain itu, interaksi magnetik di antara partikel yang berdekatan dapat mengubah SAR secara signifikan, terutama dalam agregat padat yang membentuk struktur fraktal di dalam tumor.
Karya sebelumnya telah menyoroti peran fluktuasi termal dalam SAR untuk assembly yang tidak berinteraksi dan mengidentifikasi diameter partikel optimal (~20 nm) untuk pemanasan maksimal pada frekuensi hipertemia tipikal (200–500 kHz). Namun, nilai SAR eksperimental sering kali di bawah prediksi teoretis, kemungkinan karena penguatan magneto‑dipol yang kuat dalam assembly padat. Studi ini bertujuan untuk mengukur efek tersebut menggunakan kerangka stokastik Landau–Lifshitz yang ketat.
Simulasi Numerik
Nanopartikel yang Tidak Berinteraksi
Untuk assembly superparamagnetik yang dilusi dan terorientasi secara acak, persamaan Fokker–Planck menghasilkan deskripsi kinetik perkiraan (Eq. 1) untuk populasi n1 dan n2 dari dua lubang energi. Magnetisasi tereduksi mengikuti Eq. (2). Persamaan ini, divalidasi terhadap simulasi stokastik Landau–Lifshitz langsung, memungkinkan perhitungan cepat loop histeresis untuk berbagai ukuran partikel, frekuensi, dan amplitudo medan.
Kluster Nanopartikel
Kami memodelkan dua morfologi kluster: (i) kluster acak 3‑D quasi-sferis dan (ii) kluster fraktal yang tipikal dari agregat seluler. Dalam kasus 3‑D, partikel dengan diameter D ditempatkan pada posisi acak dalam bola dengan jari-jari Rcl sambil memastikan tidak ada kontak langsung. Kepadatan pengemasan η mengontrol jarak antarpartikel rata‑rata Dav melalui Dav = (6Vcl/(πNp))1/3 dan η = (D/Dav)3. Kluster fraktal dihasilkan menggunakan algoritma Filippov dengan deskriptor Df dan kf, yang menghasilkan jumlah partikel Np = kf(2Rg/D)Df.
Persamaan stokastik Landau–Lifshitz (Eq. 3) mengatur dinamika vektor magnetisasi satuan αi dari setiap partikel, dengan medan efektif yang berasal dari anisotropi (Eq. 5), pengikat Zeeman (Eq. 6), dan interaksi dipole‑dipol (Eq. 7). Medan termal mengikuti hubungan fluktuasi‑dissipasi (Eq. 9). Skema numerik mengikuti Ref. [13, 40, 41].
Hasil dan Diskusi
Nanopartikel Besera Oksida yang Tidak Berinteraksi
Menggunakan parameter besera oksida tipikal—magnetisasi jenuh Ms = 70 Am²/kg, konstanta anisotropi K = 10⁴ J/m³, suhu 300 K, dan diameter 10–30 nm—kami menghitung SAR untuk amplitudo medan tetap H0 = 8 kA/m pada frekuensi 200–500 kHz. SAR mencapai puncak untuk diameter 20–21 nm, mencapai 350–450 kW/kg, sesuai dengan laporan eksperimental untuk partikel yang diukur optimal. Namun, banyak pengukuran in‑vivo melaporkan SAR yang lebih rendah (~20–50 kW/kg), yang menunjukkan efek interaksi tambahan.
Assembly Kluster 3‑D
Gambar 3 mengilustrasikan bagaimana peningkatan kepadatan pengemasan kluster η mengurangi area loop histeresis dan SAR. Untuk Np = 40, nilai η 0,005, 0,04, dan 0,32 sesuai dengan ras antarpartikel Dav/D sebesar 1,46, 2,92, dan 5,84, masing‑masing. Pada η ≥ 0,04, penguatan magneto‑dipol secara signifikan mengurangi SAR, dengan penurunan enam kali ketika η naik dari 0,005 ke 0,32. Diameter optimal tetap di sekitar 20 nm, meskipun puncaknya menjadi kurang jelas.
Assembly Kluster Fraktal
Agregat fraktal menunjukkan penekanan SAR yang lebih kuat. Puncak serapan bergeser ke diameter yang lebih kecil untuk semua dimensi fraktal kecuali yang tertinggi (Df = 2,7). Ketebalan cangkang non‑magnetik tSh sebesar 1 nm mengurangi penguatan dipole tetangga terdekat, secara parsial memulihkan nilai SAR dan menyempitkan ketergantungan diameter, sebagaimana digambarkan dalam Gambar 5.
Kesimpulan
Simulasi kami menunjukkan bahwa interaksi magneto‑dipol dalam kluster nanopartikel—dikwantifikasi oleh kepadatan pengemasan η—menurunkan SAR dari assembly besera oksida secara dramatis. Untuk kluster 3‑D konvensional, peningkatan η dari 0,005 ke 0,32 menurunkan SAR sekitar enam kali, menyelaraskan prediksi teoretis dengan nilai yang diamati secara eksperimental (~50–100 kW/kg). Kluster fraktal menunjukkan penekanan SAR yang serupa dan pergeseran sistematis diameter optimal ke ukuran yang lebih kecil. Memperkenalkan cangkang permukaan non‑magnetik dapat mengakomodasi penguatan dipole dan memulihkan SAR yang lebih tinggi, yang menjanjikan untuk aplikasi biomedis.
Temuan ini mempertegas pentingnya mengontrol ukuran partikel, status agregasi, dan kimia permukaan saat merancang MNP untuk hipertemia. Pekerjaan mendatang harus menyertakan distribusi ukuran yang realistis dan potensi interaksi pertukaran untuk model yang lebih komprehensif.
