CN105629461B - 一种百纳米尺度超细光针场聚焦设计方法 - Google Patents
一种百纳米尺度超细光针场聚焦设计方法 Download PDFInfo
- Publication number
- CN105629461B CN105629461B CN201610015338.3A CN201610015338A CN105629461B CN 105629461 B CN105629461 B CN 105629461B CN 201610015338 A CN201610015338 A CN 201610015338A CN 105629461 B CN105629461 B CN 105629461B
- Authority
- CN
- China
- Prior art keywords
- micro
- field
- metal film
- ultra
- structure metal
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Polarising Elements (AREA)
Abstract
本发明公开了一种百纳米尺度超细光针场聚焦设计方法,属于纳米光子学聚焦及微加工技术领域,包括四个步骤:第一、利用矢量角谱理论来积分描述矢量偏振光束照明微结构金属膜环带片时其后光场的衍射传播行为;第二、利用等价数值孔径NAeq和归一化中心遮挡因子ε约束聚焦光束横向尺度,而在轴向引入超高斯函数约束形成超细光针场,由此建立非线性约束优化模型;第三、利用配置的遗传算法和快速汉克尔变换算法编程求解上述优化模型;第四、设置优化初始参数,包括结构参数和算法参数,多次执行优化算法,优选微结构金属膜环带片。以本发明为基础灵活设计不同偏振光束照明、多种尺度微结构金属膜环带片,可应用于实现微纳光刻、纳米打印、超分辨显微成像等。
Description
技术领域
本发明属于纳米光子学聚焦及微加工技术领域,特别涉及一种百纳米尺度超细光针场聚焦设计方法。
背景技术
2008年以来,超细光针场聚焦问题引起科学界广泛的关注,成为国际前沿热点研究课题之一。超细光针场具有重要的工程应用价值,如在激光直写微光刻、光学数据存储、扫描光学纳米显微技术、原子光学、微观粒子加速、光学操控等领域。
目前,超细光针场聚焦调制方法主要分为三大类,分别基于折射光学***、反射光学***和衍射光学***。折射光学***调制方法以2008年新加坡数据存储研究院提出的多环带光瞳滤波、矢量光束调制、折射透镜聚焦为代表(参见文献H.Wang,L.Shi,B.Lukyanchuk,C.Sheppard,C.T.Chong.Creation of a needle of longitudinallypolarized light in vacuum using binary optics.Nature Photonics,2008,2:501-505),这类方法研究最为集中,优点是理论清晰、设计灵活,缺点是径向偏振矢量光束产生***复杂、环带形光瞳滤波器加工要求苛刻、高端复消色差显微物镜成本昂贵。第二类方法是基于反射光学***,主要是采用环形光瞳滤波和抛物反射镜聚焦来实现(参见文献H.Dehez,A.April,M.Piche.Needles of longitudinally polarized light:guidelinesfor minimum spot size and tunable axial extent.Optics Express,2012,20:14891-14905),优点是理论上容易实现极限数值孔径(NA=1,空气中)聚焦,缺点是聚焦光针场在入射光场一侧,且依赖于超精密非球面加工工艺,金属抛物反射镜面形精度(PV)加工要求在λ/8~λ/6量级。第三类方法直接采用单一平面片型微结构进行聚焦(参见文献T.Liu,J.Tan,J.Liu,H.Wang.Modulation of a super-Gaussian optical needle by high-NAFresnel zone plate.Optics Letters,2013,38:2742-2745),这是一种无透镜、多光束衍射干涉聚焦***,突出优点是结构轻巧、成本低廉、设计灵活,而由于菲涅尔波带片在短波长及大数值孔径情形下,外圈环带径向宽度接近且径向宽度在数十纳米,对现有微加工刻蚀工艺提出了苛刻的要求。
上述问题是目前超细光针场调制面临的实际困难,特别是现有国内外报道尚未有一种可行的方法能够实现横向尺度接近百纳米量级的超细光针场聚焦。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种百纳米尺度超细光针场聚焦设计方法,基于微结构金属膜环带片,采用短波长激光照明,从深紫外(DUV)至极紫外(EUV),具有成本低廉、***结构简单、设计灵活方便的特点,可用于微纳光刻、纳米打印、超分辨显微成像、高密度光学存储等领域。
为了实现上述目的,本发明采用的技术方案是:
一种百纳米尺度超细光针场聚焦设计方法,包括如下步骤:
步骤一,描述线偏振光和圆偏振光照明微结构金属膜环带片条件下光场的衍射传播过程,计算微结构金属膜环带片后电场各分量及总电能密度分布;
步骤二,建立非线性约束优化模型,优化决策标量为微结构金属膜环带片各环透过率函数ti和归一化中心遮挡半径ε;
步骤三,利用遗传算法和快速汉克尔变换算法求解优化模型;
步骤四,设置优化初始参数,执行优化算法,优选微结构金属膜环带片。
与现有技术相比,本发明的有益效果是:
本发明***简单、结构轻巧,仅利用单一平面微结构金属膜环带片;利用现有成熟的微加工工艺即可实现,成本低廉,可批量制作;设计方法灵活,适用于多种紫外激光波长照明、不同径向环带宽度、任意直径和工作距离设计。
附图说明
图1为本发明中的微结构金属膜环带片光场衍射聚焦示意图,图中A表示微结构金属膜环带片,B表示观察平面。
图2为本发明微结构金属膜环带片M1对应的径向振幅透过率示意图。
图3为本发明微结构金属膜环带片M1对应的轴向归一化强度(光针)场分布示意图。
图4为本发明微结构金属膜环带片M1对应的Y-Z平面内的光针场强度分布示意图。
图5为本发明微结构金属膜环带片M2对应的径向振幅透过率示意图。
图6为本发明微结构金属膜环带片M2对应的轴向归一化强度(光针)场分布示意图。
图7为本发明微结构金属膜环带片M2对应的X-Z平面内的光针场强度分布示意图。
具体实施方式
下面结合附图和实施例详细说明本发明的实施方式。
如图1所示,采用微结构金属膜环带片直接聚焦短波长紫外激光光束,在距离环带片后表面数十微米处调制产生纳米尺度超细光针场分布,利用矢量角谱理论(VectorialAngular Spectrum Theory)进行光场传播分析。
(1)矢量角谱理论光场积分表示
假设沿X轴方向振动的线偏振光(LPB)沿Z轴正向传播,如图1所示,经微结构金属膜环带片衍射后,在z>0垂轴平面内任意一点位置处电场E的直角分量根据矢量角谱理论推导得到
式中,Ex(r,z)表示x向分量,Ey(r,z)表示y向分量,表示z向分量,q(l)=(1/λ2-l2)1/2,l表示径向空间频率分量;在介质中波长为λ=λ0/μ,λ0为真空波长,μ是浸入介质折射率;J0和J1分别是第一类零阶和一阶贝塞尔函数,j是虚数单位;A0(l)表示为:
式中,t(r)表示微结构金属膜环带片对应的振幅透过率函数,不透光时为0,透光时为1;g(r)表示矢量照明光束在微结构金属膜环带片孔径平面内的光场振幅分布,这里假设为均匀平面波照射,即g(r)=1;由公式(1)得到微结构金属膜环带片后总电能密度即光强为
当照明光束是(左旋)圆偏振光(CPB)时,同理可推导得到电场E的各分量为
式中,A0(l)由公式(2)给出。光场强度分布为Ivas(r,z)=2|Ex(r,z)|2+|Ez(r,z)|2。显然对于圆偏振光,强度分布满足圆对称性,有别于线偏振光情形。
(2)建立多变量、单目标、非线性约束优化模型
根据公式(1)或公式(3),给定所需的聚焦超细光针场基本特性,建立多变量、单目标、非线性约束优化问题。必须在横向(X-Y)和轴向(Z)两方面对超细光针场进行约束,本发明采用如下方案:在横向,规定光针场中间位置zf对应的等价数值孔径(NAeq),同时对微结构金属膜环带片中心近轴区域进行遮挡,以抑制近轴低频光场,由此实现光针场横向宽度的进一步约束;在轴向,引入具有平顶分布特性的超高斯函数(Super-Gaussian Function)来限制轴向光针场强度分布,最后得到的单目标优化问题数学描述为:
Subject to:
Ical=Ivas(0,Z;T)
T={ti},Z={zn}
|zn-zf|≤DOF/2
εmin≤ε≤εmax
式中,参数NAeq=ηsin{[tan-1[D/(2zf)]},K=Ivas(0,zf);Nsg是超高斯函数阶次,取为2至5的整数,DOF是设计的光针场轴向焦深,D表示微结构金属膜环带片直径;ti表示各环的透过率,为0或1;ε是归一化中心遮挡半径,0<ε<1,且0<εmin≤εmax<1,表示对微结构金属膜环带片进行一定程度的中心遮挡,可有效实现轴向光针场拉伸和横向光束的压缩;zf是光针场的轴向居中位置,zn=zf±nΔz是参与评价的轴向采样点(n=0,1,2,…Nzp),Nzp决定取样点数,取2以上的整数,Δz=DOF/(2Nzp+1)。利用均方根误差RMSE(Root Mean Square Error)来评价实际光场分布与优化目标的逼近程度。在横向利用FWHM(半高全宽)、在轴向利用DOF(焦深)来评价超细光针场的基本特性。
(3)利用遗传算法和快递汉克尔变换算法求解优化问题
采用遗传算法(Genetic Algorithm,GA)和快速汉克尔变换算法(Fast HankelTransform Algorithm)来求解公式(4)描述的优化问题。GA是一种经典进化算法,模拟生物进化“物竞天择,适者生存”的自然法则,以进化论和遗传定律为理论依据,具有内在隐并行性和全局搜索能力。遗传算法的基本操作包括:种群的初始化配置、遗传编码、个体适应度评价、个体的选择、交叉和变异等。适应度函数即为公式(4)中的目标函数RMSE,采用{0,1}二进制代码直接编码各环带的透过率ti、中心遮挡因子ε和轴向采用间隔zn等。将GA配置成两点交叉和两点变异,同时保留精英选择策略,即将父代中直接适应度最高的个体直接遗传至子代中,将种群规模将种群规模设定为M=10~60,选取交叉概率为Pc≥0.75,变异概率为Pm≤0.15,迭代次数Ni≥10。通过大量计算实例表明,上述配置条件下,GA能够快速稳定收敛至近似最优解。
在评价任意个体的适应度时,需要执行大量的零阶和一阶汉克尔变换,因此整个优化过程中汉克尔变换的计算效率及计算精度是设计方法的关键,为了加速优化运算,编程实现一种快速汉克尔变换算法(参考文献A.E.Siegman.Quasi Fast HankelTransform.Optics Letters,1977,1:13-15),该算法具有计算速度快、精度高、极低计算机存储要求等显著优点,在标准汉克尔变换积分表达式中,利用非线性指数函数变量替换,将标准单边汉克尔变换表示为双边互相关积分,经过这样的变换后可以利用傅里叶变换计算互相关。
(4)设计结果及实施例
193nm准分子浸入式光刻***是目前100nm以下线宽尺度微纳刻蚀加工的主流方法,因此采用波长λ0=193nm,选用水浸(η=1.44),照明激光偏振态分别取LPB和CPB,微结构金属膜环带片径向基本环带宽设为Δr=200nm,各环带具体是开(透光)还是关(遮挡)状态由优化算法决定。典型环带片结构参数见表1。
表1 优化设计的微结构金属膜环带片参数
表1中,器件M1直径50μm,其结构如图2所示。采用X偏振LPB(线偏振光束)照明,工作距WD为14.5μm,等价数值孔径为NAeq=1.25,轴向焦深12.4λ0,如图3所示在13~16μm范围内的轴向归一化强度(光针)分布,图4所示为M1对应的Y-Z平面内的强度分布。在14.5μm位置处横向光斑不对称,在垂直振动方向光斑FWHMy=0.31λ0,而在X方向光斑较宽达到FWHMx=0.67λ0。
表1中,器件M2直径24μm,其结构如图5所示。采用CPB(圆偏振光束)照明,工作距WD为10.5μm,等价数值孔径为NAeq=1.08,轴向焦深6.6λ0,如图6所示在8.5~12.5μm范围内的轴向归一化强度(光针)分布,图7所示为M2对应的X-Z平面内的强度分布。由于是CPB照明,所以强度场分布满足旋转轴对称性,在10.5μm位置处横向光斑FWHMx,y=0.60λ0。
以上结合附图对本发明的具体实施方式作了说明,但这些说明不能被理解为限制了本发明的范围,本发明的保护范围由随附的权利要求书限定,任何在本发明权利要求基础上的改动都是本发明的保护范围。
Claims (8)
1.一种百纳米尺度超细光针场聚焦设计方法,其特征在于,包括如下步骤:
步骤一,描述线偏振光和圆偏振光照明微结构金属膜环带片条件下光场的衍射传播过程,计算微结构金属膜环带片后电场各分量及总电能密度分布;
步骤二,建立非线性约束优化模型,优化决策标量为微结构金属膜环带片各环透过率函数ti和归一化中心遮挡半径ε;
步骤三,利用遗传算法和快速汉克尔变换算法求解优化模型;
步骤四,设置优化初始参数,执行优化算法,优选微结构金属膜环带片。
2.根据权利要求1所述百纳米尺度超细光针场聚焦设计方法,其特征在于,所述步骤一中,光场的衍射传播过程根据矢量角谱理论来描述,所述电场各分量均表示为汉克尔变换的一维积分形式。
3.根据权利要求1所述百纳米尺度超细光针场聚焦设计方法,其特征在于,所述步骤一中,当照明光束是沿X方向振动、沿Z轴正向传输的线偏振光时,则经微结构金属膜环带片衍射后,在z>0的垂轴平面内任意一点位置处电场E的各直角分量为:
式中,Ex(r,z)、Ey(r,z)和分别是x向、y向和z向分量;q( l )=(1/λ2- l 2)1/2, l 表示径向空间频率分量;在介质中波长为λ=λ0/μ,λ0为真空波长,μ是浸入介质折射率;J0和J1分别是第一类零阶和一阶贝塞尔函数,j是虚数单位;A0(l)表示角谱,表示为式中,t(r)表示微结构金属膜环带片对应的振幅透过率函数,不透光时为0,透光时为1;g(r)表示矢量照明光束在微结构金属膜环带片孔径平面内的光场振幅分布,当为均匀平面波照射时,g(r)=0;则微结构金属膜环带片后总电能密度即光强为
当照明光束是左旋圆偏振光时,则其电场E的各分量表示为
光场强度分布为Ivas(r,z)=2|Ex(r,z)|2+|Ez(r,z)|2。
4.根据权利要求1所述百纳米尺度超细光针场聚焦设计方法,其特征在于,所述步骤三中,配置遗传算法,对ti和ε进行二值编码,采用两点交叉和两点变异,并执行精英选择策略,同时利用快速汉克尔变换算法加速积分数值计算。
5.根据权利要求1所述百纳米尺度超细光针场聚焦设计方法,其特征在于,所述步骤三中,将种群规模设定为M=10~60,选取交叉概率为Pc≥0.75,变异概率为Pm≤0.15,迭代次数Ni≥10。
6.根据权利要求1或4所述百纳米尺度超细光针场聚焦设计方法,其特征在于,所述步骤三中,所述快速汉克尔变换算法,是指在标准汉克尔变换积分表达式中,利用非线性指数函数变量替换,将标准单边汉克尔变换表示为双边互相关积分,从而实现利用傅里叶变换计算互相关。
7.根据权利要求1所述百纳米尺度超细光针场聚焦设计方法,其特征在于,所述步骤四中,优化初始参数包括***结构参数:照明波长和偏振态、浸入介质折射率μ、微结构金属膜环带片径向基本宽度Δr和环带数N、等价数值孔径NAeq、光针场轴向焦深DOF,同时设定遗传算法交叉概率Pc、变异概率Pm、种群规模M、迭代次数Ni。
8.根据权利要求1或7所述百纳米尺度超细光针场聚焦设计方法,其特征在于,所述步骤四中,照明波长,从紫外至极紫外波段,取10nm<λ0<400nm,偏振态为线偏振和圆偏振光;微结构金属膜环带片径向基本宽度取为Δr≥25nm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610015338.3A CN105629461B (zh) | 2016-01-11 | 2016-01-11 | 一种百纳米尺度超细光针场聚焦设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610015338.3A CN105629461B (zh) | 2016-01-11 | 2016-01-11 | 一种百纳米尺度超细光针场聚焦设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105629461A CN105629461A (zh) | 2016-06-01 |
CN105629461B true CN105629461B (zh) | 2017-12-26 |
Family
ID=56044558
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610015338.3A Active CN105629461B (zh) | 2016-01-11 | 2016-01-11 | 一种百纳米尺度超细光针场聚焦设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105629461B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111240011A (zh) * | 2020-01-09 | 2020-06-05 | 西安交通大学 | 一种不等环宽金属膜超振荡环带片设计方法 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106773030B (zh) * | 2016-11-24 | 2018-12-18 | 西安交通大学 | 一种均匀光针聚焦微结构及设计方法 |
CN106733573B (zh) * | 2016-12-15 | 2019-10-15 | 南京大学 | 一种三维宽带能量聚焦装置 |
CN110007460B (zh) * | 2018-01-05 | 2021-08-17 | 上海三思电子工程有限公司 | 透镜制作方法、计算机设备及存储介质 |
CN108169920B (zh) * | 2018-01-24 | 2020-03-20 | 河南科技大学 | 可传输吞吐微粒的椭圆光学传送带光束掩模板的设计方法 |
CN109613698B (zh) * | 2019-01-03 | 2020-03-17 | 西安交通大学 | 一种金属膜超振荡环带片的设计方法 |
CN111222287B (zh) * | 2020-01-09 | 2022-04-05 | 西安交通大学 | 一种不等环宽介质超振荡环带片设计方法 |
CN114137720A (zh) * | 2021-10-15 | 2022-03-04 | 中国科学院国家天文台南京天文光学技术研究所 | 一种用于高对比度成像的光瞳调制星冕仪***及工作方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100451689C (zh) * | 2007-03-21 | 2009-01-14 | 中国科学院上海光学精密机械研究所 | 偏振透镜及其制备方法 |
US7778498B2 (en) * | 2008-02-12 | 2010-08-17 | Ofs Fitel Llc | Systems and techniques for generating cylindrical vector beams |
CN101290396A (zh) * | 2008-05-14 | 2008-10-22 | 中国科学技术大学 | 差值反馈优化的衍射光学元件 |
CN104297925B (zh) * | 2014-10-11 | 2016-07-13 | 北京工业大学 | 一种实现飞秒激光长焦深的折衍混合元件的设计方法 |
-
2016
- 2016-01-11 CN CN201610015338.3A patent/CN105629461B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111240011A (zh) * | 2020-01-09 | 2020-06-05 | 西安交通大学 | 一种不等环宽金属膜超振荡环带片设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105629461A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105629461B (zh) | 一种百纳米尺度超细光针场聚焦设计方法 | |
Kazanskiy | Research and education center of diffractive optics | |
CN109613698A (zh) | 一种金属膜超振荡环带片的设计方法 | |
CN106773030B (zh) | 一种均匀光针聚焦微结构及设计方法 | |
Kotlyar et al. | Photonic nanojets generated using square-profile microsteps | |
CN109656019B (zh) | 一种电介质超振荡环带片的设计方法 | |
US10416383B2 (en) | Spatial control of the optical focusing properties of photonic nanojets | |
Khonina et al. | Analysis of the formation of a longitudinally polarized optical needle by a lens and axicon under tightly focused conditions | |
CN106547091B (zh) | 一种由轴向分辨率设计大数值孔径菲涅尔波带片的方法 | |
CN100547440C (zh) | 一种用于双光子微细加工的三维超分辨衍射光学器件及其设计方法 | |
CN102269837A (zh) | 一种实现激光束纵向偏振双焦点聚焦的位相片 | |
Liu et al. | Focusing far-field nanoscale optical needles by planar nanostructured metasurfaces | |
Savelyev | The comparison of laser radiation focusing by diffractive axicons and annular gratings with variable height using high-performance computer systems | |
Slevas et al. | Laser induced modifications in transparent materials using azimuthally modulated axicon beams | |
Piccirillo et al. | Geometric-phase waveplates for free-form dark hollow beams | |
Chen et al. | Generation of shaping nondiffracting structured caustic beams on the basis of stationary phase principle | |
Danilov | Focusing DOEs (focusators): design and investigation | |
Zhang et al. | Laser beam homogenization with randomly distributed freeform cylindrical microlens | |
Khorin et al. | Influence of 3D helical microstructure shape deviations on the properties of a vortex beam generated in the near diffraction zone | |
CN102402006A (zh) | 一种用于产生空心光束的位相型光瞳滤波器 | |
Xie et al. | Research on beam shaping based on grating bars | |
CN110967842A (zh) | 基于光镊技术的局域空心光束自由开闭*** | |
Liu et al. | A new method to produce optical nano-needle | |
Khorin et al. | Modeling diffraction of a polarized light by three-dimensional nonlinear spiral phase in the near zone | |
Wu et al. | Superoscillatory metalens for an azimuthally polarized wave with different orbital angular momentum |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |