Impacto de las interacciones magneto-dipolares en la tasa de absorción específica de ensamblajes de nanopartículas de óxido de hierro

Resumen

Calculamos la tasa de absorción específica (SAR) de ensamblajes diluidos de nanopartículas de óxido de hierro dispuestas en clusters aleatorios sometidas a un campo magnético alterno. Los cálculos emplearon la ecuación estocástica de Landau-Lifshitz, que incorpora tanto fluctuaciones térmicas de los momentos magnéticos como interacciones magneto-dipolares entre partículas. Para clusters 3-D convencionales, la fuerza de interacción se gobierna principalmente por la densidad de empaquetamiento η = N_p V / V_cl, donde N_p es el número promedio de partículas por cluster, V el volumen de partícula y V_cl el volumen del cluster. Observamos una reducción marcada en el área de histéresis a baja frecuencia y en SAR a medida que η aumentaba de 0,005 a <0,4. La dependencia del SAR en el diámetro medio persistió pero se debilitó con mayores densidades de empaquetamiento. En clusters fractales —comunes en medios biológicos— el SAR disminuyó aún más, con el pico desplazándose hacia diámetros menores. Aumentar el espesor de capas superficiales no magnéticas restauró parcialmente los valores de SAR.

Antecedentes

La hipertermia magnética es una estrategia líder para tratamiento de cáncer, basándose en calor generado por nanopartículas magnéticas (MNPs) en campo magnético alterno. La eficiencia de calentamiento, expresada como SAR, depende de tamaño de partícula, material, amplitud de campo, frecuencia y propiedades del medio circundante. Las nanopartículas de óxido de hierro son particularmente atractivas debido a su biocompatibilidad, biodegradabilidad y visibilidad en MRI.

In vivo, las MNPs a menudo se unen fuertemente a tejidos, limitando su rotación física. Consecuentemente, la relajación Browniana es insignificante, y solo la relajación Néel —reorientación del momento magnético dentro de la partícula— contribuye al calentamiento. Además, las interacciones magnéticas entre partículas cercanas pueden alterar significativamente el SAR, especialmente en agregados densos que forman estructuras fractales dentro de tumores.

Trabajos previos han destacado el papel de fluctuaciones térmicas en SAR para ensamblajes diluidos no interactuantes e identificado un diámetro de partícula óptimo (~20 nm) para calentamiento máximo a frecuencias típicas de hipertermia (200-500 kHz). Sin embargo, valores experimentales de SAR frecuentemente caen por debajo de predicciones teóricas, probablemente debido a fuerte acoplamiento magneto-dipolar en ensamblajes densos. Este estudio busca cuantificar ese efecto usando un marco estocástico riguroso de Landau-Lifshitz.

Simulación Numérica

Nanopartículas No Interactuantes

Para un ensamblaje superparamagnético diluido y orientado aleatoriamente, la ecuación de Fokker-Planck produce una descripción cinética aproximada (Ec. 1) para las poblaciones n₁ y n₂ de los dos pozos de energía. La magnetización reducida sigue Ec. (2). Estas ecuaciones, validadas contra simulaciones estocásticas directas de Landau-Lifshitz, permiten cálculo rápido de bucles de histéresis para varios diámetros de partícula, frecuencias y amplitudes de campo.

Clusters de Nanopartículas

Modelamos dos morfologías de cluster: (i) clusters 3-D aleatorios cuasi-esféricos y (ii) clusters fractales típicos de agregados intracelulares. En el caso 3-D, partículas de diámetro D se colocan en posiciones aleatorias dentro de una esfera de radio R_cl asegurando ningún contacto directo. La densidad de empaquetamiento η controla la distancia inter-partícula media D_av via D_av = (6 V_cl / (π N_p))¹ᐟ³ y η = (D / D_av)³. Clusters fractales se generan usando el algoritmo de Filippov con descriptores D_f y k_f, produciendo conteo de partículas N_p = k_f (2 R_g / D)^D_f.

La ecuación estocástica de Landau-Lifshitz (Ec. 3) gobierna la dinámica del vector de magnetización unitaria α_i de cada partícula, con campos efectivos derivados de anisotropía (Ec. 5), acoplamiento Zeeman (Ec. 6) e interacciones dipolo-dipolo (Ec. 7). Los campos térmicos obedecen la relación fluctuación-disipación (Ec. 9). El esquema numérico sigue Refs. [13, 40, 41].

Resultados y Discusión

Nanopartículas de Óxido de Hierro No Interactuantes

Usando parámetros típicos de óxido de hierro —magnetización de saturación M_s = 70 Am²/kg, constante de anisotropía K = 10⁴ J/m³, temperatura 300 K, y diámetros 10-30 nm— calculamos SAR para amplitud de campo fija H₀ = 8 kA/m a frecuencias 200-500 kHz. El SAR alcanza picos para diámetros 20-21 nm, llegando a 350-450 kW/kg, de acuerdo con reportes experimentales para partículas óptimamente dimensionadas. Sin embargo, muchas mediciones in vivo reportan SAR menor (~20-50 kW/kg), sugiriendo efectos de interacción adicionales.

Ensamblaje de Clusters 3-D

La Figura 3 ilustra cómo el aumento de densidad de empaquetamiento η disminuye el área del bucle de histéresis y SAR. Para N_p = 40, valores η de 0,005, 0,04 y 0,32 corresponden a razones inter-partícula D_av/D de 1,46, 2,92 y 5,84, respectivamente. A η ≥ 0,04, el acoplamiento magneto-dipolar reduce significativamente SAR, con caída de seis veces cuando η sube de 0,005 a 0,32. El diámetro óptimo permanece cerca de 20 nm, aunque el pico se vuelve menos pronunciado.

Ensamblaje de Clusters Fractales

Los agregados fractales muestran una supresión de SAR aún más fuerte. El pico de absorción se desplaza hacia diámetros menores para todas menos la dimensión fractal más alta (D_f = 2,7). Espesores de cáscara no magnética t_Sh de 1 nm mitigan el acoplamiento dipolar de vecinos más cercanos, restaurando parcialmente valores SAR y estrechando la dependencia del diámetro, como se muestra en Figura 5.

Conclusiones

Nuestras simulaciones demuestran que las interacciones magneto-dipolares dentro de clusters de nanopartículas —cuantificadas por la densidad de empaquetamiento η— reducen dramáticamente el SAR de ensamblajes de óxido de hierro. Para clusters 3-D convencionales, aumentar η de 0,005 a 0,32 baja SAR aproximadamente en factor seis, alineando predicciones teóricas con valores observados experimentalmente (~50-100 kW/kg). Clusters fractales exhiben supresión similar de SAR y desplazamiento sistemático del diámetro óptimo a tamaños menores. La introducción de capas superficiales no magnéticas puede contrarrestar el acoplamiento dipolar y recuperar SARs más altos, lo que es prometedor para aplicaciones biomédicas.

Estos hallazgos refuerzan la importancia de controlar tamaño de partícula, estado de agregación y química superficial al diseñar MNPs para hipertermia. Trabajos futuros deberían incorporar distribuciones de tamaño realistas e interacciones de intercambio potenciales para un modelo más completo.