使用聚焦激光快速加热材料常被用在在各种应用中,包括半导体加工行业。这篇博客,我们将研究具有周期性脉冲强度的高斯轮廓激光束,来加热沉积在硅衬底上的两种不同的半透明材料。为了建立此模型,我们将使用温度场和比尔-朗伯定律求解一个多物理场建模问题。接下来,让我们进一步探讨这个模型,看看如何设置它。
高斯轮廓激光束照亮硅晶圆
我们将以一个直径为 2 英寸(约 5cm)的硅晶圆为例,如下图所示,该晶圆的中心有两种不同的材料,每种材料厚度为 100μm,半径为 1cm。晶圆从顶部被一束高斯轮廓激光热源照射,该热源在时间上被快速脉冲化。这两种材料在 700nm 的激光波长下都是半透明的,但在更长波长的红外辐射下是不透明的。硅衬底是掺杂的并且在所有波长下都是高吸收性的。
脉冲激光在不透明晶圆上照亮了两层半透明材料。
由于所有材料都具有与入射光束垂直的平面边界,所有入射光都将沿平行于入射光束的均匀方向传播。材料之间的界面会有反射,但没有折射或衍射。两层材料的厚度都远大于波长,因此我们可以假设相干长度远小于层厚度。我们可以使用比尔-朗伯定律来解决这个问题,该定律描述了半透明介质中光的衰减。该方程使用 COMSOL Multiphysics® 软件中的吸收介质中的辐射束 接口求解。但是,由于存在反射,我们需要仔细研究一些细微差别。
了解物理场并设置模型
由于沉积层是圆形的,并且由于激光聚焦在中心上,我们可以忽略晶圆平面并将模型视为完全轴对称的。这使我们能够将模型简化为 2D 轴对称建模平面。在这个平面中,我们简单地绘制三个矩形来定义晶圆和两个沉积层,并为这三个矩形分配不同的材料属性。这样,几何形状和材料就定义好了,我们可以专注于物理场的研究。
首先,沿着光束路径穿过自由空间,从晶圆上方的激光源沿着 z 轴向下。我们有一个 40W、700nm 波长的激光器,并且光束具有标准偏差为 1.5mm 的高斯轮廓。激光器开启 75ms,然后关闭 25ms,或者激光器使用占空比为 75%,周期为 100ms 的脉冲加热。这种时间上的阶梯式加载是通过事件 接口解决的,该接口用于引入一个 离散状态 变量 ONOFF
,即时间为 0
或 1
。
我们不会明确地模拟激光源或通过自由空间的光束路径;我们将只对与材料相互作用的光进行建模。在顶层的边界处,折射率为 n_{top}=2.4 的材料会因为折射率的差异而有一些反射,如菲涅耳方程所示:
虽然这个方程适用于复值折射率,但在我们的计算中只考虑折射率的实值分量是合理的,因为折射率的虚部非常小。在界面上没有任何吸收的附加假设下(例如由于吸收材料的非常薄的涂层),透射率为 T=1-R。这样就完成了我们在吸收介质中辐射束 接口设置入射强度 功能所需的信息,如下面的屏幕截图所示。
入射强度功能的设置。
当光束穿过材料的第一层时,其强度与吸收系数成比例减小,吸收系数 \kappa 由下式确定:
其中,k 是折射率的虚部,\lambda_0 是自由空间激光波长。吸收系数可能与温度有关,但我们将从它是一个常数开始。给定光束轮廓在顶面上的强度分布,剩下整个域的光束强度通过计算获得。
在沉积材料的顶层和底层之间的介电界面,将再次存在菲涅耳方程描述的反射和透射。光束的反射分量使用已有的 吸收介质中的辐射光束 接口进行处理,只需添加第二个入射强度 功能就可以了。可以向这个界面添加任意数量的入射强度 功能;每个都将引入一个额外的变量来求解,这些变量将被命名 rbam.I1, rbam.I2, …,
依此类推。在第二个入射强度 功能中,我们可以引入基于第一个光束强度和菲涅耳反射系数的用户定义的光束轮廓。通过改变光束方向的符号,可以完全考虑光在此接口上的部分反射,如下面的屏幕截图所示。从理论上讲,在顶部边界会有一个额外的光束反射,但是这个二次反射足够小,所以我们将忽略它。
第二个 入射强度功能的屏幕截图,考虑了介电界面处的反射。
接下来,我们跟随光束穿过介电界面进入第二层半透明材料。由于跨越此边界的光强度发生了变化,因此必须添加第二个 吸收介质中的辐射束 接口,并根据菲涅耳透射率和来自第一个吸收介质中的辐射束 接口的第一束光束来定义入射强度。
吸收介质中的辐射束接口中第二个辐射束的 入射强度特征的屏幕截图,用于底部域中的强度。
最后,让我们讨论当光线到达第二层底部并击中硅晶圆衬底时会发生什么。我们将假设硅晶片是掺杂的,因此它具有高度吸收性和非反射性。由于所有到达这个边界的光都将在足够小的距离内被吸收,因此可以说光在边界处被吸收。对于这种情况,不透明表面 边界条件将在所选边界处沉积所有能量,这就完成了激光在结构中传播时的建模。通过这些功能的组合,我们已经完全模拟了入射激光束穿过模型。现在我们可以将注意力转向热模型。
模拟温度随时间的变化
晶圆最初处于 300K 的均匀温度。所有域都有传导传热,我们假设材料之间的界面没有明显的热阻,即材料界面之间没有温差,磁通量是连续的。这种情况是软件的默认假设,但如果我们确实希望覆盖它,可以添加薄层 或热接触 功能。
在 100μm 处,层厚已经足够适用经典的傅里叶传热定律,值得一提的是,纳米级传热是 COMSOL 用户研究的一个活跃领域;例如,请参阅我们的特邀博客“动力学集体模型中的流体动力热输送”。
至于热边界条件,我们将假设晶圆位于完全绝缘的底座上,并且位于一个近真空的工艺室内。这意味着不会有传导或对流传热冷却,但会向腔室壁进行辐射传热,假设保持在 300K。接着,我们假设晶圆温度只会上升几百K,因此与入射激光相比,辐射发射将处于一个更长的波长带。这意味着,从概念上讲,我们可以使用双波段模型进行辐射传热。来自激光的入射辐射已经通过吸收介质中的辐射束 接口完全处理。较长波段的发射辐射(由于晶圆相对于工艺室壁的温度升高)可以使用单波段表面对表面辐射 接口与固体传热 接口进行建模。表面到表面辐射 接口计算所有暴露表面与周围空间之间的角系数。
值得一提的是,在这种情况下,只有在晶圆上方的小内角附近存在表面对表面辐射;其他地方对环境的角系数都是统一的。如果我们想稍微简化一下,可以不使用表面对表面辐射 接口,而是在固体传热 接口中使用表面到环境辐射 边界条件。计算时间和结果的差异可以忽略不计,因此这里我们使用更准确的方法,即使用表面对表面辐射 接口计算角系数。
我们还需要特别注意这个装置的网格划分。吸收介质中的辐射束 接口求解的是一阶偏微分方程,默认情况下使用场的线性离散化。根据吸收系数,我们知道强度会随着两层的厚度发生明显变化。我们还知道,激光束轮廓在表面上的强度变化是相当渐进的。这证明了层内具有高纵横比矩形单元的映射网格是合理的。当然,随着建模复杂性的提高,我们总是希望研究网格和求解器的相对公差细化,就像我们之前的博客文章“在 COMSOL Multiphysics® 中模拟固体瞬态加热简介“所讨论的那样。
设置完成后,我们将使用随时间变化的求解器解决这个问题,并按照求解器采取的步骤保存数据。然后,我们可以绘制出温度曲线和吸收的热量,以及一段时间内中上部点的温度,如下图所示。
沿 z 轴的高度与温度的关系。
最后,为了说明,我们将介绍一种非线性材料,使底层的吸收系数随着温度的升高而上升。两种半透明材料的吸收系数比较如下图所示。随着非线性吸收系数的升高,材料的加热更大。由于这种材料的非线性,我们还需要细化具有非线性属性的层中的网格。
使用两种不同的材料模型比较温度随时间的变化。
结束语
我们介绍了一种解决半透明材料的加热问题的建模方法。准直辐射热源(激光)通过一组吸收介质中的辐射束 接口进行建模,该接口可以处理材料在激光波长下的半透明性质以及介电界面处的反射。脉冲热源通过事件 接口处理,较长波长的红外再辐射通过表面对表面辐射 接口处理。这种建模方法适用于半导体加工领域或准直光入射到半透明材料上的任何情况。
如果你对这些类型的模拟感兴趣,请随时单击下面的按钮下载文中讨论的示例模型:
评论 (13)
治 乔
2023-09-04Dear author, is it possible to solve nanosecond pulsed laser with this method?
Xiaohan Jiang
2023-09-05 COMSOL 员工It should be OK to simulate nanosecond laser by this approach. But the computation will be much larger.
Yang Liu
2023-11-27激光的波长怎么设置呢
Qingbin Yuan
2023-11-30 COMSOL 员工在这种方法中,激光波长的影响主要体现在吸收系数的变化上。因此可以将吸收系数定义为对应波长的数值(可以是单波长下的定值,也可以是宽光谱下的函数),来体现波长的影响。
鸿 刘
2023-12-11标准差1.5mm代表激光光束什么?半径为1.5mm?
Qihang Lin
2023-12-14 COMSOL 员工建议搜索高斯函数在几何上的分布进行理解
如天 乐
2024-03-28激光占空比具体怎么设置
hao huang
2024-04-02 COMSOL 员工您好
占空比是通过事件接口进行设置的,请参考博客文末所附的模型中的设置方法。
如果有进一步问题,建议您联系 COMSOL 的技术支持团队:
在线支持中心:cn.comsol.com/support
Email: support@comsol.com
谢谢!
泽宇 顾
2024-10-23Dear author,I would like to ask how the formula for the nonlinear absorption coefficient of semi transparent materials with temperature variation is obtained in your model. I really need your help.
子奇 陈
2024-10-25 COMSOL 员工In the case model above, the formula used for the nonlinear absorption coefficient with temperature is a linear variation model, meaning that the absorption rate varies linearly and linearly with temperature, which in turn leads to a nonlinear phenomenon in the variation of the absorption rate and the final heat transfer field. The form is α(T) = α0(1 + α(T-T0)), where α0 is the reference absorption coefficient at temperature T0 and α is the temperature dependence coefficient to be determined, which can be obtained either from the literature or experimentally. The purpose of using this formula in the case is only to illustrate the ability and results of calculations for nonlinearities and does not represent real values.
志键 黄
2024-11-15尊敬的作者/工程师,你好!我最近在研究分析激光诱导光学材料损伤,需要分析纳秒级脉冲激光作用后光学材料的温度场和应力场分布。对于光学材料,基本上对激光都有高的透射率,低反射率,我通过相关文献了解到了激光对金属材料加工的体热源表达式为Q=(1-R)*I(r,t)*alpha*exp(-alpha*z)*g(t),其中R代表材料对激光的反射率,I(r,t)为激光的功率密度(峰值功率密度)alpha为材料对激光的吸收系数;基于上式,我对激光加热光学材料的体热源表达式做了如下变化,Q=(1-R-T)*I(r,t)*alpha*exp(-alpha*z)*g(t),其中T为材料对激光的透射率,基于此热源表达式,我尝试着复现相关文献,但我一直无法获得正确的温升,在单脉冲后温度可能仅上升了一点点,与参考文献相差甚大,我想请问我这样对体热源表达式做出变化是否合理?或者说有其他的方式来模拟对于激光加热透明电介质(透明光学窗口)的接口?这个问题已经困恼我许久,这是关乎我毕业的一大关卡,希望能得到您的回复!谢谢!
Yuqing Ge
2024-11-22 COMSOL 员工您好,具体温升的仿真值涉及理论和模型搭建的细节,不好判断,初步推测您修改的热源Q中(1-R-T)中考虑了吸收的系数,alpha再次考虑了吸收的系数,有重复考虑的可能。
Jianxing Xiao
2024-12-06你好,我最近也在做脉冲激光对光学介质膜结构的温升模型,我看到一篇文献中对介质膜中的热源是有定义的,它与其中的归一化电场相关,(参考文献:DOI 10.1007/s00231-017-2147-8)但是我也没法复现这篇论文。如果可以的话,我们可以加个联系方式交流一下吗?谢谢!