Impact van magnetische-dipool-interacties op de specifieke absorptiesnelheid van ijzeroxide-nanopartikel-assemblages

Samenvatting

We berekenden de specifieke absorptiesnelheid (SAR) van verdunne assemblages van ijzeroxide-nanopartikels gerangschikt in willekeurige clusters onderworpen aan een wisselend magnetisch veld. De berekeningen gebruikten de stochastische Landau-Lifshitz-vergelijking, die zowel thermische fluctuaties van de magnetische momenten als magnetische-dipool-interacties tussen de deeltjes incorporeert. Voor conventionele 3-D clusters wordt de interactiestructuur primair bepaald door de pakkingsdichtheid η = N_pV/V_cl, waar N_p het gemiddelde deeltjes-aantal per cluster is, V het deeltjesvolume, en V_cl het clustervolume. We observeerden een duidelijke reductie in de lage-frequentie-hysteresis-oppervlakte en in de SAR toen η toenem van 0,005 tot <0,4. De SAR's afhankelijkheid van gemiddelde diameter bleef bestaan maar verzwakte met hogere pakkingsdichtheden. In fractale clusters – veelvoorkomend in biologische media – daalde de SAR verder, met de piek verschuivend naar kleinere diameters. Het vergroten van de dikte van niet-magnetische oppervlakte-schillen herstelde de SAR-waarden deels.

Achtergrond

Magnetische hyperthermie is een toonaangevende strategie voor kankerbehandeling, die vertrouwt op warmte gegenereerd door magnetische nanopartikels (MNPs) in een wisselend magnetisch veld. De verwarmings-efficiëntie, uitgedrukt als SAR, hangt af van deeltjesgrootte, materiaal, veldamplitude, frequentie en de omringende medium-eigenschappen. Ijzeroxide-nanopartikels zijn met name aantrekkelijk door hun biocompatibiliteit, biologische afbreekbaarheid en zichtbaarheid in MRI.

In vivo raken MNPs vaak stevig gebonden aan weefsels, wat hun fysieke rotatie beperkt. Gevolgelijk is Browniaanse relaxatie verwaarloosbaar, en alleen Néel-relaxatie – namelijk heroriëntatie van het magnetische moment binnen het deeltje – bijdraagt aan verwarming. Bovendien kunnen magnetische interacties tussen dichtbij gelegen deeltjes de SAR aanzienlijk veranderen, vooral in dichte aggregaten die fractale structuren vormen binnen tumoren.

Eerlier werk benadrukte de rol van thermische fluctuaties in SAR voor verdunne, niet-interacterende assemblages en identificeerde een optimale deeltjesdiameter (~20 nm) voor maximaal verwarming bij typische hyperthermie-frequenties (200–500 kHz). Experimentele SAR-waarden vallen echter frequent onder theoretische voorspellingen, waarschijnlijk door sterke magnetische-dipool-koppeling in dichte assemblages. Deze studie beoogt dat effect te kwantificeren met een rigoureus stochastisch Landau-Lifshitz-raamwerk.

Numerieke Simulatie

Niet-interacterende Nanopartikels

Voor een verdunne, willekeurig georiënteerde superparamagnetische assemblage levert de Fokker-Planck-vergelijking een benaderde kinetische beschrijving (Verg. 1) voor de populaties n₁ en n₂ van de twee energieputten. De gereduceerde magnetisatie volgt Verg. (2). Deze vergelijkingen, gevalideerd tegen directe stochastische Landau-Lifshitz-simulaties, maken snelle berekening van hysterese-lussen mogelijk voor diverse deeltjesdiameters, frequenties en veldamplituden.

Nanopartikel-clusters

We modelleren twee cluster-morfologieën: (i) quasi-sferische 3-D willekeurige clusters en (ii) fractale clusters typisch voor intracellulaire aggregaten. In het 3-D geval worden deeltjes van diameter D op willekeurige posities binnen een straal R_cl geplaatst terwijl direct contact vermeden wordt. De pakkingsdichtheid η controleert de gemiddelde inter-deeltjesafstand D_av via D_av = (6V_cl/(πN_p))^1/3 en η = (D/D_av)³. Fractale clusters worden gegenereerd met het Filippov-algoritme met descriptors D_f en k_f, wat een deeltjes-aantal N_p = k_f(2R_g/D)^D_f oplevert.

De stochastische Landau-Lifshitz-vergelijking (Verg. 3) regeert de dynamica van elk deeltjes eenheids-magnetisatievector α_i, met effectieve velden afgeleid van anisotropie (Verg. 5), Zeeman-koppeling (Verg. 6) en dipool-dipool-interacties (Verg. 7). Thermische velden gehoorzamen de fluctuatie-dissipatie-relatie (Verg. 9). Het numerieke schema volgt Refs. [13, 40, 41].

Resultaten en Discussie

Niet-interacterende Ijzeroxide-nanopartikels

Met typische ijzeroxide-parameters – verzadigingsmagnetisatie M_s = 70 Am²/kg, anisotropie-constante K = 10⁴ J/m³, temperatuur 300 K, en diameters 10–30 nm – berekenden we SAR voor een vaste veldamplitude H₀ = 8 kA/m bij frequenties 200–500 kHz. De SAR piekt voor diameters 20–21 nm, bereikend 350–450 kW/kg, in overeenstemming met experimentele rapporten voor optimaal gegroeide deeltjes. Echter, veel in-vivo metingen rapporteren lagere SAR (~20–50 kW/kg), wat additionele interactie-effecten suggereert.

Assemblage van 3-D Clusters

Figuur 3 illustreert hoe toenemende cluster-pakkingsdichtheid η de hysterese-lus-oppervlakte en SAR vermindert. Voor N_p = 40, corresponderen η-waarden van 0,005, 0,04 en 0,32 met inter-deeltjes-verhoudingen D_av/D van 1,46, 2,92 en 5,84 respectievelijk. Bij η ≥ 0,04 vermindert magnetische-dipool-koppeling SAR significant, met een zesvoudige daling wanneer η stijgt van 0,005 tot 0,32. De optimale diameter blijft nabij 20 nm, hoewel de piek minder uitgesproken wordt.

Assemblage van Fractale Clusters

Fractale aggregaten tonen een nog sterkere SAR-onderdrukking. De absorptiepiek verschuift naar kleinere diameters voor alle behalve de hoogste fractale dimensie (D_f = 2,7). Niet-magnetische schil-diktes t_Sh van 1 nm verzwakken nearest-neighbour dipool-koppeling, wat deels SAR-waarden herstelt en de diameter-afhankelijkheid vernauwt, zoals afgebeeld in Figuur 5.

Conclusies

Onze simulaties tonen aan dat magnetische-dipool-interacties binnen nanopartikel-clusters – gekwantificeerd door de pakkingsdichtheid η – de SAR van ijzeroxide-assemblages dramatisch verlagen. Voor conventionele 3-D clusters verlaagt het verhogen van η van 0,005 tot 0,32 de SAR met ongeveer een factor zes, wat theoretische voorspellingen in lijn brengt met experimenteel waargenomen waarden (~50–100 kW/kg). Fractale clusters vertonen vergelijkbare SAR-onderdrukking en een systematische verschuiving van de optimale diameter naar kleinere maten. Het introduceren van niet-magnetische oppervlakte-schillen kan dipool-koppeling tegengaan en hogere SARs herwinnen, wat veelbelovend is voor biomedische toepassingen.

Deze bevindingen versterken het belang van het controleren van deeltjesgrootte, aggregatiestoestand en oppervlaktechemie bij het ontwerpen van MNPs voor hyperthermie. Toekomstig werk zou realistische grootteverdelingen en potentiële uitwisselingsinteracties moeten incorporeren voor een meer compleet model.