[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN112114501A - 一种计算光刻投影物镜中朗奇剪切干涉图像的方法 - Google Patents

一种计算光刻投影物镜中朗奇剪切干涉图像的方法 Download PDF

Info

Publication number
CN112114501A
CN112114501A CN202011148935.6A CN202011148935A CN112114501A CN 112114501 A CN112114501 A CN 112114501A CN 202011148935 A CN202011148935 A CN 202011148935A CN 112114501 A CN112114501 A CN 112114501A
Authority
CN
China
Prior art keywords
electric field
expression
calculating
light
function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011148935.6A
Other languages
English (en)
Other versions
CN112114501B (zh
Inventor
牛志元
施伟杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dongfang Jingyuan Electron Ltd
Original Assignee
Shenzhen Branch Dongfang Jingyuan Microelectronics Tech Beijing Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Branch Dongfang Jingyuan Microelectronics Tech Beijing Co ltd filed Critical Shenzhen Branch Dongfang Jingyuan Microelectronics Tech Beijing Co ltd
Priority to CN202011148935.6A priority Critical patent/CN112114501B/zh
Publication of CN112114501A publication Critical patent/CN112114501A/zh
Priority to PCT/CN2021/125704 priority patent/WO2022083735A1/zh
Application granted granted Critical
Publication of CN112114501B publication Critical patent/CN112114501B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03F—PHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00—Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70—Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483—Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • G03F7/70591—Testing optical components
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
    • G01J9/0215—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods by shearing interferometric methods
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03F—PHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00—Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70—Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483—Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • G03F7/70605—Workpiece metrology
    • G03F7/70616—Monitoring the printed patterns
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)
  • Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
  • Exposure Of Semiconductors, Excluding Electron Or Ion Beam Exposure (AREA)

Abstract

本发明提供了一种计算光刻投影物镜中朗奇剪切干涉图像的方法，在该方法中，光源发出的光从物方光栅到像方光栅的过程中采用严格的光刻矢量模型进行计算，同时可以适用于任意形状的物方、像方光栅；并且光从像方光栅到接收屏的过程中采用严格的Rayleigh‑Sommerfeld衍射积分进行计算。本发明提供的计算光刻投影物镜中朗奇剪切干涉图像的方法由于采用严格的光刻矢量模型以及Rayleigh‑Sommerfeld衍射积分进行计算，所以无需考虑物方、像方光栅的近场或远场，接收屏也无需考虑傍轴近似等因素，避免了传统测试方法局限性大导致难以精确的得到剪切干涉图像的问题。

Description

一种计算光刻投影物镜中朗奇剪切干涉图像的方法
【技术领域】
本发明涉及光学测量技术领域，特别涉及一种计算光刻投影物镜中朗奇剪切干涉图像的方法。
【背景技术】
高精度、高速度提取投影物镜波像差的检测系统是高端光刻机的核心分系统。波像差是投影物镜成像质量的主要参数，直接决定光刻机的分辨率和套刻精度，高端光刻机投影物镜波像差已达到1nm以下，而随着产率提高，透镜的热像差更加显著。
在现有的基于朗奇剪切干涉图像的检测方法中，通常假定物方光栅是理想的周期线条光栅(即具有无穷多线条的光栅，并且线条长度为无穷长)，该物方光栅为一维图形，而不是复杂的二维图形，具有局限性；并且传统检测方法的计算过程中假定成像模型是标量模型，而非严格的矢量模型，并且传统检测方法中的接收屏干涉图像计算默认采用远场近似技术(即夫琅禾费衍射区域)，然而这些近似条件在实际应用当中都具有一定局限性，难以精确的得到剪切干涉图像。
【发明内容】
为了克服目前现有的基于朗奇剪切干涉图像的检测方法中计算方法局限性大导致难以精确的得到剪切干涉图像的问题，本发明提供一种计算光刻投影物镜中朗奇剪切干涉图像的方法。
本发明为解决上述技术问题，提供的技术方案如下：本发明提供了一种计算光刻投影物镜中朗奇剪切干涉图像的方法，用于计算并得出被测投影物镜的波像差，该方法包括以下步骤，步骤S1：通过计算得出射入物方光栅的光的标量电场函数及物方光栅的物方光栅频谱以得出被物方光栅衍射后的光的电场频谱空间分布函数；步骤S2：光通过被测投影物镜，计算得出被测投影物镜的光瞳函数并结合所述电场频谱空间分布函数算得光的出瞳处电场分布表达式；步骤S3：光照射在像方光栅上表面时得出其偏振信息，并通过几何旋转矩阵R3×2将出瞳处电场分布表达式从局部坐标系的sp分量转化为全局坐标系的xyz分量以计算得出入射光的像方光栅上表面电场表达式；步骤S4：计算得出RS积分核函数，并将像方光栅上表面电场表达式进行计算得出光的像方光栅下表面电场表达式，基于像方光栅下表面电场表达式并结合RS积分核函数进行计算以得到多个光在接收屏上形成的剪切干涉图形；步骤S5：定义IMM(fx,fy；f′x,f′y)变量矩阵并对其进行本征值分解以重写多个剪切干涉图形并得到分解结果，通过重复计算所述分解结果以得出光在接收屏上不同剪切位置的精确剪切干涉图像。
优选地，照射到物方光栅上的光为光源输出的多方向的入射光，所述步骤S1进一步包括以下步骤；步骤S11：通过计算得出照射在物方光栅上的光的光源频率(sx,sy)及空间频率(fx,fy)并得出射入平行光的标量电场函数为U(x,y)；步骤S12：测得物方光栅透视率函数为m1(x′,y′)并对其进行傅里叶变换以得出物方光栅频谱为M1(fx,fy)；步骤S13：进一步算出经过物方光栅衍射后的光的电场频谱空间分布函数为M1(fx-sx,fy-sy)，其表达式为：
Figure BDA0002740555090000022
在步骤S11中光源频率(sx,sy)用于描述光的入射方向，标量电场函数U(x,y)的表达式为：
Figure BDA0002740555090000021
优选地，在所述步骤S2中，所述光瞳函数为P(fx,fy)，其函数表达式为：
Figure BDA0002740555090000031
计算得出光在出瞳处成像时的辐射修正因子，其表达式为：
Figure BDA0002740555090000032
将所述电场频谱空间分布函数、所述光瞳函数及所述辐射修正因子依次相乘以得到所述出瞳处电场分布表达式：
Figure BDA0002740555090000033
优选地，所述步骤S3包括以下步骤：步骤S31：所述几何旋转矩阵R3×2为3x2维的几何旋转矩阵，并通过琼斯矢量将所述出曈处电场分布表达式扩展为矢量形式以得到出曈处电场分布矢量表达式；步骤S32：将所述出曈处电场分布矢量表达式的局部坐标系的sp分量转化为全局坐标系的xyz分量以计算得出所述像方光栅上表面电场表达式；在步骤S31中，对于x偏振所述琼斯矢量为
Figure BDA0002740555090000034
对于y偏振所述琼斯矢量为
Figure BDA0002740555090000035
利用所述几何旋转矩阵R3×2得出所述出瞳处电场分布表达式的矢量形式表达式为：
Figure BDA0002740555090000036
在步骤S32中，所述像方光栅上表面电场表达式为：
Figure BDA0002740555090000037
Figure BDA0002740555090000038
优选地，所述像方光栅下表面电场表达式用于描述光在经过像方光栅后的电场分布，所述像方光栅下表面电场表达式为：
Figure BDA0002740555090000039
优选地，所述步骤S4包括以下步骤：
步骤S41：基于所述像方光栅下表面电场表达式并结合RS积分核函数经过衍射积分计算得出接收屏电场分布函数表达式，所述接收屏电场分布函数表达式为：
Figure BDA00027405550900000310
Figure BDA00027405550900000311
步骤S42：将所述接收屏电场分布函数表达式中的(ε,η)自由度通过积分消除后进行快速傅里叶计算以得到光在接收屏上形成电场的标量电场函数；所述标量电场函数的表达式为:
Figure BDA0002740555090000041
Figure BDA0002740555090000042
优选地，所述RS积分核函数为Rayleigh-Sommerfeld积分核函数，所述RS积分核函数的表达式为：
Figure BDA0002740555090000043
在所述RS积分核函数的表达式中，r的函数表达式为：
Figure BDA0002740555090000044
优选地，所述步骤S4进一步包括以下步骤：步骤S43：通过计算得到用于描述光源在频率空间中强度分布的光源函数I(sx,sy)，将所述标量电场函数进行非相干叠加以得到多个入射方向的光在接收屏上形成的剪切干涉图形，其叠加权重为光源函数I(sx,sy)；
所述剪切干涉图形的表达式为：
Figure BDA0002740555090000045
优选地，所述步骤5进一步包括以下步骤：步骤S51：将IMM(fx,fy；f′x,f′y)变量矩阵本征值分解得到的所述分解结果储存以供多次计算使用，所述分解结果包括通过本征值分解得到的本征值及本征向量。
优选地，在所述步骤S5中，所述IMM(fx,fy；f′x,f′y)变量矩阵为厄密的矩阵且本征值为实数，其函数表达式为：IMM(fx,fy；f′x,f′y)＝∫dsxdsyI(sx,sy)M(fx-sx,fy-sy)M*(f′x-sx,f′y-sy)；对所述IMM(fx,fy；f′x,f′y)变量矩阵进行本征值分解可得：
Figure BDA0002740555090000046
Figure BDA0002740555090000047
以此将所述剪切干涉图形的表达式重写以得到所述精确剪切干涉图像的表达式：
Figure BDA0002740555090000051
与现有技术相比，本发明提供的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，具有以下优点：
(1)在本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法中，由于采用严格的光刻矢量模型以及Rayleigh-Sommerfeld衍射积分对电场分布(即光强分布)进行计算，所以无需考虑物方光栅及像方光栅的近场或远场，接收屏也无需考虑傍轴近似等因素，避免了传统测试方法局限性大导致难以精确的得到剪切干涉图像的问题。
(2)在本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法中，在计算光在像方光栅上的电场分布时，通过定义一几何旋转矩阵R3×2以表示把出曈处电场分布表达式的局部坐标系的sp分量转化为全局坐标系的xyz分量，并进一步通过琼斯矢量把出曈处电场分布表达式扩展为矢量形式以得到出曈处电场分布的矢量表达式，技术人员最终可通过该矢量表达式计算得出像方光栅上表面电场表达式；可以理解的是，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法在计算像方光栅上表面电场表达式的过程中采用了严格的矢量模型，因为光在通过投影物镜时，用sp分量对计算比较方便，而在光通过投影物镜后并照射到像方光栅的过程中就需要把电场转化为全局坐标系的xyz分量以方便计算，最终通过计算得出像方光栅上表面电场表达式，这种计算方法使得计算流程得到了简化，技术人员在进行计算时更加方便，有效地提高了技术人员的工作效率，避免了计算效率低导致费时费力的问题。
(3)本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法通过将像方光栅下表面电场表达式结合RS积分核函数并经过衍射积分计算以最终得出用于表达光在接收屏上电场分布的接收屏电场分布函数表达式；可以理解的是，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法在计算光从像方光栅到接收屏这一过程中的电场分布时，采用了严格的Rayleigh-Sommerfeld衍射积分对接收屏电场分布函数表达式进行计算，使得最终得到的光的电场分布结果更为准确，有利于最终得到精确的剪切干涉图像，解决了传统检测方法难以精确地得到剪切干涉图像的问题。
(4)在本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法中，通过在计算光在接收屏上形成的剪切干涉图像的过程中定义一IMM(fx,fy；f′x,f′y)变量矩阵以用于分解不同剪切位置的剪切干涉图像，且IMM(fx,fy；f′x,f′y)变量矩阵为光源函数I(sx,sy)和物方光栅频谱的函数，技术人员可通过IMM(fx,fy；f′x,f′y)变量矩阵研究不同波像差(即不同的有效光瞳函数P(fx,fy))、不同像方光栅及不同剪切位置的剪切干涉图像；可以理解的是，精确剪切干涉图像主要用于反映接收屏上
(5)光的光强分布，技术人员可通过精确剪切干涉图像进行反推计算以得出待测投影物镜的波像差，由于精确剪切干涉图像是通过采用严格的矢量模型及严格的衍射积分计算所得，使其精确度得到了提高，使得最终得出的波像差数据精确可靠，避免了传统检测方法检测精确度不足导致检测结果不可靠的问题；进一步地，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法通过一次性分解IMM(fx,fy；f′x,f′y)变量矩阵即可得到多个精确剪切干涉图像，技术人员通过将分解结果存入磁盘中，以方便技术人员在后续的计算中直接使用该分解结果，有效提高了技术人员的工作效率，实现了一次分解多次使用的功能。
【附图说明】
图1是本发明第一实施例提供的一种计算光刻投影物镜中朗奇剪切干涉图像的方法的流程示意框图；
图2是本发明第一实施例提供的一种计算光刻投影物镜中朗奇剪切干涉图像的方法之朗奇剪切干涉原理示意图；
图3是本发明第二实施例提供的一种具体的检测流程之物方光栅图形；
图4是本发明第二实施例提供的一种具体的检测流程之物方光栅空间成像；
图5是本发明第二实施例提供的一种具体的检测流程之剪切干涉图像。
附图标记说明：
1、物方光栅；2、被测投影物镜；3、像方光栅；4、接收屏；z、接收屏距离像方光栅的距离。
【具体实施方式】
为了使本发明的目的，技术方案及优点更加清楚明白，以下结合附图及实施实例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
需要说明的是，当元件被称为“固定于”另一个元件，它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件，它可以是直接连接到另一个元件或者可能同时存在居中元件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的。
请结合参阅图1-图2，本发明第一实施例提供了一种计算光刻投影物镜中朗奇剪切干涉图像的方法，用于计算并得出被测投影物镜2的波像差，该方法包括以下步骤，步骤S1：通过计算得出射入物方光栅1的光的标量电场函数及物方光栅1的物方光栅频谱以得出被物方光栅1衍射后的光的电场频谱空间分布函数；步骤S2：光通过被测投影物镜2，计算得出被测投影物镜2的光瞳函数并结合电场频谱空间分布函数算得光的出瞳处电场分布表达式；步骤S3：光照射在像方光栅上表面时(像方光栅上表面即光首先接触像方光栅3的一面)得出其偏振信息，并通过几何旋转矩阵R3×2将出瞳处电场分布表达式从局部坐标系的sp分量转化为全局坐标系的xyz分量以计算得出入射光的像方光栅上表面电场表达式(即出瞳处电场分布表达式扩展为矢量形式)；步骤S4：计算得出RS积分核函数，并将像方光栅上表面电场表达式进行计算得出光的像方光栅下表面电场表达式(像方光栅下表面即光离开像方光栅3时最后接触的一面)，基于像方光栅下表面电场表达式并结合RS积分核函数进行计算以得到多个光在接收屏4上形成的剪切干涉图形；步骤S5：定义IMM(fx,fy；f′x,f′y)变量矩阵并对其进行本征值分解以重写多个剪切干涉图形并得到分解结果，通过重复计算分解结果以得出光在接收屏4上不同剪切位置的精确剪切干涉图像。
具体地，照射到物方光栅1上的光为光源输出的多个方向的入射光，其中相同方向的光为相干光，而不同方向的光是非相干光；步骤S1进一步包括以下步骤；步骤S11：通过计算得出照射在物方光栅1上的光的光源频率(sx,sy)及空间频率(fx,fy)并得出射入平行光的标量电场函数为U(x,y)；步骤S12：测得物方光栅透视率函数为m1(x′,y′)并对其进行傅里叶变换以得出物方光栅频谱为M1(fx,fy)；步骤S13：进一步算出经过物方光栅1衍射后的光的电场频谱空间分布函数为M1(fx-sx,fy-sy)，其表达式为：
Figure BDA0002740555090000081
在步骤S11中光源频率(sx,sy)用于描述光的入射方向，标量电场函数U(x,y)的表达式为：
Figure BDA0002740555090000091
光源频率(sx,sy)用于描述光的方向。
请继续参阅图2，不同方向的平行光(不同方向的光是非相干的，即非相干光，如图2中箭头所示方向射入)照射到物方光栅1上，光被物方光栅1衍射后，进入被测投影物镜2；光在通过被测投影物镜2过程中，携带了投影物镜的波像差信息，在其通过被测投影物镜2后，被位于像平面的像方光栅3再次衍射，最终成像在远处的接收屏4上；当像方光栅3移动不同的距离，即接收屏4距离像方光栅3的距离z发生变化(即移相)时，接收屏4上将得到不同的剪切干涉图像。在该过程中，可以看出，接收屏4上的剪切干涉图像包含了被测投影物镜2的像差信息，技术人员通过综合不同移相的剪切干涉图像就可以提取出被测投影物镜2的波像差。
具体地，在步骤S2中，光瞳函数为P(fx,fy)，其函数表达式为：
Figure BDA0002740555090000092
计算得出光在出瞳处成像时的辐射修正因子，其表达式为：
Figure BDA0002740555090000093
将电场频谱空间分布函数、光瞳函数及辐射修正因子依次相乘以得到出瞳处电场分布表达式：
Figure BDA0002740555090000094
具体地，步骤S3包括以下步骤：步骤S31：几何旋转矩阵R3×2为3x2维的几何旋转矩阵，并通过琼斯矢量将出曈处电场分布表达式扩展为矢量形式以得到出曈处电场分布矢量表达式；步骤S32：将出曈处电场分布矢量表达式的局部坐标系的sp分量转化为全局坐标系的xyz分量以计算得出像方光栅上表面电场表达式；在步骤S31中，对于x偏振琼斯矢量为
Figure BDA0002740555090000095
对于y偏振琼斯矢量为
Figure BDA0002740555090000096
利用几何旋转矩阵R3×2得出出瞳处电场分布表达式的矢量形式表达式为：
Figure BDA0002740555090000101
在步骤S32中，像方光栅上表面电场表达式为：
Figure BDA0002740555090000102
Figure BDA0002740555090000103
进一步地，几何旋转矩阵R3×2表示把局部坐标系的sp分量转化为全局坐标系的xyz分量，光在通过被测投影物镜2时用sp分量对其进行计算比较方便，而在其通过投影物镜后，为了计算方便则需要把光的电场转化为全局坐标系的xyz分量，最终得到像方光栅上表面电场表达式。
可以理解的是，在计算光在像方光栅3上的电场分布时，通过定义一几何旋转矩阵R3×2以表示把出曈处电场分布表达式的局部坐标系的sp分量转化为全局坐标系的xyz分量，并进一步通过琼斯矢量把出曈处电场分布表达式扩展为矢量形式以得到出曈处电场分布的矢量表达式，技术人员最终可通过该矢量表达式计算得出像方光栅上表面电场表达式；可以理解的是，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法在计算像方光栅上表面电场表达式的过程中采用了严格的矢量模型，因为光在通过投影物镜时，用sp分量对计算比较方便，而在光通过投影物镜后并照射到像方光栅3的过程中就需要把电场转化为全局坐标系的xyz分量以方便计算，最终通过计算得出像方光栅上表面电场表达式，这种计算方法使得计算流程得到了简化，技术人员在进行计算时更加方便，有效地提高了技术人员的工作效率，避免了计算效率低导致费时费力的问题。
具体地，像方光栅下表面电场表达式用于描述光在经过像方光栅3后的电场分布，像方光栅下表面电场表达式为：
Figure BDA0002740555090000104
具体地，步骤S4包括以下步骤：
步骤S41：基于像方光栅下表面电场表达式并结合RS积分核函数经过衍射积分计算得出接收屏电场分布函数表达式，接收屏电场分布函数表达式为：
Figure BDA0002740555090000111
Figure BDA0002740555090000112
步骤S42：将接收屏电场分布函数表达式中的(ε,η)自由度通过积分消除后进行快速傅里叶计算以得到多个光在接收屏4上形成电场的标量电场函数；
标量电场函数的表达式为:
Figure BDA0002740555090000113
具体地，RS积分核函数为Rayleigh-Sommerfeld积分核函数，RS积分核函数的表达式为：
Figure BDA0002740555090000114
在RS积分核函数的表达式中，r的函数表达式为：
Figure BDA0002740555090000115
可以理解的是，由于本发明中提供的计算光刻投影物镜中郎奇剪切干涉图像的方法采用的是严格的光刻矢量模型以及Rayleigh-Sommerfeld衍射积分对电场分布(即光强分布)进行计算，所以无需考虑物方光栅1及像方光栅3的近场或远场，接收屏4也无需考虑傍轴近似等因素，避免了传统测试方法局限性大导致难以精确的得到剪切干涉图像的问题。
进一步地，通过将像方光栅下表面电场表达式结合RS积分核函数并经过衍射积分计算，以得出用于表达光在接收屏4上电场分布的接收屏电场分布函数表达式，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法在计算光从像方光栅3到接收屏4这一过程中的电场分布时，采用了严格的Rayleigh-Sommerfeld衍射积分对接收屏电场分布函数表达式进行计算，使得最终得到的光的电场分布结果更为准确，有利于最终得到精确的剪切干涉图像，解决了传统检测方法难以精确地得到剪切干涉图像的问题。
具体地，步骤S4进一步包括以下步骤：步骤S43：通过计算得到用于描述光源在频率空间中强度分布的光源函数I(sx,sy)，将标量电场函数进行非相干叠加以得到多个入射方向的光在接收屏4上形成的剪切干涉图形，其叠加权重为光源函数I(sx,sy)；
剪切干涉图形的表达式为：
Figure BDA0002740555090000121
具体地，步骤5进一步包括以下步骤：步骤S51：将IMM(fx,fy；f′x,f′y)变量矩阵本征值分解得到的分解结果储存以供多次使用，分解结果包括通过本征值分解得到的本征值及本征向量；可以理解的是，本发明中提供的计算光刻投影物镜中朗奇剪切干涉图像的方法通过一次性分解IMM(fx,fy；f′x,f′y)变量矩阵即可得到多个精确剪切干涉图像，技术人员可将本征值分解得出的分解结果存入磁盘中，以方便技术人员在后续的计算中直接使用该分解结果，有效提高了技术人员的工作效率，实现了一次分解多次使用的功能。
具体地，在步骤S5中，IMM(fx,fy；f′x,f′y)变量矩阵为厄密的矩阵且本征值为实数，其函数表达式为：IMM(fx,fy；f′x,f′y)＝∫dsxdsyI(sx,sy)M(fx-sx,fy-sy)M*(f′x-sx,f′y-sy)；容易证明IMM(fx,fy；f′x,f′y)＝IMM*(f′x,f′y；fx,fy)，所以可得IMM(fx,fy；f′x,f′y)变量矩阵为厄密的，且其本征值为实数；
对IMM(fx,fy；f′x,f′y)变量矩阵进行本征值分解可得到其本征值以及本征值向量，IMM(fx,fy；f′x,f′y)变量矩阵进行本征值分解的表达公式为：
Figure BDA0002740555090000122
技术人员可根据上述公式将剪切干涉图形的表达式重写以得到精确剪切干涉图像的表达式为：
Figure BDA0002740555090000131
进一步地，IMM(fx,fy；f′x,f′y)变量矩阵为光源函数I(sx,sy)和物方光栅频谱的函数，并用于研究不同波像差(即不同的有效光瞳函数P(fx,fy))、不同像方光栅3(即m2(ε,η)参数不同)及不同剪切位置的剪切干涉图像；
可以理解的是，精确剪切干涉图像主要用于反映接收屏4上光的光强分布，技术人员可通过精确剪切干涉图像进行反推计算以得出待测投影物镜的波像差；由于精确剪切干涉图像是通过采用严格的矢量模型及严格的衍射积分计算所得，使其精确度得到了提高，使得最终得出的波像差数据精确可靠，避免了传统检测方法检测精确度不足导致检测结果不可靠的问题。
请结合参阅图3-图5，本发明第二实施例提供一种具体的检测流程，该检测流程以本发明第一实施例中提供的计算光刻投影物镜中朗奇剪切干涉图像的方法为计算方法，以将其代入到具体构建的应用场景中以进一步说明本发明内容，但并不用于限定本发明，其具体实施步骤如下所示：
步骤1：对光源输入照明模式，使其照明模式为传统照明，光源在该照明模式下发出的光的相干因子为1，且其波长为632nm、偏振方向为X；
步骤S2：输入被测投影物镜2的数值孔径NA＝0.4，并输入被测投影物镜2的波像差Z9＝632nm；
步骤3：将物方光栅1参数输入，该物方光栅1的GDS图形为4线条光栅(如图3所示)，且该物方光栅1的具体参数为：物方光栅1的周期为25um，物方光栅1的线条宽度为12.5um，物方光栅1的长度为100um；输入的物方光栅1的物方掩模类型为二元掩模，并具有暗场；
步骤4：输入像方光栅3信息，这里设置为同物方掩模完全相同。
步骤5：输入像方光栅3的移相，其移相依次为0um和6.25um(即1/4周期)；
步骤6：输入接收屏距离像方光栅3的距离z为5000um，计算图像的大小为5000um*5000um，且该接收屏4的采样像素为256*256像素；
步骤7：以本发明第一实施例中提供的计算光刻投影物镜中朗奇剪切干涉图像的方法进行计算任务，当计算任务完成后，可以得到物方光栅1的空间像(如图4所示)以及剪切干涉图像的结果(如图5所示)。
与现有技术相比，本发明提供的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，具有以下优点：
(1)在本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法中，由于采用严格的光刻矢量模型以及Rayleigh-Sommerfeld衍射积分对电场分布(即光强分布)进行计算，所以无需考虑物方光栅及像方光栅的近场或远场，接收屏也无需考虑傍轴近似等因素，避免了传统测试方法局限性大导致难以精确的得到剪切干涉图像的问题。
(2)在本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法中，在计算光在像方光栅上的电场分布时，通过定义一几何旋转矩阵R3×2以表示把出曈处电场分布表达式的局部坐标系的sp分量转化为全局坐标系的xyz分量，并进一步通过琼斯矢量把出曈处电场分布表达式扩展为矢量形式以得到出曈处电场分布的矢量表达式，技术人员最终可通过该矢量表达式计算得出像方光栅上表面电场表达式；可以理解的是，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法在计算像方光栅上表面电场表达式的过程中采用了严格的矢量模型，因为光在通过投影物镜时，用sp分量对计算比较方便，而在光通过投影物镜后并照射到像方光栅的过程中就需要把电场转化为全局坐标系的xyz分量以方便计算，最终通过计算得出像方光栅上表面电场表达式，这种计算方法使得计算流程得到了简化，技术人员在进行计算时更加方便，有效地提高了技术人员的工作效率，避免了计算效率低导致费时费力的问题。
(3)本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法通过将像方光栅下表面电场表达式结合RS积分核函数并经过衍射积分计算以最终得出用于表达光在接收屏上电场分布的接收屏电场分布函数表达式；可以理解的是，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法在计算光从像方光栅到接收屏这一过程中的电场分布时，采用了严格的Rayleigh-Sommerfeld衍射积分对接收屏电场分布函数表达式进行计算，使得最终得到的光的电场分布结果更为准确，有利于最终得到精确的剪切干涉图像，解决了传统检测方法难以精确地得到剪切干涉图像的问题。
(4)在本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法中，通过在计算光在接收屏上形成的剪切干涉图像的过程中定义一IMM(fx,fy；f′x,f′y)变量矩阵以用于分解不同剪切位置的剪切干涉图像，且IMM(fx,fy；f′x,f′y)变量矩阵为光源函数I(sx,sy)和物方光栅频谱的函数，技术人员可通过IMM(fx,fy；f′x,f′y)变量矩阵研究不同波像差(即不同的有效光瞳函数P(fx,fy))、不同像方光栅及不同剪切位置的剪切干涉图像；可以理解的是，精确剪切干涉图像主要用于反映接收屏上光的光强分布，技术人员可通过精确剪切干涉图像进行反推计算以得出待测投影物镜的波像差，由于精确剪切干涉图像是通过采用严格的矢量模型及严格的衍射积分计算所得，使其精确度得到了提高，使得最终得出的波像差数据精确可靠，避免了传统检测方法检测精确度不足导致检测结果不可靠的问题；进一步地，本发明涉及的计算光刻投影物镜中朗奇剪切干涉图像的方法通过一次性分解IMM(fx,fy；f′x,f′y)变量矩阵即可得到多个精确剪切干涉图像，技术人员通过将分解结果存入磁盘中，以方便技术人员在后续的计算中直接使用该分解结果，有效提高了技术人员的工作效率，实现了一次分解多次使用的功能。

Claims (10)

1.一种计算光刻投影物镜中朗奇剪切干涉图像的方法，用于计算并得出被测投影物镜的波像差，其特征在于：该方法包括以下步骤，
步骤S1：通过计算得出射入物方光栅的光的标量电场函数及物方光栅的物方光栅频谱以得出被物方光栅衍射后的光的电场频谱空间分布函数；
步骤S2：光通过被测投影物镜，计算得出被测投影物镜的光瞳函数并结合所述电场频谱空间分布函数算得光的出瞳处电场分布表达式；
步骤S3：光照射在像方光栅上表面时得出其偏振信息，并通过几何旋转矩阵R3×2将出瞳处电场分布表达式从局部坐标系的sp分量转化为全局坐标系的xyz分量以计算得出入射光的像方光栅上表面电场表达式；
步骤S4：计算得出RS积分核函数，并将像方光栅上表面电场表达式进行计算得出光的像方光栅下表面电场表达式，基于像方光栅下表面电场表达式并结合RS积分核函数进行计算以得到光在接收屏上形成的剪切干涉图形；
步骤S5：定义IMM(fx,fy；f′x,f′y)变量矩阵并对其进行本征值分解以重写多个剪切干涉图形并得到分解结果，通过重复计算所述分解结果以得出光在接收屏上不同剪切位置的精确剪切干涉图像。
2.如权利要求1所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：照射到物方光栅上的光为光源输出的多方向的入射光，所述步骤S1进一步包括以下步骤；
步骤S11：通过计算得出照射在物方光栅上的光的光源频率(sx,sy)及空间频率(fx,fy)并得出射入平行光的标量电场函数为U(x,y)；
步骤S12：测得物方光栅透视率函数为m1(x′,y′)并对其进行傅里叶变换以得出物方光栅频谱为M1(fx,fy)；
步骤S13：进一步算出经过物方光栅衍射后的光的电场频谱空间分布函数为M1(fx-sx,fy-sy)，其表达式为：
Figure FDA0002740555080000021
在步骤S11中光源频率(sx,sy)用于描述光的入射方向，标量电场函数U(x,y)的表达式为：
Figure FDA0002740555080000022
3.如权利要求2所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：在所述步骤S2中，所述光瞳函数为P(fx,fy)，其函数表达式为：
Figure FDA0002740555080000023
计算得出光在出瞳处成像时的辐射修正因子，其表达式为：
Figure FDA0002740555080000024
将所述电场频谱空间分布函数、所述光瞳函数及所述辐射修正因子依次相乘以得到所述出瞳处电场分布表达式：
Figure FDA0002740555080000025
4.如权利要求1所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：所述步骤S3包括以下步骤：
步骤S31：所述几何旋转矩阵R3×2为3x2维的几何旋转矩阵，并通过琼斯矢量将所述出曈处电场分布表达式扩展为矢量形式以得到出曈处电场分布矢量表达式；
步骤S32：将所述出曈处电场分布矢量表达式的局部坐标系的sp分量转化为全局坐标系的xyz分量以计算得出所述像方光栅上表面电场表达式；
在步骤S31中，对于x偏振所述琼斯矢量为
Figure FDA0002740555080000026
对于y偏振所述琼斯矢量为
Figure FDA0002740555080000027
利用所述几何旋转矩阵R3×2得出所述出瞳处电场分布表达式的矢量形式表达式为：
Figure FDA0002740555080000031
在步骤S32中，所述像方光栅上表面电场表达式为：
Figure FDA0002740555080000032
5.如权利要求1所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：所述像方光栅下表面电场表达式用于描述光在经过像方光栅后的电场分布，所述像方光栅下表面电场表达式为：
Figure FDA0002740555080000033
6.如权利要求5所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：所述步骤S4包括以下步骤：
步骤S41：基于所述像方光栅下表面电场表达式并结合RS积分核函数经过衍射积分计算得出接收屏电场分布函数表达式，所述接收屏电场分布函数表达式为：
Figure FDA0002740555080000034
步骤S42：将所述接收屏电场分布函数表达式中的(ε,η)自由度通过积分消除后进行快速傅里叶计算以得到光在接收屏上形成电场的标量电场函数；
所述标量电场函数的表达式为:
Figure FDA0002740555080000035
7.如权利要求6所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：所述RS积分核函数为Rayleigh-Sommerfeld积分核函数，所述RS积分核函数的表达式为：
Figure FDA0002740555080000036
在所述RS积分核函数的表达式中，r的函数表达式为：
Figure FDA0002740555080000037
8.如权利要求7所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：所述步骤S4进一步包括以下步骤：
步骤S43：通过计算得到用于描述光源在频率空间中强度分布的光源函数I(sx,sy)，将所述标量电场函数进行非相干叠加以得到多个入射方向的光在接收屏上形成的剪切干涉图形，其叠加权重为光源函数I(sx,sy)；
所述剪切干涉图形的表达式为：
Figure FDA0002740555080000041
9.如权利要求8所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：所述步骤5进一步包括以下步骤：
步骤S51：将IMM(fx,fy；f′x,f′y)变量矩阵本征值分解得到的所述分解结果储存以供多次计算使用，所述分解结果包括通过本征值分解得到的本征值及本征向量。
10.如权利要求9所述的一种计算光刻投影物镜中朗奇剪切干涉图像的方法，其特征在于：在所述步骤S5中，所述IMM(fx,fy；f′x,f′y)变量矩阵为厄密的矩阵且本征值为实数，其函数表达式为：IMM(fx,fy；f′x,f′y)＝∫dsxdsyI(sx,sy)M(fx-sx,fy-sy)M*(f′x-sx,fy′-sy)；
对所述IMM(fx,fy；f′x,f′y)变量矩阵进行本征值分解可得：
Figure FDA0002740555080000042
以此将所述剪切干涉图形的表达式重写以得到所述精确剪切干涉图像的表达式：
Figure FDA0002740555080000043
CN202011148935.6A 2020-10-23 2020-10-23 一种计算光刻投影物镜中朗奇剪切干涉图像的方法 Active CN112114501B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202011148935.6A CN112114501B (zh) 2020-10-23 2020-10-23 一种计算光刻投影物镜中朗奇剪切干涉图像的方法
PCT/CN2021/125704 WO2022083735A1 (zh) 2020-10-23 2021-10-22 一种计算光刻投影物镜中朗奇剪切干涉图像的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011148935.6A CN112114501B (zh) 2020-10-23 2020-10-23 一种计算光刻投影物镜中朗奇剪切干涉图像的方法

Publications (2)

Publication Number Publication Date
CN112114501A true CN112114501A (zh) 2020-12-22
CN112114501B CN112114501B (zh) 2023-06-02

Family

ID=73794345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011148935.6A Active CN112114501B (zh) 2020-10-23 2020-10-23 一种计算光刻投影物镜中朗奇剪切干涉图像的方法

Country Status (2)

Country Link
CN (1) CN112114501B (zh)
WO (1) WO2022083735A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022083735A1 (zh) * 2020-10-23 2022-04-28 深圳晶源信息技术有限公司东方晶源微电子科技（北京）有限公司深圳分公司 一种计算光刻投影物镜中朗奇剪切干涉图像的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102681365A (zh) * 2012-05-18 2012-09-19 中国科学院光电技术研究所 一种投影物镜波像差检测装置及方法
US20130063805A1 (en) * 2011-05-24 2013-03-14 Craig B. Arnold Tunable acoustic gradient index of refraction lens and system
CN105723283A (zh) * 2013-11-20 2016-06-29 Asml荷兰有限公司 光刻方法和设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0240177B2 (ja) * 1982-12-24 1990-09-10 Rikagaku Kenkyusho Kogakuhikyumennokensaho
CN102889980B (zh) * 2012-10-15 2015-07-15 中国科学院光电技术研究所 一种基于光栅剪切干涉检测系统的微透镜定焦检测方法
CN104111120B (zh) * 2014-07-25 2017-05-31 中国科学院上海光学精密机械研究所 基于朗奇剪切干涉仪的相位提取方法
CN107367906A (zh) * 2017-08-16 2017-11-21 吴飞斌 一种光学系统的检焦装置和检焦方法
CN112114501B (zh) * 2020-10-23 2023-06-02 东方晶源微电子科技（北京）有限公司深圳分公司 一种计算光刻投影物镜中朗奇剪切干涉图像的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130063805A1 (en) * 2011-05-24 2013-03-14 Craig B. Arnold Tunable acoustic gradient index of refraction lens and system
CN102681365A (zh) * 2012-05-18 2012-09-19 中国科学院光电技术研究所 一种投影物镜波像差检测装置及方法
CN105723283A (zh) * 2013-11-20 2016-06-29 Asml荷兰有限公司 光刻方法和设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JOSEPH BRAAT等: "Improved Ronchi test with extended source", 《JOURNAL OF THE OPTICAL SOCIETY OF AMERICA A-OPTICS IMAGE SCIENCE AND VISION》 *
YAO, ZHENGPENG等: "Vector analysis of two-dimensional", 《INTERNATIONAL SYMPOSIUM ON PHOTOELECTRONIC DETECTION AND IMAGING 2013: MICRO/NANO OPTICAL IMAGING TECHNOLOGIES AND APPLICATIONS》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022083735A1 (zh) * 2020-10-23 2022-04-28 深圳晶源信息技术有限公司东方晶源微电子科技（北京）有限公司深圳分公司 一种计算光刻投影物镜中朗奇剪切干涉图像的方法

Also Published As

Publication number Publication date
CN112114501B (zh) 2023-06-02
WO2022083735A1 (zh) 2022-04-28

Similar Documents

Publication Publication Date Title
CN106575088B (zh) 三维测量光刻掩模的3d空间像的方法
US10162161B2 (en) Ptychography imaging systems and methods with convex relaxation
CN102681365B (zh) 一种投影物镜波像差检测装置及方法
CN110095968B (zh) 一种离轴数字全息图像重建装置和方法及显微成像系统
KR102416784B1 (ko) 임의의 입사각에서의 코히어런트 회절 영상
US20140347672A1 (en) Apparatus and method for quantitive phase tomography through linear scanning with coherent and non-coherent detection
WO2003050472A1 (en) Systems and methods for wavefront measurement
CN105765463B (zh) 用于计算结构的电磁散射性质及用于估计其几何和材料参数的方法和装置
AU2008252706A1 (en) Three dimensional imaging
CN106802233B (zh) 一种微透镜阵列测试装置及方法
CN101482503B (zh) 一种相干衍射成像方法及其处理装置
KR20220006584A (ko) 리소그래피 마스크의 에어리얼 이미지를 3차원적으로 결정하는 방법
US20220381643A1 (en) Method for determining an imaging quality of an optical system when illuminated by illumination light within an entrance pupil to be measured
JP2018521337A (ja) リソグラフィのｉｎ−ｓｉｔｕ高速高空間解像度の波面収差測定デバイス及び測定方法
CN102269926B (zh) 基于Abbe矢量成像模型的非理想光刻系统OPC的优化方法
CN102636882A (zh) 一种分析高数值孔径成像系统空间像的方法
Stockton et al. Tomographic single pixel spatial frequency projection imaging
JP6000773B2 (ja) 収差推定方法、プログラムおよび撮像装置
CN102323723B (zh) 基于Abbe矢量成像模型的光学邻近效应校正的优化方法
CN102236268B (zh) 基于空间像频谱的光刻投影物镜波像差检测方法
CN112114501A (zh) 一种计算光刻投影物镜中朗奇剪切干涉图像的方法
CN102323722B (zh) 基于Abbe矢量成像模型获取掩膜空间像的方法
CN102323721B (zh) 基于Abbe矢量成像模型获取非理想光刻系统空间像的方法
CN110927116B (zh) 一种测量标记结构的方法、装置及系统
Guzhov et al. Method of increasing the spatial resolution in digital holographic microscopy

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20240710

Address after: 100176 building 12, yard 156, Jinghai 4th Road, Beijing Economic and Technological Development Zone, Daxing District, Beijing

Patentee after: Dongfang Jingyuan Microelectronics Technology (Beijing) Co.,Ltd.

Country or region after: China

Address before: No.518000, Shenzhen Futian Free Trade Zone

Patentee before: SHENZHEN BRANCH DONGFANG JINGYUAN MICROELECTRONICS TECH (BEIJING) Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right