激光辅助近场纳米制造中的物理机制探索 下载: 1183次
Significance Laser-assisted near-field nanomanufacturing uses a near-field focused laser beam to break down the diffraction limit and heat materials to induce phase change or phase explosion to fabricate nanoscale materials and complex structures. It has an outstanding feature in both academic and industrial fields due to its highly coherent features, including continuously adjustable incident laser power, highly controllable processing position, and accessible to nanodomains. Among various near-field technologies, the tip-based near-field technique utilizes the tiny but sharp geometry of scanning probe microscope (SPM) tips to focus the incident optical field into an extremely small area in proximity to the tip apex. The highly enhanced electromagnetic field generates huge but localized heat in the surface to be manufactured and modifies its morphology through photon absorption and consequent phase change at the nanoscale.
Progress In the processing domain, in-situ information about the optical field, temperature rise, stress, and material structure evolution is critical for understanding and refining nanomanufacturing. It is helpful for in-depth understanding of the physical mechanism of multiphysics interaction and further optimization and process control. The enhanced electromagnetic field in and around SPM tips has been fully studied in recent decades; however, a knowledge gap remains relative to temperature rise and thermal stress evolution in the same region. A high-fidelity simulation on atomic force microscopy (AFM) tungsten tip under laser irradiation demonstrated that the electric field was primarily concentrated at the apex outside the tip rather than inside (
Micro/nanoparticles represent an alternative geometry employed to introduce near-field focusing and simultaneous heating at the sub-wavelength scale, especially in the application of laser-assisted nanopatterning and nanolithography on a large-area substrate. Knowledge of the temperature and thermal stress beneath particles is critical; however, acquiring this knowledge remains a significant challenge for immediate sensing. Tang et al. has conducted pioneering experimental work on direct measurement of temperature and thermal stress in a silicon substrate beneath various transparent microstructures, including a single silica particle, a glass fiber, and a layer of particles (
Although various experimental methods have been developed for nanoscale probing, the theoretical simulation, e.g., molecular dynamic simulation, is a predominant method to predict atomic-level phenomena in an ultrafast period in laser-assisted near-field nonmanufacturing processes. In addition, the previous work primarily focused on enhancing the optical field under SPM tips, the light confinement effect, and the temperature evolution of the tip and substrate only under the condition of no phase change. However, the incident laser is extremely focused and enhanced in a nanodomain in the surface to be processed under the tip; thus, the material in this domain undergoes intense heating, phase change, phase explosion, stress generation and propagation, and rapid structural evolution in a very short time. Computer simulation research has focused more on small-scale and short-time heating processes, which would lose important information on material/structure evolution. Wang established a large-scale system with hundreds of millions of atoms and studied the long-term behavior and structural evolution of heating, melting, phase transition, solidification, and defect formation in laser-assisted tip-based near-field nanomanufacturing (Figs. 9--14). In addition, the shock wave during the nanomanufacturing process has been studied extensively (Figs. 16--18). Although simulation results were based on a few assumptions, the models used in the simulations were designed according to real dimensions, and this provided cutting-edge and detailed knowledge and understanding of physical images.
Conclusions and Prospects This review primarily focuses on experimental and theoretical investigations of the optical, temperature, and stress fields, and simulation study of the structural evolution in SPM tip-based laser-assisted near-field nanomanufacturing. However, unsolved problems remain in laser-assisted near-field nanomanufacturing. In the ultrafast and dynamic process of nanostructure manufacturing, the physical phenomena in such a small domain (e.g., nano to subnano domains) are difficult to observe and detect. Although newly developed laser-assisted probe methods are nanoscale accessible, they take much more time to satisfy signal collecting than the entire time range of ultrafast dynamic process of nanomanufacturing. In addition, SPM tips are fragile and can be easily damaged and contaminated, and their geometric shape change significantly affects the near-field electromagnetic field distribution. The unexpected change in tips shape will then widen the near-field light field distribution around the needle apex and make the manufactured structure deviate from expectation. The high temperature rise in laser-assisted near-field nanomanufacturing will introduce large thermal expansion in both the tip and the processed surface, and yield unstable measurements because the distance between the tip and substrate is only several nanometers. In terms of rapid response measurement, large-scale manufacturing, and improving and expanding industrial applications, this technology still has room for significant development and improvement.
1 引言
纳米制造和加工涉及纳米尺度的材料、结构以及相关装置和系统的制造[1-3]。纳米制造过程可分为两种基本类别,即自上而下和自下而上或自组装过程,其中包含多种多样的化学和物理技术,例如化学合成、化学气相沉积、分子束外延、原子层外延、光刻、自组装以及很多新开发的技术[4-6]。在这些技术中,激光辅助技术由于激光光源能量连续可调、加工位置高度可控以及其他各种物理参数可协调性高而在纳米制造和加工领域被广泛应用。为了获得低于衍射极限的分辨率,在过去的几十年中,研究人员陆续开发出了基于激光辅助技术的近场纳米制造技术,例如,激光辅助纳米图形化技术和纳米压印光刻技术、近场扫描光学显微镜(NSOM)技术、超透镜技术、纳米光圈技术以及扫描探针显微镜(SPM)技术等。
基于激光辅助探针近场技术的制造方法是使用SPM探针针尖进行激光辅助制造的一大分支。在此技术中,SPM探针针尖很小,而且针尖的几何形状使其成为无孔NSOM的理想工具。当激光照射SPM探针针尖时,针尖可以使激光的传播方向弯曲,并使激光聚焦,增强针尖尖端附近的光场。使用扫描隧道显微镜(STM)探针可以实现大于100倍的光场增强倍率[7]。这种增强的光场会在样品中产生大量的热,并导致表面发生纳米级的相变,最终导致表面形貌发生变化。因此,可以利用这种增强光场的技术开发新的制造方法。与STM相似,原子力显微镜(AFM)探针的针尖也可以将激光束聚焦至纳米尺度,实现纳米级的表面制造。
本文重点介绍激光辅助近场纳米制造过程中纳米域中光场、温度、应力的实验探测和计算机模拟方面的工作。本文第二部分主要介绍AFM探针辅助近场聚焦激光加热时温度和光场分布的实验和模拟,第三部分主要讨论AFM探针尖端附近硅基底温度测量和光场分布的模拟研究,第四和第五部分分别介绍微米/亚微米颗粒下硅板的近场热、光、和力学现象的实验和计算机模拟研究,第六部分主要讨论近场激光加热下材料结构变化的模拟工作,包括结构演变和激波的影响。
2 基于AFM探针的近场聚焦和加热:针尖效应
已有大量工作报道了纳米级针尖内部和周围的增强光场分布,但该区域温度分布和传热过程的研究报道较少。针对AFM针尖内部和周围区域的传热问题,Chen等[8]系统地模拟了激光辐照下AFM钨探针和硅基底中的电磁场和温度场,以及钨探针中的近场热输运过程。
2.1 激光照射下AFM探针的光学和热学物理过程模拟
在Chen等[8]的工作中,AFM探针可以近似为悬于硅基底上方的一个圆锥形针尖,如
图 1. 激光辐照下探针针尖的近场光、热和应力模拟[8]。(a)探针针尖-基底模型结构示意图;(b)探针针尖的几何模型;(c)激光强度、针尖温度和热膨胀随时间的变化;(d)针尖周围的坡印廷矢量分布
Fig. 1. Simulation of near field optical, thermal, and stress fields of the tip under laser irradiation[8]. (a) Configuration of the tip-substrate system under modeling; (b) geometric configuration of the tip;(c) evolution of the laser intensity, apex temperature, and elongation (thermal expansion) with time; (d) Poynting vector distribution around the tip
坡印廷矢量DE模拟结果如
以下三个因素会显著影响针尖电场的增强:激光的偏振方向、探针的几何形状以及尖端与基底之间的距离。由于入射光中只有沿探针轴向的分量对尖端附近的电场增强有贡献,因此电场强度的增强取决于cos2φ(φ是入射激光偏振方向相对于探针轴向的极化角)。研究发现,探针的半锥角和直径都会影响电场增强。当半锥角θ为10°时,电场增强的最佳半径r约为9 nm。对于固定半径为30 nm的尖端,当θ小于25°时,光场增强随θ的增大而增加;当θ在25°~30°范围内时,光场的增强几乎保持恒定;当θ大于35°时,光场增强随着θ的进一步增大而减小。此外,随着尖端与基底距离的增加,光场增强的峰值呈指数下降,峰值位置发生偏移,并逐渐远离尖端所在轴,当尖端接触基底时,该偏移消失。
Chen等[8]通过FLUENT软件模拟研究了室温(300 K)下施加单个脉冲的探针尖端的温度演变和分布。相较于激光脉冲的峰值能量时间,探针中的热传导过程在更大程度上延迟了探针尖端最高温度的出现,如
2.2 激光照射下AFM探针温度的实验研究
为了验证理论模拟的结果,Chen等[9]在具有微米级空间分辨率的激光辐照下,使用AFM硅探针进行了相应的测温实验。在他们的工作中,锑掺杂单晶硅探针(CSG10)被固定在空间分辨率可达20 nm的三维压电驱动平台上,激光从探针侧面入射。AFM硅探针的电镜图如
图 2. 激光辐照下硅探针针尖上微米量级的热响应[9]。(a)AFM硅探针正面的电镜图;(b)AFM硅探针侧面的电镜图;(c)激光照射探针侧面时,焦平面处激光光斑与探针的位置关系;(d)实验中激光光斑相对于探针的移动方向;(e)探针在焦点附近沿z方向移动时的拉曼位移、温度和拉曼强度;(f)激光强度分别为5.6×108 W/m2和2.5×108 W/m2时的拉曼强度以及位移沿图2 (d)y方向扫描线1的变化
Fig. 2. Microscale spatially resolved thermal response of Si tip to laser irradiation[9]. (a) SEM image of the front plane of the AFM silicon tip; (b) SEM image of the side plane of the AFM silicon; (c) positional relationship between the laser beam at the focal spot and the tip when the laser illuminates the tip side; (d) the moving directions of laser spot with respect to the tip during the experiment; (e) Raman shift, temperature, and Raman intensity v
为研究近场激光的加热效果,Chen等[9]分别采用4,7,12.7,15.9 mW的激光照射硅探针尖端,分析了硅探针拉曼特征峰的位移,以确定激光辐照下探针的温升(当拉曼探测区域的温度升高时,特征峰的拉曼位移变小),然后根据McCarthy等[10]的理论计算出了准确的温升。在实验测量过程中,Chen等[9]通过监测探针在激光照射下的衍射图像来监控聚焦的水平,从而消除了散焦对实验结果的影响。所测光谱中的硅拉曼特征峰均利用高斯函数进行拟合,以提高拉曼位移测量的准确性。
Chen等[9]首先研究了探针沿z方向移动时探针内部的温升[在整个z方向上移动探针的行程中,激光光斑完全覆盖探针,如
Chen等[9]进一步研究了探针上的近场加热机理,他们使用2.1节所述有限元方法来计算探针尖端周围以及内部的电场分布。计算结果显示,探针尖端外侧和内部的最高电场增强分别为5.19和2.86。他们使用一维传热模型计算了探针轴向的温度分布[9],结果表明,在半球顶点获得的最高温度为175 ℃,略低于实际测得的温度(约327 ℃)。两者之间的差异可能是因为硅探针中的掺杂剂降低了探针的热导率,以及低估了激光在探针上的散射或吸收效率。
3 基于AFM探针的微纳尺度近场聚焦和加热:基底的影响
Chen等[9]的研究结果显示,当外部激光照射AFM探针时,探针的近场作用会使其尖端附近的电磁场增强,同时在靠近探针的基底中也会形成局部加热,并产生较大的温升。Yue等[11]首次实现了基底中温度的直接测温,他们同样采用拉曼位移法测量了AFM探针顶端下方硅基底表面以下10 nm区域内的温度分布。
3.1 实验方法
在测温之前,Yue等[11]首先校准了硅样品从室温到230 ℃之间的特征峰拉曼位移。他们在热校准实验中使用功率可调的电加热器加热硅样品,并通过T型热电偶测量硅样品的实际温度。同时,他们采用拉曼光谱仪及785 nm激发光激发并探测硅的拉曼峰,然后对拉曼峰进行高斯函数拟合,以精确获取硅的拉曼位移。该位移与温度之间的线性关系将用于下一步实验中确定硅基底的温度。
实验装置如
图 3. 近场激光加热下非接触式温度测量[11]。(a)无孔NSOM测温实验装置光学系统示意图;(b)硅基底表面温度与入射光斑照射针尖位置之间的关系(横坐标中激光散射强度弱表示激光光斑照射在针尖尖端,而激光散射强度强则表示激光光斑照射在针尖上方靠近悬臂梁的位置)
Fig. 3. Noncontact temperature measurement in near-field laser heating[11]. (a) Schematic of experimental setup for thermal probing using the apertureless NSOM; (b) relationship between the measured silicon surface temperature and laser spot position on the tip (a weak laser scattering intensity in the x axis indicates the laser spot is on tip apex, while a strong scattering intensity indicates the laser spot is on the upper part of the tip near t
3.2 近场光、热效应的模拟研究
针对上述近场加热实验,Yue等[11]对物理模型进行了光电场增强的有限元分析,并解释了上述实验的物理机制。他们应用HFSS V13软件计算了
图 4. 硅探针针尖上近场光场的模拟[11]。(a)探针-基底模型示意图;(b)探针周围的电场在y=0平面分布的前视图;(c)探针周围的电场在探针尖端下方分布的俯视图和横截面图;(d)A-A'平面中的电场分布图;(e)横截面中的电场分布图; (f)探针和硅基底热阻计算示意图
Fig. 4. Simulation of near-field optical field in silicon tip[11].(a) Tip-substrate domain in the simulation; (b) front view of the electric field around the tip distributed in the y=0 plane; (c) top view and cross-sectional view of the electric field around the tip under the tip apex; (d)electric field distribution in the A-A' plane; (e) electric field distribution in the cross-section; (f) schematic of the thermal resistance calculation o
3.3 近场聚焦作用下基底中的非线性光学吸收和弹道热传导
在近场加热区域中,当增强光场产生的热量与从加热区域向外传输的热量达到平衡时,可以通过建立传热模型来计算近场加热区域中的温升。传热模型的表达式为
式中:ΔV是近场加热域的体积,ΔV≈2π
式中:keff是局域有效热导率。但是这个极小的加热区域尺寸(约10 nm)小于室温下硅晶片中声子的平均自由程(约43 nm),由于声子的散射和碰撞不足,无法实现热平衡[13]。局域有效热导率keff可以通过弹道传热效应来确定,其表达式为
式中:k是硅晶片块材的热导率;τ为无量纲径向距离,τ=r/Λ,Λ是硅晶片块材的声子平均自由程。因此,在本例计算域中,有效热导率keff约为2.5 W/(m∙K)。在偏振角为15°和30°的情况下,模拟计算得到了该10 nm区域内的温升分别为235 ℃和105 ℃,这一结果与实验结果非常吻合。在偏振角为75°的情况下,光场增强较小,温度升高可以忽略不计。
4 微纳米球辅助近场聚焦技术:实验研究
微纳米球辅助近场聚焦技术是另一种有效的、可突破光学衍射极限的技术。微纳米球可在亚波长范围内聚焦入射激光并产生局域加热。该技术已被广泛应用于大面积基板上的激光辅助纳米图案化和纳米光刻。目前,探测纳米尺度极小区域(等效直径约为100 nm或更小)中的温度场、应力场和光场仍然是一个巨大的挑战,这是因为加热区域太小,传统的温度传感器无法直接进行测温。Tang等[14-16]首次使用拉曼光谱法分别测量了单个氧化硅微米球、单根玻璃纤维和单层氧化硅微米球阵列中硅片的温度场和应力场分布,通过实验开创性地研究了基于微米球的近场聚焦对硅片的影响。
4.1 实验方法
Tang等[14]采用倾斜技术在干净的硅片上制备了单层自组装氧化硅纳米颗粒薄膜,氧化硅纳米颗粒的直径分别为200,400,800,1210 nm。测量温度场、应力场和光场的实验方案如
图 5. 1210 nm直径氧化硅微米球在辅助近场加热下的物理场的纳米量级扫描[14]。(a)不同激光能量下,拉曼强度I沿扫描方向的变化;(b)激光强度为3.1×109 W/m2时硅的拉曼强度I、拉曼位移ω和线宽Γ沿x方向的变化;(c)激光入射轴与拉曼散射收集轴之间的位置偏移;(d)氧化硅微米球相对于激光入射轴和拉曼散射收集轴的位置
Fig. 5. Nanoscale mapping of physics fields under 1210 nm silica micro-particles assisted near-field heating[14]. (a) Under different laser irradiation energy, the variation of Raman intensity I along the mapping direction; (b) variations of Raman intensity I, Raman shift ω, and linewidth Γ of silicon in the x direction with laser intensity of 3.1×109 W/m2; (c) deviation between the laser beam axis
4.2 解耦温度、应力和光场
针对获取的硅拉曼光谱特征峰,首先使用高斯函数进行拟合,精确获取其拉曼特征峰强度、拉曼位移和拉曼线宽。拉曼特征峰强度、拉曼位移和拉曼线宽会随温度而变化:当温度升高,拉曼特征峰强度和拉曼位移减小,而拉曼线宽变大。其中,拉曼线宽与声子寿命密切相关,它取决于测量区域内的温度,而受应力的影响较小。因此,基于线宽的测温方法可用于直接测量温度。拉曼位移同时受到局部应力和温度的影响,而拉曼强度主要受入射激光强度和温度的影响。因此,通过分别分析拉曼强度和拉曼位移可以同时得到光场强度的变化和应力分布。Tang等[14]利用该方法分析了周期性结构下方硅片的拉曼特征峰强度、拉曼位移和线宽在空间中的变化,分别测量了近场光加热下温度、局部应力和光场的周期性变化。
系统中的光学准直和激光能量变化会影响拉曼强度测量的准确性,但与样品无关;因此,可在相同的测量条件下,用氧化硅微米球下方的硅晶片拉曼强度
知道硅晶片的温度分布后,可以利用拉曼位移ω和拉曼线宽Γ对温度响应灵敏度的差异来确定氧化硅微米球下方硅晶片内部的热应力σ。线宽Γ对应力的响应灵敏度低,故可忽略应力对线宽的影响。因此,实验测得的拉曼线宽和拉曼位移可以分别表示为ΔΓexp = ΔΓΔT + ΔΓfocus和Δωexp = ΔωΔT + Δωfocus + Δωσ。其中下标ΔT表示温度升高的贡献值,下标focus表示微米球引起的离焦效应使线宽产生的变化量;下标σ表示由应力引起的变化量。依据热校准获得的线宽温度系数以及根据拉曼强度计算得到的温升ΔT,可以得到ΔΓΔT,并由ΔΓexp表达式获得离焦对线宽的影响ΔΓfocus。而离焦过程中由散焦效应引起的拉曼位移Δωfocus和线宽ΔΓfocus之间的关系校准为Δωfocus=-0.210-0.767·ΔΓfocus。ΔωΔT由拉曼位移温度系数和温升确定,再由Δωexp关系式可得到源于应力作用的Δωσ。再依据文献[ 17]报道的拉曼位移/应力比值(-3.6 cm-1/GPa),就可计算得到硅晶片内部的应力σ。随着激光能量从25%增加到100%,直径为1210 nm的氧化硅微米球下方硅晶片的内部应力从140 MPa增大到370 MPa,不确定性为40 MPa。
Tang等[14]还测量了直径为800,400,200 nm等较小尺寸微米球下方硅晶片内部的温度场和应力场分布,结果发现,最大拉曼强度比(Imax/Imin)随着微米球直径的减小而减小,如
图 6. 氧化硅微米球辅助近场检测的空间极限[14]。(a)Imax/Imin随氧化硅微米球直径的变化趋势;(b)200 nm氧化硅微米球下方硅晶片的拉曼强度I沿x方向的变化
Fig. 6. Spatial resolution limit of silica microsphere assisted near-field detection[14].(a) Trend of Imax/Imin against the diameter of the silica particle; (b) variation of Raman intensity I along the x direction for a silicon wafer under a 200 nm diameter silica particle
5 近场热、光和应力场:模拟
Tang等[14-16]还采用有限元方法研究了1210 nm氧化硅微米球及其基底的光电场分布。为了优化网格密度和计算量,他们的模型使用了四分之一的氧化硅微米球和硅晶片基底。入射光(波长λ=532 nm)设为平面波,从微米球顶部入射。实验中的激光光斑小于氧化硅微米球的直径。仿照具体的实验过程,他们的模型中仅保留了部分激光照射区域内的微米球,对称平面设为完美磁场和电场对称边界,其他边界设为吸收边界。
图 7. 氧化硅微米球和硅基底中光电场和温度场模拟以及误差分析[14]。(a)当激光照射在单颗氧化硅微米球上和两颗氧化硅微米球之间时的电场分布模拟结果;(b)在近场聚焦激光加热下,1210 nm氧化硅微米球下方硅晶片的内部温度沿z方向的分布,插图显示了基底表面的温度分布;(c)物镜聚焦光轴与微米球1、2中的拉曼散射光路示意图;(d)z方向离焦效应对硅峰拉曼强度的影响
Fig. 7. Simulation and error analysis of electric and temperature fields in silica particles and the silicon substrate[14]. (a) Modeling results of the electric field for two situations: laser spot on a single particle and between two particles; (b) under the near-field focused laser heating, silicon wafer temperature variation along the thickness (z) direction directly under a 1210 nm silica particle, and the inset shows the temperature distribut
Tang等[15]将模型及温度分布从FLUENT导入至ANSYS V14中,进一步模拟了单个氧化硅微米球下方硅基底中的应力分布。模型的外侧和底部为固定边界,顶部边界不受限,两个垂直横截面都为对称边界。模拟得到的硅晶片基底内部的热应力分布如
图 8. 硅基底表面和内部横截面的光场、温度场和热应力场[15]。(a)光场;(b)温度场;(c)由温度升高引起的热应力场
Fig. 8. Optical field, temperature field, and thermal stress field at the surface and in cross-section of the silicon substrate[15]. (a) Optical field; (b) temperature field; (c) thermal stress field introduced by the temperature rise
6 近场激光加热下的材料结构演变
由于纳米级针尖的近场聚焦增强了入射激光的强度,因此针尖尖端下方材料的纳米区域在极短时间内经历了剧烈的加热、相变、相爆炸、应力产生和传播,以及快速的结构演变,最终实现了激光辅助STM探针的纳米结构制造。前人的研究主要集中在STM探针针尖下方的光场增强和光限制效应以及仅在无相变情况下样品表面的温度演变。计算机模拟研究侧重于小规模和短时间内的加热过程。除此之外,在表面纳米结构制造中,表面纳米结构的最终形貌还取决于熔融材料的固化。而固化过程则由材料中的热传导来驱动,并且比加热和熔化过程慢很多。Tang等[14-16]建立了亿级原子数的大型系统,研究了基于激光辅助STM探针技术的纳米结构制造中的加热、熔化、相变、凝固和缺陷形成等长期行为和结构演变。尽管模拟结果基于少量假设,但模拟中使用的模型是按真实尺寸设计的,为相关研究人员提供了非常前沿和细致的对物理图像的认识和理解。
6.1 大尺度分子动力学并行模拟
在激光纳米加工中,人们无法对纳米尺度上的原子建立热平衡,原子的热运动在空间上表现出强烈的统计规律。分子动力学(MD)模拟可以直接跟踪原子的运动,可以研究原子位置的变化,探索原子水平的物理现象。Wang等[18-19]利用MD模拟了温度为50 K时近场聚焦激光束照射悬浮氩固体薄膜中的原子运动。入射激光束在空间和时间上的光场分布为
式中:I0是激光束强度常数;r是激光束的位置;r0是样品表面的中心位置;rg是表征激光束长度的常数;t是时间;t0是表征激光强度在时域中峰值位置的时间常数;tg是表征激光束脉冲宽度的时间常数。
氩薄膜中原子的运动可以用牛顿第二定律描述为
式中:mi和ri分别是原子i的质量和位置;Fij是原子i和j之间的相互作用力,可由公式Fij=-∂φij/∂rij得到,其中L-J势(Lennard-Jones potential)φij的表达式为
式中:ε是L-J阱深参数,氩原子的对应值为1.653×10-21 J;rij=ri-rj;σe是平衡分离参数,氩原子的对应值为0.3406 nm。Wang等在MD模拟中使用了半步蛙跳算法。对于氩原子,截止距离rc采用公认的2.5σe。
在研究大型系统的长期行为时,Wang等[18-19]将物理域划分为N个大小相同的子域,如
图 9. MD模型及材料中的激光吸收示意图[19]。(a)MD模型结构示意图;(b)材料中的激光吸收示意图
Fig. 9. Schematics of MD model and laser beam absorption in the material[19]. (a) Schematic of MD model; (b) schematic of laser beam absorption in the material
样品中激光吸收的数值处理方法如
式中:E1是在单位时间步长δt内到达单元1的激光能量,E1=∫AIδt·dA,其中A是单元1的顶表面积;δz是单元1在z方向的单位长度。由于加热、相爆炸及其他物理过程中材料的物性会发生变化,因此τ1是考虑物性变化后得到的修正光学吸收深度,该吸收深度旨在反映样品对激光束的体积吸收,而不是代表实际的实验条件。当激光到达单元2时,能量为E2=E1-δE1。单元2吸收的激光能量可以用(4)式计算得到。激光在后续单元中持续传输,直至能量被完全吸收为止。材料中原子的动能由激光传递的能量激发,而动量在激光能量吸收过程中守恒。
6.2 材料结构演变:加热、熔化、相变、凝固和缺陷的形成
基于上述大尺度分子动力学并行模拟方法,Wang[19]建立了一个由2亿多个原子组成的系统,探索了激光辅助STM探针技术进行纳米结构制造时,纳米孔加工过程中的温度、应力和结构变化,并观察到了纳米孔边缘的不规则扰动。他们研究的样品在x和y方向上包含648个面心立方单位晶格(350.83 nm),在z方向上包含60个面心立方单位晶格(32.48 nm)。单个脉冲能量设为1.6×10-14 J,相当于100 MW/cm2水平的能量密度。由于样品吸收了激光能量会产生温升,因此样品亚表面会产生巨大的压应力。样品中的温升和应力随时间的演变过程分别如
图 11. 正应力σrr (MPa)在x-y平面内随时间的演变和传播[19]
Fig. 11. Normal stress σrr(MPa) development and propagation in x-y plane[19]
图 12. x-z截面内不同时刻的氩原子位置[19]
Fig. 12. Snapshots of argon atom positions in x-z cross-section at different moments[19]
Wang等[20]在上述工作的基础上延长了MD模拟时间,进一步研究了氩晶体薄膜在近场增强脉冲激光辐照后2 ns内液相的凝固和结晶过程。在MD模拟中,激光单个脉冲的能量取8×10-15 J,相当于50 MW/cm2水平的激光功率密度。在进行激光加热前,系统先运行100 ps调整原子速度,使原子速度逐渐达到50 K时的统计平衡态。在此之后,再进行100 ps的模拟,以消除速度调整引入的干扰。同时,晶体薄膜上下方应留有自由空间,以捕获逃逸的原子。
图 13. 表面纳米结构的凝固和外延再生长[20]。(a)长时间凝固过程中,x-z截面内不同时刻氩原子的位置;(b)2 ns时x-z截面的放大图,以显示外延生长和原子位错
Fig. 13. Solidification and epitaxial regrowth in surface nanostructuring[20]. (a) Snapshots of argon atom positions in x-z cross-section during a long time solidification; (b) magnification in the x-z cross-section at t=2 ns to show the epitaxial growth and atomic dislocation
纳米孔中液相凝固过程的液相固化速率对于了解纳米结构化过程非常重要,Wang等[20]基于X射线衍射理论建立了描述结晶度的函数Φ(ri,z),该函数值可以反映不同空间位置处材料的结晶状态。该函数的计算公式为
式中:ri,z是原子i的z坐标;N是在计算域内的总原子数;λ是晶格常数a。固相的结晶度较高,Φ(ri,z)通常大于0.8,晶体结构的Φ(ri,z)为1;相较于固相,液相的结晶度低,函数Φ(ri,z)的值一般小于0.3。
图 14. 不同时刻x-z 截面上函数Φ(ri,z)的分布[20]
Fig. 14. Distribution of function Φ(ri,z) at different times in x-z cross-section[20]
激光能量密度能够强烈地影响纳米结构的直径、深度和边缘形貌。Feng等[21]通过MD模拟研究了不同能量密度的脉冲激光加热后纳米孔的凝固过程。他们采用的激光能量分别为0.05,0.1,0.2,0.35,0.5,0.75,1.0,1.1 fJ(对应的能量密度约为12,24,48,84,120,180,240,260 MW/cm2)。结果如
图 15. 激光通量对熔体凝固的影响[21]。(a)熔融物层厚和最低位置随时间的变化趋势;(b)不同激光脉冲能量下表面完全凝固后的最终形貌
Fig. 15. Effects of laser fluence on solidification of melts[21]. (a) Evolution of thickness and lowest melting position of the molten material with time; (b) final surface profiles under different laser pulse energies after complete solidification
表 1. 文献[ 21]中三种计算域的具体尺寸信息
Table 1. Detailed information about the three computational domains studied in Ref.[21]
|
6.3 激波
在基于激光辅助STM探针技术的纳米结构制造中,高能量密度的激光加热可以直接将固相材料转化为气相,并向外发射具有很高能量和运动速度的原子,形成相爆炸。当这种高速运动的原子与环境气体在的原子相遇时,会在极短的时间内压缩环境气体原子,产生强烈的激波。尽管激波在纳米制造中较少被关注,但这是一个非常重要的现象。因为纳米加工域中激波的形成会显著影响激光材料相互作用区域以及纳米结构化过程。Zhang等[22-23]首次完成了皮秒激光与材料相互作用时环境气体中激波形成和传播的一维建模工作,系统地研究了长时间相变过程,包括固液界面速度、温度、最大熔化深度和材料剥离速率。在Zhang等的模拟中,整个模型分为两个子域,其中:子域I直接受纳秒激光照射,其在y方向的长度远大于激光熔化深度;子域II仅经历了热传导,其长度大于材料中的热扩散尺寸。子域I和II分别用有限差分法和MD方法进行模拟,界面处的温度和热通量是连续的。通过模拟观察到固液界面和熔池表面存在强烈的波动,但在固液界面上没有观察到过热现象。较高的熔化峰值速度伴随较高的激光通量一起出现,材料溅射量和激光通量之间呈近似的线性关系,这表明有弱相爆炸发生。
Feng等[24]建立了激光加热样品的三维模型,研究了近场激光-材料加工过程中纳米域中的激波。在纳米域中,从固体靶材溅射出的熔融颗粒高速运动导致气体压缩,同时熔融颗粒与气体原子之间的相互作用压缩气体并形成了陡峭的激波阵面,该阵面以超音速运动。气体快速压缩的同时形成了密度、温度和压力分布陡峭的界面,如
图 16. 200 ps时,系统在厚度方向截面上的物理特性分布[24]。(a)密度(nm-3);(b)温度(K);(c)压力(MPa)
Fig. 16. Physical properties distribution in a cross-section of the system at 200 ps[24]. (a) Density (nm-3); (b) temperature (K); (c) pressure (MPa)
图 17. 激波波前温度、位置和马赫数随时间的演变[24]
Fig. 17. Evolution of shock wave front temperature, position and Mach number with time[24]
在激光加工时,样品所处环境的气体性质(包括相对分子质量和密度)也会影响激波的产生及性质[25-26]。Guo等[25]和Li等[26]建立了准三维模型并采用了并行计算,研究了气体环境中的氩晶体薄膜。Guo等[25]分别研究了8种不同环境气体下激波的形成和传播过程,研究结果如
图 18. 在不同气体环境下的不同瞬间(β=Mgas/Msolid, λ=pambient/pref),激波形成和传播过程中原子位置的演变[25](红色:靶材原子;蓝色:环境气体原子)
Fig. 18. Evolution of atom positions during shock wave formation and propagation for different cases (β=Mgas/Msolid, λ=pambient/pref) at different instants[25] (red dots: target atoms; blue dots: ambient gas atoms)
7 结束语
本文主要总结回顾了激光辅助SPM探针和基底相互作用的实验和模拟工作,包括采用拉曼热测量方法测量探针和基底的温度,以及近场聚焦效应对针尖内部和下方光场的增强。大型并行MD模拟揭示了对靶材施加强激光脉冲后其表面的长时间演变,包括材料结构演变,如加热、熔化、相变、凝固和缺陷形成。此外还探讨了激光加工中重要的激波现象,以及环境对激波的影响。
纳米结构制造是一个动态过程,并且加工区域非常小(低至纳米级),在该区域发生的物理现象很难观察和探测。在SPM探针辅助的纳米加工中,由于针尖易损,近场效应很容易改变。针尖的损坏或污染都会使尖端直径变宽,导致针尖周围的近场光场分布发生改变。此外,在探针-基底结构中,探针和基底非常接近(纳米级),可以近似认为彼此接触。由局部激光加热引起的探针和基底的任何热膨胀行为都会导致测量不稳定。目前,与纳米制造的快速动态过程相比,基于激光辅助探针技术的纳米制造中的测量需要更长的时间,在快速响应测量、规模化制造及提高扩展工业化应用等方面,该技术仍有较大的发展空间。
[7] Lu Y F, Hu B, Mai Z H, et al. Laser-scanning probe microscope based nanoprocessing of electronics materials[J]. Japanese Journal of Applied Physics, 2001, 40(6B): 4395-4398.
[13] Chen G. Nonlocal and nonequilibrium heat conduction in the vicinity of nanoparticles[J]. Journal of Heat Transfer, 1996, 118(3): 539-545.
[18] Wang XW, LawrenceC. Molecular dynamics simulation of thermal conductivity of silicon films using environment-dependent interatomic potential[C] //Proceedings of ASME 2004 Heat Transfer/Fluids Engineering Summer Conference, July 11-15, 2004, Charlotte, North Carolina, USA. New York: American Society of Mechanical Engineers2009: 619- 625.
[19] Wang X W. Large-scale molecular dynamics simulation of surface nanostructuring with a laser-assisted scanning tunnelling microscope[J]. Journal of Physics D Applied Physics, 2005, 38(11): 1805-1823.
Article Outline
徐屾, 朱红阳, 张莉君, 岳亚楠, 王信伟. 激光辅助近场纳米制造中的物理机制探索[J]. 中国激光, 2021, 48(6): 0600001. Shen Xu, Hongyang Zhu, Lijun Zhang, Yanan Yue, Xinwei Wang. Exploration of Physical Mechanism in Laser Assisted Near-field Nanomanufacturing[J]. Chinese Journal of Lasers, 2021, 48(6): 0600001.