当金属纳米结构的尺寸和入射光波长相当时，入射光能够在近场范围内穿过金属并诱导传导电子极化。在金属纳米结构的内部，入射光的电场使导电电子相对于固定在晶格中的离子发生位移，粒子内部会形成偶极矩，在电场力的作用下，内部电子就像弹簧振子一样相对平衡位置发生振荡，进而产生一种震荡波，这就是等离激元共振效应 [ 7 ] [ 8 ] 。纳米结构在入射光的照射下，光能可以在非辐射跃迁下转变为热能。等离激元纳米结构有着优异的吸光性能和光热转换特性 [ 9 ] [ 10 ] 。光热材料和光吸收材料可作为热介导光学操作中目标物体的基底和热源。在光照下，金属纳米结构的局部表面等离激元被激发，增加导电电子和晶格原子之间的碰撞频率，从而在金属纳米结构的激元表面建立温度梯度场 [ 9 ] [ 10 ] [ 11 ] [ 12 ] 。

图1是光热汽泡去除微塑料的原理示意图，激光束照射在等离激元基底上，纳米结构阵列发生等离子体共振，吸光集热并在聚光点处达到汽泡成核温度，形成汽泡，汽泡诱发局域Marangoni涡流，水中的微塑料通过循环对流进入高温汽泡中熔融成大块塑料，再通过常规沉降、过滤等方式有效分离和去除。

图1. 光热汽泡去除微塑料示意图

表1. 模型参数设置

光热转换的几何模型是如图2(a-c)所示的三种不同阵列单元的金纳米结构阵列，球形直径和三棱锥高为30 nm、长方体高为20 nm，阵列间距均为10 nm。模拟涡流的反应水池的直径0.6 cm，高0.2 cm。其他参数设置，见表1。

模拟中用实验室中常见的520 nm激光垂直照射三种不同的金纳米结构阵列，得到相对于入射场强的稳态相对电场强度和中心截面温度分布云图，如图2(d-f)和(g-i)所示。由于纳米结构中自由电子和光场相互作用产生的等离激元共振效应，光学近场范围内有局域电场极大增强的现象，局域电场的极值出现在纳米阵列的间隙和阵列结构的棱角处。这是因为纳米结构彼此靠近，建立了强等离子体耦合，使每个阵列单元暴露在入射光电场以及相邻多个单元的近场辐射中，每对阵列单元之间的电场相互叠加 [ 13 ] 。具有表面纳米结构的金薄膜类似于宽频带高效光吸收器，将吸收的光能转换为热能。此外，金属纳米结构在热管理方面也表现出优异的性能，可以用于控制热量的产生和传播 [ 9 ] [ 23 ] [ 24 ] [ 25 ] 。通过对入射电场放大倍数和产热聚热效果比较分析，得出以棱锥为单元的金纳米结构阵列的光热性能优于立方形、大幅度领先于球形，棱锥结构中叠加电场出现极值的相应位置出现局域高温，热量聚集不易扩散，表现出优良的吸光集热性能。

图2. 以球形、长方体和棱锥为单元的金纳米结构阵列的吸光集热特性：(a-c)为纳米阵列示意图；(d-f)为相对电场强度云图；(g-i)为截面温度分布云图

分别以1、5、10、20、100个太阳(1357 W/m²)折算的功率密度照射三种不同的纳米结构阵列，结果如图3所示，三种结构的最大温度T_max和最大相对电场强度(E_max/E_in)都随功率通量增加而变大，并且不同结构的最大温度/相对场强随功率增大差异越明显。棱锥单元的最高温度和场强放大倍数的极值均高于长方体和球形单元，在10个太阳照射下，其最大电场放大倍数达到23.56，局部最高温度达到447 K。

图3. 三种结构阵列的最大相对电场强度和最大温度随入射功率的变化

将10个太阳照射棱锥金纳米结构阵列得到的光热转换最大温度447 K作为热点的边界条件，考虑汽泡成核后，耦合传热和流体流动模拟得到图4所示的结果，图4(a)左半部分是池内温度分布，热汽泡表面和底面接触部分均认为达到纳米结构阵列的最高温度，图右半部分是涡量云图，分布的圆锥箭头为对应位置的速度矢量方向，箭头大小和颜色梯度反应速度大小，图4(b)对比在相同过热度下汽泡大小对池内涡流的影响。

池底面纳米结构为汽泡成核生长提供大量潜在位点，可以在金膜上通过激光照射产生汽泡，与常规的汽泡不同，这种汽泡在生成后可以保持与表面的接触，以截断球形的形状增大到一定大小 [ 26 ] 。由于汽泡可以包覆在纳米结构上，纳米结构周围的汽泡形成可以显著提高结构温度，汽泡消散的速度比较缓慢 [ 13 ] 。这使得纳米结构表面的汽泡可以成为一个稳定的高温“熔炉”。当液体中存在微小汽泡或者其他异物时，这些微小汽泡或异物可以作为汽泡核，促使液体沸腾并形成汽泡，汽泡大小取决于过热度、液体性质、汽泡核大小和数量等多种因素。

气液交界面处存在较大温度梯度，会导致界面上的应力不连续，产生表面张力梯度，进而引起传质现象，相比于自然对流即在流体密度梯度作用下的传质流动，Marangoni对流引起的涡旋较强，可以操控纳米级物体 [ 11 ] 。如图5(a)所示，汽泡附近的微颗粒在流体域中受到附加质量力F_M、曳力F_D、浮力F_B和剪切引起的升力F_L，颗粒位于气液界面时除了受到曳力外，还会受到表面张力F_S，另外压力F_P也会对颗粒运动产生影响，这三种力沿汽泡径向平衡使得粒子被束缚在汽泡表面 [ 27 ] [ 28 ] 。图5(a)中最大速度位于汽泡表面两侧和中下部位置，可以观察到流域内流体流动成漩涡状，流速较大的流体微团向热汽泡四周自上而下裹挟带动周围低速流体微团，并将其卷入汽泡下部，流体中微塑料颗粒在高温条件下融化，并在被卷吸进入汽泡的过程中不断熔合聚集成块，最后在汽泡内部汇聚。

改变汽泡直径，选取同一时间断面，从图4(b)中折线走势可以看出，直径较小的汽泡附近产生的涡量较大，这意味着，流体速度矢量旋转的强度越大，涡旋强度的增大可以使水体微塑料颗粒的加速度增大，强化颗粒在流体中运动能力。直径较小的汽泡能够提高流速，提高对微颗粒的控制能力，这对于捕获流体域中粒径稍大的微塑料颗粒是有利的，此外由于较小汽泡比表面积较大，对微颗粒的吸附性能更好，诱导阵列结构不同位点产生多个小汽泡，有利于提高对水体微颗粒的富集效率。

图4. 汽泡直径对速度/涡量大小及分布的影响

图5. (a) 粒子受力图；(b) 粒子的位置分布和轨迹图；(c) 粒子数和粒径对热汽泡收集颗粒效率的影响

用粒子追踪模块研究粒子的运动特点，在初始时刻流体域内随机位置释放粒子，前文模拟得到的流体瞬态速度场作为初始条件，粒子的位置和速度由拉格朗日方程确定，两侧壁条件为漫散射，本模拟中，粒子示踪是单向耦合流体流动，忽略流体与颗粒以及颗粒之间的相互作用，仅考虑Marangoni对流驱动下的颗粒运动。模拟中认为汽泡捕获颗粒为理想情况，给汽泡外边界设置粒子连续性条件即接触到热汽泡的粒子会被吞噬，并且给热汽泡内边界设置黏附条件即进入热汽泡内的粒子会附着在汽泡内壁上，以此来模拟汽泡富集微颗粒的效果。图5(b)显示的是反应池内释放的200个粒子在初始0时刻、第0.5 s和第3.8 s的位置分布图，热汽泡附近的颗粒在初始较短时间内就被热汽泡捕获，第3.8 s发现绝大多数粒子(占总数约90%)已经附着在热汽泡中下部，离热汽泡较远的部分粒子还未被捕获，随模拟计算时间推移，这些粒子最终汇集在热汽泡表面。在第3.8 s位置分布图的背景上绘制了各粒子从第0 s到第3.8 s的运动轨迹，轨迹图与利用Marangoni效应和热毛细对流进行微操作的相关研究中报道的现象相似，呈现涡旋流动特征 [ 17 ] [ 29 ] 。

图5(c)研究了颗粒浓度和粒径对收集效率的影响，当热汽泡内统计的颗粒数占总数的90%时可认为近乎收集完成，记录相应时刻，收集时间越长代表效率越低。模拟结果显示在一定浓度范围内收集效率受浓度变化影响较小，在粒子数的取值范围内，收集时间上下浮动0.05 s，占总时长1%左右；随着粒径增大，收集时间减少，粒径较大的颗粒质量大，在重力影响下，做涡旋运动的轨迹半径较小，更容易被卷吸进入热汽泡。