磁偶极相互作用对氧化铁纳米颗粒组装比吸收率的影响

摘要

我们计算了交替磁场作用下稀释氧化铁纳米颗粒随机簇团组装的比吸收率（SAR）。计算采用朗道-利夫希茨随机方程，同时包含磁矩热涨落和颗粒间磁偶极相互作用。对于常规三维簇团，相互作用强度主要由堆积密度η=N_p V/V_cl决定（N_p为簇团平均颗粒数，V为颗粒体积，V_cl为簇团体积）。观察到低频滞回面积和SAR随η从0.005增至<0.4显著下降。SAR对平均直径依赖性持存但随堆积密度增加而减弱。分形簇团（生物介质常见）中SAR进一步降低，峰值向小直径偏移。增加非磁性表面壳层厚度可部分恢复SAR值。

研究背景

磁流体热疗是癌症治疗领先策略，依赖交变磁场中磁性纳米颗粒（MNPs）产热。加热效率用SAR表达，取决于颗粒尺寸、材料、场幅、频率及周围介质性质。氧化铁纳米颗粒因生物相容性、可降解性、MRI可见性备受青睐。

体内MNPs常紧密结合组织，物理旋转受限，布朗弛豫可忽略，仅内禀弛豫——即颗粒内磁矩重取向——贡献发热。此外，紧密间隔颗粒间磁相互作用可显著改变SAR，尤其肿瘤内形成分形结构的致密聚集体中。

前人研究强调热涨落在稀释非相互作用组装SAR中作用，确定典型热疗频率（200-500 kHz）下最优颗粒直径~20 nm。但实验SAR常低于理论预测，可能源于致密组装中强磁偶极耦合。本研究旨在用严格随机朗道-利夫希茨框架量化该效应。

数值模拟

非相互作用纳米颗粒

稀释随机取向超顺磁组装，福克-普朗克方程产出近似动力学描述（式1）用于双能阱布居数n₁、n₂。归一化磁化强度遵循式(2)。这些方程经直接随机朗道-利夫希茨模拟验证，可快速计算不同颗粒直径、频率、场幅下滞回环。

纳米颗粒簇团

建模两类簇团形貌：(i) 准球形三维随机簇团，(ii) 细胞内聚集体典型分形簇团。三维情况下，直径D颗粒随机置于半径R_cl球体内且无直接接触。堆积密度η控制平均颗粒间距D_av：D_av=(6V_cl/(πN_p))¹ᐟ³，η=(D/D_av)³。分形簇团用Filippov算法生成，描述符D_f和k_f，得颗粒数N_p=k_f(2R_g/D)^D_f。

随机朗道-利夫希茨方程（式3）控制各颗粒单位磁化矢量α_i动力学，有效场源自各向异性（式5）、塞曼耦合（式6）、偶极-偶极相互作用（式7）。热场满足涨落-耗散关系（式9）。数值格式遵循文献[13,40,41]。

结果与讨论

非相互作用氧化铁纳米颗粒

采用典型氧化铁参数——饱和磁化强度M_s=70 Am²/kg、各向异性常数K=10⁴ J/m³、温度300 K、直径10-30 nm——固定场幅H₀=8 kA/m、频率200-500 kHz计算SAR。SAR在20-21 nm直径峰值达350-450 kW/kg，与最优尺寸颗粒实验报道一致。但许多体内测得SAR较低（~20-50 kW/kg），提示额外相互作用效应。

三维簇团组装

图3展示增加簇团堆积密度η如何减小滞回环面积和SAR。N_p=40时，η=0.005、0.04、0.32对应颗粒间距比D_av/D=1.46、2.92、5.84。η≥0.04时磁偶极耦合显著降低SAR，η从0.005增至0.32导致SAR六倍下降。最优直径仍近20 nm，但峰值变钝化。

分形簇团组装

分形聚集体显示更强SAR抑制。除最高分形维数（D_f=2.7）外，吸收峰系统性向小直径偏移。非磁性壳层厚度t_Sh=1 nm缓解最近邻偶极耦合，部分恢复SAR值并收窄直径依赖性，如图5所示。

结论

模拟证明纳米颗粒簇团内磁偶极相互作用——以堆积密度η量化——剧烈降低氧化铁组装SAR。常规三维簇团中η从0.005增至0.32使SAR下降约六倍，使理论预测与实验观测值（~50-100 kW/kg）吻合。分形簇团呈现类似SAR抑制及最优直径向小尺寸系统偏移。引入非磁性表面壳层可对抗偶极耦合恢复较高SAR，对生物医学应用极具前景。

这些发现强化了设计热疗用MNPs时控制颗粒尺寸、聚集态、表面化学的重要性。未来工作应纳入真实尺寸分布和潜在交换相互作用构建更全面模型。